Next Article in Journal
Lenition in L2 Spanish: The Impact of Study Abroad on Phonological Acquisition
Next Article in Special Issue
Repetitive Transcranial Magnetic Stimulation for the Treatment of Depression in a Real-World Setting: Findings from a Cohort Study
Previous Article in Journal
Neuronal Mechanisms of Reading Informational Texts in People with Different Levels of Mental Resilience
Previous Article in Special Issue
Intermittent Theta Burst Stimulation Combined with Cognitive Training to Improve Negative Symptoms and Cognitive Impairment in Schizophrenia: A Pilot Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Brain Sciences
Volume 14
Issue 9
10.3390/brainsci14090945
brainsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcutaneous Auricular Vagus Nerve Stimulation Modulating the Brain Topological Architecture of Functional Network in Major Depressive Disorder: An fMRI Study

by
Zhi-Peng Guo
1,†,
Dan Liao
2,†,
Lei Chen
1,
Cong Wang
3,
Miao Qu
4,
Xue-Yu Lv
5,
Ji-Liang Fang
5 and
Chun-Hong Liu
1,*
1
Beijing Hospital of Traditional Chinese Medicine, Capital Medical University, Beijing 100010, China
2
Department of Radiology, Guizhou Provincial People’s Hospital, Guiyang 550002, China
3
Kerfun Medical (Suzhou) Co., Ltd., Suzhou 215000, China
4
Department of Neurology, Xuanwu Hospital, Capital Medical University, Beijing 100053, China
5
Guang’anmen Hospital, China Academy of Chinese Medical Sciences, Beijing 100053, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Brain Sci. 2024, 14(9), 945; https://doi.org/10.3390/brainsci14090945
Submission received: 20 May 2024 / Revised: 20 September 2024 / Accepted: 20 September 2024 / Published: 21 September 2024
(This article belongs to the Special Issue Brain Stimulation for Psychiatric Disorders: Emerging Evidence and New Perspectives)
Browse Figures
Versions Notes

Abstract

:
Background: Transcutaneous auricular vagus nerve stimulation (taVNS) is effective in regulating mood and high-level cognition in patients with major depressive disorder (MDD). This study aimed to investigate the efficacy of taVNS treatment in patients with MDD and an altered brain topological organization of functional networks. Methods: Nineteen patients with MDD were enrolled in this study. Patients with MDD underwent 4 weeks of taVNS treatments; resting-state functional magnetic resonance imaging (rs-fMRI) data of the patients were collected before and after taVNS treatment. The graph theory method and network-based statistics (NBS) analysis were used to detect abnormal topological organizations of functional networks in patients with MDD before and after taVNS treatment. A correlation analysis was performed to characterize the relationship between altered network properties and neuropsychological scores. Results: After 4 weeks of taVNS treatment, patients with MDD had increased global efficiency and decreased characteristic path length (Lp). Additionally, patients with MDD exhibited increased nodal efficiency (NE) and degree centrality (DC) in the left angular gyrus. NBS results showed that patients with MDD exhibited reduced connectivity between default mode network (DMN)–frontoparietal network (FPN), DMN–cingulo-opercular network (CON), and FPN–CON. Furthermore, changes in Lp and DC were correlated with changes in Hamilton depression scores. Conclusions: These findings demonstrated that taVNS may be an effective method for reducing the severity of depressive symptoms in patients with MDD, mainly through modulating the brain’s topological organization. Our study may offer insights into the underlying neural mechanism of taVNS treatment in patients with MDD.

1. Introduction

Major depressive disorder (MDD) is a common psychiatric disorder that affects over 185 million people globally and is expected to be the leading cause of the global burden of mental illness by 2030 [1]. However, current treatments for MDD are unsatisfactory. The most common treatment of MDD is pharmacotherapy using antidepressants, and only approximately one-third of patients with MDD are responsive to drug therapy. Moreover, most patients with MDD are at risk of drug dependence after long-term use of antidepressants [2,3]. Other non-drug therapies, such as transcranial magnetic stimulation (TMS) [4], transcranial direct current stimulation (tDCS) [5], deep brain stimulation (DBS) [6], and electroconvulsive therapy (ECT) [7], have been increasingly used for the treatment of MDD. However, these treatment options are difficult to implement widely in clinical practice, mainly due to their poorly understood side effects, complex procedures, and high costs [8].
The vagus nerve is the longest cranial nerve in the human body and has an important role in the regulation of several body systems. It plays a key role in maintaining the stability of the internal environment, including physiological and biochemical processes [9]. The auricular branch of the vagus nerve is its only cutaneous branch [10]. Researchers have suggested that electric acupuncture stimulating the external ear canal, which is innervated by the cutaneous branch of the vagus nerve, produces a “vagus nerve effect”, which can improve the clinical symptoms of MDD [11,12]. Transcutaneous auricular VNS (taVNS), as proposed by Ventureyra in 2000 [13], is effective in treating disorders such as epilepsy [14], migraines [3], insomnia [15], and MDD [16]. The basis of taVNS as a symptomatic treatment in various neuropsychiatric disorders may lie in its capacity to modulate diffuse neuromodulatory systems, including the noradrenergic, cholinergic, and serotonergic systems [17]. Furthermore, taVNS does not only reduce the shortcomings (e.g., invasive operation and high cost) of VNS but is also simple and easy to perform [10]. Similar to taVNS, TMS is noninvasive and has been used to improve symptoms of depression [18] by employing a pulsed magnetic field to induce an electrical current and modulate the activities of some brain regions. Previous studies have also shown that taVNS can significantly reduce the severity of depressive symptoms in patients with MDD [11,12,19]. However, although taVNS may be a promising and effective method for the management of MDD, its potential mechanism remains unknown.
Resting-state functional magnetic resonance imaging (rs-fMRI) is a classic tool for assessing the neural activity of the brain using blood oxygen level-dependent (BOLD) signals. These approaches have been widely utilized to investigate the underlying neural mechanism of insomnia [20], MDD [21], and the treatment efficiency of taVNS [3,8]. Sun et al. found that after 30 min of taVNS treatments, patients with MDD displayed improvements in symptoms of depression, mainly reflected by the decreased amplitude of low-frequency fluctuations (ALFFs) in the right precuneus and decreased functional connectivity (FC) between the default mode network (DMN) and the frontoparietal network (FPN) [8]. Yi et al. found that after 4 weeks of taVNS treatments, patients with MDD who had a history of suicidal attempts experienced improvement in symptoms of depression, reflected by a decrease in regional homogeneity (ReHo) in the left supplementary motor area and the right median cingulate cortex, and a decrease in ReHo values that were correlated with changes in Hamilton depression (HAMD) scores [3]. Therefore, taVNS could reduce the severity of depressive symptoms by regulating spontaneous brain activity.
As a complex hierarchical network, the human brain is capable of efficient information segmentation and integration [22]. In recent years, the advent of network science in the field of neuroscience has offered a novel approach for the topological characterization of the functional bases that underpin brain activities [23]. Numerous neuroimaging studies have emphasized the combination of rs-fMRI and graph theory analysis [24]. The graph theory method mathematically defines the entire brain as a graph, consisting of nodes (brain regions) and edges (the connectivity between any pair of brain regions) [25], and can be used to detect the transmission of information across the entire brain. Network-based statistics (NBS) analysis is another graph theory tool that provides a complementary approach for recognizing the connected subnetworks and the underlying feature of these networks [26].
Previous studies have indicated that the topological properties of the brain were altered in patients with MDD. Zhang et al. found that patients with MDD showed reduced global and local efficiency compared with controls [27]. Jacob et al. found that patients with MDD showed decreased levels of centrality (DC) in the right precuneus [28]. Furthermore, Dai et al. found that after 8 weeks of treatment using antidepressants, patients with MDD exhibited decreased topological properties (characteristic path length and the clustering coefficient) and decreased characteristic path length that correlated with the subscale of HAMD scores [29]. Wu et al. revealed that after treatment using ECT, patients with MDD exhibited increased DC in the bilateral angular cortex (AG), precuneus, and right superior frontal gyrus [30]. However, to the best of our knowledge, the topological reorganization in patients with MDD after taVNS treatment has not been investigated, and this may limit the understanding and interpretation of the neural mechanisms in patients with MDD after taVNS treatment.
In this study, we hypothesized that taVNS treatment is effective in reducing the severity of depressive symptoms in patients with MDD through the modulation of the topological organization of the brain. To test this hypothesis, the graph theory method and NBS analysis were performed in patients with MDD before and after treatment using taVNS. We also explored the relationships between changed network topological properties and a changed neuropsychological scale. Our study may provide new insights into the underlying brain mechanism for taVNS treatment in patients with MDD.

