- Open Access
Different modulation effects of 1 Hz and 20 Hz transcutaneous auricular vagus nerve stimulation on the functional connectivity of the periaqueductal gray in patients with migraine
Journal of Translational Medicine volume 19, Article number: 354 (2021)
A growing body of evidence suggests that transcutaneous auricular vagus nerve stimulation (taVNS) may relieve symptoms of migraineurs. Frequency is one of the key stimulation parameters. The aim of this study is to investigate the modulation effect of taVNS frequency on the descending pain modulation system (DPMS) in patients with migraine.
Twenty-four episodic migraineurs without aura (21 females) were recruited for the single-blind, crossover, functional magnetic resonance imaging (fMRI) study. Each participant attended two separate fMRI scan sessions, one for 1 Hz and another for 20 Hz taVNS, in a random order. Seed-based functional connectivity analysis was applied using the ventrolateral periaqueductal gray (PAG) as the region of interest.
Compared with the pre-taVNS resting state, continuous 1 Hz taVNS (during) produced a significant increase in functional connectivity between the PAG and the bilateral middle cingulate cortex (MCC), right precuneus, left middle frontal gyrus (MFG), and left cuneus. Compared with 20 Hz taVNS, 1 Hz taVNS produced greater PAG connectivity increases with the MCC, right precuneus/posterior cingulate cortex, left insula, and anterior cingulate cortex (ACC). A significant negative correlation was observed between the number of migraine attacks in the previous 4 weeks and the PAG-MCC functional connectivity in the pre-taVNS resting-state before 1 Hz taVNS.
Our findings suggest that taVNS with different frequencies may produce different modulation effects on the descending pain modulation system, demonstrating the important role of stimulation frequency in taVNS treatment.
The vagus nerve consists of a complex system that may regulate pain, mood, and the neuro-endocrine-immune axis [1,2,3,4,5,6,7]. Thus, stimulating the vagus nerve to modulate the function of the nerve and related organs has drawn the attention of clinicians and investigators for a long time. Anatomical studies found peripheral branches of the vagus nerve distributed on the ear [8, 9], and according to the bottom-up mechanism of the central nervous system, the propagation of electrical stimuli may follow an afferent path from the peripheral nerves towards the brain stem and central structures [10, 11]. Thus, direct stimulation of the nerve fibers on the ear may produce an effect similar to classic vagus nerve stimulation. This plausibility has led to the development of transcutaneous auricular vagus nerve stimulation (taVNS), a non-invasive, low-cost, and easily implementable alternative to classic vagus nerve stimulation [12,13,14,15]. A growing body of evidence suggests that taVNS can induce antinociception, which may affect peripheral and central nociception, inflammatory responses, autonomic activity, and pain-related behavior [1, 16,17,18].
While taVNS has demonstrated its potentials, the optimal parameters for taVNS, such as frequency, remain unclear . Accumulating evidence suggests different frequencies may be associated with different physiological and treatment effects. For instance, investigators compared the effect of 2, 10, and 20 Hz stimulation on heart rate in healthy subjects, and they found that both 10 and 20 Hz could decrease heart rate . Furthermore, studies suggest that the optimal taVNS frequency may vary across different disorders. For example, a recent clinical research study on taVNS treatment of drug-resistant epilepsy showed a significant reduction in seizure frequency in patients of the 25 Hz group compared to the 1 Hz group . However, in another clinical study of migraine patients, investigators found that although both 1 Hz and 25 Hz taVNS improved clinical outcomes in patients with chronic migraine, 1 Hz taVNS produced greater improvement . Nevertheless, the underlying mechanism of different taVNS frequencies remains unclear.
Recently, brain imaging has been widely used to investigate the central mechanism of taVNS, and these studies demonstrate that intermittent taVNS can modulate activity of certain brain regions consistent with the vagus nerve central projections [22,23,24,25,26,27,28]. For instance, investigators have assessed brainstem fMRI response to 2, 10, 25, and 100 Hz taVNS in healthy individuals, and found that the strongest brainstem response was evoked by 100 Hz stimulation . In recent studies, we also applied the resting-state functional connectivity method to investigate the functional connectivity alteration during “the continuous taVNS” (20 Hz) and found that taVNS can modulate the functional connectivity of the ventral striatum and hypothalamus [30, 31].
Nevertheless, the neural substrates underlying frequency have rarely been investigated in a patient population such as migraine; elucidating how different frequencies can modulate pathways associated with migraine may further facilitate the development of this promising neuromodulation method.
Recently, the role of descending pain modulatory system (DPMS) in pain modulation and the physiopathology of chronic pain has drawn more and more attention [32,33,34]. Yet, investigating the functional status of the DPMS in humans remains a challenge. In an earlier study , we investigated the resting state functional connectivity (rsFC) of the periaqueductal grey (PAG), a key region in the DPMS in healthy subjects and found significant rsFC between the PAG and central regions of the DPMS, such as the anterior cingulate cortex (ACC), rostroventral medulla (RVM) and anterior insula, demonstrating the feasibility of using functional connectivity methods to non-invasively investigate the DPMS in humans.
Following the study, the PAG functional connectivity has been applied to investigate the physiopathology of chronic pain disorders including migraine [35,36,37,38,39,40,41,42,43], menstrual pain [44,45,46], postherpetic neuralgia , fibromyalgia , myofascial pain , visceral pain , low back pain , and neck pain . Further, studies have also shown that effective treatment can significantly modulate the PAG functional connectivity in patients with migraine , chronic low back pain , and knee osteoarthritis . We also found that continuous electroacupuncture stimulation alters PAG functional connectivity . Taken together, these findings demonstrate the important role of PAG functional connectivity in pain research.
Thus, in this study, we investigate how continuous taVNS at 1 Hz versus 20 Hz (a relatively low frequency versus a moderate frequency) that are widely applied in taVNS studies  can modulate the PAG functional connectivity in patients with migraine without aura, using a cross-over design. We hypothesize that taVNS at 1 Hz versus 20 Hz may produce greater PAG functional connectivity changes due to its greater improvement in patients with migraine .
The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Ethics Committee of the Second Affiliated Hospital of Guangzhou University of Chinese Medicine (Z2016-079-01). All participants provided written informed consent before starting the study.
Twenty-four episodic migraineurs without aura were recruited in the present study from outpatient neurology clinics of the Second Affiliated Hospital of Guangzhou University of Chinese Medicine. Similar to our previous studies [55, 56], the diagnosis of migraine was based on the International Classification of Headache Disorders, 2nd Edition (ICHD-II), as diagnosed by a specialist working at the neurology outpatient service.
Patients were eligible for participation if they: (1) were 18 to 45 years of age, (2) self-reported being right-handed, (3) have at least 6 months of migraine duration, (4) have at least one headache attack per month, (5) have not taken any prophylactic headache medications during the past 4 weeks, (6) have not taken any psychoactive or vasoactive drugs during the past 3 months. Patients were excluded if there was a/an: (1) headache induced by other diseases, (2) headache attack within 48 h prior to the experiment or during the experiment, (3) pregnant or lactating, (4) any other chronic pain conditions, (5) severe head deformity or intracranial lesions, (6) score on the Self-Rating Depression Scale  or Self-Rating Anxiety Scale  > 50, and (7) inability to provide informed consent for oneself.
A single-blind, crossover functional magnetic resonance imaging (fMRI) trial design was applied in the present study to investigate the modulation effects of 1 Hz and 20 Hz taVNS in patients with migraine without aura. Specifically, each participant attended two taVNS fMRI scan sessions with identical parameters, one for 1 Hz and another for 20 Hz taVNS in a random order (Fig. 1A). Each session was separated by at least 7 days to avoid sensitization to the stimuli. All scans were applied during an interictal period when the participants were free from headache symptoms.
In the current study, we applied taVNS on the participant’s left concha (cymba and cavum, Fig. 1B, Additional file 1: Fig. S1) . The electrical stimulation was performed using the Electronic Acupuncture Treatment Instrument (SDZ IIB, Huatuo, Suzhou, China) with the self-made MRI compatible electrode to deliver electric current at 1 Hz or 20 Hz with a continuous wave (width: ~ 0.2 ms). The 1 Hz/20 Hz taVNS stimulation lasted about 8 min. Similar to our previous studies, stimulation current intensity was adjusted to the strongest nonpainful sensation that participants could tolerate (approximately 4 mA) [31, 55, 60,61,62].
Migraine duration, migraine attacks during the past 4 weeks, and average migraine intensity of the past 4 weeks on the 0 (“not at all”) to 100 (“extremely”) visual analog scale (VAS) were assessed preceding the first MRI scan session. Participants were also asked to complete the Migraine Specific Quality-of Life Questionnaire  to measure the impact of migraine on health‐related quality of life.
MRI data acquisition
All imaging data was acquired at the Second Affiliated Hospital of Guangzhou University of Chinese Medicine using a 3 T MRI System (Siemens MAGNETOM Verio 3.0 T, Erlangen, Germany) with a 24-channel phased-array head coil. Each scan session included a pre-taVNS resting-state fMRI (8 min), the 1 Hz or 20 Hz continuous taVNS (8 min) fMRI, and a post-taVNS resting-state fMRI (8 min).
fMRI scans were acquired with the following parameters: time repetition = 2000 ms, time echo = 30 ms, flip angle = 90°, field of view = 224 mm × 224 mm, matrix size = 64 × 64, slice thickness = 3.5 mm with 0.7 mm inter-slice gap, 31 axial slices paralleled and 240 time points. During the fMRI scans in resting-state and continuous taVNS, participants were asked to stay awake, keep their heads still, eyes closed, and ears plugged and to not think about any particular thing. A T1-weighted structural image was acquired by an isotropic multi-echo magnetization-prepared rapid gradient-echo pulse sequence for anatomic localization of significant signal changes: time repetition = 1900 ms, time echo = 2.27 ms, flip angle = 9°, field of view = 256 mm × 256 mm, data matrix = 256 × 256, and slice thickness = 1.0 mm.
Functional connectivity analysis
Data and calculations of functional connectivity were conducted using the CONN toolbox version 18.b (http://www.nitrc.org/projects/conn) . We used the default preprocessing pipeline for seed-to-voxel functional connectivity analysis. The specific steps were as follows: functional realignment and unwarping, slice timing correction, head motion correction, co-registration of the anatomical image to the mean functional image, segmentation of the anatomical gray matter, white matter, and CSF, normalization to Montreal Neurological Institute (MNI) 152 standard template and smoothing with a 6-mm full width at half maximum (FWHM) kernel. A default frequency window of 0.008 to 0.09 Hz was used for band-pass filtering.
To eliminate correlations caused by head motion and artifacts, we identified outlier time points in the motion parameters and global signal intensity using ART (http://www.nitrc.org/projects/artifact_detect). Images whose composite movement exceeded 0.5 mm or whose global mean intensity was greater than three standard deviations from the mean image intensity were treated as outliers. The time series of the head motion matrix of outliers was also entered as first-level covariates.
Similar to our previous studies [35, 36, 41], we selected the right ventrolateral periaqueductal gray (vlPAG) with a 2 mm radius sphere (MNI coordinates x = 4, y = − 26, z = − 14) as the region of interest (ROI). In addition, we also chose seeds with a 2 mm radius in the fourth ventricle (MNI coordinates: x = 4, y = 10, z = 12; x = − 4, y = 10, z = 12) as a control. Seeds were created using the SPM Wake Forest University Pickatlas toolbox (http://fmri.wfubmc.edu/software/pickatlas) .
In the first-level analysis, we produced a correlation map for each participant by extracting the blood oxygenation level dependent time course separately from the vlPAG and the control seeds and computing Pearson’s correlation coefficients between the time course in the vlPAG/control seeds and every voxel of the whole brain. Correlation coefficients were Fisher transformed into “z” scores to increase normality.
In seed-to-voxel functional connectivity analyses, we first used a pairwise t-test to compare the vlPAG-based functional connectivity between the pre-taVNS resting-state and during continuous taVNS (1 Hz and 20 Hz taVNS, respectively). Next, we compared the difference of vlPAG-based functional connectivity change (during continuous taVNS minus pre-taVNS resting-state) between 1 and 20 Hz taVNS. Finally, we compared the vlPAG-based functional connectivity difference between the pre-taVNS and post-taVNS resting-state between the 1 Hz and 20 Hz taVNS.
For whole brain analysis, a threshold of voxel-wise p < 0.005, and pFDR < 0.05 at cluster level was applied. Also, given the important role of the anterior cingulate cortex (ACC), medial prefrontal cortex (mPFC), insula, amygdala, and thalamus in the DPMS [32, 36, 53, 66] and pathophysiology of migraine [41, 67,68,69,70,71,72], we pre-defined these areas as regions of interest (ROIs), and derived masks of each region from the Automated Anatomical Labeling brain atlas using the Wake Forest University Pickatlas toolbox as ROIs. A threshold of voxel-wise p < 0.005 was used in data analysis. Similar to previous studies [73,74,75], Monte Carlo simulations using the 3dFWHMx and 3dClustSim (as part of the Analysis of Functional NeuroImages program [http://afni.nimh.nih.gov] released in July 2017) were applied for the p value correction for pre-defined ROIs. For each region, the minimum voxel size required for p < 0.05 cluster level p value correction is indicated as the k value in the results presented below.
To explore the association between the initial clinical assessments and the vlPAG-based pre-taVNS resting-state functional connectivity for 1 Hz and 20 Hz respectively, we also selected significantly altered vlPAG-based connectivity clusters (during continuous taVNS minus pre-taVNS resting-state) and extracted the average z-score values of peak MNI of clusters above significance in vlPAG-based pre-taVNS resting-state. Correlation analyses were conducted using the R program in JASP open-source statistical software (Version 0.8.1, http://www.jasp-stats.org), and p values were Bonferroni corrected (see “Results” for details).
Demographic and clinical assessments
Twenty-four participants completed the study and were included in the data analysis [21 females; age 31.33 ± 1.55 years, mean ± standard error (SE)]. No participant reported administration of acute migraine medication or having an attack 48 h prior to the MRI sessions. Detailed results for demographic and clinical assessments are shown in Table 1. All participants reported acceptable stimulation intensity underneath the electrodes during the continuous taVNS, with no adverse effects reported. The interval period of the two taVNS/fMRI scan sessions was 8.79 ± 0.74 (mean ± SE) days.
vlPAG-based functional connectivity analysis results
Compared with pre-taVNS resting state, 1 Hz continuous taVNS (during) produced significant functional connectivity increases between the vlPAG and the bilateral middle cingulate cortex (MCC), the right precuneus, the left middle frontal gyrus (MFG), and the left cuneus (Table 2, Fig. 2A). There was no significant finding detected when we applied the same analysis on the 20 Hz taVNS data set.
In addition, we compared the vlPAG-based connectivity difference in 1 Hz vs. 20 Hz taVNS [(during 1 Hz taVNS minus 1 Hz pre-taVNS resting-state) vs. (during 20 Hz taVNS minus 20 Hz pre-taVNS resting-state)], and found that compared to 20 Hz, 1 Hz taVNS produced greater vlPAG-based connectivity increases with the MCC, the right precuneus/posterior cingulate cortex (PCC), the left insula (k = 18), and the anterior cingulate cortex (ACC) (k = 41) (Table 3, Fig. 2C). No significant decrease in vlPAG-based functional connectivity was detected.
With the threshold we set, no significant result has been found in the comparison of vlPAG-based pre- and post-taVNS resting-state functional connectivity differences between the 1 Hz and 20 Hz taVNS.
We found that 1 Hz taVNS increased vlPAG resting state functional connectivity with the MCC, precuneus, MFG and cuneus compared with the pre-taVNS resting state. To explore the potential clinical meaning of these functional connectivity increases, we performed correlation analyses between the vlPAG resting-state functional connectivity with these regions during the pre-taVNS (1 Hz) and the clinical measures (migraine attacks in the past 4 weeks and VAS). Results showed a significant negative correlation between the number of preceding migraine attacks and the vlPAG-MCC functional connectivity in the pre-taVNS resting-state preceding the 1 Hz taVNS (r = − 0.52, p = 0.01, significant after Bonferroni correction p < 0.05/4 = 0.0125 because four significant clusters were identified, please see Table 2 and Fig. 2B for details). No other significant vlPAG-based functional connectivity finding was detected.
We also performed the above analysis using bilateral seeds from the fourth ventricle. No result was found at the threshold we set in functional connectivity analysis.
In the present study, we compared the vlPAG connectivity changes evoked by 1 Hz and 20 Hz taVNS in migraine patients. Results showed that compared to pre-taVNS resting-state, continuous 1 Hz taVNS produced increased connectivity in the MCC, MFG, precuneus and cuneus. Compared to 20 Hz, 1 Hz taVNS produced greater connectivity increases in the MCC, ACC, precuneus and left insula. There is a significant negative association between migraine attacks in the past 4 weeks and the vlPAG-MCC connectivity during resting-state. Our findings suggest that taVNS with different stimulation frequencies may produce different modulation effects on the descending pain modulation system.
As a non-invasive and safe peripheral neuromodulation method, taVNS has been applied in a wide range of disorders such as depression, epilepsy, tinnitus, migraine, as well as cognitive and behavioral disorders [12, 21, 76,77,78]. Nevertheless, one challenge for the development of taVNS is to elucidate the modulation effect of taVNS with different parameters so that we can optimize its effects for different disorders.
As a key parameter of taNVS, frequency is a continuous measurement. Thus, it is not possible to test/compare the effects of different frequencies in one study. As a start of this line of work, we have chosen 1 Hz as a representative of low frequency. The 20 Hz frequency has been used to treat depression, and previous studies have found that 20 Hz taVNS can significantly modulate the multiple brain networks [27, 60,61,62], particularly the functional connectivity of the amygdala , default mode network , hypothalamus , and ventral striatum , all of which are associated with pathophysiology of migraine [80, 81]. Further, investigators found that 20 Hz taVNS in healthy subjects could decrease heart rate , and 20 Hz is also close to the higher frequency used in a previous study in which the authors have compared the treatment effect of 1 Hz and 25 Hz in . Thus, we have chosen 20 Hz to represent a moderate frequency in this study.
We found that continuous 1 Hz taVNS can significantly increase vlPAG-MCC connectivity. In addition, the vlPAG-MCC connectivity during resting-state before 1 Hz taVNS was negatively associated with participants’ migraine attacks. Literature suggests that the MCC is involved in the affective, cognitive, attention, and orienting aspects of pain [82,83,84]. A previous study found that migraine is associated with decreased grey matter at the MCC  and increased activation during experimental heat pain (compared to healthy controls). The pain-induced MCC activation is associated with migraine attacks in migraineurs . Interestingly, we found that the vlPAG-MCC connectivity increased during 1 Hz taVNS, but not during 20 Hz taVNS, which may provide a neural mechanistic support to a previous clinical trial , in which researchers investigated the therapeutic effects of daily 1 Hz and 25 Hz taVNS on chronic migraineurs over 3 months, and demonstrated that 1 Hz taVNS was more prominent in migraine alleviation.
Furthermore, we observed that continuous 1 Hz taVNS can produce vlPAG-rACC connectivity increases compared to 20 Hz. In addition, we also detected an increase in vlPAG-rACC connectivity (compared to resting-state) at a less conservative threshold (p = 0.01, cluster size = 14). The rACC is a key region of the DPMS [35, 36], and contains numerous opioid receptors . Previous studies have suggested that the rACC plays an important role in the pathophysiology of migraine [88, 89]. Findings from the current study are consistent with our prior study, in which migraine patients are associated with reduced connectivity of the PAG-rACC, compared to healthy subjects, and effective acupuncture treatment can normalize the decreased connectivity in PAG-rACC correspondingly . Further, the study demonstrates that a DPMS abnormality might be an underlying pathological mechanism of migraine, and such an abnormality can be normalized by effective treatment.
In addition, we found that 1 Hz taVNS can increase vlPAG connectivity with the precuneus and cuneus. The precuneus is a key region in the default mode network. Studies suggest that the default mode network (DMN) is a pivotal network affected by migraine [90,91,92]. We found that migraineurs showed decreased functional connectivity between the right frontoparietal network and precuneus compared with healthy controls, and the connectivity significantly increased after effective treatment . In a more recent study, we found abnormal posterior thalamus (pulvinar nucleus) dynamic network functional connectivity with the precuneus, and the changes were significantly correlated with the headache frequency of migraine .
The cuneus is a key region of the visual network. In a recent longitudinal study on grey matter volume of migraineurs, researchers found that migraineurs developed a decreased grey matter volume of visual regions, including the cuneus. The decreased volume was associated with the level of migraine severity, in terms of disease duration, pain intensity, and attack frequency . We found migraine is associated with altered posterior thalamus dynamic network functional connectivity with the visual cortex , and the abnormal functional connectivity within the visual, default mode, sensorimotor, and frontal-parietal networks, which could discriminate migraineurs from healthy controls, with 93% sensitivity and 89% specificity . More recently, we found that 4-week taVNS at 1 Hz can decrease the connectivity between the occipital cortex-related thalamus subregion and the postcentral gyrus/precuneus . Taken together, these studies demonstrate the important role of the precuneus and cuneus in the pathology of migraine. Our study further suggests that 1 Hz taVNS may modulate the connectivity between the descending pain modulation system, the default mode network, and the visual network.
Nevertheless, the question whether the effective frequency of taVNS that influences migraine is different from other diseases remains open [19, 20, 97]. Further studies are needed to determine the optimal frequency of taVNS for different diseases. Additionally, as a brain imaging study, the aim of this study was to investigate and compare if 1 Hz and 20 Hz taVNS can modulate the vlPAG functional connectivity in a migraine population rather than assessing the efficacy/clinical effects of 1 Hz and 20 Hz taVNS. In addition, we used seeds in the ventricle as a control ROI, and the lack of significant results further validated our findings.
Potential limitations of this work include a relatively small sample size of migraine participants with low-frequency migraine attacks. Future studies are needed to investigate if the findings can be replicated in migraineurs with high attack frequencies in a larger sample size. Also, there are only three male participants (of 24 in total) included in this study. This ratio is partly consistent with epidemiology studies showing the prevalence rate of female migraineurs is much higher than male migraineurs . Nevertheless, we have applied a cross-over design, which should have controlled the potential gender effects in this study. This study is not designed to answer the question of gender differences. A future study is needed to elucidate if male and female migraineurs are associated with same taVNS response.
Furthermore, our MRI scans were applied when participants were migraine-free, so we could not assess the acute effects of taVNS on headache intensity. Also, clinical trials on migraine usually assess the clinical improvement (migraine attack time or pain intensity) in the past month, thus, we could not investigate/compare the clinical improvement produced by single 1 Hz/20 Hz taVNS treatment, as well as the association between functional connectivity changes (evoked by 1 Hz and 20 Hz taVNS) and clinical improvement. Moreover, although still under investigation, some studies suggest that different stimulation frequencies of taVNS may induce different changes in heart rate, which can be considered as a confounding factor of functional connectivity [19, 99]. Nevertheless, the heart rate changes evoked by taVNS are relatively small, and studies also show no significant change on blood pressure values after taVNS . Future study should consider measuring this confounding factor and adjust for it during data analysis.
In summary, we found continuous 1 Hz taVNS can significantly modulate functional connectivity between the vlPAG and key regions of the DPMS in patients with migraine. Our findings demonstrate the important role of stimulation parameters (particularly the frequency) in taVNS treatment of different disorders.
Availability of data and materials
Data supporting the findings of this study are available from the corresponding author, upon reasonable request.
Transcutaneous auricular vagus nerve stimulation
Descending pain modulation system
Default mode network
Magnetic resonance imaging
Functional magnetic resonance imaging
Ventrolateral periaqueductal gray
Anterior cingulate cortex
Rostral anterior cingulate cortex
Middle cingulate cortex
Posterior cingulate cortex
Middle frontal gyrus
Medial prefrontal cortex
Regions of interest
Full width at half maximum
Visual analog scale
Yuan H, Silberstein SD. Vagus nerve and vagus nerve stimulation, a comprehensive review: part III. Headache. 2016;56:479–90. https://doi.org/10.1111/head.12649.
Pereira MR, Leite PEC. The involvement of parasympathetic and sympathetic nerve in the inflammatory reflex. J Cell Physiol. 2016;231:1862–9. https://doi.org/10.1002/jcp.25307.
Bellavance MA, Rivest S. The HPA—immune axis and the immunomodulatory actions of glucocorticoids in the brain. Front Immunol. 2014;5:136.
Tracey KJ. Reflex control of immunity. Nat Rev Immunol. 2009;9:418–28.
Willemze RA, Luyer MD, Buurman WA, de Jonge WJ. Neural reflex pathways in intestinal inflammation: hypotheses to viable therapy. Nat Rev Gastroenterol Hepatol. 2015;12:353–62. https://doi.org/10.1038/nrgastro.2015.56.
Borovikova LV, Ivanova S, Zhang M, et al. Vagus nerve stimulation attenuates the systemic inflammatory response to endotoxin. Nature. 2000;405:458–62. https://doi.org/10.1038/35013070.
Browning KN, Verheijden S, Boeckxstaens GE. The vagus nerve in appetite regulation, mood, and intestinal inflammation. Gastroenterology. 2017;152:730–44. https://doi.org/10.1053/j.gastro.2016.10.046.
Henry TR. Therapeutic mechanisms of vagus nerve stimulation. Neurology. 2002. https://doi.org/10.1212/wnl.59.6_suppl_4.s3.
Peuker ET, Filler TJ. The nerve supply of the human auricle. Clin Anat. 2002;15:35–7. https://doi.org/10.1002/ca.1089.
Shiozawa P, da Silva ME, de Carvalho TC, et al. Transcutaneous vagus and trigeminal nerve stimulation for neuropsychiatric disorders: a systematic review. Arq Neuropsiquiatr. 2014;72:542–7. https://doi.org/10.1590/0004-282x20140061.
Mercante B, Deriu F, Rangon C-M. Auricular neuromodulation: the emerging concept beyond the stimulation of vagus and trigeminal nerves. Medicines. 2018;5:10. https://doi.org/10.3390/medicines5010010.
Kong J, Fang J, Park J, et al. Treating depression with transcutaneous auricular vagus nerve stimulation: state of the art and future perspectives. Front Psychiatry. 2018;9:20. https://doi.org/10.3389/fpsyt.2018.00020.
Ben-Menachem E, Revesz D, Simon BJ, Silberstein S. Surgically implanted and non-invasive vagus nerve stimulation: a review of efficacy, safety and tolerability. Eur J Neurol. 2015;22:1260–8. https://doi.org/10.1111/ene.12629.
Ellrich J. Transcutaneous auricular vagus nerve stimulation. J Clin Neurophysiol. 2019;36:437–42.
Butt MF, Albusoda A, Farmer AD, Aziz Q. The anatomical basis for transcutaneous auricular vagus nerve stimulation. J Anat. 2020;236:588–611.
Busch V, Zeman F, Heckel A, et al. The effect of transcutaneous vagus nerve stimulation on pain perception—an experimental study. Brain Stimul. 2013;6:202–9. https://doi.org/10.1016/j.brs.2012.04.006.
Mercante B, Ginatempo F, Manca A, et al. Anatomo-physiologic basis for auricular stimulation. Med Acupunct. 2018;30:141–50.
Kaniusas E, Kampusch S, Tittgemeyer M, et al. Current directions in the auricular vagus nerve stimulation I—a physiological perspective. Front Neurosci. 2019;13:854.
Badran BW, Mithoefer OJ, Summer CE, et al. Short trains of transcutaneous auricular vagus nerve stimulation (taVNS) have parameter-specific effects on heart rate. Brain Stimul. 2018;11:699–708. https://doi.org/10.1016/j.brs.2018.04.004.
Bauer S, Baier H, Baumgartner C, et al. Transcutaneous vagus nerve stimulation (tVNS) for treatment of drug-resistant epilepsy: a randomized, double-blind clinical trial (cMPsE02). Brain Stimul. 2016;9:356–63. https://doi.org/10.1016/j.brs.2015.11.003.
Straube A, Ellrich J, Eren O, et al. Treatment of chronic migraine with transcutaneous stimulation of the auricular branch of the vagal nerve (auricular t-VNS): a randomized, monocentric clinical trial. J Headache Pain. 2015. https://doi.org/10.1186/s10194-015-0543-3.
Frangos E, Ellrich J, Komisaruk BR. Non-invasive access to the vagus nerve central projections via electrical stimulation of the external ear: FMRI evidence in humans. Brain Stimul. 2015;8:624–36. https://doi.org/10.1016/j.brs.2014.11.018.
Sclocco R, Garcia RG, Kettner NW, et al. The influence of respiration on brainstem and cardiovagal response to auricular vagus nerve stimulation: a multimodal ultrahigh-field (7T) fMRI study. Brain Stimul. 2019;12:911–21. https://doi.org/10.1016/j.brs.2019.02.003.
Zhao B, Bi Y, Li L, et al. The instant spontaneous neuronal activity modulation of transcutaneous auricular vagus nerve stimulation on patients with primary insomnia. Front Neurosci. 2020;14:205. https://doi.org/10.3389/fnins.2020.00205.
Badran BW, Dowdle LT, Mithoefer OJ, et al. Neurophysiologic effects of transcutaneous auricular vagus nerve stimulation (taVNS) via electrical stimulation of the tragus: a concurrent taVNS/fMRI study and review. Brain Stimul. 2018;11:492–500. https://doi.org/10.1016/j.brs.2017.12.009.
Kraus T, Hösl K, Kiess O, et al. BOLD fMRI deactivation of limbic and temporal brain structures and mood enhancing effect by transcutaneous vagus nerve stimulation. J Neural Transm. 2007;114:1485–93. https://doi.org/10.1007/s00702-007-0755-z.
Fang J, Egorova N, Rong P, et al. Early cortical biomarkers of longitudinal transcutaneous vagus nerve stimulation treatment success in depression. NeuroImage Clin. 2017;14:105–11. https://doi.org/10.1016/j.nicl.2016.12.016.
Kraus T, Kiess O, Hösl K, et al. CNS BOLD fMRI effects of sham-controlled transcutaneous electrical nerve stimulation in the left outer auditory canal—a pilot study. Brain Stimul. 2013;6:798–804. https://doi.org/10.1016/j.brs.2013.01.011.
Sclocco R, Garcia RG, Kettner NW, et al. Stimulus frequency modulates brainstem response to respiratory-gated transcutaneous auricular vagus nerve stimulation. Brain Stimul. 2020;13:970–8. https://doi.org/10.1016/j.brs.2020.03.011.
Wang Z, Fang J, Liu J, et al. Frequency-dependent functional connectivity of the nucleus accumbens during continuous transcutaneous vagus nerve stimulation in major depressive disorder. J Psychiatr Res. 2018;102:123–31. https://doi.org/10.1016/j.jpsychires.2017.12.018.
Tu Y, Fang J, Cao J, et al. A distinct biomarker of continuous transcutaneous vagus nerve stimulation treatment in major depressive disorder. Brain Stimul. 2018;11:501–8. https://doi.org/10.1016/j.brs.2018.01.006.
Fields H. State-dependent opioid control of pain. Nat Rev Neurosci. 2004;5:565–75. https://doi.org/10.1038/nrn1431.
Borsook D, Burstein R. The enigma of the dorsolateral pons as a migraine generator. Cephalalgia. 2012;32:803–12.
Ossipov MH, Morimura K, Porreca F. Descending pain modulation and chronification of pain. Curr Opin Support Palliat Care. 2014;8:143–51.
Kong J, Tu P, Zyloney C, Su T. Intrinsic functional connectivity of the periaqueductal gray, a resting fMRI study. Behav Brain Res. 2010;211:215–9. https://doi.org/10.1016/j.bbr.2010.03.042.
Yu R, Gollub RL, Spaeth R, et al. Disrupted functional connectivity of the periaqueductal gray in chronic low back pain. NeuroImage Clin. 2014;6:100–8. https://doi.org/10.1016/j.nicl.2014.08.019.
Aurora SK, Wilkinson F. The brain is hyperexcitable in migraine. Cephalalgia. 2007;27:1442–53.
Aurora SK, Barrodale PM, Tipton RL, Khodavirdi A. Brainstem dysfunction in chronic migraine as evidenced by neurophysiological and positron emission tomography studies. Headache. 2007;47:996–1003. https://doi.org/10.1111/j.1526-4610.2007.00853.x (discussion 1004–7).
Welch KMA, Nagesh V, Aurora SK, Gelman N. Periaqueductal gray matter dysfunction in migraine: cause or the burden of illness? Headache. 2001;41:629–37. https://doi.org/10.1046/j.1526-4610.2001.041007629.x.
Akerman S, Holland PR, Goadsby PJ. Diencephalic and brainstem mechanisms in migraine. Nat Rev Neurosci. 2011;12:570–84.
Li Z, Liu M, Lan L, et al. Altered periaqueductal gray resting state functional connectivity in migraine and the modulation effect of treatment. Sci Rep. 2016;6:20298. https://doi.org/10.1038/srep20298.
Mainero C, Boshyan J, Hadjikhani N. Altered functional magnetic resonance imaging resting-state connectivity in periaqueductal gray networks in migraine. Ann Neurol. 2011;70:838–45. https://doi.org/10.1002/ana.22537.
Chen Z, Chen X, Liu M, et al. Disrupted functional connectivity of periaqueductal gray subregions in episodic migraine. J Headache Pain. 2017. https://doi.org/10.1186/s10194-017-0747-9.
Wei SY, Chen LF, Lin MW, et al. The OPRM1 A118G polymorphism modulates the descending pain modulatory system for individual pain experience in young women with primary dysmenorrhea. Sci Rep. 2017. https://doi.org/10.1038/srep39906.
Wei SY, Chao HT, Tu CH, et al. The BDNF Val66Met polymorphism is associated with the functional connectivity dynamics of pain modulatory systems in primary dysmenorrhea. Sci Rep. 2016. https://doi.org/10.1038/srep23639.
Wei SY, Chao HT, Tu CH, et al. Changes in functional connectivity of pain modulatory systems in women with primary dysmenorrhea. Pain. 2016;157:92–102. https://doi.org/10.1097/j.pain.0000000000000340.
Li H, Li X, Feng Y, et al. Deficits in ascending and descending pain modulation pathways in patients with postherpetic neuralgia. Neuroimage. 2020;221: 117186. https://doi.org/10.1016/j.neuroimage.2020.117186.
Cifre I, Sitges C, Fraiman D, et al. Disrupted functional connectivity of the pain network in fibromyalgia. Psychosom Med. 2012;74:55–62. https://doi.org/10.1097/PSY.0b013e3182408f04.
Michels L, Christidi F, Steiger VR, et al. Pain modulation is affected differently in medication-overuse headache and chronic myofascial pain—a multimodal MRI study. Cephalalgia. 2017;37:764–79. https://doi.org/10.1177/0333102416652625.
Fenske SJ, Bierer D, Chelimsky G, et al. Sensitivity of functional connectivity to periaqueductal gray localization, with implications for identifying disease-related changes in chronic visceral pain: a MAPP Research Network neuroimaging study. NeuroImage Clin. 2020;28:102443. https://doi.org/10.1016/j.nicl.2020.102443.
Yu CX, Li B, Xu YK, et al. Altered functional connectivity of the periaqueductal gray in chronic neck and shoulder pain. NeuroReport. 2017;28:720–5. https://doi.org/10.1097/WNR.0000000000000819.
Yu S, Ortiz A, Gollub RL, et al. Acupuncture treatment modulates the connectivity of key regions of the descending pain modulation and reward systems in patients with chronic low back pain. J Clin Med. 2020;9:1719. https://doi.org/10.3390/jcm9061719.
Egorova N, Gollub RL, Kong J. Repeated verum but not placebo acupuncture normalizes connectivity in brain regions dysregulated in chronic pain. NeuroImage Clin. 2015;9:430–5. https://doi.org/10.1016/j.nicl.2015.09.012.
Zyloney CE, Jensen K, Polich G, et al. Imaging the functional connectivity of the periaqueductal gray during genuine and sham electroacupuncture treatment. Mol Pain. 2010. https://doi.org/10.1186/1744-8069-6-80.
Zhang Y, Liu J, Li H, et al. Transcutaneous auricular vagus nerve stimulation at 1 Hz modulates locus coeruleus activity and resting state functional connectivity in patients with migraine: an fMRI study. NeuroImage Clin. 2019;24: 101971. https://doi.org/10.1016/j.nicl.2019.101971.
Li Z, Zeng F, Yin T, et al. Acupuncture modulates the abnormal brainstem activity in migraine without aura patients. NeuroImage Clin. 2017;15:367–75. https://doi.org/10.1016/j.nicl.2017.05.013.
Zung WWK. A self-rating depression scale. Arch Gen Psychiatry. 1965;12:63–70. https://doi.org/10.1001/archpsyc.1965.01720310065008.
Zung WWK. A rating instrument for anxiety disorders. Psychosomatics. 1971;12:371–9. https://doi.org/10.1016/S0033-3182(71)71479-0.
Badran BW, Brown JC, Dowdle LT, et al. Tragus or cymba conchae? Investigating the anatomical foundation of transcutaneous auricular vagus nerve stimulation (taVNS). Brain Stimul. 2018;11:947–8.
Fang J, Rong P, Hong Y, et al. Transcutaneous vagus nerve stimulation modulates default mode network in major depressive disorder. Biol Psychiatry. 2016;79:266–73. https://doi.org/10.1016/j.biopsych.2015.03.025.
Rong P, Liu J, Wang L, et al. Effect of transcutaneous auricular vagus nerve stimulation on major depressive disorder: a nonrandomized controlled pilot study. J Affect Disord. 2016;195:172–9. https://doi.org/10.1016/j.jad.2016.02.031.
Liu J, Fang J, Wang Z, et al. Transcutaneous vagus nerve stimulation modulates amygdala functional connectivity in patients with depression. J Affect Disord. 2016;205:319–26. https://doi.org/10.1016/j.jad.2016.08.003.
Bagley CL, Rendas-Baum R, Maglinte GA, et al. Validating migraine-specific quality of life questionnaire v2.1 in episodic and chronic migraine. Headache. 2012;52:409–21. https://doi.org/10.1111/j.1526-4610.2011.01997.x.
Whitfield-Gabrieli S, Nieto-Castanon A. Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect. 2012;2:125–41. https://doi.org/10.1089/brain.2012.0073.
Maldjian JA, Laurienti PJ, Kraft RA, Burdette JH. An automated method for neuroanatomic and cytoarchitectonic atlas-based interrogation of fMRI data sets. Neuroimage. 2003;19:1233–9.
Kwon M, Altin M, Duenas H, Alev L. The role of descending inhibitory pathways on chronic pain modulation and clinical implications. Pain Pract. 2014;14:656–67. https://doi.org/10.1111/papr.12145.
Brennan KC, Pietrobon D. A systems neuroscience approach to migraine. Neuron. 2018;97:1004–21.
Amin FM, Hougaard A, Magon S, et al. Altered thalamic connectivity during spontaneous attacks of migraine without aura: a resting-state fMRI study. Cephalalgia. 2018;38:1237–44. https://doi.org/10.1177/0333102417729113.
Chong CD, Dumkrieger GM, Schwedt TJ. Structural co-variance patterns in migraine: a cross-sectional study exploring the role of the hippocampus. Headache. 2017;57:1522–31. https://doi.org/10.1111/head.13193.
Coppola G, Petolicchio B, Di Renzo A, et al. Cerebral gray matter volume in patients with chronic migraine: correlations with clinical features. J Headache Pain. 2017;18:115. https://doi.org/10.1186/s10194-017-0825-z.
Liu H-Y, Chou K-H, Lee P-L, et al. Hippocampus and amygdala volume in relation to migraine frequency and prognosis. Cephalalgia. 2017;37:1329–36. https://doi.org/10.1177/0333102416678624.
Qin ZX, Su JJ, He XW, et al. Altered resting-state functional connectivity between subregions in the thalamus and cortex in migraine without aura. Eur J Neurol. 2020;27:2233–41. https://doi.org/10.1111/ene.14411.
Eippert F, Bingel U, Schoell ED, et al. Activation of the opioidergic descending pain control system underlies placebo analgesia. Neuron. 2009;63:533–43. https://doi.org/10.1016/j.neuron.2009.07.014.
Kong J, Jensen K, Loiotile R, et al. Functional connectivity of the frontoparietal network predicts cognitive modulation of pain. Pain. 2013;154:459–67. https://doi.org/10.1016/j.pain.2012.12.004.
Gollub RL, Kirsch I, Maleki N, et al. A Functional neuroimaging study of expectancy effects on pain response in patients with knee osteoarthritis. J Pain. 2018;19:515–27. https://doi.org/10.1016/j.jpain.2017.12.260.
Hyvarinen P, Yrttiaho S, Lehtimaki J, et al. Transcutaneous vagus nerve stimulation modulates tinnitus-related beta-and gamma-band activity. Ear Hear. 2015;36:e76–85. https://doi.org/10.1097/AUD.0000000000000123.
Jacobs HIL, Riphagen JM, Razat CM, et al. Transcutaneous vagus nerve stimulation boosts associative memory in older individuals. Neurobiol Aging. 2015;36:1860–7. https://doi.org/10.1016/j.neurobiolaging.2015.02.023.
Beste C, Steenbergen L, Sellaro R, et al. Effects of concomitant stimulation of the GABAergic and norepinephrine system on inhibitory control—a study using transcutaneous vagus nerve stimulation. Brain Stimul. 2016;9:811–8. https://doi.org/10.1016/j.brs.2016.07.004.
Tu Y, Fu Z, Zeng F, et al. Abnormal thalamocortical network dynamics in migraine. Neurology. 2019;92:e2706–16. https://doi.org/10.1212/WNL.0000000000007607.
Goadsby PJ, Holland PR, Martins-Oliveira M, et al. Pathophysiology of migraine: a disorder of sensory processing. Physiol Rev. 2017;97:553. https://doi.org/10.1152/PHYSREV.00034.2015.
Schwedt TJ, Chiang CC, Chong CD, Dodick DW. Functional MRI of migraine. Lancet Neurol. 2015;14:81–91. https://doi.org/10.1016/S1474-4422(14)70193-0.
Frot M, Mauguière F, Magnin M, Garcia-Larrea L. Parallel processing of nociceptive A-δ inputs in SII and midcingulate cortex in humans. J Neurosci. 2008;28:944–52. https://doi.org/10.1523/JNEUROSCI.2934-07.2008.
Shackman AJ, Salomons TV, Slagter HA, et al. The integration of negative affect, pain and cognitive control in the cingulate cortex. Nat Rev Neurosci. 2011;12:154–67.
Vogt BA. Pain and emotion interactions in subregions of the cingulate gyrus. Nat Rev Neurosci. 2005;6:533–44. https://doi.org/10.1038/nrn1704.
Hougaard A, Amin FM, Arngrim N, et al. Sensory migraine aura is not associated with structural grey matter abnormalities. NeuroImage Clin. 2016;11:322–7. https://doi.org/10.1016/j.nicl.2016.02.007.
Schwedt TJ, Chong CD, Chiang CC, et al. Enhanced pain-induced activity of pain-processing regions in a case–control study of episodic migraine. Cephalalgia. 2014;34:947–58. https://doi.org/10.1177/0333102414526069.
Baumgärtner U, Buchholz HG, Bellosevich A, et al. High opiate receptor binding potential in the human lateral pain system. Neuroimage. 2006;30:692–9. https://doi.org/10.1016/j.neuroimage.2005.10.033.
Magis D, D’Ostilio K, Thibaut A, et al. Cerebral metabolism before and after external trigeminal nerve stimulation in episodic migraine. Cephalalgia. 2017;37:881–91. https://doi.org/10.1177/0333102416656118.
Matharu MS, Bartsch T, Ward N, et al. Central neuromodulation in chronic migraine patients with suboccipital stimulators: a PET study. Brain. 2004;127:220–30. https://doi.org/10.1093/brain/awh022.
Tessitore A, Russo A, Giordano A, et al. Disrupted default mode network connectivity in migraine without aura. J Headache Pain. 2013;14:89. https://doi.org/10.1186/1129-2377-14-89.
Coppola G, Di Renzo A, Tinelli E, et al. Resting state connectivity between default mode network and insula encodes acute migraine headache. Cephalalgia. 2018;38:846–54. https://doi.org/10.1177/0333102417715230.
Zhang J, Su J, Wang M, et al. Increased default mode network connectivity and increased regional homogeneity in migraineurs without aura. J Headache Pain. 2016;17:98. https://doi.org/10.1186/s10194-016-0692-z.
Li Z, Lan L, Zeng F, et al. The altered right frontoparietal network functional connectivity in migraine and the modulation effect of treatment. Cephalalgia. 2017;37:161–76. https://doi.org/10.1177/0333102416641665.
Messina R, Rocca MA, Colombo B, et al. Gray matter volume modifications in migraine: a cross-sectional and longitudinal study. Neurology. 2018;91:e280–92. https://doi.org/10.1212/WNL.0000000000005819.
Tu Y, Zeng F, Lan L, et al. An fMRI-based neural marker for migraine without aura. Neurology. 2020;94:e741–51. https://doi.org/10.1212/wnl.0000000000008962.
Zhang Y, Huang Y, Li H, et al. Transcutaneous auricular vagus nerve stimulation (taVNS) for migraine: an fMRI study. Reg Anesth Pain Med. 2021;46:145–50. https://doi.org/10.1136/rapm-2020-102088.
Milby AH, Halpern CH, Baltuch GH. Vagus nerve stimulation for epilepsy and depression. Neurotherapeutics. 2008;5:75–85. https://doi.org/10.1016/j.nurt.2007.10.071.
Vetvik KG, MacGregor EA. Sex differences in the epidemiology, clinical features, and pathophysiology of migraine. Lancet Neurol. 2017;16:76.
Staley R, Garcia RG, Stowell J, et al. Modulatory effects of respiratory-gated auricular vagal nerve stimulation on cardiovagal activity in hypertension. In: Proceedings of the annual international conference of the IEEE Engineering in Medicine and Biology Society, EMBS. Institute of Electrical and Electronics Engineers Inc.; 2020. p. 2581–4.
Hein E, Nowak M, Kiess O, et al. Auricular transcutaneous electrical nerve stimulation in depressed patients: a randomized controlled pilot study. J Neural Transm. 2013;120:821–7. https://doi.org/10.1007/s00702-012-0908-6.
This work was supported by Medical Scientific Research Foundation of Guangdong Province of China (A2017234) and the Administration of Traditional Chinese Medicine of Guangdong Province of China (20182047) to Yue Zhang.
Ethics approval and consent to participate
Study protocol was approved by the Ethics Committee of the Second Affiliated Hospital of Guangzhou University of Chinese Medicine (Z2016-079-01), and a written informed consent was obtained from all participants in accordance with the Declaration of Helsinki prior to study enrollment.
Consent for publication
Jian Kong has a disclosure to report (holding equity in a startup company, MNT, and a patent to develop new neuromodulation devices) but declares no conflict of interest. All other authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cao, J., Zhang, Y., Li, H. et al. Different modulation effects of 1 Hz and 20 Hz transcutaneous auricular vagus nerve stimulation on the functional connectivity of the periaqueductal gray in patients with migraine. J Transl Med 19, 354 (2021). https://doi.org/10.1186/s12967-021-03024-9
- Transcutaneous auricular vagus nerve stimulation
- Functional connectivity
- Periaqueductal gray
- Descending pain modulation network