Front. Neurosci.Frontiers in NeuroscienceFront. Neurosci.1662-453XFrontiers Media S.A.10.3389/fnins.2022.855483NeuroscienceOriginal ResearchAltered Brain Function and Causal Connectivity Induced by Repetitive Transcranial Magnetic Stimulation Treatment for Major Depressive DisorderGuanMuzhen12*†WangZhongheng2†ShiYanru2XieYuanjun3MaZhujing4LiuZirong5LiuJunchang2GaoXinyu4TanQingrong2WangHuaning2*1Department of Mental Health, Xi’an Medical University, Xi’an, China2Deptartment of Psychiatry, Xijing Hospital, Air Force Medical University, Xi’an, China3Department of Radiology, Xijing Hospital, Fourth Military Medical University, Xi’an, China4Deptartment of Psychology, Air Force Medical University, Xi’an, China5Deptartment of Psychiatry, Yulin Fifth Hospital, Yuling, China
Edited by: Guolin Ma, China-Japan Friendship Hospital, China
Reviewed by: Weihua Yue, Peking University Sixth Hospital, China; Xiaowei Ma, Central South University, China
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Objective
Repetitive transcranial magnetic stimulation (rTMS) can effectively improve depression symptoms in patients with major depressive disorder (MDD); however, its mechanism of action remains obscure. This study explored the neuralimaging mechanisms of rTMS in improving depression symptoms in patients with MDD.
Methods
In this study, MDD patients with first-episode, drug-naive (n = 29) and healthy controls (n = 33) were enrolled. Depression symptoms before and after rTMS treatment were assessed using the Hamilton Depression Rating Scale (HAMD-17). Resting-state functional magnetic resonance imaging (rs-fMRI) data were collected both before and after the treatment. Changes in the brain function after the treatment were compared using the following two indices: the amplitude of the low-frequency fluctuation (ALFF) and regional homogeneity (ReHo), which are sensitive for evaluating spontaneous neuronal activity. The brain region with synchronous changes was selected as the seed point, and the differences in the causal connectivity between the seed point and whole brain before and after rTMS treatment were investigated via Granger causality analysis (GCA).
Results
Before treatment, patients with MDD had significantly lower ALFF in the left superior frontal gyrus (p < 0.01), higher ALFF in the left middle frontal gyrus and left precuneus (p < 0.01), and lower ReHo in the left middle frontal and left middle occipital gyri (p < 0.01) than the values observed in healthy controls. After the rTMS treatment, the ALFF was significantly increased in the left superior frontal gyrus (p < 0.01) and decreased in the left middle frontal gyrus and left precuneus (p < 0.01). Furthermore, ReHo was significantly increased in the left middle frontal and left middle occipital gyri (p < 0.01) in patients with MDD. Before treatment, GCA using the left middle frontal gyrus (the brain region with synchronous changes) as the seed point revealed a weak bidirectional causal connectivity between the middle and superior frontal gyri as well as a weak causal connectivity from the inferior temporal to the middle frontal gyri. After treatment, these causal connectivities were strengthened. Moreover, the causal connectivity from the inferior temporal gyrus to the middle frontal gyri negatively correlated with the total HAMD-17 score (r = −0.443, p = 0.021).
Conclusion
rTMS treatment not only improves the local neural activity in the middle frontal gyrus, superior frontal gyrus, and precuneus but also strengthens the bidirectional causal connectivity between the middle and superior frontal gyri and the causal connectivity from the inferior temporal to the middle frontal gyri. Changes in these neuroimaging indices may represent the neural mechanisms underlying rTMS treatment in MDD.
Clinical Trial Registration
This study was registered in the Chinese Clinical Trial Registry (Registration number: ChiCTR1800019761).
major depressive disorderamplitude of the low-frequency fluctuationregional homogeneityGranger causality analysis (GCA)repetitive transcranial magnetic stimulationIntroduction
Major depressive disorder (MDD) is a mental illness characterized by mood disorders and cognitive dysfunction. By the year 2020, MDD had become the second most disabling illness and a major social and public health concern (Cheng et al., 2018). Blood-oxygen-level-dependent magnetic resonance imaging works based on the principle of detecting changes in the magnetic field properties due to the changes in local metabolism and blood oxygen level after neuronal activity (Tian et al., 2020). This imaging method has been widely used in exploring the pathophysiology of and therapeutic mechanism for MDD (Tateishi et al., 2020). The amplitude of the low-frequency fluctuation (ALFF) in the blood-oxygen-level-dependent signal is used to describe the amplitude of the spontaneous activity of the brain neurons, and larger ALFF indicates greater activity in the brain region. Regional homogeneity (ReHo) refers to the correlation of the synchronous activity in the blood-oxygen-level-dependent signal between a given voxel and its neighboring voxels (Dichter et al., 2015). The ALFF and ReHo can evaluate the spontaneous activity of neurons from different perspectives, and these two indices are complementary to each other (Kong et al., 2017).
Patients with MDD have been reported to exhibit abnormal ALFF and ReHo in the brain function (Yao et al., 2014; Li et al., 2018), A recent study investigated the MDD patients with negative bias and showed significantly higher neuronal activation in the frontal lobe (Gollan et al., 2015; Tadayonnejad et al., 2015). Studies have also indicated increased ReHo in the right postcentral gyrus and increased ALFF in the left triangular part of the inferior frontal gyrus (Xia et al., 2019; Rosenbaum et al., 2020; Liu et al., 2021). Liang et al. (2020) demonstrated that compared with healthy controls, patients with MDD had lower ReHo in the calcarine gyrus, insula, right superior temporal gyrus, left cuneus, right precuneus and right postcentral gyrus as well as lower ALFF in the calcarine gyrus, left cuneus, right inferior parietal cortex, right precuneus, and right postcentral gyrus (Kaiser et al., 2015), indicating that brain dysfunction in these patients affects their emotion regulation and social function.
Medication is the primary treatment modality used for patients with MDD. Antidepressants usually begin to show effects after 2 weeks of intake; hence, such medications (e.g., venlafaxine) should be taken for at least 4 weeks before determining the drug efficacy (Gex-Fabry et al., 2004; Otte et al., 2016). However, the first-line medication and psychotherapy often fail to show any improvement in the condition of 20–30% of patients with MDD (Loerinc et al., 2015; Maneeton et al., 2020; Cohen et al., 2021), repetitive transcranial magnetic stimulation—a non-invasive and well-tolerated physical therapy (McGirr et al., 2015) in the domain of antidepressant therapy—has been an important treatment choice (Antonenko et al., 2019; Yin et al., 2020). Zheng et al. (2020) collected the functional magnetic resonance imaging (fMRI) data of 15 patients with MDD both before and after repetitive transcranial magnetic stimulation (rTMS) treatment and found an increase in the ALFF in the left dorsolateral prefrontal cortex and left superior frontal gyrus after 2 weeks of treatment; however, they found no correlation between the change in the ALFF and symptoms of depression.
Granger causality analysis (GCA) is a method used for identifying the effective connectivity of the directed functional interactions from time-series data, and it has been widely used in resting-state fMRI studies (Berlim et al., 2014; Gaynes et al., 2014; Yang et al., 2021). Feng et al. (2016) analyzed the data of 23 drug-naive patients with first-onset MDD using GCA. They revealed that these patients had a stronger causal connectivity from the right insula, right putamen, and right caudate nucleus to the cingulate gyrus seed point as well as a weaker causal connectivity from the bilateral dorsolateral prefrontal and left orbitofrontal cortices to the cingulate gyrus seed point than those of the healthy controls (Feng et al., 2016); this indicated the dysfunction of the prefrontal cortex–limbic system in patients with MDD (Ishida et al., 2020).
However, the mechanism through which rTMS improves depression symptoms remains unclear, and few studies have explored the changes in both the ALFF and ReHo of patients with MDD receiving rTMS treatment. In this study, resting-state fMRI (rs-fMRI) was used to evaluate the changes in the abovementioned indices in patients with MDD after rTMS treatment. The brain region with synchronous changes was selected as the seed point, and the causal connectivity between the seed point and whole brain was analyzed via GCA to investigate the mechanism through which rTMS improves the symptoms of depression in patients with MDD.
Materials and MethodsParticipants
This study included drug-naive patients with first-onset MDD aged 18–45 years who visited the psychiatry outpatient clinic of the First Affiliated Hospital of Air Force Military Medical University, China, from May 2019 to October 2021. The inclusion criteria were as follows: (1) Patients who fulfilled the MDD diagnostic criteria of the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition, (2) patients with a total score of ≥ 18 points on the Hamilton Depression Rating Scale (HAMD-17), and (3) patients who met the enrollment requirements of the study and provided their informed consent. The exclusion criteria were as follows: (1) patients with severe physical diseases, (2) patients with a history of traumatic brain injury or brain surgery, (3) patients with a history of alcohol or substance abuse, and patients with a history of other mental or nervous system disorders.
Healthy controls matching the characteristics of patients with MDD were also recruited in this study. All participants read and understood the experimental procedures and precautions and signed informed consent forms before the commencement of the study.
Clinical Evaluation
In this study, the depression symptoms of patients with MDD were evaluated using the Hamilton Depression Rating Scale (HAMD-17) both before and after rTMS treatment (Hamilton, 1967). Higher HAMD-17 scores indicated more severe depression.
Repetitive Transcranial Magnetic Stimulation Procedure
The repetitive transcranial magnetic stimulator used in this study is a magnetic field stimulator from YIRYIDE Medical (MagPro R30, Dantec Medtronic, Denmark, CCY-IA), and the stimulation coil is a 100 mm figure of eight shaped coil. The treatment was administered at 10 Hz and 110% of the resting motor threshold (RMT) coil. The RMT is the minimum stimulation intensity that can elicit at least 5 motor evoked potential with an amplitude > 50 μV with 10 consecutive stimuli to the patient (Bashir et al., 2019). According to the international 10–20 system, the rTMS treatment target is located at the F3 point of the left dorsolateral prefrontal region, one pulse per second for 10 s, with 10 s interval, consisted of 1,000 pluses. rTMS treatment was provided for 15 min per day on 15 successive days.
Research Procedures
In the current study, we conducted a single-blind, randomized, controlled study at a general hospital. The randomization program was created using a computer and executed by an investigator who is not involved in the treatment and recruitment of patients. The allocation of patients was screened, applying numbered in the sealed and opaque envelopes. The antidepressant taken by the patients was venlafaxine (H32022135). The doses of venlafaxine prescribed were tailored based on clinical considerations determined by the physicians. The patients taken oral venlafaxine were initiated at a dose of 75 mg/d. Based on response and tolerability, the dose could be titrated upward to 150 mg/d at a 2-week interval. The patient with MDD used HAMD to assess the severity of depressive symptoms before and after the treatment rTMS treatment and completed the fMRI scan.
Imaging Data Acquire and Preprocessing
The patients underwent scanning within 48 h before the commencement of rTMS treatment and on the day following the end of the treatment course. The healthy controls were only scanned at baseline. Imaging data were acquired on a 3.0 Tesla MRI system with a standard 8-channel head coil (GE Medical Systems, Milwaukee, WI, United States). Functional images were acquired using a gradient echo-planar imaging (EPI) sequence (repetition time, 2,000 ms; echo time, 30 ms; FOV = 240 mm × 240 mm; FA = 90; matrix = 64 × 64; slice thickness, 4.5 mm; 45 axial slices no gap. A total of 210 volumes were collected for a total scan time of 420 s.
We performed image preprocessing using DPABI1 and SPM12.2 Briefly steps are as follow: (1) The first ten scan volumes were discarded for steady-state magnetization; (2) subsequent images were corrected for temporal differences by slice timing and head motion by alignment; (3) the resulting functional images were spatially normalized to the standard space of the Montreal Neurological Institute (MNI) using an optimum affine transformation and non-linear deformations, and then resampled to 3 mm × 3 mm × 3 mm isotropic voxels. Nuisance signals, including those from Fristo; (4) then all the functional images were smoothed with a 6-mm full-width at half-maximum (FWHM) Gaussian filter; (5) time series linear detrending was conducted to remove low-frequency drifts and high-frequency physiological noise.
Amplitude of the Low-Frequency Fluctuation Analysis
The whole-brain ALFF was calculated using the DPABI software. First, the time series of a given voxel was extracted; next, the amplitude of all frequencies in a frequency range (which was set as 0.01–0.1 Hz in this study) was calculated via Fourier transform; then, it was converted to the power spectrum for radication to determine the ALFF. The ALFF of each voxel was standardized (subtracted from the mean value of the whole-brain signal and then divided by the standard deviation) to obtain the standardized ALFF for each participant, which was used in subsequent statistical analyses.
Regional Homogeneity Analysis
By calculating Kendall’s coefficient of concordance of each voxel with a total of 26 voxel time series adjacent to the point, line, and plane of the given voxel, ReHo was determined and then divided by the mean value of the whole-brain signal to get the standardized ReHo; next, the standardized ReHo was subjected to Gaussian smoothing with the full width at a half-maximum of 6 × 6 × 6 mm3 and then used for subsequent statistical analyses.
Granger Causality Model
In this study, the left middle frontal gyrus with synchronous changes was used as the seed point, and dual-coefficient GCA was used to evaluate the causal connectivity between the time series describing the seed point and that of each voxel in the whole brain. GCA based on the whole-brain voxel level was used in the representational state transfer toolkit; the mean time series of the left middle frontal gyrus was defined as the seed-based time series x, whereas the time series of each voxel in the whole brain was represented by the time series y. A positive value of the causal connectivity from x to y indicated the existence of a causal connectivity between the activity in the middle frontal gyrus and that in the brain in the same direction, whereas a negative value indicated the opposite direction. The seed point was analyzed on the causal connectivity from the seed point to the other voxels in the whole brain (x to y) and that from the other voxels in the whole brain to the seed point (y to x). Afterward, the Fisher R-to-Z transformation was performed to obtain the voxel-level connectivity Z-map for subsequent statistical analyses.
Statistical Analysis
Statistical analysis of the clinical data was performed using SPSS 22.0. The paired-samples t-test was used to compare the ALFF and ReHo of patients with MDD before and after the treatment, whereas the independent-samples t-test was used to compare these neuroimaging indices of patients with MDD before and after the treatment with those of the healthy controls. The results of multiple comparisons were corrected using a Gaussian random field. The cluster level p < 0.01 (false discovery rate correction), voxel level p < 0.001, and voxel number > 70 were used to identify the brain region that differed significantly between groups. The imaging indices of different brain regions were extracted and their correlations with HAMD-17 scores were analyzed. p < 0.05 was considered statistically significant.
ResultsParticipant Demographics
A total of 29 patients with MDD, with a mean age of 28.44 ± 7.909 years and mean MDD duration of 4.05 ± 3.77 months were included in the treatment group, whereas a total of 33 healthy controls with a mean age of 26.53 ± 5.563 years were included in the control group. There were no significant differences between the two groups in terms of age, sex, and education level (p > 0.05; Table 1).
Group demographics and clinical measures (mean ± SD).
MDD (n = 29)
Healthy controls (n = 33)
p
Age (years)
28.44 ± 7.91
26.53 ± 5.56
0.52
Female/male
20/9
22/11
0.88
Education (years)
14.71 ± 2.01
15.18 ± 1.99
0.36
Duration of illness (months)
4.05 ± 3.77
-
Dose of venlafaxine (mg/d)
98.27 ± 22.16
-
HAMD-17 Score
The HAMD-17 total score was significantly higher in patients with MDD than in healthy controls, and the total score after treatment was significantly lower in patients with MDD than that before treatment in (t = 11.313, p < 0.001; Table 2).
Comparisons of the score of HAMD-17 between MDD and healthy controls (mean ± SD).
MDD
Healthy controls
t
Pre-rTMS
21.00 ± 5.382
2.85 ± 2.279
18.571***
Post-rTMS
9.52 ± 4.106
-
7.961***
*** means p < 0.001.
Amplitude of the Low-Frequency Fluctuation Analysis
Patients with MDD had lower ALFF in the left superior frontal, right inferior temporal, left middle occipital, and right inferior occipital gyri as well as higher ALFF in the left precuneus and left middle frontal gyrus than the values observed in healthy controls (Table 3 and Figure 1).
Brain regions in ALFF between MDD and healthy controls.
Brain regions
Brodmann areas
Side
MNI coordinates
Cluster size
t-values
X
Y
Z
Superior frontal gyrus
11/10
Left
−12
60
−21
133
−8.651
Inferior temporal gyrus
20
Right
69
−27
−27
76
−6.483
Middle occipital gyrus
17
Left
−27
−102
3
207
−7.827
Inferior occipital gyrus
18
Right
30
−99
−15
79
−7.421
Precuneus
40
Left
−42
−42
36
370
6.711
Middle frontal gyrus
8
Left
−36
12
60
313
6.619
(A) Brain regions showing significant differences of amplitude of low-frequency fluctuation (ALFF) between MDD and healthy controls. (B) Brain regions showing significant differences of ALFF between pre- and post-rTMS. The warm color denoted the region where ALFF is higher, and the cool color denotes the region where ALFF is lower.
After 15 days of rTMS treatment, the ALFF was increased in the left superior frontal gyrus and decreased in the left precuneus and left middle frontal gyrus of patients with MDD (Table 4 and Figure 1).
Brain regions in ALFF between post- and pre-rTMS.
Brain regions
Brodmann areas
Side
MNI coordinates
Cluster size
t-values
X
Y
Z
Superior frontal gyrus
11
Left
−12
57
−21
115
8.5434
Precuneus
23
Left
−6
−57
24
141
−5.859
Middle frontal gyrus
44
Left
−48
12
36
99
−5.747
Regional Homogeneity Analysis
Patients with MDD had higher resting-state ReHo in the right inferior temporal gyrus and lower ReHo in the left middle occipital, left middle frontal, and left postcentral gyri than the values obtained in healthy controls (Table 5 and Figure 2).
Brain regions in ReHo between MDD and healthy controls.
Brain regions
Brodmann areas
Side
MNI coordinates
Cluster size
t-values
X
Y
Z
Inferior temporal gyrus
20
Right
54
−18
−24
164
4.851
Middle occipital gyrus
17
Left
0
−90
−3
488
−7.469
Middle frontal gyrus
24
Left
−34
39
15
87
−5.282
Postcentrol gyrus
3
Left
−57
−21
45
62
−5.8024
(A) Brain regions showing significant differences of regional homogeneity (ReHo) between MDD and healthy controls. (B) Brain regions showing significant differences of ReHo between pre- and post-rTMS. The warm color denoted the region where ReHo is higher, and the cool color denotes the region where ReHo is lower.
After 15 days of rTMS treatment, ReHo was increased in the left middle occipital and frontal gyri (Table 6 and Figure 2).
Brain regions in ReHo between post- and pre-rTMS.
Brain regions
Brodmann areas
Side
MNI coordinates
Cluster size
t-values
X
Y
Z
Middle occipital gyrus
18
Left
−33
−96
0
74
4.189
Middle frontal gyrus
47
Left
−39
39
−3
57
4.237
Granger Causality Analysis Results
The ALFF and ReHo analyses identified the left middle frontal gyrus as a brain region with synchronous changes; therefore, it was used as the seed point for GCA. The results showed that the causal connectivity from the seed point to the left superior frontal gyrus, right middle temporal gyrus, and right caudate nucleus was significantly weak (p < 0.05; Table 7 and Figure 3). Moreover, the causal connectivity from the left superior frontal, left inferior temporal, and left middle occipital gyri to the seed point was significantly weak (p < 0.05; Table 7 and Figure 3). However, the causal connectivity from the right caudate nucleus to the seed point was significantly strong (p < 0.05; Table 7 and Figure 3).
Significant differences in Granger causality analysis between MDD and healthy controls.
Brain regions
Brodmann areas
Side
MNI coordinates
Cluster size
t-values
X
Y
Z
x to y
MDD VS. healthy controls
Superior frontal gyrus
46
Left
−30
18
39
345
−4.251
Middle temporal gyrus
21
Right
66
−24
0
95
−3.558
Caudate
48
Right
21
6
12
149
−3.964
y to x
MDD VS. healthy controls
Superior frontal gyrus
10
Left
−33
54
6
98
−3.258
Inferior temporal gyrus
20
Left
−48
−18
−30
145
−3.707
Middle occipital gyrus
19
Left
−33
−81
15
88
−3.284
Caudate
48
Right
24
0
27
98
3.514
Granger causality analysis between patients with MDD and healthy controls. (A) Significant different brain regions from seed point (LMFG) to whole brain regions; (B) significant different brain regions from whole brain regions to seed point (LMFG). Red areas show brain regions where patients with MDD had increased causal connectivity than healthy controls; Blue areas show brain regions where patients with MDD had decreased causal connectivity than healthy controls. The color bar represents t-values.
After rTMS treatment, the causal connectivity from the seed point to the superior frontal gyrus of patients with MDD was strengthened significantly (p < 0.05; Table 8 and Figure 4) and the causal connectivity from the left superior frontal and left inferior temporal gyri of patients with MDD to the seed point was also strengthened significantly (p < 0.05; Table 8 and Figure 4) compared with the causal connectivity before treatment.
Significant differences in Granger causality analysis between pre- and post- rTMS.
Brain regions
Brodmann areas
Side
MNI coordinates
Cluster size
t-values
X
Y
Z
x to y
Post-rTMS vs. pre- rTMS
Superior frontal gyrus
46
Left
−33
18
39
129
4.081
y to x
Post-rTMS vs. pre- rTMS
Superior frontal gyrus
10
Left
−24
60
9
76
3.678
Inferior temporal gyrus
48
Left
−30
6
−12
246
3.583
Granger causality analysis pre- and post- rTMS in MDD. (A) Brain regions showing group differences in causal connectivity from the seed point (LMFG) to whole brain in pre- and post-rTMS comparison; (B) Brain regions showing group differences in causal connectivity from the whole brain to seed point (LMFG) in pre- and post-rTMS comparison. Red areas show brain regions where post- rTMS patients with MDD had increased causal connectivity than pre- rTMS; Blue areas show brain regions where post- rTMS patients with MDD had decreased causal connectivity than pre-rTMS. The color bar represents t-values.
Correlation of Granger Causality Analysis Outcomes With HAMD-17 Score in Patients With Major Depressive Disorder After Repetitive Transcranial Magnetic Stimulation Treatment
The bidirectional causal connectivity between the left middle frontal gyri and left superior frontal gyri showed no significant negative correlation with the total HAMD-17 score (p > 0.05), whereas the causal connectivity from the left inferior temporal to the left middle frontal gyri showed a significant negative correlation with the total HAMD-17 score (r = −0.443, p = 0.021; Figure 5). The receiver operator characteristic (ROC) analysis was used to assess whether causal connectivity from inferior temporal gyrus to seed point can distinguish patients (pretreatment and posttreatment) and controls. The results showed that causal connectivity from inferior temporal gyrus to seed point successfully distinguished the patient group from the healthy control group and also pretreatment patients from the posttreatment patients (Figure 5).
Scatter plot between causal connectivity from left inferior temporal gyrus to left middle frontal gyrus and score of HAMD-17 in post- rTMS of patients with MDD (left); ROC curve for causal connectivity from inferior temporal gyrus to seed point with patients between patients (pretreatment and posttreatment) and control (right).
Discussion
In this study, the left dorsolateral prefrontal cortex was selected as the therapeutic target where the depression symptoms of patients with MDD were improved effectively. On the basis of ALFF and ReHo analyses, the ALFF and ReHo in the left superior frontal gyrus, left precuneus, and left middle frontal gyrus were significantly increased after rTMS treatment. GCA with the left middle frontal gyrus as the seed point revealed that after rTMS treatment, the bidirectional causal connectivity between the seed point and superior frontal gyrus was significantly improved and the causal connectivity from the inferior temporal gyrus to the seed point was also significantly improved and showed a remarkable negative correlation with the total HAMD-17 score. These findings demonstrate that the changes in the brain activity and causal connectivity in the frontal, parietal, and temporal regions may constitute the neural mechanisms underlying rTMS treatment in MDD.
As an important part of the prefrontal cortex (Li et al., 2013), the left superior frontal gyrus plays a vital role in self-awareness (Niendam et al., 2012) and negative emotion regulation (Chen et al., 2022), and it is involved in attention control and emotion regulation in patients with MDD (Dutta et al., 2014; Hamilton et al., 2015). In this study, we found that patients with MDD had significantly lower ALFF in the left superior frontal gyrus than healthy controls, which is consistent with the findings of previous studies (Liu et al., 2020; Zhang et al., 2021). The decreased spontaneous activities of brain neurons in the superior frontal gyrus may lead to low mood, anhedonia, and excessive negative self-evaluation in patients with MDD (Guo et al., 2013).
The left precuneus is a part of the superior parietal lobule (Zhu et al., 2018), which is involved in high-level cognitive functions such as episodic memory and self-reflection (Li et al., 2018). Studies have found that the abnormal neuronal activity in the precuneus of patients with MDD (Zhou et al., 2017) weakens the regulation of negative emotions (Zhong et al., 2016). The number of depressive episodes has a significant positive correlation with the ALFF in the precuneus (Jing et al., 2013). The present study demonstrated that patients with MDD had significantly higher ALFF in the left precuneus than healthy controls, indicating an excessive increase and abnormal activation of neuronal activities in this brain region of these patients; this finding is consistent with the that of previous studies (Guo et al., 2013; Yao et al., 2014; Zhang et al., 2021) and may be associated with a stronger negative self-awareness and sensory dysfunction (Gong et al., 2020; Wang et al., 2020).
The middle frontal gyrus is an important part of the ventromedial prefrontal cortex, which is involved in functions such as advanced cognition, autobiographical memory (van Heeringen et al., 2017), and emotion processing (Whalley et al., 2012) and cognition control dysfunction (Han et al., 2014). Evidence shows that an abnormal middle frontal gyrus leads to the dysregulation of top-down emotions and cognitive control disorder in patients with MDD (Lu et al., 2012). The ALFF in the middle frontal gyrus can not only specifically distinguish between the depressive episode and convalescence in MDD but also objectively reflect the severity of depression symptoms to a certain extent (Zhao et al., 2020). In the present study, the ALFF in the resting-state middle frontal gyrus was increased, further supporting the abnormal function of the middle frontal gyrus in patients with MDD. Abnormalities in functions such as abnormal emotional processing and cognitive memory may result from the compensatory increase in local neuronal activities in the middle frontal gyrus of patients with MDD (Guo et al., 2013). ReHo is also a sensitivity index of neuronal activity. The present study found that ReHo in the left middle frontal gyrus of patients with MDD was decreased, which is consistent with the finding of previous studies (Liu et al., 2021). The weakened synchronous and consistent activity of neurons in the middle frontal gyrus may lead to a decline in the ability to regulate emotions, thereby causing low mood with pessimistic and negative emotions in patients with MDD. In essence, abnormal ALFF and ReHo in the middle frontal gyrus reflect the disorders of local neuronal activity in patients with MDD in the resting state, which lead to the dysfunction of the middle frontal gyrus and dysregulation of emotional responses; these may represent important pathophysiological mechanisms of MDD (Geng et al., 2019).
High-frequency rTMS can regulate the excitability of the cerebral cortex (Xie et al., 2021). In this study, the left dorsolateral prefrontal cortex was targeted by high-frequency stimulation at 10 HZ, which directly enhanced the neural activity in the target-related brain region (Du et al., 2018). Therefore, an increase in the ALFF in the left superior frontal gyrus was observed in patients with MDD. The left superior frontal gyrus is the key brain region of the left dorsolateral prefrontal cortex. If it is damaged, the surrounding non-core organization brain regions compensate for the core brain region (Wang et al., 2019). The left middle frontal gyrus and left precuneus have a strong structural and functional connectivity with the left superior frontal gyrus (van den Heuvel and Sporns, 2011). If the function of the superior frontal gyrus is impaired, the middle frontal gyrus and precuneus may compensate for this impairment. Thus, an excessive increase in activity is observed in the left middle frontal gyrus and left precuneus. However, after treatment with high-frequency repetitive transcranial magnetic stimulation (HF-rTMS), the neural activity in the left superior frontal gyrus—the key brain region—was restored and the functions of the left middle frontal gyrus and left precuneus showed a tendency to be normal.
Our findings suggest that rTMS directly improves the neural activity in the target brain region stimulated indirectly improves the activity in local brain regions surrounding the target. The middle frontal gyrus is anatomically connected to the superior frontal gyrus through the uncinate fasciculus (Briggs et al., 2021), and this connection is involved in decision-making and emotional control. The functional connectivity between the middle and superior frontal gyri has been reported to be strengthened after antidepressant therapy (Porta-Casteràs et al., 2020). The findings of the present study are not only consistent with the findings of previous studies but also identify the direction of the causal connectivity between the middle and superior frontal gyri. The causal connectivity between the middle and superior frontal gyri was weak before treatment but strengthened after 15 days of continuous rTMS treatment; this suggests that the functional brain connectivity in patients with MDD tends to be normal and the negative thoughts can be improved, thus alleviating depression (Gröhn et al., 2019). Furthermore, the inferior temporal and middle frontal gyri are anatomically connected through the arcuate fasciculus, and these two also have a close functional connectivity, which is associated with emotional processing and episodic memory (Briggs et al., 2021). After rTMS treatment, the causal connectivity from the inferior temporal to the middle frontal gyri was strengthened and found to be negatively correlated with the total HAMD-17 score, indicating that the strengthening of such causal connectivity leads to the remission of depression symptoms.
Conclusion
In this study, the ALFF, ReHo, and GCA were used to explore the mechanisms through which HF-rTMS improves depression symptoms in patients with MDD. Our findings indicate that the rTMS treatment was efficacious for reducing depressive symptoms. And the rTMS treatment could regulate the neuronal activity in frontal gyrus of patients with MDD; It further enhanced the connectivity among the frontal—temporal—parietal regions of the brain after rTMS treatment. These changes may represent the neural mechanism through HF-rTMS treatment and may be used as a clinical biomarker for the objective evaluation of the therapeutic efficacy of rTMS.
Limitations
This study has some limitations. The sample size was small; this might have influenced the results and conclusion of this study. Moreover, Venlafaxine is the first of the SNRIs that provides dose-dependent norepinephrine reuptake inhibition; a dosage of 150 mg/day or higher is sufficient to produce noradrenergic activity, and it has low affinity for the postsynaptic receptors (Strawn et al., 2018; Zhou et al., 2021). A latest review has demonstrated that the brain function changes at least 4 weeks’ antidepressant pharmacotherapy in patients with MDD (Cattarinussi et al., 2021). In current study, venlafaxine, which is added from 75 to 150 mg/d within 15 days during rTMS treatment, may marginally affect brain function, but post-hoc test is difficult to control the confusion between drugs and rTMS treatment results. Future studies should ensure a larger sample size and enhance the experimental design to further discuss the therapeutic effects of rTMS in MDD.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
This study was reviewed and approved by the Medical Ethics Committee of Xijing Hospital (Approval document Number: KY20202055-F-1). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
MG, HW, and QT conceived and designed the experiments. MG, YS, ZM, YX, and ZW performed the experiments. MG, ZL, JL, and XG analyzed the data. HW and QT conceived the project and modified the manuscript. All authors read and approved the final manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Funding
This work was supported by the National Nature Science Foundation of China (No. 81974215 to HW), the National Nature Science Founda tion of China (No. 81630032 to QT), and the Project of Natural Science of Shaanxi Science and Technology Department (No. 2021JQ-783).
ReferencesAntonenkoD.ThielscherA.SaturninoG. B.AydinS.IttermannB.GrittnerU. (2019). Towards precise brain stimulation: is electric field simulation related to neuromodulation?Brain Stimul.121159–1168. 10.1016/j.brs.2019.03.07230930209BashirS.AhmadS.AlatefiM.HamzaA.SharafM.FecteauS. (2019). Effects of anodal transcranial direct current stimulation on motor evoked potentials variability in humans.Physiol. Rep.7:e14087. 10.14814/phy2.1408731301123BerlimM. T.van den, EyndeF.Tovar-PerdomoS.DaskalakisZ. J. (2014). Response, remission and drop-out rates following high-frequency repetitive transcranial magnetic stimulation (rTMS) for treating major depression: a systematic review and meta-analysis of randomized, double-blind and sham-controlled trials.Psychol. Med.44225–239. 10.1017/S003329171300051223507264BriggsR. G.LinY.-H.DadarioN. B.KimS. J.YoungI. M.BaiM. Y. (2021). Anatomy and white matter connections of the middle frontal gyrus.World Neurosurg.150e520–e529. 10.1016/j.wneu.2021.03.04533744423CattarinussiG.DelvecchioG.MaggioniE.BressiC.BrambillaP. (2021). Ultra-high field imaging in major depressive disorder: a review of structural and functional studies. J. Affect Disord. 290, 65–73. 10.1016/j.jad.2021.04.056ChenZ.ZhaoS.TianS.YanR.WangH.WangX. (2022). Diurnal mood variation symptoms in major depressive disorder associated with evening chronotype: evidence from a neuroimaging study.J. Affect. Disord.298151–159. 10.1016/j.jad.2021.10.08734715183ChengW.RollsE. T.RuanH.FengJ. (2018). Functional connectivities in the brain that mediate the association between depressive problems and sleep quality.JAMA Psychiatry751052–1061. 10.1001/jamapsychiatry.2018.194130046833CohenS. E.ZantvoordJ. B.WezenbergB. N.BocktingC. L. H.Van WingenG. A. (2021). Magnetic resonance imaging for individual prediction of treatment response in major depressive disorder: a systematic review and meta-analysis.Transl. Psychiat.111–10. 10.1038/s41398-021-01286-x33723229DichterG. S.GibbsD.SmoskiM. J. (2015). A systematic review of relations between resting-state functional-MRI and treatment response in major depressive disorder.J. Affect. Disord.1728–17. 10.1016/j.jad.2014.09.02825451389DuL.LiuH.DuW.ChaoF.ZhangL.WangK. (2018). Stimulated left DLPFC-nucleus accumbens functional connectivity predicts the anti-depression and anti-anxiety effects of rTMS for depression.Transl. Psychiat.71–10. 10.1038/s41398-017-0005-629520002DuttaA.MckieS.DeakinJ. F. W. (2014). Resting state networks in major depressive disorder.Psychiatry Res.224139–151. 10.1016/j.pscychresns.2014.10.003FengZ.XuS.HuangM.ShiY.XiongB.YangH. (2016). Disrupted causal connectivity anchored on the anterior cingulate cortex in first-episode medication-naive major depressive disorder.Prog. Neuro-Psychoph.64124–130. 10.1016/j.pnpbp.2015.07.00826234517GaynesB. N.LloydS. W.LuxL.GartlehnerG.HansenR. A.BrodeS. (2014). Repetitive transcranial magnetic stimulation for treatment-resistant depression: a systematic review and meta-analysis.J. Clin. Psychiatry75477–489.GengJ.YanR.ShiJ.ChenY.MoZ.ShaoJ. (2019). Altered regional homogeneity in patients with somatic depression: A resting-state fMRI study.J. Affect Disord.246498–505. 10.1016/j.jad.2018.12.06630599374Gex-FabryM.Balant-GorgiaA. E.BalantL. P.RudazS.VeutheyJ.-L.BertschyG. (2004). Time course of clinical response to venlafaxine: relevance of plasma level and chirality.Eur. J. Clin. Pharmacol.59883–891. 10.1007/s00228-003-0710-314704834GollanJ. K.ConnollyM.BuchananA.HoxhaD.RosebrockL.CacioppoJ. (2015). Neural substrates of negativity bias in women with and without major depression.Biol. Psychol.109184–191. 10.1016/j.biopsycho.2015.06.00326073417GongJ.WangJ.QiuS.ChenP.LuoZ.WangJ. (2020). Common and distinct patterns of intrinsic brain activity alterations in major depression and bipolar disorder: voxel-based meta-analysis.Transl. Psychiat.101–13. 10.1038/s41398-020-01036-533077728GröhnH.GillickB. T.TkáčI.BednaříkP.MascaliD.DeelchandD. K. (2019). Influence of repetitive transcranial magnetic stimulation on human neurochemistry and functional connectivity: a Pilot MRI/MRS study at 7 T.Front. Neurosci.13:1260. 10.3389/fnins.2019.0126031827419GuoW.-B.LiuF.XunG.-L.HuM.-R.GuoX.-F.XiaoC.-Q. (2013). Reversal alterations of amplitude of low-frequency fluctuations in early and late onset, first-episode, drug-naive depression.Prog. Neuro-Psychoph.40153–159.HamiltonJ. P.FarmerM.FogelmanP.GotlibI. H. (2015). Depressive rumination, the default-mode network, and the dark matter of clinical neuroscience.Biol. Psychiatry78224–230. 10.1016/j.biopsych.2015.02.02025861700HamiltonM. A. X. (1967). Development of a rating scale for primary depressive illness.Br. J. Soc. Clin. Psychol.6278–296. 10.1111/j.2044-8260.1967.tb00530.xHanK.-M.ChoiS.JungJ.NaK.-S.YoonH.-K.LeeM.-S. (2014). Cortical thickness, cortical and subcortical volume, and white matter integrity in patients with their first episode of major depression.J. Affect. Disord.15542–48. 10.1016/j.jad.2013.10.02124210630IshidaT.DierksT.StrikW.MorishimaY. (2020). Converging Resting State Networks Unravels Potential Remote Effects of Transcranial Magnetic Stimulation for Major Depression.Front. Psychia.111–13. 10.3389/fpsyt.2020.0083632973580JingB.LiuC.-H.MaX.YanH.-G.ZhuoZ.-Z.ZhangY. (2013). Difference in amplitude of low-frequency fluctuation between currently depressed and remitted females with major depressive disorder.Brain Res.154074–83. 10.1016/j.brainres.2013.09.03924121137KaiserR. H.Andrews-HannaJ. R.WagerT. D.PizzagalliD. A. (2015). Large-scale network dysfunction in major depressive disorder: a meta-analysis of resting-state functional connectivity.JAMA Psychiatry72603–611. 10.1001/jamapsychiatry.2015.007125785575KongX.-M.XuS.-X.SunY.WangK.-Y.WangC.ZhangJ. (2017). Electroconvulsive therapy changes the regional resting state function measured by regional homogeneity (ReHo) and amplitude of low frequency fluctuations (ALFF) in elderly major depressive disorder patients: an exploratory study.Psychiat. Res. Neuroim.26413–21. 10.1016/j.pscychresns.2017.04.00128412557LiG.RossbachK.ZhangA.LiuP.ZhangK. (2018). Resting-state functional changes in the precuneus within first-episode drug-naive patients with MDD.Neuropsych. Dis. Treat.14:1991. 10.2147/NDT.S16806030122932LiW.QinW.LiuH.FanL.WangJ.JiangT. (2013). Subregions of the human superior frontal gyrus and their connections.Neuroimage7846–58. 10.1016/j.neuroimage.2013.04.01123587692LiangS.XueK.WangW.YuW.MaX.LuoS. (2020). Altered brain function and clinical features in patients with first-episode, drug naïve major depressive disorder: A resting-state fMRI study.Psychiat. Res. Neuroim.303:111134. 10.1016/j.pscychresns.2020.11113432652482LiuP.TuH.ZhangA.YangC.LiuZ.LeiL. (2021). Brain functional alterations in MDD patients with somatic symptoms: A resting-state fMRI study.J. Affect. Disord.295788–796. 10.1016/j.jad.2021.08.14334517253LiuX.HouZ.YinY.XieC.ZhangH.ZhangH. (2020). Dopamine multilocus genetic profile, spontaneous activity of left superior temporal gyrus, and early therapeutic effect in major depressive disorder.Front. Psychiatry11:1432. 10.3389/fpsyt.2020.59140733414733LoerincA. G.MeuretA. E.TwohigM. P.RosenfieldD.BluettE. J.CraskeM. G. (2015). Response rates for CBT for anxiety disorders: Need for standardized criteria.Clin. Psychol. Rev.4272–82. 10.1016/j.cpr.2015.08.00426319194LuQ.LiH.LuoG.WangY.TangH.HanL. (2012). Impaired prefrontal–amygdala effective connectivity is responsible for the dysfunction of emotion process in major depressive disorder: a dynamic causal modeling study on MEG.Neurosci. Lett.523125–130. 10.1016/j.neulet.2012.06.05822750155ManeetonB.ManeetonN.WoottilukP.LikhitsathianS. (2020). Repetitive transcranial magnetic stimulation combined with anti-depressants for the first episode of major depressive disorder.Curr. Neuropharmacol.18852–860. 10.2174/1570159X1866620022111313432091338McGirrA.Van Den EyndeF.Tovar-PerdomoS.FleckM. P. A.BerlimM. T. (2015). Effectiveness and acceptability of accelerated repetitive transcranial magnetic stimulation (rTMS) for treatment-resistant major depressive disorder: an open label trial.J. Affect. Disord.173216–220. 10.1016/j.jad.2014.10.06825462419NiendamT. A.LairdA. R.RayK. L.DeanY. M.GlahnD. C.CarterC. S. (2012). Meta-analytic evidence for a superordinate cognitive control network subserving diverse executive functions.Cogn. Affect Behav. Neurosci.12241–268.OtteC.GoldS. M.PenninxB. W.ParianteC. M.EtkinA.FavaM. (2016). Major depressive disorder.Nat. Rev. Dis. Prim.21–20. 10.1093/med/9780198810940.003.0001Porta-CasteràsD.CanoM.CamprodonJ. A.LooC.PalaoD.Soriano-MasC. (2020). A multimetric systematic review of fMRI findings in patients with MDD receiving ECT.Prog. Neuro-Psychoph2020:110178. 10.1016/j.pnpbp.2020.11017833197507RosenbaumD.Int-VeenI.KroczekA.HilsendegenP.Velten-SchurianK.BihlmaierI. (2020). Amplitude of low frequency fluctuations (ALFF) of spontaneous and induced rumination in major depression: an fNIRS study.Sci. Rep. UK101–13. 10.1038/s41598-020-78317-y33299001StrawnJ. R.GeraciotiL.RajdevN.ClemenzaK.LevineA. (2018). Pharmacotherapy for generalized anxiety disorder in adult and pediatric patients: an evidencebased treatment review.Expert Opin. Pharmacother.191057–1070. 10.1080/14656566.2018.1491966TadayonnejadR.YangS. L.KumarA.AjiloreO. (2015). Clinical, cognitive, and functional connectivity correlations of resting-state intrinsic brain activity alterations in unmedicated depression.J. Affect. Disord.172241–250. 10.1016/j.jad.2014.10.01725451423TateishiH.MizoguchiY.KawaguchiA.ImamuraY.MatsushimaJ.KunitakeH. (2020). Changes in interleukin-1 beta induced by rTMS are significantly correlated with partial improvement of cognitive dysfunction in treatment-resistant depression: a pilot study.Psychiat. Res.289:112995. 10.1016/j.psychres.2020.11299532371273TianL.SunS. S.CuiL. B.WangS. Q.PengZ. W.TanQ. R. (2020). Repetitive transcranial magnetic stimulation elicits antidepressant-and anxiolytic-like effect via nuclear factor-E2-related factor 2-mediated anti-inflammation mechanism in rats.Neuroscience429119–133. 10.1016/j.neuroscience.2019.12.02531918011van den HeuvelM. P.SpornsO. (2011). Rich-club organization of the human connectome.J. Neurosci.3115775–15786. 10.1523/JNEUROSCI.3539-11.201122049421van HeeringenK.WuG.-R.VervaetM.VanderhasseltM.-A.BaekenC. (2017). Decreased resting state metabolic activity in frontopolar and parietal brain regions is associated with suicide plans in depressed individuals.J. Psychiatric. Res.84243–248. 10.1016/j.jpsychires.2016.10.01127771590WangM.JuY.LuX.SunJ.DongQ.LiuJ. (2020). Longitudinal changes of amplitude of low-frequency fluctuations in MDD patients: A 6-month follow-up resting-state functional magnetic resonance imaging study.J. Affect. Disord.276411–417. 10.1016/j.jad.2020.07.06732871671WangX.QinJ.ZhuJ.BiK.ZhangS.YanR. (2019). Rehabilitative compensatory mechanism of hierarchical subnetworks in major depressive disorder: a longitudinal study across multi-sites.Eur. Psychiat.5854–62. 10.1016/j.eurpsy.2019.02.00430822739WhalleyM. G.RuggM. D.BrewinC. R. (2012). Autobiographical memory in depression: an fMRI study.Psychiatry Res.20198–106.XiaM.SiT.SunX.MaQ.LiuB.WangL. (2019). Reproducibility of functional brain alterations in major depressive disorder: Evidence from a multisite resting-state functional MRI study with 1,434 individuals.Neuroimage189700–714. 10.1016/j.neuroimage.2019.01.07430716456XieY.GuanM.WangZ.MaZ.WangH.FangP. (2021). rTMS induces brain functional and structural alternations in schizophrenia patient with auditory verbal hallucination.Front. Neurosci.2021:15. 10.3389/fnins.2021.72289434539338YangL.WeiA. H.OuyangT. T.CaoZ. Z.DuanA. W.ZhangH. H. (2021). Functional plasticity abnormalities over the lifespan of first-episode patients with major depressive disorder: a resting state fMRI study.Ann. Transl. Med.9:349. 10.21037/atm-21-36733708976YaoZ.YanR.WeiM.TangH.QinJ.LuQ. (2014). Gender differences in brain activity and the relationship between brain activity and differences in prevalence rates between male and female major depressive disorder patients: a resting-state fMRI study.Clin. Neurophysiol.1252232–2239. 10.1016/j.clinph.2014.03.00624746685YinM.LiuY.ZhangL.ZhengH.PengL.AiY. (2020). Effects of rTMS Treatment on cognitive impairment and resting-state brain activity in stroke patients: a randomized clinical trial.Front. Neural. Circuits14:62. 10.3389/fncir.2020.56377733117131ZhangB.QiS.LiuS.LiuX.WeiX.MingD. (2021). Altered spontaneous neural activity in the precuneus, middle and superior frontal gyri, and hippocampus in college students with subclinical depression.BMC Psychiatry211–10. 10.1186/s12888-021-03292-134074266ZhaoP.YanR.WangX.GengJ.ChattunM. R.WangQ. (2020). Reduced resting state neural activity in the right orbital part of middle frontal gyrus in anxious depression.Front. Psychiatry10:994. 10.3389/fpsyt.2019.0099432038329ZhengA.YuR.DuW.LiuH.ZhangZ.XuZ. (2020). Two-week rTMS-induced neuroimaging changes measured with fMRI in depression.J. Affect. Disord.27015–21. 10.1016/j.jad.2020.03.03832275215ZhongX.PuW.YaoS. (2016). Functional alterations of fronto-limbic circuit and default mode network systems in first-episode, drug-naïve patients with major depressive disorder: a meta-analysis of resting-state fMRI data.J. Affect. Disord.206280–286. 10.1016/j.jad.2016.09.00527639862ZhouM.HuX.LuL.ZhangL.ChenL.GongQ. (2017). Intrinsic cerebral activity at resting state in adults with major depressive disorder: a meta-analysis.Prog. Neuropsychopharmacol. Biol. Psychiatry75157–164. 10.1016/j.pnpbp.2017.02.00128174129ZhouJ. J.WangX.FengL.XiaoL.YangR.ZhuX. Q. (2021). Venlafaxine vs. fluoxetine in postmenopausal women with major depressive disorder: an 8-week, randomized, single-blind, active-controlled study.BMC. psychia.211–10. 10.1186/s12888-021-03253-8ZhuJ.LinX.LinC.ZhuoC.YuY. (2018). Selective functional dysconnectivity of the dorsal-anterior subregion of the precuneus in drug-naive major depressive disorder.J. Affect. Disord.225676–683. 10.1016/j.jad.2017.08.08428917194