ORIGINAL RESEARCH article
Longitudinal Changes in Depressive Circuitry in Response to Neuromodulation Therapy
- 1Department of Biomedical Engineering, Marquette University, Milwaukee, WI, USA
- 2Department of Psychiatry, Medical College of Wisconsin, Milwaukee, WI, USA
- 3McConnell Brain Imaging Centre, Montreal Neurological Institute, McGill University, Montreal, QC, Canada
- 4Department of Neurology, Medical College of Wisconsin, Milwaukee, WI, USA
- 5Department of Bioengineering, Scientific Computing and Imaging (SCI) Institute, University of Utah, Salt Lake City, UT, USA
Background: Major depressive disorder (MDD) is a public health problem worldwide. There is increasing interest in using non-invasive therapies such as repetitive transcranial magnetic stimulation (rTMS) to treat MDD. However, the changes induced by rTMS on neural circuits remain poorly characterized. The present study aims to test whether the brain regions previously targeted by deep brain stimulation (DBS) in the treatment of MDD respond to rTMS, and whether functional connectivity (FC) measures can predict clinical response.
Methods: rTMS (20 sessions) was administered to five MDD patients at the left-dorsolateral prefrontal cortex (L-DLPFC) over 4 weeks. Magnetoencephalography (MEG) recordings and Montgomery-Asberg depression rating scale (MADRS) assessments were acquired before, during and after treatment. Our primary measures, obtained with MEG source imaging, were changes in power spectral density (PSD) and changes in FC as measured using coherence.
Results: Of the five patients, four met the clinical response criterion (40% or greater decrease in MADRS) after 4 weeks of treatment. An increase in gamma power at the L-DLPFC was correlated with improvement in symptoms. We also found that increases in delta band connectivity between L-DLPFC/amygdala and L-DLPFC/pregenual anterior cingulate cortex (pACC), and decreases in gamma band connectivity between L-DLPFC/subgenual anterior cingulate cortex (sACC), were correlated with improvements in depressive symptoms.
Conclusions: Our results suggest that non-invasive intervention techniques, such as rTMS, modulate the ongoing activity of depressive circuits targeted for DBS, and that MEG can capture these changes. Gamma oscillations may originate from GABA-mediated inhibition, which increases synchronization of large neuronal populations, possibly leading to increased long-range FC. We postulate that responses to rTMS could provide valuable insights into early evaluation of patient candidates for DBS surgery.
Major depressive disorder (MDD) is characterized by depressed mood, anhedonia, irritability, poor concentration, feelings of hopelessness, thoughts of self-harm, inappropriate guilt, and abnormal appetite and sleep patterns. About a sixth of the American population will suffer from clinical depression in their lifetime. Because of its significant societal and economic impact on the general population, both the underlying mechanism(s) of MDD and markers of individual response to treatment are important areas of ongoing research.
Previous studies have proposed that MDD is a circuit-based disorder in which several regions of the brain (prefrontal cortex, amygdala, striatum, pallidum and medial thalamus) are functionally aberrant (Drevets et al., 1992). Activity in the limbic regions, specifically the amygdala and anterior cingulate cortex (ACC; Brodmann’s area 25), is associated with depressive symptoms. Several hypotheses have been proposed to characterize MDD. The monoamine hypothesis (Krishnan and Nestler, 2008; Pittenger and Duman, 2008) suggests that a decrease in monoamine function can lead to depression. Among the monoamine projections from mid-brain to brainstem nuclei are: (1) dopamine from the ventral tegmental area (VT area); (2) serotonin from dorsal raphe in the periaqueductal gray area (PAG); and (3) noradrenaline from locus coeruleus. Secondary neuroplasticity on a cellular level is also induced by some antidepressants that acutely increase synaptic monoamines. In addition, the brain-derived neurotrophic factor (BDNF) hypothesis (Nestler and Carlezon, 2006; Krishnan and Nestler, 2008; Pittenger and Duman, 2008) implies that an increase in BDNF at the VT area and nucleus accumbens (NAc) may increase neurovegetative symptoms associated with depression.
Due to the heterogeneity of MDD, about a third of the patients suffering from the disorder are resistant to pharmacological antidepressants and/or psychotherapy. Neuromodulation therapies have emerged as potential treatments for depression. Dorsolateral prefrontal cortex (DLPFC) and ventromedial prefrontal cortex (VMPFC) are two of the main targets of neuromodulation therapy. These two targets exhibit opposite metabolic activity in the pathological depressed state (VMPFC: hyper, DLPFC: hypo; Ressler and Mayberg, 2007; Koenigs and Grafman, 2009; Kopell et al., 2011; Pathak et al., 2013). Deep brain stimulation (DBS) at the subgenual ACC (sACC) of the VMPFC network has shown promising results (Mayberg et al., 2005; Lozano et al., 2008; Mayberg, 2009; Hamani et al., 2010, 2011; Holtzheimer et al., 2012). The mechanisms of action in DBS are likely complex; however, a suppression of activity at the sACC has generally been reported to be correlated with improvement in depressive symptoms. Conversely, excitatory, noninvasive stimulation at the DLPFC via repetitive transcranial magnetic stimulation (rTMS) has efficaciously ameliorated depressive symptoms for patients (O’Reardon et al., 2007; Fitzgerald et al., 2009; Herbsman et al., 2009; Downar and Daskalakis, 2013). Lesion studies and advances in neuroimaging have provided evidence of the role of PFC in depression. In this study, we applied excitatory rTMS at the DLPFC because of its superficial location and because it plays a major role in higher cognitive functions and affective regulation.
rTMS alters neuronal activity directly beneath the excitation coil and at distant associated areas, as mediated by functional connections (Walsh and Cowey, 2000; Goodwin and Butson, 2015). Slow rTMS (<1 Hz) has been shown to be inhibitory, whereas fast rTMS (>5 Hz) is considered excitatory (Chen et al., 1997; Nakamura et al., 1997; Burt et al., 2002; Lefaucheur, 2008). Recent studies have investigated TMS as both a priming method and a treatment option for depression (Nahas et al., 2001; O’Reardon et al., 2007; Herbsman et al., 2009). The stimulation location and intensity, which are determined based on resting motor threshold (RMT), are important parameters in rTMS therapy. Lastly, magnetoencephalography (MEG) or electroencephalography (EEG) have been used to investigate local responses and remote interactions induced by rTMS (Siebner et al., 2009; Thut and Pascual-Leone, 2010).
In terms of possible biomarkers of MDD, several studies have investigated electrophysiological correlates as predictors of clinical response. Pre-treatment increase in theta power at the ACC as computed with EEG source imaging may reflect increased metabolism in ACC and is correlated with anti-depressive response (Pizzagalli, 2011; Spronk et al., 2011; Olbrich and Arns, 2013). Greater alpha power and hemispheric asymmetry have also been reported as potential predictors of response (Ulrich et al., 1984; Bruder et al., 2008; Breitenstein et al., 2014). In addition to resting-state studies, event-related potentials (ERPs) have been assessed in context of depression. Strong pre-treatment loudness dependent auditory evoked potentials (LDAEPs; Gallinat et al., 2000; Juckel et al., 2007) and P300 (Kalayam and Alexopoulos, 1999) have been associated with antidepressive response after drug treatment. These biomarkers in response to an auditory stimulus are modulated by the serotonergic system (Breitenstein et al., 2014).
We used resting-state MEG to assess possible longitudinal changes induced by high-frequency rTMS therapy. MEG source imaging is a functional neuroimaging method that senses magnetic fields generated by neural activity (Baillet et al., 2001). Specifically, scalp MEG recordings are generated essentially by synchronous postsynaptic currents driven by assemblies of pyramidal neurons oscillating over a relatively wide frequency spectrum. MEG imaging has excellent temporal resolution (millisecond) and magnetic fields are less distorted by skull and scalp tissues than EEG signals. This enables more robust modeling of neural generators in terms of mapping local activity and interregional coupling, based on time-dependent measures (Baillet et al., 2001; Gross et al., 2013).
Choosing left-dorsolateral prefrontal cortex (L-DLPFC) as a target for stimulation remains ambiguous due its broad anatomical definition. In a resting-state functional magnetic resonance imaging (fMRI) study (Fox et al., 2012a), measures from healthy participants were used to develop a connectivity-based targeting strategy to identify optimized L-DLPFC coordinates for stimulation in depressed patients. The approach included a paired comparison of the functional connectivity (FC) between two sites that previously showed clinical significance, and a correlation analysis between connectivity measures and clinical efficacy based on a model developed by Herbsman et al. (2009). The results indicate that better clinical efficacy was obtained when L-DLPFC coordinates were negatively correlated with sACC, which accounted for 70% of the variance in clinical efficacy. This is in line with the view of MDD as a “network disease” because it affects multiple, anatomically distributed brain regions.
FC is a rapidly developing subfield of neuroimaging that has contributed to better definition and testing of network models in clinical and cognitive neuroscience (Schoffelen and Gross, 2009; Sakkalis, 2011). In line with the theory of communication through coherence as a vehicle of FC in the brain, we computed magnitude-squared coherence between regions of interest (Fries, 2005). Coherence identifies significant frequency-dependent correlations between two time variant signals. It is sensitive to power and phase of neural oscillations and measures relationship stability.
We emphasize that the efficacy of rTMS therapy in treating MDD has previously been established and is not the focus of this study. Rather, the current study aims to identify objective indicators of antidepressive response to rTMS using functional neuroimaging. The long-term goal is to predict which patients will respond to anti-depressive treatment and to customize neuromodulation therapy accordingly. We are interested in two primary questions: (1) Which depressive circuits previously targeted for DBS are modulated by successful rTMS therapy for MDD? (2) Can changes in FC predict clinical response? Our hypothesis is that modulation of the L-DLPFC will be reflected through longitudinal changes in the power spectrum of local neural oscillations, along with changes in FC with other regions of the depressive circuit.
The major contribution of this study is to provide a bridge between physiological markers of MDD (neuroimaging) and clinical response to treatment. In the “Materials and Methods” Section, we describe the acquisition of high-temporal resolution MEG before, during and after delivery of high-frequency rTMS. Next, we present the approach of combining MEG source imaging with neuronavigated-rTMS to quantify how changes in depressive symptoms correlate with modulating specific circuits. In the results section, we report the significance of engaging specific functional connections in the fronto-limbic circuit for effective MDD improvement demonstrated through frequency-dependent connectivity analysis. Lastly, in the discussion section, we highlight the advances in source-imaging analysis techniques that have enabled the current study along with the limitations of those methods. We also discuss the relevance of our findings in the context of the field of neuromodulation and MDD.
Materials and Methods
High frequency (10 Hz) rTMS was administered to the L-DLPFC (Brodmann areas 9 and 46) in five MDD patients (n = 5). Patient selection and the decision to treat using rTMS was made through standard care, independent of this study. In tandem with rTMS treatment, we enrolled patients into a research study approved by the Institutional Review Board (IRB) at the Medical College of Wisconsin. During treatment, a total of 1000 biphasic pulses were delivered in 25-s cycles (5 s on, 20 s off) over 20 trains per session using a Magstim Rapid stimulator (Magstim Company Limited, Whitland, UK). Treatment was delivered for five 15-min sessions per week over 4 weeks. Stimulation was applied at 110% of RMT, which was determined by stimulating the hand-knob of motor cortex. Periodic behavioral assessments were performed with the Montgomery-Asberg depression rating scale (MADRS).
High-resolution (3T) structural magnetic resonance imaging (MRI) was acquired before treatment to assist with targeting and navigation. It was used for definition of TMS targets and registration with MEG imaging using FreeSurfer (Dale et al., 1999; Martinos Center for Biomedical Imaging) and Brainstorm (Tadel et al., 2011). A Visor neuronavigation system (ANT Neuro, Enschede, Netherlands) was used for anatomical targeting, and to record the position and orientation of the TMS coil during each pulse (Figure 1).
Figure 1. Neuronavigation for repetitive transcranial magnetic stimulation (rTMS) therapy. Screenshots from ANT software outline the pipeline used to target and administer neuro-navigated rTMS. Subject-specific magnetic resonance imaging (MRI) was used for 3D reconstruction of the head volume. High-frequency rTMS (10 Hz) at 110% of the resting motor threshold (RMT) was then applied at the left-dorsolateral prefrontal cortex (L-DLPFC) with real-time targeting using the ANT neuro-navigation system. (+x towards the nose, +y towards the left ear, and +z in the superior direction with the origin at the center of the head).
MEG was recorded at 2000 Hz using a 306-channel VectorView system (Elekta Neuromag Ltd, Helsinki, Finland). EOG and EKG were also recorded simultaneous to MEG for post hoc artifact removal. MEG was acquired at four time points: before treatment, halfway through the treatment course (2 weeks), at the end of the treatment course (4 weeks), and 2 months after the end of treatment. Empty-room recordings were also collected at each acquisition for noise modeling in source imaging. During MEG, the subject was supine under the MEG sensor array. A total of 24 min of data was collected in four 6-min runs. Each condition (eyes-open and eyes-closed) was repeated once.
MEG data for the eyes-closed condition were analyzed for cortical source modeling, power spectral density (PSD), and FC using Brainstorm, with default imaging parameters (Tadel et al., 2011). The results were assessed to identify longitudinal changes in PSD over standard frequency bands (delta = 2–4 Hz; theta = 5–7 Hz; alpha = 8–12 Hz; beta = 15–29 Hz; gamma = 30–59 Hz) of ongoing oscillatory neural activity. Continuous signals were cleaned from physiological and other artifacts with the signal-space separation (SSS) method (Taulu et al., 2005) and further corrected with specific signal-space projectors (SSP), a method based on principal components analysis to eliminate artifacts caused by eye-blinks and heartbeats; all were performed using default settings in Brainstorm. The data was then processed with a 60-Hz notch filter and a 2-Hz high-pass filter to remove residual noise due to the power line and slow drifts of sensor DC offset.
The weighted minimum-norm estimate (wMNE), a subtype of distributed source models, was used to generate the source image maps, with default imaging parameters in Brainstorm. Elementary dipole current sources were constrained perpendicularly to the individual cortex; additive sensor noise was modeled with a Gaussian distribution with known covariance statistics estimated from empty-room recordings. The Welch periodogram was used to obtain the PSD of each individual source time series on the cortex. Signals from the L-DLPFC were averaged over the stimulation area to determine the power at that region (Goodwin and Butson, 2015). For this analysis, L-DLPFC was defined based on the Desikan-Killiany atlas (Desikan et al., 2006) from FreeSurfer. A linear regression analysis was performed in MATLAB (Mathworks Inc., Natick, MA, USA) to determine if power within a particular frequency band was correlated with clinical improvement.
Coherence measures computed on 10-s epochs were used to determine FC between L-DLPFC and several regions of interest (ROIs) previously investigated for DBS therapy: pregenual anterior cingulate cortex (pACC) and sACC, NAc, and amygdala. The Desikan-Killiany atlas (Desikan et al., 2006) from FreeSurfer was used to define the ROIs. In addition, a cluster analysis was conducted in MATLAB and visualized in SCIRun (Center for Integrative Biomedical Computing, University of Utah). For this analysis, each node on the L-DLPFC surface was assigned a mean connectivity score based on the average of its coherence with signals at each of the nodes in the other ROIs (pACC, sACC or NAc). For example, if the L-DLPFC was parcellated into N1 nodes and sACC was parcellated into N2 nodes, a coherence matrix of N1 × N2 would be computed to reflect the connectivity between L-DLPFC and sACC. Coherence would then be averaged for all nodes in the sACC to assign one connectivity score to each node of the L-DLPFC, collapsing the dimensions of the initial matrix to N1 × 1. Mean connectivity was represented with a colormap to visualize the differential effects of sub-regions in the L-DLPFC.
We used a multiple regression model (Equation 1) to investigate the relationship between change in behavior (MADRS) and change in brain activity, controlling for time and subject effects which were treated as categorical variables. We performed this regression analysis for each frequency band in the source analysis (L-DLPFC) and the connectivity analysis (Amygdala, sACC, pACC and NAc). We did not correct for multiple comparisons given the exploratory nature of this study and its small sample size. All tests were two-tailed and statistical significance was assessed at the 0.05 level.
Data from five MDD patients were analyzed in this study. Of them, four were responders and one was a non-responder using the clinical response criterion of a 40% reduction in MADRS scores from baseline. Stimulation was applied to the L-DLPFC using a predefined target coordinate on an anatomical MRI. By midpoint, an improvement in depressive symptoms was observed in responders (Table 1). These participants continued to improve during treatment, and were still in remission at the 2-month follow-up (Figure 2).
Table 1. Summary of Montgomery-Asberg depression rating scale (MADRS) scores for all subjects over 4 weeks of repetitive transcranial magnetic stimulation (rTMS) therapy and at 12-week follow-up (R, Responder; NR, Non-responder).
Figure 2. Relationship between Montgomery-Asberg depression rating scale (MADRS) scores and time. Drastic improvement occurred in responders during the course of treatment, which persisted at least for another 8 weeks, as demonstrated by the scores at follow-up.
Increase in High-Frequency Power at the L-DLPFC Correlates with an Improvement in Depressive Symptoms
The results from the PSD analysis indicated that an increase in ongoing gamma (γ; 30–59 Hz) power (p(γ) = 0.01) and beta (12–29 Hz; β) power (p(β) = 0.04) at the L-DLPFC was significantly correlated with an improvement in depressive symptoms (Table 2). A power distribution over cortical maps (Figure 3) showed a marked increase of high-frequency power bilaterally in the frontal regions of responders throughout the course of treatment. These observations were corroborated by a multiple regression analysis (Figure 4).
Figure 3. Changes in cortical power spectral density (PSD). These cortical maps display the normalized difference in PSD from baseline in (A) beta (13–29 Hz) and (B) gamma (30–59 Hz) frequency bands. For the responders, there is a trend towards increase in high-frequency power bilaterally, which corresponds to improvements in depressive symptoms, as determined by the change in MADRS (right column for each patient).
Figure 4. Increase in high-frequency power at the L-DLPFC is correlated with an improvement in depressive symptoms. Multiple regression models of MADRS vs. PSD controlled for time and subject, show that an increase in high-frequency power is indicative of lower MADRS or better clinical outcome. Each symbol represents a different subject and the trend line is a reflection of the model to display the relationship between (A) MADRS and beta power, and (B) MADRS and gamma power.
Opposite Effects of Functional Connectivity in High and Low Frequency Bands Between L-DLPFC and Remote Regions
Changes in FC of L-DLPFC and several ROI in the depressive circuit were assessed. An increase in mean coherence in the delta band (2–4 Hz) correlated with a better clinical outcome for depression in the connectivity analysis between L-DLPFC and amygdala (p(δ) = 0.04; Table 2), as demonstrated by a multiple regression model of MADRS with coherence (Figure 5). A qualitative cluster analysis also demonstrated an increase in DLPFC-amygdala coherence from baseline to endpoint in delta band over the entire area of stimulation for the responders when mapped onto the cortex. In contrast, a decrease in DLPFC-sACC coherence in the gamma band was observed. Furthermore, the nodes involved in modulating connectivity were common for amygdala and sACC (Figure 6).
Figure 5. Increase in low-frequency connectivity between L-DLPFC and amygdala. Multiple regression models of MADRS vs. coherence show that there is a significant linear relationship between connectivity and an improvement in depressive symptoms. Each symbol represents a different subject and the trend line is a reflection of the model to display the relationship between MADRS and coherence in the delta band.
Figure 6. rTMS modulates connectivity over the L-DLPFC in responders. L-DLPFC/amygdala coherence in delta band (bottom) increased from baseline to endpoint over the entire area of stimulation for responders. An opposite effect was observed with a decrease in L-DLPFC/subgenual anterior cingulate cortex (sACC) coherence in gamma band (top) from baseline to endpoint over the entire area of stimulation for responders. The nodes in L-DLPFC contributing to change in connectivity are common for both amygdala and sACC for all patients.
The target of rTMS stimulation was the L-DLPFC as this region tends to show diminished function in individuals with MDD. This area is implicated in working-memory, attention, and long-term memory storage and consolidation and previous studies have shown improvements following stimulation in frontal cortices (Curtis and D’Esposito, 2003; Fregni et al., 2005). We did not formally test working memory, but questions relating to recall, attention and orientation were asked daily to ensure there were no negative consequences of stimulation. Our results support the hypothesis that longitudinal changes from rTMS therapy are reflected in the PSD and connectivity analyses of ongoing neural oscillations in the resting-state. We emphasize that the cohort for this study consisted of treatment-resistant patients who failed previously at least two treatment strategies; this implies that a placebo response is unlikely and cannot fully explain the improvement in MDD symptoms. Furthermore, we show that correlation of such improvement with the neurophysiological changes captured by MEG is an objective way of assessing rTMS therapy.
A local increase in beta and gamma power at the L-DLPFC was correlated with an improvement in depressive symptoms. A comparison of PSD in the right occipital cortex as a control region did not result in a trend similar to the one observed at the L-DLPFC. We posit that this increased power may reflect a rise in metabolic activity of the DLPFC, which was shown to contribute to improvement of symptoms. High-frequency oscillations in the beta and gamma range have previously been linked to increased attentional processes (Benchenane et al., 2011). One mechanism of gamma modulation at a pharmacological level involves GABAergic interneurons. These cells contribute to synchronization of the pyramidal neurons in the DLPFC during tasks that require increased cognitive effort. Furthermore, GABA-mediated inhibition was shown to be sufficient in generating gamma oscillations (Barr et al., 2009).
Previous evidence also points to GABAergic involvement in antidepressant response (Bajbouj et al., 2005; Yue et al., 2009). Animal models of depression have shown that GABA agonists have antidepressant activity (Petty, 1995; Yue et al., 2009). If inhibition is facilitated by fast acting GABAA receptors, oscillation frequency will correspond to the gamma band (Buzsaki, 2006). Additionally, increased cortical beta is known to enhance the effect of GABA. Consequently, it is plausible that increased cortical beta activity, which was also observed in this study, boosted the effect of GABA leading to inhibition that increased gamma power. Another consequence of GABA-mediated inhibition is synchronization of large neuronal populations, leading to increased long-range FC mediated by interregional coherence (Fingelkurts et al., 2004). Alternatively, post-synaptic targets are also known to up-regulate expression patterns of GABAergic receptors or change their expression pattern of the GABAergic receptor subtypes that could lead to changes in behavioral outcome (Lloyd et al., 1989; Petty, 1995; Brambilla et al., 2003). Though GABA may be involved, it is difficult to determine exactly how it is exerting its effects in this case.
Changes in synaptic plasticity are directly related to changes measured in FC and the strength of a synapse can further dictate the effect of inhibition (Buzsaki, 2006). The depressive circuit targeted for DBS includes amygdala, pACC, sACC, and the NAc (Ressler and Mayberg, 2007; Krishnan and Nestler, 2008; Holtzheimer and Mayberg, 2011). Recent developments in source imaging methods enable a more accurate representation of neural activity at deeper sources. Attal and Schwartz (2013) demonstrated that realistic anatomical and electrophysiological models of deep brain activity (DBA) allow for the estimation of contributions from deep sources such as thalamus and amygdala at MEG sensors. In the present study, the measures extracted are based on integration over a large amount of time samples (power spectrum density and coherence). Therefore, these are equivalent to a form of averaging effect (enhancement of SNR), which indicates that signal detection was substantially improved. Additionally, source imaging was obtained using the MRI obtained from each participant to account for individual anatomy and noise statistics were captured using an empty-room recording, and embedded into the source estimator. The most effective targets for DBS are sACC and amygdala. This evidence further lends credence to our results that indicate the significance of modulating connectivity between these specific regions.
Activity in sACC contributes substantially to limbic-cortical dysregulation (Mayberg, 1997; Greicius et al., 2007). Depression severity has been shown to correlate positively with FC in sACC. In addition, increased connectivity of the default mode network (DMN) in depression was also driven by this region (Greicius et al., 2007). Fox et al. (2012b) also showed that the magnitude of relative disconnectivity (anticorrelation of fMRI signals) between L-DLPFC and ACC was indicative of depression severity. The distinct role of subregion FC is corroborated by differential anatomical connectivity of the pACC and sACC, suggesting a strong connection of pACC to frontal regions vs. sACC to limbic regions (Johansen-Berg et al., 2008). These studies imply that when an individual is improving, anti-correlation between sACC and DLPFC weakens, which would in turn manifest in a coherence decrease, supporting our results for the cluster analysis that clearly show the nodes involved in decreased gamma coherence between L-DLPFC and sACC (Figure 6). This further suggests that part of the antidepressant mechanism enabled by rTMS may be in the remote suppression of the sACC, which has been shown to be hypermetabolic in depressed individuals.
Another finding is that an increase in coherence was observed in the delta band for amygdala (Figures 5, 6). Different frequency bands may enable the multiplexing of channels of communication for different functional brain networks (Lakatos et al., 2005). Communication through coherence is thought to be regulated by inhibitory windows (Benchenane et al., 2011). FC in lower frequency bands has been associated with higher cognitive tasks. Previous research has also demonstrated that anatomically distributed EEG coherence increased in lower frequencies, particularly the delta band during creative thought (Petsche, 1996; Gruzelier, 2009). Furthermore, neuroimaging and neurophysiological studies have reported that long-distance integration is coordinated by lower frequencies, whereas short-distance networks are coordinated by higher frequencies (Buzsaki, 2006; Leuchter et al., 2012). In particular for mood disorders, an fMRI study (Anand et al., 2005) showed an increase in low-frequency blood-oxygen-level dependent (BOLD) fluctuation (LFBF) correlation between ACC and limbic regions in patients after treatment for depression. Another study (Chepenik et al., 2010) identified a negative correlation in low-frequency resting-state fMRI activity between the ventral PFC and amygdala in healthy subjects compared to those with bipolar disorder. These reports corroborate our current findings of increased connectivity between L-DLPFC and amygdala in the delta band correlating with better clinical outcome.
According to an MEG study of resting-state FC (Hillebrand et al., 2012), there is a positive correlation between patterns of source power and FC. A modeling study (Chawla et al., 1999, 2000) also showed increased phase locking between two connected neuronal populations resulting in increases in both power and FC. Though we did not note a statistical significance in the connectivity analysis between L-DLPFC and sACC in our multiple regression analysis, we would allude to the results of our cluster analysis that align with these previous findings and demonstrate decreased gamma connectivity between the two regions.
Lastly, the modulation in coherence between frontal and limbic regions was observed over the entire area of stimulation for responders, when mean connectivity was mapped onto the L-DLPFC. This finding implies that stimulation of responders engaged functional connections that were relevant for clinical improvement and nonresponders were effectively lacking in engaging those essential connections. In a previous study investigating chronic epidural cortical stimulation (EpCS) for depression (Kopell et al., 2011; Pathak et al., 2013), stimulation using the mixed mode polarity, in which the configuration of the active electrode alternated between anode and cathode in a pseudorandom manner, was effective in achieving an antidepressant response. This response was likely due to recruitment of a wider range of neurons owing to characteristic differences in anodic and cathodic stimulations (Wongsarnpigoon and Grill, 2008; Pathak et al., 2013). Another study (Riva-Posse et al., 2014) investigating fiber tracts concluded that the involvement of bilateral forceps minor, the cingulate bundle, and short descending midline fibers was predictive of clinical response. Riva-Posse et al. (2014) emphasized that the anatomical location of active electrodes was not the source of variability, but instead it was the hub location where these three fibers intersected that determined potential of response.
In terms of study design, a sham group was not included because the study was not testing actual treatment efficacy of TMS. Instead, our goal was to conduct a longitudinal study to evaluate neuroimaging changes in patient volunteers showing MDD clinical improvements following TMS treatment. While this still does not entirely rule out a possible placebo effect, it remains highly unlikely since these patients had failed at least two previous treatment options that are otherwise effective. Further, the STAR * D report elegantly explains that there is a reduction in remission rates as the number of treatment steps increase (Rush et al., 2006).
When working with resting state signals, we cannot disregard the possibility of negative self-talk or rumination. Future studies should attempt to control for this. However, this is a pilot study. Further, resting state here only implies task-free; how an individual controls their thoughts are beyond the scope of this study. In fact, this state of rumination might reflect the resting state of an MDD patient. Additionally, our data has been averaged over several segments of time for the analysis. If these occurrences are spontaneous, it is likely that their effect will be eliminated.
There are a few limitations that should be addressed in future studies as the field advances. First, interaction of TMS with medication is not clearly understood. There is a possibility that changes in medication can affect cortical excitability which might necessitate a re-evaluation of the TMS “dose”. Secondly, we analyzed data from a relatively small sample size. As this is a small patient cohort, we did not correct for multiple comparisons which could increase the probability of a Type I error. However, the longitudinal nature of our study design provides the opportunity to observe changes at an individual level over the course of neuromodulation therapy. With a push towards precision medicine, it becomes important to measure improvement in a patient-specific manner. Additionally, our findings regarding neuroimaging correlates of improvement in depression take into account response and non-response, and the results are selective for specific anatomical regions.
The primary finding from this study is that in addition to local changes, high-frequency rTMS modulates remote nodes of the depressive circuit during an effective clinical response. This conclusion was reached by combining anatomical targeting (from rTMS navigation) with source analysis and FC as determined from MEG. We postulate that this approach could fill two important gaps in our current knowledge. First, it could improve identification of specific regions and functional connections to target in order to maximize the therapeutic response. Second, this approach could provide insights about common mechanisms for different types of neuromodulation therapy including rTMS, cortical stimulation and DBS.
YP, OS, SB, ZL, and CRB participated in manuscript preparation and editing; YP and CRB conducted the analysis to determine the effect of repetitive transcranial magnetic stimulation on changes in the power spectral density and functional connectivity; OS was the psychiatrist for the clinical study and assisted with the administration of therapy; ZL assisted with the acquisition of neuroimaging data; SB guided the analysis and interpretation of the neuroimaging data.
Conflict of Interest Statement
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.
The reviewer CP and the handling Editor declared their shared affiliation, and the handling Editor states that the process nevertheless met the standards of a fair and objective review.
This project was supported by the Department of Neurology at Medical College of Wisconsin. Statistical analysis was supported, in part, by grant 1UL1RR031973 from the Clinical and Translational Science Institute (CTSI), and the Division of Biostatistics. SCIRun was provided by a grant from the NIH/NCRR Center for Integrative Biomedical Computing, P41GM103545. SB was supported by the Killam Foundation, a Senior-Researcher grant from the Fonds de Recherche du Québec-Santé, a Discovery Grant from the National Science and Engineering Research Council of Canada and the NIH (2R01EB009048-05).
Anand, A., Li, Y., Wang, Y., Wu, J., Gao, S., Bukhari, L., et al. (2005). Antidepressant effect on connectivity of the mood-regulating circuit: an FMRI study. Neuropsychopharmacology, 30, 1334–1344. doi: 10.1038/sj.npp.1300725
Attal, Y., and Schwartz, D. (2013). Assessment of subcortical source localization using deep brain activity imaging model with minimum norm operators: a MEG study. PLoS One 8:e59856. doi: 10.1371/journal.pone.0059856
Bajbouj, M., Brakemeier, E.-L., Schubert, F., Lang, U. E., Neu, P., Schindowski, C., et al. (2005). Repetitive transcranial magnetic stimulation of the dorsolateral prefrontal cortex and cortical excitability in patients with major depressive disorder. Exp. Neurol. 196, 332–338. doi: 10.1016/j.expneurol.2005.08.008
Barr, M. S., Farzan, F., Rusjan, P. M., Chen, R., Fitzgerald, P. B., and Daskalakis, Z. J. (2009). Potentiation of gamma oscillatory activity through repetitive transcranial magnetic stimulation of the dorsolateral prefrontal cortex. Neuropsychopharmacology, 34, 2359–2367. doi: 10.1038/npp.2009.79
Benchenane, K., Tiesinga, P. H., and Battaglia, F. P. (2011). Oscillations in the prefrontal cortex: a gateway to memory and attention. Curr. Opin. Neurobiol. 21, 475–485. doi: 10.1016/j.conb.2011.01.004
Bruder, G. E., Sedoruk, J. P., Stewart, J. W., McGrath, P. J., Quitkin, F. M., and Tenke, C. E. (2008). Electroencephalographic alpha measures predict therapeutic response to a selective serotonin reuptake inhibitor antidepressant: pre- and post-treatment findings. Biol. Psychiatry 63, 1171–1177. doi: 10.1016/j.biopsych.2007.10.009
Burt, T., Lisanby, S. H., and Sackeim, H. A. (2002). Neuropsychiatric applications of transcranial magnetic stimulation: a meta analysis. Int. J. Neuropsychopharmacol. 5, 73–103. doi: 10.1017/s1461145702002791
Chawla, D., Lumer, E. D., and Friston, K. J. (1999). The relationship between synchronization among neuronal populations and their mean activity levels. Neural Comput. 11, 1389–1411. doi: 10.1162/089976699300016287
Chawla, D., Lumer, E. D., and Friston, K. J. (2000). Relating macroscopic measures of brain activity to fast, dynamic neuronal interactions. Neural Comput. 12, 2805–2821. doi: 10.1162/089976600300014737
Chepenik, L. G., Raffo, M., Hampson, M., Lacadie, C., Wang, F., Jones, M. M., et al. (2010). Functional connectivity between ventral prefrontal cortex and amygdala at low frequency in the resting state in bipolar disorder. Psychiatry Res. 182, 207–210. doi: 10.1016/j.pscychresns.2010.04.002
Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D. et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021
Fingelkurts, A. A. A., Kivisaari, R., Pekkonen, E., Ilmoniemi, R. J., and Ka, S. (2004). Enhancement of GABA-related signalling is associated with increase of functional connectivity in human cortex. Hum. Brain Mapp. 22, 27–39. doi: 10.1002/hbm.20014
Fitzgerald, P. B., Hoy, K., McQueen, S., Maller, J. J., Herring, S., Segrave, R., et al. (2009). A randomized trial of rTMS targeted with MRI based neuro-navigation in treatment-resistant depression. Neuropsychopharmacology, 34, 1255–1262. doi: 10.1038/npp.2008.233
Fox, M. D., Buckner, R. L., White, M. P., Greicius, M. D., and Pascual-Leone, A. (2012a). Efficacy of transcranial magnetic stimulation targets for depression is related to intrinsic functional connectivity with the subgenual cingulate. Biol. Psychiatry 72, 595–603. doi: 10.1016/j.biopsych.2012.04.028
Fox, M. D., Halko, M., Eldaief, M., and Pascual-Leone, A. (2012b). Measuring and manipulating brain connectivity with resting state functional connectivity magnetic resonance imaging (fcMRI) and transcranial magnetic stimulation (TMS). Neuroimage 62, 2232–2243. doi: 10.1016/j.neuroimage.2012.03.035
Fregni, F., Boggio, P. S., Nitsche, M., Bermpohl, F., Antal, A., Feredoes, E., et al. (2005). Anodal transcranial direct current stimulation of prefrontal cortex enhances working memory. Exp. Brain Res. 166, 23–30. doi: 10.1007/s00221-005-2334-6
Gallinat, J., Bottlender, R., Juckel, G., Munke-Puchner, A., Stotz, G., Kuss, H.-J., et al. (2000). The loudness dependency of the auditory evoked N1/P2-component as a predictor of the acute SSRI response in depression. Psychopharmacology (Berl) 148, 404–411. doi: 10.1007/s002130050070
Greicius, M. D., Flores, B. H., Menon, V., Glover, G. H., Solvason, H. B., Kenna, H., et al. (2007). Resting-state functional connectivity in major depression: abnormally increased contributions from subgenual cingulate cortex and thalamus. Biol. Psychiatry 62, 429–437. doi: 10.1016/j.biopsych.2006.09.020
Gross, J., Baillet, S., Barnes, G. R., Henson, R. N., Hillebrand, A., Jensen, O., et al. (2013). Good practice for conducting and reporting MEG research. Neuroimage 65, 349–363. doi: 10.1016/j.neuroimage.2012.10.001
Gruzelier, J. (2009). A theory of alpha/theta neurofeedback, creative performance enhancement, long distance functional connectivity and psychological integration. Cogn. Process. 10, 101–109. doi: 10.1007/s10339-008-0248-5
Hamani, C., Diwan, M., Isabella, S., Lozano, A. M., and Nobrega, J. N. (2010). Effects of different stimulation parameters on the antidepressant-like response of medial prefrontal cortex deep brain stimulation in rats. J. Psychiatr. Res. 44, 683–687. doi: 10.1016/j.jpsychires.2009.12.010
Hamani, C., Mayberg, H. S., Stone, S., Laxton, A., Haber, S. N., and Lozano, A. M. (2011). The subcallosal cingulate gyrus in the context of major depression. Biol. Psychiatry 69, 301–308. doi: 10.1016/j.biopsych.2010.09.034
Herbsman, T., Avery, D., Ramsey, D., Holtzheimer, P. E., Wadjik, C., Hardaway, F., et al. (2009). More lateral and anterior prefrontal coil location is associated with better repetitive transcranial magnetic stimulation antidepressant response. Biol. Psychiatry 66, 509–515. doi: 10.1016/j.biopsych.2009.04.034
Hillebrand, A., Barnes, G. R., Bosboom, J. L., Berendse, H. W., and Stam, C. J. (2012). Frequency-dependent functional connectivity within resting-state networks: an atlas-based MEG beamformer solution. NeuroImage 59, 3909–3921. doi: 10.1016/j.neuroimage.2011.11.005
Holtzheimer, P. E., Kelley, M. E., Gross, R. E., Filkowski, M. M., Garlow, S. J., Barrocas, A., et al. (2012). Subcallosal cingulate deep brain stimulation for treatment-resistant unipolar and bipolar depression. Arch. Gen. Psychiatry 69, 150–158. doi: 10.1001/archgenpsychiatry.2011.1456
Johansen-Berg, H., Gutman, D. A., Behrens, T. E. J., Matthews, P. M., Rushworth, M. F. S., Katz, E., et al. (2008). Anatomical connectivity of the subgenual cingulate region targeted with deep brain stimulation for treatment-resistant depression. Cereb. Cortex 18, 1374–1383. doi: 10.1093/cercor/bhm167
Juckel, G., Pogarell, O., Augustin, H., Mulert, C., Müller-Siecheneder, F., Frodl, T., et al. (2007). Differential prediction of first clinical response to serotonergic and noradrenergic antidepressants using the loudness dependence of auditory evoked potentials in patients with major depressive disorder. J. Clin. Psychiatry 68, 1206–1212. doi: 10.4088/jcp.v68n0806
Koenigs, M., and Grafman, J. (2009). The functional neuroanatomy of depression: distinct roles for ventromedial and dorsolateral prefrontal cortex. Behav. Brain Res. 201, 239–243. doi: 10.1016/j.bbr.2009.03.004
Kopell, B. H., Halverson, J., Butson, C. R., Dickinson, M., Bobholz, J., Harsch, H., et al. (2011). Epidural cortical stimulation of the left dorsolateral prefrontal cortex for refractory major depressive disorder. Neurosurgery 69, 1015–1029; discussion 1029. doi: 10.1227/NEU.0b013e318229cfcd
Lakatos, P., Shah, A. S., Knuth, K. H., Ulbert, I., Karmos, G., and Schroeder, C. E. (2005). An oscillatory hierarchy controlling neuronal excitability and stimulus processing in the auditory cortex an oscillatory hierarchy controlling neuronal excitability and stimulus processing in the auditory cortex. J. Neurophysiol. 94, 1904–1911. doi: 10.1152/jn.00263.2005
Leuchter, A. F., Cook, I. A., Hunter, A. M., Cai, C., and Horvath, S. (2012). Resting-state quantitative electroencephalography reveals increased neurophysiologic connectivity in depression. PLoS One 7:e32508. doi: 10.1371/journal.pone.0032508
Lloyd, K. G., Zivkovic, B., Scatton, B., Morselli, P. L., and Bartholimi, G. (1989). The gabaergic hypothesis of depression. Prog. Neuropsychopharmacol. Biol. Psychiatry 13, 341–351. doi: 10.1016/0278-5846(89)90123-1
Lozano, A. M., Mayberg, H. S., Giacobbe, P., Hamani, C., Craddock, R. C., and Kennedy, S. H. S. H. (2008). Subcallosal cingulate gyrus deep brain stimulation for treatment-resistant depression. Biol. Psychiatry 64, 461–467. doi: 10.1016/j.biopsych.2008.05.034
Mayberg, H. S., Lozano, A. M., Voon, V., McNeely, H. E., Seminowicz, D., Hamani, C., et al. (2005). Deep brain stimulation for treatment-resistant depression. Neuron 45, 651–660. doi: 10.1016/j.neuron.2005.02.014
Nahas, Z., Teneback, C. C., Kozel, F. A., Speer, A. M., DeBrux, C., Molloy, M., et al. (2001). Brain effects of TMS delivered over prefrontal cortex in depressed adults: role of stimulation frequency and coil-cortex distance. J. Neuropsychiatry Clin. Neurosci. 13, 459–470. doi: 10.1176/appi.neuropsych.13.4.459
Nakamura, H., Kitagawa, H., Kawaguchi, Y., and Tsuji, H. (1997). Intracortical facilitation and inhibition after transcranial magnetic stimulation in conscious humans. J. Physiol. 498, 817–823. doi: 10.1113/jphysiol.1997.sp021905
Olbrich, S., and Arns, M. (2013). EEG biomarkers in major depressive disorder: discriminative power and prediction of treatment response. Int. Rev. Psychiatry 25, 604–618. doi: 10.3109/09540261.2013.816269
O’Reardon, J. P., Solvason, H. B., Janicak, P. G., Sampson, S. M., Isenberg, K. E., Nahas, Z., et al. (2007). Efficacy and safety of transcranial magnetic stimulation in the acute treatment of major depression: a multisite randomized controlled trial. Biol. Psychiatry 62, 1208–1216. doi: 10.1016/j.biopsych.2007.01.018
Pathak, Y., Kopell, B. H., Szabo, A., Rainey, C., Harsch, H., and Butson, C. R. (2013). The role of electrode location and stimulation polarity in patient response to cortical stimulation for major depressive disorder. Brain Stimul. 6, 254–260. doi: 10.1016/j.brs.2012.07.001
Riva-Posse, P., Choi, K. S., Holtzheimer, P. E., McIntyre, C. C., Gross, R. E., Chaturvedi, A., et al. (2014). Defining critical white matter pathways mediating successful subcallosal cingulate deep brain stimulation for treatment-resistant depression. Biol. Psychiatry 76, 963–969. doi: 10.1016/j.biopsych.2014.03.029
Rush, A. J., Trivedi, M. H., Wisniewski, S. R., Nierenberg, A. A., Stewart, J. W., Warden, D., et al. (2006). Acute and longer-term outcomes in depressed outpatients requiring one or several treatment steps: a STAR * D report. Am. J. Psychiatry 163, 1905–1917. doi: 10.1176/appi.ajp.163.11.1905
Siebner, H. R., Bergmann, T. O., Bestmann, S., Massimini, M., Johansen-Berg, H., Mochizuki, H., et al. (2009). Consensus paper: combining transcranial stimulation with neuroimaging. Brain Stimul. 2, 58–80. doi: 10.1016/j.brs.2008.11.002
Spronk, D., Arns, M., Barnett, K. J., Cooper, N. J., and Gordon, E. (2011). An investigation of EEG, genetic and cognitive markers of treatment response to antidepressant medication in patients with major depressive disorder: a pilot study. J. Affect. Disord. 128, 41–48. doi: 10.1016/j.jad.2010.06.021
Tadel, F., Baillet, S., Mosher, J. C., Pantazis, D., and Leahy, R. M. (2011). Brainstorm: a user-friendly application for MEG/EEG analysis. Comput. Intell. Neurosci. 2011:879716. doi: 10.1155/2011/879716
Ulrich, G., Renfordt, E., Zeller, G., and Frick, K. (1984). Interrelation between changes in the EEG and psychopathology under pharmacotherapy for endogenous depression. A contribution to the predictor question. Pharmacopsychiatry 17, 178–183. doi: 10.1055/s-2007-1017433
Keywords: repetitive transcranial magnetic stimulation (rTMS), neuromodulation, dorsolateral prefrontal cortex (DLPFC), depression, deep brain stimulation (DBS), functional connectivity
Citation: Pathak Y, Salami O, Baillet S, Li Z and Butson CR (2016) Longitudinal Changes in Depressive Circuitry in Response to Neuromodulation Therapy. Front. Neural Circuits 10:50. doi: 10.3389/fncir.2016.00050
Received: 16 July 2015; Accepted: 29 June 2016;
Published: 29 July 2016.
Edited by:Alexander Rotenberg, Boston Children’s Hospital, Harvard Medical School, USA
Reviewed by:Christos Papadelis, Boston Children’s Hospital, Harvard Medical School, USA
Alik Sunil Widge, Massachusetts General Hospital, USA
Copyright © 2016 Pathak, Salami, Baillet, Li and Butson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution and reproduction in other forums is permitted, provided the original author(s) or licensor 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.
*Correspondence: Christopher R. Butson, email@example.com