Key PointsQuestion
Is repetitive transcranial magnetic stimulation associated with changes in subgenual cingulate cortex (SGC) activity in patients with major depressive disorder (MDD)?
Findings
This diagnostic study, which compared 30 patients with MDD and 30 healthy controls, found that using transcranial magnetic stimulation combined with electroencephalography, SGC activity in patients with MDD was significantly higher compared with healthy controls. After active repetitive transcranial magnetic stimulation, SGC hyperactivity in patients with MDD was attenuated toward the levels of healthy controls.
Meaning
These SGC-localized findings support SGC hyperactivity as a central construct in the pathophysiology of MDD, which future work might develop into a clinically significant biological target.
Importance
Hyperactivity in the subgenual cingulate cortex (SGC) is associated with major depressive disorder (MDD) and anticorrelated with activity in the dorsolateral prefrontal cortex (DLPFC). This association was found to be predictive of responsiveness to repetitive transcranial magnetic stimulation (rTMS) treatment. Such findings suggest that DLPFC-SGC connectivity is important for understanding both the therapeutic mechanism of rTMS in patients with MDD and the underlying pathophysiology of MDD.
Objective
To evaluate SGC hyperactivity in patients with MDD before and after rTMS treatment.
Design, Setting, and Participants
In this diagnostic study, among participants recruited from the adult and geriatric mood and anxiety services at the Centre for Addiction and Mental Health, Toronto, Ontario, Canada, who had participated in a randomized clinical trial, baseline SGC activity of patients with MDD was compared with healthy controls. In patients with MDD, SGC activity was compared before and after active or sham high-frequency rTMS treatment. Data collection started in July 2008 and concluded in March 2012. Neurophysiological data analysis started in January 2017 and ended in May 2018.
Main Outcomes and Measures
Hyperactivity in the SGC before and after rTMS treatment was measured. Subgenual cingulate cortex hyperactivity activity was quantified using significant current density (SCD), and effective connectivity between the left DLPFC and SGC was computed using significant current scattering (SCS). Both measures were computed around TMS evoked potentials standard peak latencies prior to rTMS and after rTMS treatment, comparing patients with MMD treated with active and sham rTMS. Patients with MDD were assessed with the 17-item Hamilton Rating Scale for Depression.
Results
Of 121 patients with MDD in the initial trial, 30 (15 [50.0%] women) were compared with 30 healthy controls (15 [50.0%] women) at rTMS treatment baseline. The mean (SD) age of the cohort with MDD was 39.1 (10.9) years, and the mean (SD) age of healthy controls was 37.0 (11.0) years. Following rTMS treatment, 26 patients with MDD who had active rTMS treatment (21.5%) were compared with 17 patients with MDD who had sham treatment (14.0%). At baseline, the SGC mean (SD) SCD and mean (SD) SCS at 200 milliseconds after TMS pulse were higher in participants with MDD compared with healthy controls (SCD: 1.04 × 10−6 [1.41 × 10−6] μA/mm2 vs 3.8 × 10−7 [7.8 × 10−7] μA/mm2; z = –2.95; P = .004; SCS: 0.87 [0.86] mm vs 0.54 [0.87] mm; z = –2.27; P = .02). Baseline source current density was able to classify MDD with 77% accuracy. Scores on the 17-item Hamilton Rating Scale for Depression were correlated with current density at the SGC (ρ = 0.41; P = .03). After rTMS treatment, SGC mean (SD) SCD and mean (SD) SCS at 200 milliseconds after rTMS pulse were attenuated to approximately the standard TMS-evoked potential latencies in the active rTMS group compared with the sham rTMS group (SCD: 1.57 × 10−7 [3.67 × 10−7] μA/mm2 vs 7.00 × 10−7 [7.51 × 10−7] μA/mm2; z = −2.91; P = .004; SCS: 0.20 [0.44] mm vs 0.74 [0.73] mm; z = −2.78; P = .006). Additionally, the SGC SCS change was correlated with symptom improvement on the 17-item Hamilton Rating Scale for Depression in the active rTMS group (ρ = 0.58; P = .047).
Conclusions and Relevance
The findings of this study further implicate left DLPFC-SGC effective connectivity and SGC excitability in the pathophysiology of MDD and treatment with rTMS. These findings suggest that DLPFC-SGC connectivity may be a marker of rTMS treatment responsiveness.
Trial Registration
ClinicalTrials.gov identifier: NCT01515215
Major depressive disorder (MDD) is a debilitating psychiatric condition with a 16% lifelong prevalence.1 Antidepressant medication is widely used and has been researched for more than 4 decades. However, the mechanism of antidepressant action remains unclear, and nearly 40% of patients with MDD experience persistent depression even after 2 antidepressant treatments.2 To improve outcomes, the neurophysiology of MDD and underlying therapeutic mechanisms of treatments need to be better understood.
Hyperactivity of the subgenual cingulate cortex (SGC) is associated with MDD pathophysiology. Based on positron emission tomography studies,3 the SGC was shown to be overactive in the depressed state and during transient sadness in healthy controls. Furthermore, higher SGC activity was also found in patients with MDD in both depressed and remitted states.3 The SGC has also been a primary target for deep brain stimulation in the treatment of MDD.4
Evidence also implicates dorsolateral prefrontal cortex (DLPFC)–SGC connectivity in MDD pathophysiology. An inverse association of the SGC with right DLPFC activity was found in patients with MDD and healthy patients,3 implying that there is a functional association of the SGC with the DLPFC. Moreover, functional magnetic resonance imaging showed that DLPFC-SGC activity is anticorrelated,5 and the magnitude of this anticorrelation can predict the antidepressant efficacy of repetitive transcranial magnetic stimulation (rTMS) targeting the DLPFC. Finally, both positron emission tomography and functional magnetic resonance imaging studies have shown that baseline SGC activation abnormalities were attenuated toward healthy levels using a variety of antidepressant treatments, including conventional antidepressant medications,6-9 rTMS,10,11 electroconvulsive therapy,12 and deep brain stimulation.4
Combining TMS with electroencephalography (TMS-EEG) is a noninvasive technique that is used to assess neuronal activity,13,14 connectivity,15,16 and plasticity.17 Compared with other neurophysiological methods, TMS-EEG enables a more reliable and causal assessment of the neurophysiology of brain function. This is owing to the fact that TMS-EEG stimulates the cortex directly without dependence on prior activation of upstream brain functional regions (eg, sensory constructs).16 Significant current density (SCD) is a standardization measure that sums all of the significant TMS-responsive current sources (compared with baseline). Significant current scattering (SCS) is a standardization measure for inferring activation propagation; SCS may also measure the effective connectivity of a source relative to the stimulation site.13,18 Previously, SCD has been computed to differentiate between patients with Alzheimer disease and healthy controls.19 The 2 measures have also been used to capture changes in brain network connectivity during loss of consciousness,18 saccadic movement,20 and task performance.21 Significant current scattering has also been effective in distinguishing patients with schizophrenia from healthy controls.22
We used TMS-EEG measures of SCD to evaluate SGC excitability and TMS-EEG measures of SCS to evaluate DLPFC-SGC effective connectivity in patients with MDD. These measures were taken at baseline and after rTMS treatment. We hypothesized that SCD would demonstrate higher excitability in the SGC of patients with MDD. We also hypothesized that SCS would show a stronger effective connectivity between the left DLPFC and the SGC in patients with MDD. Finally, we hypothesized that increased excitability in (ie, SCD) and connectivity to (ie, SCS) the SGC would be attenuated after applying rTMS over the DLPFC in patients with MDD.
Recruitment and Treatment
Overall, 121 participants with MDD (77 [63.6%] women) were recruited at the Centre for Addiction and Mental Health in Toronto, Ontario, Canada. These patients participated in a randomized clinical trial, as previously described.23 The demographic and clinical characteristics are described in the eTable in the Supplement. The adverse effects of rTMS on study participants have been previously described.23 A random subset of participants in the randomized clinical trial (n = 30) underwent a pretreatment neurophysiological assessment. We recruited an additional 30 healthy controls separately for a neurophysiological measures comparison. Of 101 participants who completed week 6 of the trial, 43 underwent posttreatment neurophysiological assessments, had usable EEG recordings, and were analyzed for TMS-EEG neurophysiological effects. The participants with MDD were randomized to 3 different treatment arms: unilateral rTMS stimulation (n = 40), bilateral rTMS stimulation (n = 40), and sham stimulation (n = 41) (Figure 1). The rTMS protocol and parameters were conducted per the article by Blumberger et al.23 Neurophysiological measurements by TMS-EEG were performed 1 week prior to rTMS treatment and within 48 hours of the last rTMS treatment. The 17-item Hamilton Rating Scale for Depression (HRSD-17)24 was given to patients before and after rTMS treatment. Participants gave written informed consent, and the protocol was approved by the Centre for Addiction and Mental Health ethics board in accordance with Declaration of Helsinki.25 This report follows the Consolidated Standards of Reporting Trials (CONSORT) reporting guideline.
Acquisition and Preprocessing of TMS-EEG Data
In all groups, a figure-of-eight coil connected to Magstim 200 stimulators (Magstim) was used to stimulate the brain over the left DLPFC. Locating the motor cortex was done by eliciting motor-evoked potentials at the right abductor pollicis brevis. The left DLPFC was identified using the miniBIRD neuronavigation system (Ascension Technologies). Before each experiment, resting motor threshold was determined, as previously described.26 Each patient’s stimulus intensity was determined as the intensity eliciting peak-to-peak amplitude of 1 mV averaged over 20 trials. Electroencephalography was performed using a 64-channel Synamps 2 EEG system (Compumedics Neuroscan). It was recorded applying a 200-Hz low-pass filter at a 20-kHz sampling rate. All electrodes (silver/silver chloride ring) impedances were kept less than 5 kΩ throughout the session, and the ground electrode was positioned posteriorly to the Cz electrode. Electroencephalography analysis was performed using MATLAB version r2017b (The MathWorks) EEGLAB27 and FieldTrip28 toolboxes. Noisy channels were removed (with high-powered 60-Hz amplitudes or with extremely high and variable amplitudes). The EEG recording was epoched around the TMS pulse (1 second before and after). Transcranial magnetic stimulation pulses were removed from each epoch (−2 milliseconds to 20 milliseconds) and linearly interpolated. Noisy epochs were removed. The data were baseline-corrected (500 milliseconds to 200 milliseconds prepulse) and average-referenced. Data was resampled to 1 kHz. A first round of independent component analysis was used to remove large muscle artifacts. The data were bandpass filtered (1-100 Hz), with a 58- to 62-Hz notch. A second round of independent component analysis was used to remove eye blinks, eye movements, and additional muscle artifacts. Finally, missing electrodes were interpolated. Further details regarding our TMS-EEG preprocessing can be found in the article by Rogasch et al.29
Source Localization Procedure
Source analysis of the TMS-evoked potential (TEP) was done using the MATLAB Brainstorm toolbox.30 A generic brain with 15 002 voxels, based on the Montreal Neurological Institute International Consortium for Brain Mapping 152–averaged magnetic resonance imaging for the extracted cortex surface,31 was used. The EEG cap used in the experiment, Neuroscan 64-channel quik-cap (Compumedics Neuroscan), was coregistered to the generic head model. The forward model was computed using the OpenMEEG approach with a solution space limited to the cortex surface.32 The prestimulus period of individual trials was used to calculate the noise covariance. Finally, the inverse solution was computed based on the standardized low-resolution brain electromagnetic tomography algorithm,33 with dipoles constrained normally to the cortex surface. For each patient, the source localization procedure generated a 15 002-voxel current density map in brain space for every point of the TEP.
Significant current scatter was calculated based on the following equation adapted from methods previously published in the article by Casali et al13: SCS = SS(x, t) × d(x – xstim), in which SS(x, t) is a binary matrix of significant sources across the brain at each point of the TEP and d(x – xstim) is calculated as the distance of every source voxel from site of stimulation (ie, F5). To calculate SS(x, t), each poststimulus point trial distribution was compared by a paired t test with a surrogate distribution taken from a point in a prestimulus time segment. A current dipole at a specific point was considered significant if its 2-sided α was less than .05 in a paired t test compared with the surrogate responses at that time point.
To compare the overall scatter values from the TMS target (the DLPFC) to the SGC, SCS magnitudes were summed across the TEP peak standard time periods (approximately 30 milliseconds, 60 milliseconds, 100 milliseconds, and 200 milliseconds) for the right and left SGC regions as defined by the Destrieux atlas.34
Both bilateral and unilateral rTMS treatment groups were pooled together into the active rTMS treatment group owing to statistical power considerations. To avoid normality assumptions, the Wilcoxon rank sum test was used to examine the differences between current density, SCD, and SCS in patients with MDD and healthy controls and in patients with MDD in the different rTMS treatment arms. Statistical significance was set at P < .05, and all tests were 2-tailed. No correction for multiple comparisons was applied. A Spearman correlation was used to quantify the association of SCS change with HRSD-17 score change before and after treatment. Finally, a receiver operating characteristic analysis was computed for the pre-rTMS sample; a source current density amplitude window from 15 milliseconds to 350 milliseconds was used as the criterion to classify healthy controls vs patients with MDD. This time segment was chosen because significant sources dropped drastically 350 milliseconds after stimulation. The probability of correct prediction was quantified by the area under the receiver operating characteristic curve,35 while the optimal threshold was determined as the source current density value associated with the maximum Youden index, or height above the diagonal line of no discrimination.36 All statistical analysis was done using MATLAB version r2017b.
Overall, 30 of 121 trial participants with MDD (15 [50.0%] women) pre-rTMS treatment were compared with 30 healthy controls (15 [50.0%] women). The mean (SD) age of the cohort with MDD was 39.1 (10.9) years, and the mean (SD) HRSD-17 score was 24.8 (3.5). The mean (SD) age of the healthy controls was 37.0 (11.0) years (eTable in the Supplement). There was no age difference between the active and sham rTMS treatment groups (t = −0.47; P = .46).
Overall, 26 patients with MDD (21.5%; 17 [65.4%] women) who received active rTMS treatment were analyzed. They had a mean (SD) age of 47.3 (14.0) years and mean (SD) HRSD-17 score of 25.6 (2.9). We analyzed 17 patients with MDD (14.0%; 8 [47.1%] women) who received sham rTMS. They had a mean (SD) age of 45.7 (12.5) years and a mean (SD) HRSD-17 score of 25.3 (3.2) (eTable in the Supplement). There were no age differences between the active and sham rTMS treatment groups (t = 0.17; P = .70). Patients with MDD demonstrated significant differences on SCD and SCS computations compared with healthy controls (Figure 2A). These differences were localized in voxels in the region of the SGC, and their timings aligned with standard TEP temporal components (P30 milliseconds, N100 milliseconds, and P200 milliseconds). Mean (SD) current density in the SGC was higher at 30 milliseconds for patients with MDD compared with healthy controls (1.51 × 10−7 [7.18 × 10−8] μA/mm2 vs 9.42 × 10−8 [4.23 × 10−8] μA/mm2; z = −3.42; P < .001), at 100 milliseconds (2.05 × 10−7 [1.28 × 10−7] μA/mm2 vs 1.02 × 10−7 [4.7 × 10−8] μA/mm2; z = –4.13; P < .001), and at 200 milliseconds (2.1 × 10−7 [1.37 × 10−7] μA/mm2 vs 1.01 × 10−7 [7.45 × 10−8] μA/mm2; z = –4.26; P < .001) after the TMS pulse. Mean (SD) SCD in the SGC in patients with MDD compared with healthy controls was higher at 30 milliseconds (4.26 × 10−7 [5.04 × 10−7] μA/mm2 vs 1.12 × 10−7 [1.67 × 10−7] μA/mm2; z = –2.57; P = .01), at 100 milliseconds (9.43 × 10−7 [1.41 × 10−6] μA/mm2 vs 2.83 × 10−7 [4.47 × 10−7] μA/mm2; z = –2.55; P = .01), and at 200 milliseconds (1.04 × 10−6 [1.41 × 10−6] μA/mm2 vs 3.8 × 10−7 [7.8 × 10−7] μA/mm2; z = –2.95; P = .004) after the pulse. Additionally, the mean (SD) SCS between the stimulation site (left DLPFC) and SGC was higher in patients with MDD compared with healthy controls at 100 milliseconds (0.93 [0.99] mm vs 0.47 [0.66] mm; z = –1.97; P = .048) and at 200 milliseconds (0.87 [0.86] mm vs 0.54 [0.87] mm; z = –2.27; P = .02) after the TMS pulse. A receiver operating characteristic curve analysis was used to identify the predictive power of source current density in patients with MDD compared with healthy controls. The model differentiated patients with MDD from the healthy controls with 77% accuracy (70% sensitivity and 83% specificity) (Figure 3).
We then compared patients with MDD after active vs sham rTMS treatment (Figure 2B). Differences in mean (SD) current density between unilateral and bilateral rTMS treatment groups were not statistically significant (1.77 × 10–7 [4.63 × 10–7] μA/mm2 vs 1.38 × 10–7 [2.43 × 10–7] μA/mm2; z = −1.00; P = .31). Therefore, we pooled these treatment groups to achieve better averaged neurophysiological signals and greater statistical power. Marked differences of current density, SCD, and SCS were found between the 2 experimental groups and were also localized in SGC-related voxels. These timing differences were also associated with the known TEP components (P60 milliseconds and P200 milliseconds). Mean (SD) current density in the SGC was higher for the sham group compared with the active rTMS group at 50 milliseconds (9.81 × 10−8 [4.26 × 10−8] μA/mm2 vs 1.87 × 10−7 [2.05 × 10−7] μA/mm2; z = 2.53; P = .01) and at 200 milliseconds (1.80 × 10−7 [1.10 × 10−7] μA/mm2 vs 1.10 × 10−7 [5.60 × 10−8] μA/mm2; z = −2.02; P = .04) after the TMS pulse. This comparison is also presented in Figure 4 as a region of interest activation around 200 milliseconds after stimulation. Mean (SD) SCD in the sham group was higher at 200 milliseconds compared with the active group (7.00 × 10−7 [7.51 × 10−7] μA/mm2 vs 1.57 × 10−7 [3.67−7 × 10−7] μA/mm2; z = −2.91; P = .004). Additionally, the SCS between the stimulation site (left DLPFC) and SGC was higher in the sham group than the active group at 200 milliseconds after the TMS pulse (0.74 [0.73] mm vs 0.20 [0.44] mm; z = −2.78; P = .006). Figure 5A shows the association of SGC source current density with HRSD-17 score in participants with MDD before rTMS treatment (ρ = 0.41; P = .03). The correlation shows that symptoms of depression were more pronounced when SGC current density was higher at 100 milliseconds after the TMS pulse. Figure 5B shows a correlation of the change in HRSD-17 score from baseline to post-rTMS treatment with change in SCS from baseline to post-rTMS treatment at approximately 100 milliseconds after the pulse. Patients with MDD showed a significant correlation between these measures after active treatment (ρ = 0.58; P = .047). The correlation between HRSD-17 change and SCS change in patients with MDD after sham treatment was weaker and not significant (ρ = 0.22; P = .54). The HRSD-17 score change demonstrated a stronger correlation with SCS change after active treatment compared with the correlation of SCS change after sham treatment. However, these 2 correlations were not significantly different (z = 0.877; P = .19).
This study’s results show that patients with MDD had higher-amplitude SGC-localized source measures and higher DLPFC-SGC effective connectivity compared with healthy controls. We found that SCD indexing had a high accuracy in discriminating patients with MDD from healthy controls and that the SGC-localized signal at approximately 100 milliseconds post-TMS pulse was correlated with depression severity before rTMS treatment. Moreover, the physiological signature of higher SGC-localized source measures in patients with MDD was attenuated after active rTMS treatment compared with patients with MDD who received sham rTMS. Finally, SGC-localized source measures (ie, SCS) were significantly correlated with improvement in depressive symptoms in the active treatment group but not in the sham group.
The higher SGC source activation of patients with MDD compared with healthy controls is in agreement with the limbic cortical dysregulation model,37 which suggests that the SGC is metabolically overactive in the depressed state38 and attenuated to lower levels with successful treatment.6 The attenuation of the SGC-localized SCD and SCS by rTMS treatment to values similar to those of untreated healthy controls, and the correlation of MDD symptom improvement with the change in left DLPFC-SGC SCS after treatment (Figure 5B) is supported by several converging lines of evidence that demonstrate reduction in SGC hyperactivity after responding to treatment, such as pharmacotherapy,6,39 deep brain stimulation,4 vagus nerve stimulation,40 and rTMS.5,41-44
The SGC-DLPFC interconnectivity may, in part, be governed by γ-aminobutyric acid (GABA)–ergic neurotransmission.45 Because rTMS is capable of modulating deep brain structure activity transsynaptically,26 the observed reduction of the effective connectivity signal (ie, SCS) between those 2 regions might be associated with the reconfiguration of these long-range connections. The SGC current density at approximately the 100 millisecond peak, its correlation with symptom severity at baseline (Figure 5A), and the correlation between symptom change and SCS signal change at approximately 100 milliseconds (Figure 5B) may also be associated with GABAergic neurotransmission because the evoked potential at approximately 100 milliseconds post-TMS pulse is associated with neuronal inhibitory processes.46-49 Using the EEG sensor distribution for localizing electrical generators in the cortex yields different possible spatial arrangements of source current generators. Incorporating some assumptions regarding brain anatomy and electrophysiological energetic constraints into the source localizing computation will restrict the multiplicity of possible solutions, but it will remain nonabsolute. Further, the immense spatial and temporal complexity of brain activity and the noisy nature of the EEG signal increase uncertainties when trying to solve this inverse problem.50 However, when TMS-EEG is applied, the timing and location of brain activation is constrained by the experimenter, reducing substantially the short latency inaccuracies of the source estimations. Additionally, as TMS was applied over the DLPFC, we can rely on prior knowledge regarding DLPFC connectivity and its association with our main region of interest, the SGC. This supports the potential of SGC activation after longer latencies following the TMS induction.
Apart from the advantages of using TMS-EEG to assess activation of localized generators in the brain mentioned earlier, we also implemented SCD and SCS statistical indexing over our data, which further mitigates erroneous source evaluations. On its own, simple signal averaging is a strong tool for reducing noise in our electrophysiological recordings. When applied over the averaged signal, SCD is far more statistically stringent by taking into account the whole distribution of the activation—making this indexing robust when dealing with signal outliers.13 The SCS takes the SCD approach a step further by accumulating the distances between significant sources in the SGC and the site of stimulation. Hence, the SCS computation estimates the effective connectivity between those 2 sites after stimulation. The methodological and computational steps we took to ensure our data quality and the fact that our a priori hypotheses were supported by our results were key in making our source estimations findings more deterministic and valid.
This study had some limitations. Our study design allowed only a between-participant rather than within-participant statistical inference, which is not ideal in cases that look at the potential effects of treatment. Despite this issue, each group in this study had an adequate number of participants. Moreover, the measured associations produced highly comparable values between the independent groups, enabling valid statistical inferences. It is important to note that to achieve adequate neurophysiological signals, the active treatment group was pooled from the bilateral rTMS treatment arm and the unilateral rTMS treatment arm (as depicted in Figure 1). However, the unilateral and bilateral treatment groups did not show statistical differences in EEG activations. As another limitation, patients with MDD in this study were receiving heterogeneous antidepressant pharmacotherapy. This fact may have confounded our clinical symptom assessments and EEG signal observations. However, the use of sham control treatment and the finding that only the active rTMS group demonstrated a correlation of the EEG signal with symptom improvement suggests that these changes may stem from the rTMS treatment itself rather than from concomitant pharmacotherapy. Moreover, the depression change score had a stronger correlation with SCS change after active treatment compared with the correlation after sham treatment. However, our sample size was likely not large enough to find a significant difference in these 2 correlations.
In conclusion, this study demonstrated the usefulness of TMS-EEG and the SCD/SCS computation when investigating the association of SGC activation and DLPFC-SGC effective connectivity with MDD pathophysiology and clinical improvement. Left DLPFC rTMS may have improved MDD symptoms by altering the connectivity between the DLPFC and the SGC, likely via GABAergic neurotransmission. These findings support the hypothesis of SGC involvement in the pathophysiology and clinical improvement of MDD and also provide an objective biological target to differentiate mood states in MDD and to differentiate patients with MDD from healthy controls. Further research in larger patient populations and across different treatment modalities are warranted to further assess the diagnostic reliability and clinical usefulness of TMS-EEG in MDD.
Accepted for Publication: April 25, 2019.
Published: June 5, 2019. doi:10.1001/jamanetworkopen.2019.5578
Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2019 Hadas I et al. JAMA Network Open.
Corresponding Author: Zafiris J. Daskalakis, MD, PhD, FRCPC, Temerty Centre for Therapeutic Brain Intervention, Centre for Addiction and Mental Health, 1001 Queen St W, Unit 4-1, Toronto, ON M6J 1H4, Canada (jeff.daskalakis@camh.ca).
Author Contributions: Dr Hadas had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Concept and design: Hadas, Sun, Lioumis, Blumberger, Daskalakis.
Acquisition, analysis, or interpretation of data: All authors.
Drafting of the manuscript: Hadas, Lioumis, Jones, Voineskos, Downar, Daskalakis.
Critical revision of the manuscript for important intellectual content: Hadas, Sun, Lioumis, Zomorrodi, Voineskos, Downar, Fitzgerald, Blumberger, Daskalakis.
Statistical analysis: Hadas, Zomorrodi, Jones, Daskalakis.
Obtained funding: Daskalakis.
Administrative, technical, or material support: Sun, Blumberger, Daskalakis.
Supervision: Lioumis, Downar, Fitzgerald, Blumberger, Daskalakis.
Conflict of Interest Disclosures: Dr Voineskos reported receiving grants from the Ontario Mental Health Foundation, the Eli Lilly Research Fellowship Fund, and the Alternative Funding Plan of the Academic Health Sciences Centres of Ontario during the conduct of the study. Dr Downar reported serving as an advisor to TMS Neuro Health and to Restorative Brain Clinics outside the submitted work. Dr Fitzgerald reported receiving equipment for research from Magventure, Neurosoft, and Medtronic; being the founder and medical director of TMS Australia; receiving personal fees from LivaNova and Bionomics; and holding patent 2013266017 outside the submitted work. Dr Blumberger reported receiving in-kind equipment support from Magventure; serving as primary site investigator for 3 investigator-initiated trials for Brainsway; receiving grants from the Canadian Institutes of Health Research, Brain Canada, and Temerty Centre for Therapeutic Brain Intervention during the conduct of the study; attending an advisory board meeting for Janssen Pharmaceuticals; and receiving grants from the Canadian Institutes of Health Research, the National Institutes of Health, the Patient-Centered Outcomes Research Institute, and the Weston Foundation outside the submitted work. Dr Daskalakis reported receiving grants from Brainsway and Magventure during the conduct of the study. No other disclosures were reported.
Meeting Presentation: This article was presented at the World Congress of Biological Psychiatry; June 5, 2019; Vancouver, British Columbia, Canada.
1.Kessler
RC, Berglund
P, Demler
O,
et al; National Comorbidity Survey Replication. The epidemiology of major depressive disorder: results from the National Comorbidity Survey Replication (NCS-R).
JAMA. 2003;289(23):3095-3105. doi:
10.1001/jama.289.23.3095PubMedGoogle ScholarCrossref 3.Mayberg
HS, Liotti
M, Brannan
SK,
et al. Reciprocal limbic-cortical function and negative mood: converging PET findings in depression and normal sadness.
Am J Psychiatry. 1999;156(5):675-682. doi:
10.1176/ajp.156.5.675PubMedGoogle Scholar 5.Fox
MD, Buckner
RL, White
MP, Greicius
MD, Pascual-Leone
A. Efficacy of transcranial magnetic stimulation targets for depression is related to intrinsic functional connectivity with the subgenual cingulate.
Biol Psychiatry. 2012;72(7):595-603. doi:
10.1016/j.biopsych.2012.04.028PubMedGoogle ScholarCrossref 9.Keedwell
P, Drapier
D, Surguladze
S, Giampietro
V, Brammer
M, Phillips
M. Neural markers of symptomatic improvement during antidepressant therapy in severe depression: subgenual cingulate and visual cortical responses to sad, but not happy, facial stimuli are correlated with changes in symptom score.
J Psychopharmacol. 2009;23(7):775-788. doi:
10.1177/0269881108093589PubMedGoogle ScholarCrossref 20.Railo
H, Tuominen
J, Kaasinen
V, Pesonen
H. Dynamic changes in cortical effective connectivity underlie transsaccadic integration in humans.
Cereb Cortex. 2017;27(7):3609-3617. doi:
10.1093/cercor/bhw182PubMedGoogle Scholar 22.Ferrarelli
F, Riedner
BA, Peterson
MJ, Tononi
G. Altered prefrontal activity and connectivity predict different cognitive deficits in schizophrenia.
Hum Brain Mapp. 2015;36(11):4539-4552. doi:
10.1002/hbm.22935PubMedGoogle ScholarCrossref 23.Blumberger
DM, Maller
JJ, Thomson
L,
et al. Unilateral and bilateral MRI-targeted repetitive transcranial magnetic stimulation for treatment-resistant depression: a randomized controlled study.
J Psychiatry Neurosci. 2016;41(4):E58-E66. doi:
10.1503/jpn.150265PubMedGoogle ScholarCrossref 26.Rossini
PM, Burke
D, Chen
R,
et al. Non-invasive electrical and magnetic stimulation of the brain, spinal cord, roots and peripheral nerves: basic principles and procedures for routine clinical and research application: an updated report from an IFCN Committee.
Clin Neurophysiol. 2015;126(6):1071-1107. doi:
10.1016/j.clinph.2015.02.001PubMedGoogle ScholarCrossref 28.Oostenveld
R, Fries
P, Maris
E, Schoffelen
J-M. FieldTrip: Open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data.
Comput Intell Neurosci. 2011;2011:156869.
PubMedGoogle ScholarCrossref 31.Mazziotta
J, Toga
A, Evans
A,
et al. A probabilistic atlas and reference system for the human brain: International Consortium for Brain Mapping (ICBM).
Philos Trans R Soc Lond B Biol Sci. 2001;356(1412):1293-1322. doi:
10.1098/rstb.2001.0915PubMedGoogle ScholarCrossref 33.Pascual-Marqui
RD. Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details.
Methods Find Exp Clin Pharmacol. 2002;24(suppl D):5-12.
PubMedGoogle Scholar 41.Teneback
CC, Nahas
Z, Speer
AM,
et al. Changes in prefrontal cortex and paralimbic activity in depression following two weeks of daily left prefrontal TMS.
J Neuropsychiatry Clin Neurosci. 1999;11(4):426-435. doi:
10.1176/jnp.11.4.426PubMedGoogle Scholar 42.Nahas
Z, Teneback
CC, Kozel
A,
et al. Brain effects of TMS delivered over prefrontal cortex in depressed adults: role of stimulation frequency and coil-cortex distance.
J Neuropsychiatry Clin Neurosci. 2001;13(4):459-470. doi:
10.1176/jnp.13.4.459PubMedGoogle ScholarCrossref 43.Kito
S, Fujita
K, Koga
Y. Regional cerebral blood flow changes after low-frequency transcranial magnetic stimulation of the right dorsolateral prefrontal cortex in treatment-resistant depression.
Neuropsychobiology. 2008;58(1):29-36. doi:
10.1159/000154477PubMedGoogle ScholarCrossref