Brain-based biomarker of treatment recovery for deep brain stimulation for treatment-resistant depression
Brain-based biomarker of treatment recovery for deep brain stimulation for treatment-resistant depression
Sankaraleengam Alagapan1, Stephen Heisig2, Ki Sueng Choi2, Allison C. Waters2, Ashan Veerakumar3, Vineet Tiruvadi1,3, Mosadoluwa Obatusin2, Tanya Nauvel2, Jungho Cha2, Andrea Crowell3, Martijn Figee2, Patricio Riva-Posse3, Robert Butera1, Helen Mayberg2, Christopher Rozell1
Georgia Institute of Technology1, Icahn School of Medicine at Mount Sinai2, Emory University3
Deep brain stimulation (DBS) of the subcallosal cingulate cortex (SCC) has shown effectiveness in treating patients with treatment-resistant depression (TRD). However, inter-individual variability in recovery and intra-individual variability in hard-to-measure symptoms hinder clinical decision-making, including adjustment of stimulation parameters. To overcome this obstacle, which also limits the scalability, we aim to identify electrophysiological markers of recovery from local field potential (LFP) recordings.
10 participants with TRD were recruited for this study. LFPs were acquired weekly from 6 participants during 24 weeks of SCC DBS therapy using implanted pulse generators (Activa PC+S, Medtronic Inc, USA). 4 participants were excluded for missing LFP data and LFP data contaminated with artifacts. Spectral features were extracted from the LFPs, and a neural network classifier was used to differentiate the beginning and the end of the 24 weeks. A novel explainable artificial intelligence (xAI) tool was used to identify a biomarker that explained the feature changes that contributed to classifier performance. Depression severity was measured using the Hamilton Depression Rating Scale (HDRS). Instead of tracking weekly HDRS scores, the HDRS was binarized into ‘sick’ and ‘stable response’ states and compared to the states derived from the biomarker.
The neural network classifier was able to differentiate the beginning and the end of 24 weeks (Area under ROC curve: 0.87 ± 0.09) from the LFP features. The features that contributed to classifier performance were identified to be alpha (8 – 12 Hz), beta (12 – 30 Hz) and gamma (30 – 40 Hz) bands. The binarized states (‘sick’ and ‘stable response’) derived from the biomarker predicted HDRS-derived states with high accuracy (Area under ROC curve: 0.94 ± 0.04). In addition, the biomarker predicted relapse in a participant whose LFP data was not used for training the different machine learning models. Notably, the biomarker exhibited changes in the week after the stimulation dose was adjusted.
The results indicate that the biomarker derived from LFP satisfies two essential requirements - tracking relevant changes in disease state and responding to intervention. Thus, the biomarker can potentially be utilized to inform clinical decision-making thereby improving the scalability of this promising approach beyond a few highly specialized centers and reducing variability in outcomes.