2. Materials and Methods

2.1. Participants

Nineteen patients with MDD were recruited from the Psychiatry Department of Guang’anmen Hospital. Due to safety and ethical concerns, we only recruited patients with mild or moderate MDD. No participant had suicidal ideation before or after treatment. Patients with MDD were diagnosed according to the American Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) by two qualified psychiatrists using a clinically structured interview. This study was approved by the Ethics Committee of the Guang’anmen Hospital, Beijing Hospital of Traditional Chinese Medicine. Written informed consent was obtained from each participant before being enrolled in the study.
The G*Power (version 3.1.9.7) was used to calculate the sample size. A previous study reported an effect size of 0.57 [11]. Based on a modest power (1—beta = 0.8) and the standard alpha level (0.05, two-tailed), 15 patients with MDD were needed. Considering the loss of visits and refusals (20%), at least 19 patients with MDD were needed.
The inclusion criteria for patients with MDD were as follows: (1) patients who were willing to receive taVNS treatment; (2) HAMD scores > 17; (3) age between 18 and 60 years; (4) right-handedness; (5) patients willing to discontinue the use of antidepressants or other psychiatric drugs two weeks before beginning the treatment; and (6) patients who had exhibited depressive symptoms for at least two weeks.
The exclusion criteria for patients with MDD were as follows: (1) patients with other psychiatric disorders (such as bipolar disorder, etc.); (2) pregnancy; (3) contraindications to MRI; (4) patients with a history of drug and/or alcohol abuse; and (5) history of severe medical diseases.

2.2. TaVNS Treatment

In this study, all patients with MDD received 4 weeks of taVNS treatment, and a taVNS device (Hwato, SDZ-IIB, Suzhou, China) was used to stimulate the auricular concha area where there is an abundant branch distribution of the vagus nerve (Figure 1).
The stimulation parameters of taVNS were as follows: (1) the intensity is 4–6 milliamperes (adjusted according to the tolerance of patients); (2) the frequency is 20 Hz, with a wave width of <1 millisecond; and (3) each treatment lasted for 30 min and was performed twice a day (morning and evening) for a minimum of five days per week.

2.3. Neuropsychological Scales

The participants were evaluated using the Hamilton Anxiety Rating Scale (HAMA), the 24-item Hamilton Depression Rating Scale, the Self-Rating Anxiety Scale (SAS), and the Self-Rating Depression Scale (SDS) before and after 4 weeks of taVNS treatments by qualified psychiatrists. The severity of depressive symptoms was assessed by the HAMD score [31]. The severity of anxiety symptoms was detected using the HAMA score [32]. The SDS and SAS scores were used to measure the depressive and anxiety symptoms [33].

2.4. MRI Data Acquisition

The functional MR data were acquired from a GE Signa 1.5 Tesla scanner (GE Healthcare, Buckinghamshire, UK) at the Guang’anmen Hospital. All participants underwent rs-fMRI scans before and after the 4-week treatment of taVNS. During the MRI scans, the participants were instructed to keep their eyes closed without falling asleep and to have relaxed thoughts while refraining from engaging in any specific thought.
The rs-fMRI data were acquired with the single-shot gradient-echo EPI sequence as follows: echo time (TE) = 30 ms, repetition time (TR) =2500 ms, field of view (FOV) = 240 × 240 mm2, fip angle (FA) = 90°, gap = 0.5 mm, matrix = 64 × 64, slice thickness = 3.0 mm, and number of slices = 41. The parameters of T1-weighted images were as follows: TR = 8.1 ms, TE = 3.7 ms, FA = 8°, number of slices = 170, FOV = 256 × 256 mm2, slice thickness = 1 mm, and matrix = 256 × 256. The rs-fMRI sequence scan took 8 min and 6 s.

2.5. Data Preprocessing

