Non-linear EEG analyses predict non-response to rTMS treatment in major depressive disorder
Introduction
The application of repetitive transcranial magnetic stimulation (rTMS) for major depressive disorder (MDD) has been investigated intensively over the last 15 years. Several meta-analyses have demonstrated that compared to placebo, the effects of rTMS applied to the left or right dorsolateral prefrontal cortex (DLPFC) have antidepressant effects (Berlim et al., 2013, Schutter, 2009, Schutter, 2010). With the establishment of the efficacy of rTMS, there has been increased interest in finding predictors of clinical response. The value of clinical features in predicting treatment outcome in MDD is very limited (Bagby et al., 2002, Simon and Perlis, 2010) and a shift towards biomarkers is evident, evidenced by new initiatives such as the NIMH research domain criteria (RDoC) and precision medicine.
Since the first description by Lemere in 1936 who described that the capability to ‘…produce “good” alpha waves…’ was associated with the ‘…affective capacity of the individual…” (Lemere, 1936), increased EEG alpha power is an often replicated finding in MDD (Olbrich and Arns, 2013). Recent studies further suggest that alpha power can be considered an endophenotype, mediating the pathway between the Brain Derived Neurotrophic Factor (BDNF) Val66Met polymorphism and depressive mood (Gatt et al., 2008, Zoon et al., 2013). Increased alpha power has also been associated with a favourable treatment outcome to antidepressant drug-treatments (Bruder et al., 2008, Bruder et al., 2001, Tenke et al., 2011, Ulrich et al., 1986), and a slow alpha peak frequency (APF) has been associated with less favourable treatment response to rTMS and antidepressants (Arns et al., 2012, Arns et al., 2010, Conca et al., 2000). Conceptually, the increased posterior alpha power can be regarded as a hyperstable vigilance regulation resulting in MDD behaviour, and it is known that treatments that destabilize vigilance (e.g., sleep deprivation) have strong antidepressant effects (Hegerl and Hensch, 2012). Given this supposed effect of destabilization, we hypothesised whether nonlinear measures would be more appropriate to assess the predictive power of these two dimensions of alpha (its amplitude fluctuations and frequency domain).
Nonlinear measures such as global and local scaling exponents and quantification of phase space dynamics have been successfully used to distinguish between healthy and pathological time series of human physiology and performance (e.g., Attention deficit disorders Gilden and Hancock, 2007; Ageing and disease Goldberger et al., 2002, Little et al., 2006; Developmental dyslexia Wijnants et al., 2012).
Fitting within a more general framework of health, wellbeing and complexity (Stam, 2010, Van Orden et al., 2009), analyses of neurophysiological time series and functional and structural networks of the brain also reveal dynamics such as long-range temporal correlations (LRTC) and topological structure associated with self-organised critical states of complex systems (Berthouze et al., 2010, Freeman et al., 2009, Orlandi et al., 2013, Palva et al., 2013, Rubinov and Sporns, 2010, Tagliazucchi et al., 2012, Yu et al., 2013). Nonlinear analyses of neurophysiological data have been used to characterize sleep, coma, anaesthesia, epilepsy, drug effects, schizophrenia, Alzheimer’s disease and dementia as deviations from such optimal critical states or “small-world” network topology (Stam, 2005 for a review).
Nonlinear measures have been used to study MDD, Linkenkaer-Hansen et al. (2005) found a breakdown of LRTCs in the theta band (3–7 Hz) range in MDD patients during eyes-closed wakeful rest. The theta band LRCT magnitude over the left temporocentral region was negatively associated with severity of depression (r = −0.79), whereas the LRCT magnitude in the alpha band over the occipitoparietal region was positively correlated with severity (r = 0.59). These results seem to be in line with the reported increased alpha in MDD mentioned earlier, however, Linkenkaer-Hansen et al. (2005) did not find any correlations between amplitude and severity in different frequency bands which warrants further study of the relationship between these measures in MDD. In addition, few studies exist where nonlinear analyses have been applied to investigate efficacy of intervention in general and to the best of the authors knowledge, no studies exist where they have been used to predict treatment outcome in MDD to rTMS or antidepressant medication. Two studies have investigated non-linear measures and their association to treatment outcome to Electroconvulsive Treatment (ECT), in which a smaller post-seizure fractal dimension (Gangadhar et al., 1999) or smaller Largest Lyapunov Exponent (LLE) were interpreted to indicate a more predictable pattern of EEG seizure activity (Krystal et al., 1997) that was associated with a more favourable treatment outcomes.
The goal of many intervention studies using nonlinear analyses is to evaluate the added value of these methods above more traditional linear methods (Granic and Hollenstein, 2003, Hayes et al., 2007, Lichtwarck-Aschoff et al., 2012). Especially in the context of physiological measurements, the potential for clinical applications of these techniques receives much attention (Bravi et al., 2011, Harbourne and Stergiou, 2009, Huikuri and Stein, 2012). Although the potential of nonlinear tools is widely recognised, Bravi et al. (2011) note in their evaluation of 70 different variability analyses that the challenges for the field are to develop a shared vocabulary and increase coherence between results of different studies and the techniques that were used. With these challenges in mind and given the results described earlier, in this study we compared three non-linear EEG metrics and their added value to linear metrics to predict treatment outcomes. For this study the dataset from Arns and colleagues (2012) was used which consisted of 90 MDD patients, all treated with rTMS and psychotherapy. In this previous study the slow alpha peak frequency (APF) was the strongest measure in the discriminant analysis, hence in the present study we investigated the non-linear dynamics of the alpha band, and evaluated whether such measures would add to the prediction of treatment outcome, using exactly the same EEG channels and frequency bands as reported in the previous study. The alpha band employed by Arns et al. (2012) was 7–13 Hz for posterior APF and 6–13 Hz for anterior APF. These bandwidths differ from traditionally used frequencies that have studied alpha power. However, it is well known that the APF can range from 4–6 Hz in very young children (Niedermeyer and Da Silva, 1999), to 10 Hz in adults and 7–9 Hz in chronic pain patients (Boord et al., 2008) and Alzheimer patients (Rodriguez et al., 1999). Given that specifically the slow APF values have yielded predictive power (for overview see Arns and Olbrich, 2013) and LRCTs indicate power laws in which slow frequencies dominate the signal (Babloyantz, 1991, Linkenkaer-Hansen et al., 2005, Linkenkaer-Hansen et al., 2001), we used the above frequency band. A second reason was to keep the results comparable to the previous report and contribute to improve coherence between different studies and their interpretation and evaluation of these nonlinear measures (cf. Bravi et al., 2011). Furthermore, in previous analyses as well as in the present study, the focus will be on predicting non-response to treatment. This excludes the effects of placebo response and differential effects of rTMS and psychotherapy. Characterizing non-response to treatment is not only clinically meaningful, but could potentially also result in the development of new treatments based on such biomarkers.
Section snippets
Design
This study was a multi-site open-label study and inclusion criteria were: (1) a primary diagnosis of Depression or Dysthymic disorder according to DSM-IV criteria rated using the MINI (MINI Plus Dutch version 5.0.0) and (2) a Becks Depression Inventory (BDI) score of 14 or higher at enrolment. Exclusion criteria were: previously treated with ECT, epilepsy, wearing a cardiac pacemaker, metal parts in the head, pregnancy and the presence of paroxysmal EEG. All patients signed an informed consent
Results
A total of 90 patients were enrolled (average age: 42.9 yrs, range 19–69 yrs; 49 females and 41 males). There were no differences in any of the clinical outcome measures and demographics between the HF and LF TMS groups nor between the MDD-R and MDD-NR groups.
From the 70 responders, 58 (83%) achieved remission and the remainder demonstrated a more than 50% improvement on the BDI, resulting in 20 MDD-NR and 70 MDD-R and 29 of the 90 patients were unmedicated. For the whole group, the average BDI
Discussion
This study demonstrated that there is no difference in non-linear metrics between depressed patients and healthy participants (no between group difference, nor a correlation with baseline MDD severity), in line with other findings using linear approaches, which have revealed EEG has no diagnostic value for MDD (reviewed in Arns, 2012) or for ADHD in (Arns et al., 2013). However, this study did demonstrate there is utility for non-linear metrics, specifically the Lempel–Ziv Complexity measure,
Acknowledgements
We wish to acknowledge Rosalinde van Ruth, Aukje Bootsma, Anne Schreiber, Vera Kruiver, Irene Giesbers, Rik van Dinteren, Hanneke Friesen, Niels Veth and Desiree Spronk for their help and support in the treatment and data collection of the reported data. Finally we would like to thank our colleagues at Psychologenpraktijk Timmers for sharing the rTMS treatment data, especially Dagmar Timmers.
References (59)
- et al.
Neurophysiological predictors of non-response to rtms in depression
Brain Stimul
(2012) - et al.
Potential differential effects of 9 Hz rTMS and 10 Hz rTMS in the treatment of depression
Brain Stimul
(2010) Evidence for slow brain waves: a dynamical approach
Electroencephalogr Clin Neurophysiol
(1991)- et al.
Human EEG shows long-range temporal correlations of oscillation amplitude in Theta, Alpha and Beta bands across a wide age range
Clin Neurophysiol
(2010) - et al.
Electroencephalographic alpha measures predict therapeutic response to a selective serotonin reuptake inhibitor antidepressant: pre- and post-treatment findings
Biol Psychiatry
(2008) - et al.
Electroencephalographic and perceptual asymmetry differences between responders and nonresponders to an SSRI antidepressant
Biol Psychiatry
(2001) - et al.
Combining fMRI with EEG and MEG in order to relate patterns of brain activity to cognition
Int J Psychophysiol
(2009) - et al.
Post-seizure EEG fractal dimension of first ECT predicts antidepressant response at two weeks
J Affect Disord
(1999) - et al.
Association between BDNF val66met polymorphism and trait depression is mediated via resting EEG alpha band activity
Biol Psychol
(2008) - et al.
What is physiologic complexity and how does it change with aging and disease?
Neurobiol Aging
(2002)
A new method for off-line removal of ocular artifact
Electroencephalography and Clinical Neurophysiology
Change is not always linear: the study of nonlinear and discontinuous patterns of change in psychotherapy
Clin Psychol Rev
The largest lyapunov exponent of the EEG during ECT seizures as a measure of ECT seizure adequacy
Electroencephalogr Clin Neurophysiol
Response and remission criteria in major depression–a validation of current practice
J Psychiatr Res
EEG spectral profile to stage alzheimer’s disease
Clin Neurophysiol
Complex network measures of brain connectivity: uses and interpretations
Neuroimage
Nonlinear dynamical analysis of EEG and MEG: review of an emerging field
Clin Neurophysiol
Characterization of anatomical and functional connectivity in the brain: a complex networks perspective
Int J Psychophysiol
Current source density measures of electroencephalographic alpha predict antidepressant treatment response
Biol Psychiatry
Temporal order of nonlinear dynamics in human brain
Brain Res Brain Res Rev
Interpretation of the Lempel–Ziv complexity measure in the context of biomedical signal analysis
IEEE Transactions Bio-med Engin
EEG-based personalized medicine in ADHD: individual alpha peak frequency as an endophenotype associated with nonresponse
J Neurother
A decade of EEG theta/beta ratio research in ADHD: a meta-analysis
J Atten Disord
Neurofeedback in ADHD and insomnia: vigilance stabilization through sleep spindles and circadian networks
Neurosci Biobehav Rev
Two EEG channels do not make a quantitative EEG (QEEG): a response to Widge, Avery and Zarkowski
Brain Stim
Psychosocial and clinical predictors of response to pharmacotherapy for depression
J Psychiatry Neurosci
Clinically meaningful efficacy and acceptability of low-frequency repetitive transcranial magnetic stimulation (rtms) for treating primary major depression: a meta-analysis of randomized, double-blind and sham-controlled trials
Neuropsychopharmacology
Electroencephalographic slowing and reduced reactivity in neuropathic pain following spinal cord injury
Spinal cord
Review and classification of variability analysis techniques with clinical applications
Biomed Eng Online
Cited by (44)
Non-linear processing and reinforcement learning to predict rTMS treatment response in depression
2024, Psychiatry Research - NeuroimagingA polygenic-informed approach to a predictive EEG signature empowers antidepressant treatment prediction: A proof-of-concept study
2022, European NeuropsychopharmacologyAberrant EEG signal variability at a specific temporal scale in major depressive disorder
2021, Clinical NeurophysiologyRemission of depression is associated with asymmetric hemispheric variation in EEG complexity before antidepressant treatment
2021, Journal of Affective Disorders
- 1
Both authors contributed equally to this manuscript.