The rs-fMRI data were preprocessed by using the graph theoretical network analysis toolbox (GRETNA, http://www.nitrc.org/projects/gretna/ (accessed on 25 September 2021)) and Statistical Parametric Mapping, version 12 (SPM12, http://www.fil.ion.ucl.ac.uk/spm/ (accessed on 13 August 2021)) [34]. The steps of data preprocessing were as follows: (1) the first 10 volumes were removed because of possible equilibration of the magnetic field and movement; (2) slice timing and realignment of head motion were performed in the remaining 133 time points; participants with a head motion (<2 mm or <2°) were excluded; (3) diffeomorphic anatomical registration through the exponentiated lie algebra (DARTEL) were used to normalize the generated images into the Montreal Neurological Institute (MNI) space and then resampled to a 3 × 3 × 3 mm isotropic voxel size; (4) to reduce the influence of MRI equipment, linear detrending was performed; (5) to reduce the physiological high-frequency respiratory and cardiac noise, low-frequency band-pass filtering (0.01–0.10 Hz) was performed; (6) a Gaussian kernel with a full width at half maximum of 6 mm was used for spatial smoothing; and (7) the framewise displacement values were also calculated to minimize the possible effects of micromovements of the head.

2.6. Functional Network Construction

As the basic elements of a network, the definition of the nodes and edges was key to the construction of the network. In the current study, the brain network nodes were defined by the Dosenbach atlas, which includes 160 cortical or subcortical regions [35]. Afterwards, Pearson correlation was employed to determine the edges of the network, which was calculated based on the mean time series of each node pair. Subsequently, a 160 × 160 functional connectivity (FC) matrix was obtained for each participant and then transformed into a binary matrix through a set of thresholds.
Consistent with the methods used in previous studies, our study utilized a sparsity range (0.05–0.4) with a step of 0.01 to construct the brain network [36]. The ratio of edges in a network compared with the maximum possible number of edges was the sparsity value. The minimum sparsity was the lower bound when all the nodes were connected in a network. Finally, we performed Fisher’s z-transformation on the individual matrix.

2.7. Network Metrics Analysis

Global and nodal network properties of the brain network were calculated at a serial of the sparsity using the GRETNA toolbox. Global network properties included the following: clustering coefficient (Cp), normalized Cp (γ), characteristic path length (Lp), normalized Lp (λ), small-worldness (σ), global efficiency (Eglob), and local efficiency (Eloc). The nodal network properties included the betweenness centrality (BC), degree centrality (DC), and nodal efficiency (NE) [36].
In order to offer a summarized scalar independent of single-threshold selection, the area under the curve (AUC) of topological properties were calculated. Detailed information is presented in Table S1.

2.8. Statistical Analysis

SPSS version 25.0 software (IBM Corp., Armonk, NY, USA) was used to compare demographic and clinical data, and the normality of the distribution was detected using the Kolmogorov–Smirnov test. A paired t-test was used to detect the alternation of the neuropsychological scale in patients with MDD before and after taVNS treatment, and p < 0.05 was used as the threshold of statistical significance.
Global and nodal properties were analyzed using the GRETNA toolbox. Paired t-tests were used to detect differences in patients with MDD before and after taVNS treatment, with covariates (age, sex, and mean FD values). The false discovery rate (FDR corrected, p < 0.05) was used for multiple comparisons.

2.8.1. NBS Analysis

NBS analysis was used to further investigate the edge-based FC of the network in patients with MDD before and after taVNS treatment. This analysis was performed using the NBS method (version 1.2, https://www.nitrc.org/projects/nbs (accessed on 13 August 2021)) [37], which offers a higher statistical power [38]. The steps involved were as follows: (1) the t-statistic was used to determine a range of suprathreshold links for each network edge; (2) connected components were detected using a breadth-first search, where any two regions were linked by suprathreshold connections; (3) random permutation testing (5000 times) was conducted with the p-value controlled for a family-wise error, which was attributed to the number of links in each connected component; and (4) a corrected p-value was obtained by computing the proportion of permutations, as the null distribution for the size of the largest connected component was generated during the permutation process.

2.8.2. Correlation Analysis

We also conducted correlation analyses between altered topological properties and the change in neuropsychological scale (HAMD, HAMA, SAD, and SDS) scores. The threshold of statistical significance was p < 0.05.

3. Results

3.1. Demographic and Neuropsychological Scores

The demographic and neuropsychological scores of participants are depicted in Table 1. After 4 weeks of taVNS treatment, there were significant differences between pre- and post-treatment HAMD, HAMA, SAS, and SDS scores of patients with MDD (all p < 0.001). These changes in neuropsychological scales are shown in Table 2 and Figure S1.

3.2. Global Topological Properties

Patients with MDD (before and after the taVNS treatment) all demonstrated small-worldness (sigma > 1), with gamma > 1 and lambda approximately equal to one across the sparsity range of 0.05–0.40 (step = 0.01) (Figure 2). After 4 weeks of taVNS treatment, patients with MDD demonstrated increased Eglob and decreased Lp compared to baseline values (p < 0.001, FDR corrected). No significant difference was observed in other global properties. Detailed information is presented in Table 3 and Figure 3.

3.3. Nodal Topological Properties

After 4 weeks of taVNS treatment, the patients with MDD showed significantly increased DC values in the left angular gyrus (AG) (p < 0.001, FDR corrected). In addition, the patients with MDD showed significantly increased NE values of the left AG (p = 0.048, FDR corrected). The changes in the participants’ global network topological properties are presented in Table 3 and Figure 4.

3.4. NBS Analysis

Compared with baseline values, after 4 weeks of taVNS treatment, patients with MDD showed decreased FC in several regions, which contained 16 nodes and 10 edges (p < 0.05, FDR corrected), mainly relating to the default mode network (DMN), frontoparietal network (FPN), and cingulo-opercular network (CON). These edges connected the left thalamus and the right ventromedial prefrontal cortex (vmPFC), the right thalamus to the right superior frontal cortex, the right inferior parietal lobe and left middle insula to the left occipital cortex, the left inferior parietal lobe to the right precuneus, the right dorsolateral prefrontal cortex (dlPFC) to the right precuneus, the right superior temporal cortex to the right occipital cortex, the right inferior parietal lobe to the right dlPFC, the left middle insula to the left anterior cingulate cortex, and the left middle insula to the left parietal lobe. Specifically, after taVNS treatment, patients with MDD exhibited reduced connectivity between the DMN–FPN, DMN–CON, and FPN–CON (all p < 0.05, FDR corrected). The visualization of these alterations in FC is shown in Figure 5.

3.5. Correlation Analysis

The correlation analysis indicated that the change in Lp was associated with changes in HAMD scores in patients with MDD (r = 0.616, p = 0.005). Furthermore, changes in DC in the left AG were negatively associated with changes in HAMD scores (r = −0.703, p = 0.001) (Figure 6).

4. Discussion

By combining the graph theory method and NBS analysis, our study investigated topological metrics and FC changes in patients with MDD before and after taVNS treatment. After 4 weeks of taVNS treatment, patients with MDD had increased Eglob and decreased Lp. Additionally, patients with MDD exhibited increased DC and NE in the left AG and reduced FC between the DMN–FPN, DMN–CON, and FPN–CON. The alternations in Lp and DC were correlated with changes in HAMD scores. Our study may offer insights into the underlying neural mechanisms of taVNS treatment in patients with MDD.

4.1. TaVNS Treatment Efficiency

In this study, after 4 weeks of treatment using taVNS, patients with MDD exhibited significantly decreased HAMD scores, which indicated improvements in depressive symptoms. Our study’s results are consistent with those of previous studies that concluded that depressive symptoms were relieved in patients with MDD after treatment with taVNS [11,39]. Liu et al. found that after 4 weeks of taVNS treatments, patients with MDD showed significantly decreased HAMD scores [39]. Rong et al. also found that after 4 weeks of taVNS treatment, patients with MDD demonstrated significantly reduced HAMD scores [11]. To summarize, our findings, together with those of the previous literature, indicate that the taVNS method may be an effective tool to relieve the depressive symptoms in patients with MDD.

4.2. Global Topological Properties’ Alteration after taVNS Treatment

A brain with small-worldness is associated with a high Cp and a low Lp, supporting the balance of two fundamental brain organization principles, namely segregation and integration [40]. These properties facilitate the efficient transfer of information throughout the brain’s network with minimal wiring and energy [41]. In the current study, patients with MDD (both before and after taVNS treatment) exhibited small-worldness, indicating that the small-world model may not be easily changed after taVNS treatment in order to ensure an accurate transmission of information throughout the entire brain [42].
The Eglob and Lp represent the capacity of global information transmission in a brain network, and a high Eglob and low Lp indicate stronger information integration and faster information communication [22,29]. Daws et al. found that after psilocybin therapy, patients with MDD experienced enhanced global integration [43]. Dai et al. revealed that after 9 weeks of treatment using antidepressants, the Eglob was increased while the Lp was decreased in patients with MDD [29]. In the current study, the findings showed that patients with MDD had increased Eglob and decreased Lp after 4 weeks of taVNS treatment; these findings are consistent with those of previous studies [42,44,45].
Furthermore, changes in Lp were associated with changes in HAMD scores in patients with MDD, suggesting that the change in Lp was related to improvement in the severity of depressive symptoms in patients with MDD. Zhang et al. observed that after 8 weeks of treatment using antidepressants, changes in Lp were associated with a reduction in HAMD scores in patients with MDD [46]. This suggests that taVNS could improve the depressive symptoms in patients with MDD by normalizing the disrupted network global topological organization.

4.3. Nodal Topological Properties’ Alteration after taVNS Treatment

Within the mathematical framework of graph theory, DC is often used to quantify the functional connectivity of a node [47]. Higher DC values represent a stronger influence on other nodes and a greater capacity to communicate information in the network, making it potentially a good target for intervention at the nodal level [48]. In addition, NE is used to describe the ability of information to transmit from one node to other nodes [49]. Higher NE values means a more efficient information transfer between network nodes [50].
After 4 weeks of taVNS treatment, our study found that patients with MDD demonstrated increased DC and NE in the left AG, which was a key node of the DMN. Dysfunction of the DMN has been widely reported in patients with MDD, which is associated with the processing of self-reference, emotional appraisals, and rumination [51]. Previous studies have reported that patients with MDD had decreased nodal topological properties within the DMN [28,41,52], suggesting that the ability of information transmission was impaired in MDD. Li et al. also found that ECT can improve the depressive symptoms in patients with MDD by enhancing the NE of the DMN [22].
Wu et al. observed that after ECT intervention, patients with MDD exhibited increased DC of the DMN (precuneus and bilateral AG) [30]. Our findings were consistent with prior studies, indicating the function recovery of DMN (e.g., the AG) from the aspect of topological organization after taVNS treatment. In addition, our study found that the alternation of DC in the left AG was associated with the alternation of HAMD scores in patients with MDD after taVNS treatment. Thus, we may speculate that the taVNS may regulate the disrupted DC and NE of the DMN to achieve an antidepressant effect, which could be the potential intervention target for the treatment of patients with MDD.

4.4. Functional Connectivity Alterations after taVNS Treatment

In our study, after treatment with taVNS, patients with MDD exhibited reduced connectivity between the DMN–FPN, DMN–CON, and FPN–CON. The abnormality of these brain networks has been documented in previous research [53,54,55], which suggests that patients with MDD may have disruptions in multiple networks. Hence, our findings are consistent with those of the previous literature.
The DMN, a large-scale distributed network, is involved in the processing of self-reference and the regulation of affectivity [56]. The FPN, a flexible cognitive control center, is a part of the top–down control system and is essential in a wide range of cognitive processes [57]. A meta-analysis identified an imbalance in FC between the DMN and FPN in patients with MDD [58]. Previous studies have shown that ECT, antidepressant medications, and taVNS can relieve depressive symptoms by reducing FC between the DMN and FPN in patients with MDD [8,59,60]. In this study, patients with MDD had decreased FC between the DMN and FPN after 4 weeks of taVNS treatment, which is consistent with the reported findings of previous studies.
The CON, a key part of the top–down control system, is involved in the flexible control of goal-directed performance [61]. Recent neuroimaging research indicated that patients with MDD have stronger FC between the DMN and CON [62]. Wu et al. found that patients with MDD showed abnormal FC between the FPN and CON [63], and the FC was associated with the HAMD scores. Argyelan et al. indicated that after ECT intervention, patients with MDD exhibited reduced FC between the DMN and CON [64]. Moreover, Kang et al. observed that after 2 weeks of rTMS treatment, patients with MDD showed decreased FC between the FPN and CON [65]. Our findings were consistent with those of these studies, further highlighting the importance of these high-level cognitive networks (e.g., DMN, FPN, CON) in MDD. Therefore, it is our opinion that taVNS may improve the clinical symptoms of depression by downregulating the FC between DMN–FPN, DMN–CON, and FPN–CON in patients with MDD.

4.5. Limitations

This study has several limitations. First, the sample size was small, which may have constrained the generalizability of the results. Second, data on the number of episodes and the history of medication use were not collected. Third, taVNS is not yet widely accepted clinically; therefore, further studies are necessary to confirm the result of the present study and to verify the practicability of taVNS in clinical settings. Finally, our current study was conducted over a 4-week period, and the results may have been influenced by the relatively short washout period of 2 weeks. We intend to assess the long-term effects of taVNS therapy in a future study.

5. Conclusions

By combining the graph theory and NBS analysis, our study demonstrated that taVNS can modulate disrupted global (Eglob and Lp) and nodal properties (DC and NE), and that changes in network properties (Lp and DC) were correlated with changes in HAMD scores. Furthermore, patients with MDD exhibited reduced connectivity between the DMN–FPN, DMN–CON, and FPN–CON. These findings demonstrate that taVNS treatment may improve depressive symptoms by modulating the disrupted network topological organization in patients with MDD, which may provide new insights into the underlying brain mechanism for patients with MDD.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/brainsci14090945/s1, Figure S1: Paired t-test of neuropsychological scales in patients with MDD before and after taVNS treatment.; Table S1: Description of topological properties.

Author Contributions

Conceptualization: C.-H.L.; formal analysis: Z.-P.G. and D.L.; writing—original draft: Z.-P.G., D.L. and L.C.; validation and visualization: Z.-P.G.; data curation: J.-L.F. and X.-Y.L.; supervision: C.-H.L., J.-L.F., C.W. and M.Q.; writing—review and editing: Z.-P.G. and D.L.; funding acquisition: C.-H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Beijing Natural Science Foundation of China (7212200) and the National Natural Science Foundation of China (81871507).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the Guang’anmen Hospital, Beijing Hospital of Traditional Chinese Medicine, with approval number 2021020443.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author due to the data not being easily accessible.

Acknowledgments

We thank all individuals who participated in this study and gratefully acknowledge the invaluable contributions of all authors.

Conflicts of Interest

Author Cong Wang was employed by the company Kerfun Medical (Suzhou) Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Marx, W.; Penninx, B.; Solmi, M.; Furukawa, T.A.; Firth, J.; Carvalho, A.F.; Berk, M. Major depressive disorder. Nat. Rev. Dis. Primers 2023, 9, 44. [Google Scholar] [CrossRef] [PubMed]
  2. Herrman, H.; Patel, V.; Kieling, C.; Berk, M.; Buchweitz, C.; Cuijpers, P.; Furukawa, T.A.; Kessler, R.C.; Kohrt, B.A.; Maj, M.; et al. Time for united action on depression: A Lancet-World Psychiatric Association Commission. Lancet 2022, 399, 957–1022. [Google Scholar] [CrossRef] [PubMed]
  3. Yi, S.; Wang, Z.; Yang, W.; Huang, C.; Liu, P.; Chen, Y.; Zhang, H.; Zhao, G.; Li, W.; Fang, J.; et al. Neural activity changes in first-episode, drug-naive patients with major depressive disorder after transcutaneous auricular vagus nerve stimulation treatment: A resting-state fMRI study. Front. Neurosci. 2022, 16, 1018387. [Google Scholar] [CrossRef] [PubMed]
  4. van Rooij, S.J.H.; Arulpragasam, A.R.; McDonald, W.M.; Philip, N.S. Accelerated TMS—Moving quickly into the future of depression treatment. Neuropsychopharmacology 2024, 49, 128–137. [Google Scholar] [CrossRef]
  5. Burkhardt, G.; Kumpf, U.; Crispin, A.; Goerigk, S.; Andre, E.; Plewnia, C.; Brendel, B.; Fallgatter, A.; Langguth, B.; Abdelnaim, M.; et al. Transcranial direct current stimulation as an additional treatment to selective serotonin reuptake inhibitors in adults with major depressive disorder in Germany (DepressionDC): A triple-blind, randomised, sham-controlled, multicentre trial. Lancet 2023, 402, 545–554. [Google Scholar] [CrossRef]
  6. Johnson, K.A.; Okun, M.S.; Scangos, K.W.; Mayberg, H.S.; de Hemptinne, C. Deep brain stimulation for refractory major depressive disorder: A comprehensive review. Mol. Psychiatry 2024, 29, 1075–1087. [Google Scholar] [CrossRef]
  7. Guo, Q.; Wang, Y.; Guo, L.; Li, X.; Ma, X.; He, X.; Li, J.; Zhang, X.; Shang, S. Long-term cognitive effects of electroconvulsive therapy in major depressive disorder: A systematic review and meta-analysis. Psychiatry Res. 2024, 331, 115611. [Google Scholar] [CrossRef]
  8. Sun, J.; Guo, C.; Ma, Y.; Gao, S.; Luo, Y.; Chen, Q.; Hong, Y.; Hou, X.; Xiao, X.; Yu, X.; et al. Immediate modulatory effects of transcutaneous auricular vagus nerve stimulation on the resting state of major depressive disorder. J. Affect. Disord. 2023, 325, 513–521. [Google Scholar] [CrossRef]
  9. Kong, J.; Fang, J.; Park, J.; Li, S.; Rong, P. Treating Depression with Transcutaneous Auricular Vagus Nerve Stimulation: State of the Art and Future Perspectives. Front. Psychiatry 2018, 9, 20. [Google Scholar] [CrossRef]
  10. Wang, Y.; Li, S.Y.; Wang, D.; Wu, M.Z.; He, J.K.; Zhang, J.L.; Zhao, B.; Hou, L.W.; Wang, J.Y.; Wang, L.; et al. Transcutaneous Auricular Vagus Nerve Stimulation: From Concept to Application. Neurosci. Bull. 2021, 37, 853–862. [Google Scholar] [CrossRef]
  11. Rong, P.; Liu, J.; Wang, L.; Liu, R.; Fang, J.; Zhao, J.; Zhao, Y.; Wang, H.; Vangel, M.; Sun, S.; et al. Effect of transcutaneous auricular vagus nerve stimulation on major depressive disorder: A nonrandomized controlled pilot study. J. Affect. Disord. 2016, 195, 172–179. [Google Scholar] [CrossRef] [PubMed]
  12. Hein, E.; Nowak, M.; Kiess, O.; Biermann, T.; Bayerlein, K.; Kornhuber, J.; Kraus, T. Auricular transcutaneous electrical nerve stimulation in depressed patients: A randomized controlled pilot study. J. Neural Transm. (Vienna) 2013, 120, 821–827. [Google Scholar] [CrossRef] [PubMed]
  13. Ventureyra, E.C. Transcutaneous vagus nerve stimulation for partial onset seizure therapy. A new concept. Childs Nerv. Syst. 2000, 16, 101–102. [Google Scholar] [CrossRef] [PubMed]
  14. Yang, H.; Shi, W.; Fan, J.; Wang, X.; Song, Y.; Lian, Y.; Shan, W.; Wang, Q. Transcutaneous Auricular Vagus Nerve Stimulation (ta-VNS) for Treatment of Drug-Resistant Epilepsy: A Randomized, Double-Blind Clinical Trial. Neurotherapeutics 2023, 20, 870–880. [Google Scholar] [CrossRef]
  15. Mao, Y.; Chen, C.; Falahpour, M.; MacNiven, K.H.; Heit, G.; Sharma, V.; Alataris, K.; Liu, T.T. Effects of Sub-threshold Transcutaneous Auricular Vagus Nerve Stimulation on Cingulate Cortex and Insula Resting-state Functional Connectivity. Front. Hum. Neurosci. 2022, 16, 862443. [Google Scholar] [CrossRef]
  16. Fang, J.; Rong, P.; Hong, Y.; Fan, Y.; Liu, J.; Wang, H.; Zhang, G.; Chen, X.; Shi, S.; Wang, L.; et al. Transcutaneous Vagus Nerve Stimulation Modulates Default Mode Network in Major Depressive Disorder. Biol. Psychiatry 2016, 79, 266–273. [Google Scholar] [CrossRef]
  17. van Midden, V.M.; Pirtosek, Z.; Kojovic, M. The Effect of taVNS on the Cerebello-Thalamo-Cortical Pathway: A TMS Study. Cerebellum 2024, 23, 1013–1019. [Google Scholar] [CrossRef]
  18. Philip, N.S.; Barredo, J.; van ’t Wout-Frank, M.; Tyrka, A.R.; Price, L.H.; Carpenter, L.L. Network Mechanisms of Clinical Response to Transcranial Magnetic Stimulation in Posttraumatic Stress Disorder and Major Depressive Disorder. Biol. Psychiatry 2018, 83, 263–272. [Google Scholar] [CrossRef] [PubMed]
  19. Trevizol, A.P.; Shiozawa, P.; Taiar, I.; Soares, A.; Gomes, J.S.; Barros, M.D.; Liquidato, B.M.; Cordeiro, Q. Transcutaneous Vagus Nerve Stimulation (taVNS) for Major Depressive Disorder: An Open Label Proof-of-Concept Trial. Brain Stimul. 2016, 9, 453–454. [Google Scholar] [CrossRef]
  20. Chellappa, S.L.; Aeschbach, D. Sleep and anxiety: From mechanisms to interventions. Sleep. Med. Rev. 2022, 61, 101583. [Google Scholar] [CrossRef]
  21. Kang, S.G.; Cho, S.E. Neuroimaging Biomarkers for Predicting Treatment Response and Recurrence of Major Depressive Disorder. Int. J. Mol. Sci. 2020, 21, 2148. [Google Scholar] [CrossRef] [PubMed]
  22. Li, Y.; Li, Y.; Wei, Q.; Bai, T.; Wang, K.; Wang, J.; Tian, Y. Mapping intrinsic functional network topological architecture in major depression disorder after electroconvulsive therapy. J. Affect. Disord. 2022, 311, 103–109. [Google Scholar] [CrossRef] [PubMed]
  23. Korhonen, O.; Zanin, M.; Papo, D. Principles and open questions in functional brain network reconstruction. Hum. Brain Mapp. 2021, 42, 3680–3711. [Google Scholar] [CrossRef] [PubMed]
  24. Yun, J.Y.; Kim, Y.K. Graph theory approach for the structural-functional brain connectome of depression. Prog. Neuropsychopharmacol. Biol. Psychiatry 2021, 111, 110401. [Google Scholar] [CrossRef]
  25. Zhou, Y.; Zhu, Y.; Ye, H.; Jiang, W.; Zhang, Y.; Kong, Y.; Yuan, Y.; Consortium, R.E.-m.-M. Abnormal changes of dynamic topological characteristics in patients with major depressive disorder. J. Affect. Disord. 2024, 345, 349–357. [Google Scholar] [CrossRef]
  26. Li, Y.; Chu, T.; Che, K.; Dong, F.; Shi, Y.; Ma, H.; Zhao, F.; Mao, N.; Xie, H. Altered gray matter structural covariance networks in postpartum depression: A graph theoretical analysis. J. Affect. Disord. 2021, 293, 159–167. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, R.; Kranz, G.S.; Zou, W.; Deng, Y.; Huang, X.; Lin, K.; Lee, T.M.C. Rumination network dysfunction in major depression: A brain connectome study. Prog. Neuropsychopharmacol. Biol. Psychiatry 2020, 98, 109819. [Google Scholar] [CrossRef]
  28. Jacob, Y.; Morris, L.S.; Huang, K.H.; Schneider, M.; Rutter, S.; Verma, G.; Murrough, J.W.; Balchandani, P. Neural correlates of rumination in major depressive disorder: A brain network analysis. Neuroimage Clin. 2020, 25, 102142. [Google Scholar] [CrossRef]
  29. Dai, Y.R.; Wu, Y.K.; Chen, X.; Zeng, Y.W.; Li, K.; Li, J.T.; Su, Y.A.; Zhu, L.L.; Yan, C.G.; Si, T.M. Eight-week antidepressant treatment changes intrinsic functional brain topology in first-episode drug-naive patients with major depressive disorder. J. Affect. Disord. 2023, 329, 225–234. [Google Scholar] [CrossRef]
  30. Wu, Y.; Ji, Y.; Bai, T.; Wei, Q.; Zu, M.; Guo, Y.; Lv, H.; Zhang, A.; Qiu, B.; Wang, K.; et al. Nodal degree changes induced by electroconvulsive therapy in major depressive disorder: Evidence in two independent cohorts. J. Affect. Disord. 2022, 307, 46–52. [Google Scholar] [CrossRef]
  31. Croarkin, P.E.; Elmaadawi, A.Z.; Aaronson, S.T.; Schrodt, G.R., Jr.; Holbert, R.C.; Verdoliva, S.; Heart, K.L.; Demitrack, M.A.; Strawn, J.R. Left prefrontal transcranial magnetic stimulation for treatment-resistant depression in adolescents: A double-blind, randomized, sham-controlled trial. Neuropsychopharmacology 2021, 46, 462–469. [Google Scholar] [CrossRef] [PubMed]
  32. Thompson, E. Hamilton Rating Scale for Anxiety (HAM-A). Occup. Med. (Lond.) 2015, 65, 601. [Google Scholar] [CrossRef] [PubMed]
  33. Shi, D.; Li, Z.; Li, Y.; Jiang, Q. Variables associated with self-reported anxiety and depression symptoms in patients with chronic myeloid leukemia receiving tyrosine kinase inhibitor therapy. Leuk. Lymphoma 2021, 62, 640–648. [Google Scholar] [CrossRef] [PubMed]
  34. 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. [Google Scholar] [CrossRef]
  35. Dosenbach, N.U.; Nardos, B.; Cohen, A.L.; Fair, D.A.; Power, J.D.; Church, J.A.; Nelson, S.M.; Wig, G.S.; Vogel, A.C.; Lessov-Schlaggar, C.N.; et al. Prediction of individual brain maturity using fMRI. Science 2010, 329, 1358–1361. [Google Scholar] [CrossRef] [PubMed]
  36. Amiri, S.; Arbabi, M.; Kazemi, K.; Parvaresh-Rizi, M.; Mirbagheri, M.M. Characterization of brain functional connectivity in treatment-resistant depression. Prog. Neuropsychopharmacol. Biol. Psychiatry 2021, 111, 110346. [Google Scholar] [CrossRef]
  37. Korgaonkar, M.S.; Fornito, A.; Williams, L.M.; Grieve, S.M. Abnormal structural networks characterize major depressive disorder: A connectome analysis. Biol. Psychiatry 2014, 76, 567–574. [Google Scholar] [CrossRef]
  38. Li, L.; Su, Y.A.; Wu, Y.K.; Castellanos, F.X.; Li, K.; Li, J.T.; Si, T.M.; Yan, C.G. Eight-week antidepressant treatment reduces functional connectivity in first-episode drug-naive patients with major depressive disorder. Hum. Brain Mapp. 2021, 42, 2593–2605. [Google Scholar] [CrossRef]
  39. Liu, J.; Fang, J.; Wang, Z.; Rong, P.; Hong, Y.; Fan, Y.; Wang, X.; Park, J.; Jin, Y.; Liu, C.; et al. Transcutaneous vagus nerve stimulation modulates amygdala functional connectivity in patients with depression. J. Affect. Disord. 2016, 205, 319–326. [Google Scholar] [CrossRef]
  40. Wang, X.; Xia, Y.; Yan, R.; Wang, H.; Sun, H.; Huang, Y.; Hua, L.; Tang, H.; Yao, Z.; Lu, Q. The relationship between disrupted anhedonia-related circuitry and suicidal ideation in major depressive disorder: A network-based analysis. Neuroimage Clin. 2023, 40, 103512. [Google Scholar] [CrossRef]
  41. Yang, H.; Chen, X.; Chen, Z.B.; Li, L.; Li, X.Y.; Castellanos, F.X.; Bai, T.J.; Bo, Q.J.; Cao, J.; Chang, Z.K.; et al. Disrupted intrinsic functional brain topology in patients with major depressive disorder. Mol. Psychiatry 2021, 26, 7363–7371. [Google Scholar] [CrossRef] [PubMed]
  42. Lv, X.; Lu, F.; Zhang, J.; Chen, H.; Zhang, L.; Wang, X.; Fan, Y.; Fang, J.; Hong, L.; Wang, J.; et al. Effects of TIP treatment on brain network topology of frontolimbic circuit in first-episode, treatment-naive major depressive disorder. J. Affect. Disord. 2021, 279, 122–130. [Google Scholar] [CrossRef] [PubMed]
  43. Daws, R.E.; Timmermann, C.; Giribaldi, B.; Sexton, J.D.; Wall, M.B.; Erritzoe, D.; Roseman, L.; Nutt, D.; Carhart-Harris, R. Increased global integration in the brain after psilocybin therapy for depression. Nat. Med. 2022, 28, 844–851. [Google Scholar] [CrossRef]
  44. He, M.; Shen, Z.; Ping, L.; Zhou, C.; Cheng, Y.; Xu, X. Age-related heterogeneity revealed by disruption of white matter structural networks in patients with first-episode untreated major depressive disorder: WM Network In OA-MDD. J. Affect. Disord. 2022, 303, 286–296. [Google Scholar] [CrossRef] [PubMed]
  45. Manelis, A.; Almeida, J.R.; Stiffler, R.; Lockovich, J.C.; Aslam, H.A.; Phillips, M.L. Anticipation-related brain connectivity in bipolar and unipolar depression: A graph theory approach. Brain 2016, 139, 2554–2566. [Google Scholar] [CrossRef]
  46. Zhang, Y.; Liu, X.; Hou, Z.; Yin, Y.; Xie, C.; Zhang, H.; Zhang, H.; Kong, Y.; Gao, S.; Zhang, Z.; et al. Global topology alteration of the brain functional network affects the 8-week antidepressant response in major depressive disorder. J. Affect. Disord. 2021, 294, 491–496. [Google Scholar] [CrossRef] [PubMed]
  47. Ghanbari, M.; Soussia, M.; Jiang, W.; Wei, D.; Yap, P.T.; Shen, D.; Zhang, H. Alterations of dynamic redundancy of functional brain subnetworks in Alzheimer’s disease and major depression disorders. Neuroimage Clin. 2022, 33, 102917. [Google Scholar] [CrossRef] [PubMed]
  48. Jacob, Y.; Morris, L.S.; Verma, G.; Rutter, S.B.; Balchandani, P.; Murrough, J.W. Altered hippocampus and amygdala subregion connectome hierarchy in major depressive disorder. Transl. Psychiatry 2022, 12, 209. [Google Scholar] [CrossRef]
  49. Zhang, A.; Qiao, D.; Wang, Y.; Yang, C.; Wang, Y.; Sun, N.; Hu, X.; Liu, Z.; Zhang, K. Distinguishing between bipolar depression and unipolar depression based on the reward circuit activities and clinical characteristics: A machine learning analysis. J. Affect. Disord. 2023, 327, 46–53. [Google Scholar] [CrossRef]
  50. Li, X.; Steffens, D.C.; Potter, G.G.; Guo, H.; Song, S.; Wang, L. Decreased between-hemisphere connectivity strength and network efficiency in geriatric depression. Hum. Brain Mapp. 2017, 38, 53–67. [Google Scholar] [CrossRef]
  51. Mo, Y.; Wei, Q.; Bai, T.; Zhang, T.; Lv, H.; Zhang, L.; Ji, G.; Yu, F.; Tian, Y.; Wang, K. Bifrontal electroconvulsive therapy changed regional homogeneity and functional connectivity of left angular gyrus in major depressive disorder. Psychiatry Res. 2020, 294, 113461. [Google Scholar] [CrossRef] [PubMed]
  52. Tan, W.; Ouyang, X.; Huang, D.; Wu, Z.; Liu, Z.; He, Z.; Long, Y.; Consortium, R.E.-m.-M. Disrupted intrinsic functional brain network in patients with late-life depression: Evidence from a multi-site dataset. J. Affect. Disord. 2023, 323, 631–639. [Google Scholar] [CrossRef] [PubMed]
  53. Liu, J.; Xu, X.; Zhu, C.; Luo, L.; Wang, Q.; Xiao, B.; Feng, B.; Hu, L.; Liu, L. Disrupted Structural Brain Network Organization Behind Depressive Symptoms in Major Depressive Disorder. Front. Psychiatry 2020, 11, 565890. [Google Scholar] [CrossRef] [PubMed]
  54. Brakowski, J.; Spinelli, S.; Dorig, N.; Bosch, O.G.; Manoliu, A.; Holtforth, M.G.; Seifritz, E. Resting state brain network function in major depression—Depression symptomatology, antidepressant treatment effects, future research. J. Psychiatr. Res. 2017, 92, 147–159. [Google Scholar] [CrossRef]
  55. Ye, M.; Qing, P.; Zhang, K.; Liu, G. Altered network efficiency in major depressive disorder. BMC Psychiatry 2016, 16, 450. [Google Scholar] [CrossRef] [PubMed]
  56. Lu, F.; Chen, Y.; Cui, Q.; Guo, Y.; Pang, Y.; Luo, W.; Yu, Y.; Chen, J.; Gao, J.; Sheng, W.; et al. Shared and distinct patterns of dynamic functional connectivity variability of thalamo-cortical circuit in bipolar depression and major depressive disorder. Cereb. Cortex 2023, 33, 6681–6692. [Google Scholar] [CrossRef]
  57. Zanto, T.P.; Gazzaley, A. Fronto-parietal network: Flexible hub of cognitive control. Trends Cogn. Sci. 2013, 17, 602–603. [Google Scholar] [CrossRef]
  58. Kaiser, R.H.; Andrews-Hanna, J.R.; Wager, T.D.; Pizzagalli, D.A. Large-Scale Network Dysfunction in Major Depressive Disorder: A Meta-analysis of Resting-State Functional Connectivity. JAMA Psychiatry 2015, 72, 603–611. [Google Scholar] [CrossRef]
  59. Rubin-Falcone, H.; Weber, J.; Kishon, R.; Ochsner, K.; Delaparte, L.; Dore, B.; Raman, S.; Denny, B.T.; Oquendo, M.A.; Mann, J.J.; et al. Neural predictors and effects of cognitive behavioral therapy for depression: The role of emotional reactivity and regulation. Psychol. Med. 2020, 50, 146–160. [Google Scholar] [CrossRef]
  60. Yan, C.G.; Chen, X.; Li, L.; Castellanos, F.X.; Bai, T.J.; Bo, Q.J.; Cao, J.; Chen, G.M.; Chen, N.X.; Chen, W.; et al. Reduced default mode network functional connectivity in patients with recurrent major depressive disorder. Proc. Natl. Acad. Sci. USA 2019, 116, 9078–9083. [Google Scholar] [CrossRef]
  61. Hawkey, E.J.; Tillman, R.; Luby, J.L.; Barch, D.M. Preschool Executive Function Predicts Childhood Resting-State Functional Connectivity and Attention-Deficit/Hyperactivity Disorder and Depression. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 2018, 3, 927–936. [Google Scholar] [CrossRef] [PubMed]
  62. Shao, J.; Meng, C.; Tahmasian, M.; Brandl, F.; Yang, Q.; Luo, G.; Luo, C.; Yao, D.; Gao, L.; Riedl, V.; et al. Common and distinct changes of default mode and salience network in schizophrenia and major depression. Brain Imaging Behav. 2018, 12, 1708–1719. [Google Scholar] [CrossRef] [PubMed]
  63. Wu, X.; Lin, P.; Yang, J.; Song, H.; Yang, R.; Yang, J. Dysfunction of the cingulo-opercular network in first-episode medication-naive patients with major depressive disorder. J. Affect. Disord. 2016, 200, 275–283. [Google Scholar] [CrossRef] [PubMed]
  64. Argyelan, M.; Lencz, T.; Kaliora, S.; Sarpal, D.K.; Weissman, N.; Kingsley, P.B.; Malhotra, A.K.; Petrides, G. Subgenual cingulate cortical activity predicts the efficacy of electroconvulsive therapy. Transl. Psychiatry 2016, 6, e789. [Google Scholar] [CrossRef]
  65. Kang, J.I.; Lee, H.; Jhung, K.; Kim, K.R.; An, S.K.; Yoon, K.J.; Kim, S.I.; Namkoong, K.; Lee, E. Frontostriatal Connectivity Changes in Major Depressive Disorder After Repetitive Transcranial Magnetic Stimulation: A Randomized Sham-Controlled Study. J. Clin. Psychiatry 2016, 77, e1137–e1143. [Google Scholar] [CrossRef]
Brainsci 14 00945 g001
Figure 1. Study workflow. (A) taVNS intervention and bold signal extraction. (B) Graph theory analysis. (C) Statistical analysis. Abbreviations: taVNS, transcutaneous auricular vagus nerve stimulation. * p < 0.05, *** p < 0.001.
Figure 1. Study workflow. (A) taVNS intervention and bold signal extraction. (B) Graph theory analysis. (C) Statistical analysis. Abbreviations: taVNS, transcutaneous auricular vagus nerve stimulation. * p < 0.05, *** p < 0.001.
Brainsci 14 00945 g001
Brainsci 14 00945 g002
Figure 2. Differences in global network topological properties of patients with MDD across the sparsity range (0.05–0.4). (A) Gamma in patients with MDD before and after taVNS treatment; (B) lambda in patients with MDD; (C) sigma in patients with MDD; (D) Lp in patients with MDD; (E) Eglob in patients with MDD; (F) Eloc in patients with MDD. Abbreviations: Lp, characteristic path length; MDD, major depressive disorder; Eglob, global efficiency; taVNS, transcutaneous auricular vagus nerve stimulation; Eloc, local efficiency.
Figure 2. Differences in global network topological properties of patients with MDD across the sparsity range (0.05–0.4). (A) Gamma in patients with MDD before and after taVNS treatment; (B) lambda in patients with MDD; (C) sigma in patients with MDD; (D) Lp in patients with MDD; (E) Eglob in patients with MDD; (F) Eloc in patients with MDD. Abbreviations: Lp, characteristic path length; MDD, major depressive disorder; Eglob, global efficiency; taVNS, transcutaneous auricular vagus nerve stimulation; Eloc, local efficiency.
Brainsci 14 00945 g002
Brainsci 14 00945 g003
Figure 3. Differences in global properties of patients with MDD before and after taVNS treatment based on AUC values. (A) Eglob in patients with MDD; (B) Lp in patients with MDD. Abbreviations: Eglob, global efficiency; taVNS, transcutaneous auricular vagus nerve stimulation; Lp, characteristic path length; MDD, major depressive disorder; AUC, area under the curve. *** p < 0.001.
Figure 3. Differences in global properties of patients with MDD before and after taVNS treatment based on AUC values. (A) Eglob in patients with MDD; (B) Lp in patients with MDD. Abbreviations: Eglob, global efficiency; taVNS, transcutaneous auricular vagus nerve stimulation; Lp, characteristic path length; MDD, major depressive disorder; AUC, area under the curve. *** p < 0.001.
Brainsci 14 00945 g003
Brainsci 14 00945 g004
Figure 4. Differences in nodal properties of patients with MDD before and after taVNS treatment based on the AUC values. (A) The left AG; (B) the DC of the left AG in patients with MDD; (C) the NE of the left AG in patients with MDD. Abbreviations: AG, angular gyrus; AUC, area under the curve; DC, degree centrality; taVNS, transcutaneous auricular vagus nerve stimulation; NE, nodal efficiency; MDD, major depressive disorder. * p < 0.05, *** p < 0.001.
Figure 4. Differences in nodal properties of patients with MDD before and after taVNS treatment based on the AUC values. (A) The left AG; (B) the DC of the left AG in patients with MDD; (C) the NE of the left AG in patients with MDD. Abbreviations: AG, angular gyrus; AUC, area under the curve; DC, degree centrality; taVNS, transcutaneous auricular vagus nerve stimulation; NE, nodal efficiency; MDD, major depressive disorder. * p < 0.05, *** p < 0.001.
Brainsci 14 00945 g004
Brainsci 14 00945 g005
Figure 5. Differences in network connectivity of patients with MDD before and after taVNS treatment. The blue lines indicate the decreased functional connections between the DMN and FPN, DMN and CON, and FPN and CON. Abbreviations: CON, cingulo-opercular network; CN, cerebellar network; DMN, default mode network; MDD, major depressive disorder; FPN, frontoparietal network; SMN, sensorimotor network; taVNS, transcutaneous auricular vagus nerve stimulation; VN, visual network.
Figure 5. Differences in network connectivity of patients with MDD before and after taVNS treatment. The blue lines indicate the decreased functional connections between the DMN and FPN, DMN and CON, and FPN and CON. Abbreviations: CON, cingulo-opercular network; CN, cerebellar network; DMN, default mode network; MDD, major depressive disorder; FPN, frontoparietal network; SMN, sensorimotor network; taVNS, transcutaneous auricular vagus nerve stimulation; VN, visual network.
Brainsci 14 00945 g005
Brainsci 14 00945 g006
Figure 6. Relationships between changes in global and nodal network properties and changes in HAMD scores of patients with MDD after taVNS treatment. (A) Change in Lp was positively associated with alternation in HAMD scores (r = 0.616, p = 0.005). (B) Changes in DC in the left AG were negatively associated with changes in HAMD scores (r = −0.703, p = 0.001). Abbreviations: AG, angular gyrus; DC, degree centrality; MDD, major depressive disorder; HAMD, 24-Hamilton Depression Rating Scale; Lp, characteristic path length; taVNS, transcutaneous auricular vagus nerve stimulation.
Figure 6. Relationships between changes in global and nodal network properties and changes in HAMD scores of patients with MDD after taVNS treatment. (A) Change in Lp was positively associated with alternation in HAMD scores (r = 0.616, p = 0.005). (B) Changes in DC in the left AG were negatively associated with changes in HAMD scores (r = −0.703, p = 0.001). Abbreviations: AG, angular gyrus; DC, degree centrality; MDD, major depressive disorder; HAMD, 24-Hamilton Depression Rating Scale; Lp, characteristic path length; taVNS, transcutaneous auricular vagus nerve stimulation.
Brainsci 14 00945 g006
Table 1. Demographic and clinical characteristics of participants.
Table 1. Demographic and clinical characteristics of participants.
MDD (n = 19)
Sex (male/female)13/6
Age (years)38.89 ± 14.48
HAMD29.11 ± 5.59
HAMA21.00 ± 7.10
SAS57.53 ± 9.01
SDS66.89 ± 10.70
Abbreviations: HAMA, Hamilton Anxiety Rating Scale; HAMD, 24-Hamilton Depression Rating Scale; MDD, major depressive disorder; SAS, Self-Rating Anxiety Scale; SDS, Self-Rating Depression Scale.
Table 2. Changes in neuropsychological scales in patients with MDD after taVNS treatment.
Table 2. Changes in neuropsychological scales in patients with MDD after taVNS treatment.
PrePosttp
HAMD29.11 ± 5.5915.05 ± 4.479.692<0.001
HAMA21.00 ± 7.1010.74 ± 4.047.446<0.001
SAS57.53 ± 9.0142.89 ± 9.677.664<0.001
SDS66.89 ± 10.7050.84 ± 11.787.160<0.001
Abbreviations: HAMA, Hamilton Anxiety Rating Scale; HAMD, 24-Hamilton Depression Rating Scale; MDD, major depressive disorder; SAS, Self-Rating Anxiety Scale; SDS, Self-Rating Depression Scale; taVNS, transcutaneous vagus nerve stimulation.
Table 3. Changes in topological properties in patients with MDD after taVNS treatment.
Table 3. Changes in topological properties in patients with MDD after taVNS treatment.
PrePosttp
Eglob0.269 ± 0.0040.274 ± 0.004−4.537<0.001
Lp0.790 ± 0.0250.760 ± 0.0156.353<0.001
DC of left AG18.43 ± 5.7225.34 ± 3.71−7.003<0.001
Ne of left AG0.274 ± 0.0220.291 ± 0.025−2.1230.048
Abbreviations: MDD, major depressive disorder; taVNS, transcutaneous vagus nerve stimulation; Eglob, global efficiency; Lp, characteristic path length; DC, degree centrality; AG, angular gyrus; Ne, nodal efficiency.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Z.-P.; Liao, D.; Chen, L.; Wang, C.; Qu, M.; Lv, X.-Y.; Fang, J.-L.; Liu, C.-H. Transcutaneous Auricular Vagus Nerve Stimulation Modulating the Brain Topological Architecture of Functional Network in Major Depressive Disorder: An fMRI Study. Brain Sci. 2024, 14, 945. https://doi.org/10.3390/brainsci14090945

AMA Style

Guo Z-P, Liao D, Chen L, Wang C, Qu M, Lv X-Y, Fang J-L, Liu C-H. Transcutaneous Auricular Vagus Nerve Stimulation Modulating the Brain Topological Architecture of Functional Network in Major Depressive Disorder: An fMRI Study. Brain Sciences. 2024; 14(9):945. https://doi.org/10.3390/brainsci14090945

Chicago/Turabian Style

Guo, Zhi-Peng, Dan Liao, Lei Chen, Cong Wang, Miao Qu, Xue-Yu Lv, Ji-Liang Fang, and Chun-Hong Liu. 2024. "Transcutaneous Auricular Vagus Nerve Stimulation Modulating the Brain Topological Architecture of Functional Network in Major Depressive Disorder: An fMRI Study" Brain Sciences 14, no. 9: 945. https://doi.org/10.3390/brainsci14090945

APA Style

Guo, Z. -P., Liao, D., Chen, L., Wang, C., Qu, M., Lv, X. -Y., Fang, J. -L., & Liu, C. -H. (2024). Transcutaneous Auricular Vagus Nerve Stimulation Modulating the Brain Topological Architecture of Functional Network in Major Depressive Disorder: An fMRI Study. Brain Sciences, 14(9), 945. https://doi.org/10.3390/brainsci14090945

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop