Effect of non-invasive spinal cord stimulation in unmedicated adults with major depressive disorder: a pilot randomized controlled trial and induced current flow pattern.

Converging theoretical frameworks suggest a role and a therapeutic potential for spinal interoceptive pathways in major depressive disorder (MDD). Here, we aimed to evaluate the antidepressant effects and tolerability of transcutaneous spinal direct current stimulation (tsDCS) in MDD. This was a double-blind, randomized, sham-controlled, parallel group, pilot clinical trial in unmedicated adults with moderate MDD. Twenty participants were randomly allocated (1:1 ratio) to receive "active" 2.5 mA or "sham" anodal tsDCS sessions with a thoracic (anode; T10)/right shoulder (cathode) electrode montage 3 times/week for 8 weeks. Change in depression severity (MADRS) scores (prespecified primary outcome) and secondary clinical outcomes were analyzed with ANOVA models. An E-Field model was generated using the active tsDCS parameters. Compared to sham (n = 9), the active tsDCS group (n = 10) showed a greater baseline to endpoint decrease in MADRS score with a large effect size (-14.6 ± 2.5 vs. -21.7 ± 2.3, p = 0.040, d = 0.86). Additionally, compared to sham, active tsDCS induced a greater decrease in MADRS "reported sadness" item (-1.8 ± 0.4 vs. -3.2 ± 0.4, p = 0.012), and a greater cumulative decrease in pre/post tsDCS session diastolic blood pressure change from baseline to endpoint (group difference: 7.9 ± 3.7 mmHg, p = 0.039). Statistical trends in the same direction were observed for MADRS "pessimistic thoughts" item and week-8 CGI-I scores. No group differences were observed in adverse events (AEs) and no serious AEs occurred. The current flow simulation showed electric field at strength within the neuromodulation range (max. ~0.45 V/m) reaching the thoracic spinal gray matter. The results from this pilot study suggest that tsDCS is feasible, well-tolerated, and shows therapeutic potential in MDD. This work also provides the initial framework for the cautious exploration of non-invasive spinal cord neuromodulation in the context of mental health research and therapeutics. The underlying mechanisms warrant further investigation. Clinicaltrials.gov registration: NCT03433339 URL: https://clinicaltrials.gov/ct2/show/NCT03433339 .

CLoSES: A platform for closed-loop intracranial stimulation in humans.

Targeted interrogation of brain networks through invasive brain stimulation has become an increasingly important research tool as well as therapeutic modality. The majority of work with this emerging capability has been focused on open-loop approaches. Closed-loop techniques, however, could improve neuromodulatory therapies and research investigations by optimizing stimulation approaches using neurally informed, personalized targets. Implementing closed-loop systems is challenging particularly with regard to applying consistent strategies considering inter-individual variability. In particular, during intracranial epilepsy monitoring, where much of this research is currently progressing, electrodes are implanted exclusively for clinical reasons. Thus, detection and stimulation sites must be participant- and task-specific. The system must run in parallel with clinical systems, integrate seamlessly with existing setups, and ensure safety features are in place. In other words, a robust, yet flexible platform is required to perform different tests with a single participant and to comply with clinical requirements. In order to investigate closed-loop stimulation for research and therapeutic use, we developed a Closed-Loop System for Electrical Stimulation (CLoSES) that computes neural features which are then used in a decision algorithm to trigger stimulation in near real-time. To summarize CLoSES, intracranial electroencephalography (iEEG) signals are acquired, band-pass filtered, and local and network features are continuously computed. If target features are detected (e.g. above a preset threshold for a certain duration), stimulation is triggered. Not only could the system trigger stimulation while detecting real-time neural features, but we incorporated a pipeline wherein we used an encoder/decoder model to estimate a hidden cognitive state from the neural features. CLoSES provides a flexible platform to implement a variety of closed-loop experimental paradigms in humans. CLoSES has been successfully used with twelve patients implanted with depth electrodes in the epilepsy monitoring unit. During cognitive tasks (N=5), stimulation in closed loop modified a cognitive hidden state on a trial by trial basis. Sleep spindle oscillations (N=6) and sharp transient epileptic activity (N=9) were detected in near real-time, and stimulation was applied during the event or at specified delays (N=3). In addition, we measured the capabilities of the CLoSES system. Total latency was related to the characteristics of the event being detected, with tens of milliseconds for epileptic activity and hundreds of milliseconds for spindle detection. Stepwise latency, the actual duration of each continuous step, was within the specified fixed-step duration and increased linearly with the number of channels and features. We anticipate that probing neural dynamics and interaction between brain states and stimulation responses with CLoSES will lead to novel insights into the mechanism of normal and pathological brain activity, the discovery and evaluation of potential electrographic biomarkers of neurological and psychiatric disorders, and the development and testing of patient-specific stimulation targets and control signals before implanting a therapeutic device.

Decoding Hidden Cognitive States From Behavior and Physiology Using a Bayesian Approach.

Cognitive processes, such as learning and cognitive flexibility, are both difficult to measure and to sample continuously using objective tools because cognitive processes arise from distributed, high-dimensional neural activity. For both research and clinical applications, that dimensionality must be reduced. To reduce dimensionality and measure underlying cognitive processes, we propose a modeling framework in which a cognitive process is defined as a low-dimensional dynamical latent variable-called a cognitive state, which links high-dimensional neural recordings and multidimensional behavioral readouts. This framework allows us to decompose the hard problem of modeling the relationship between neural and behavioral data into separable encoding-decoding approaches. We first use a state-space modeling framework, the behavioral decoder, to articulate the relationship between an objective behavioral readout (e.g., response times) and cognitive state. The second step, the neural encoder, involves using a generalized linear model (GLM) to identify the relationship between the cognitive state and neural signals, such as local field potential (LFP). We then use the neural encoder model and a Bayesian filter to estimate cognitive state using neural data (LFP power) to generate the neural decoder. We provide goodness-of-fit analysis and model selection criteria in support of the encoding-decoding result. We apply this framework to estimate an underlying cognitive state from neural data in human participants (N=8) performing a cognitive conflict task. We successfully estimated the cognitive state within the 95% confidence intervals of that estimated using behavior readout for an average of 90% of task trials across participants. In contrast to previous encoder-decoder models, our proposed modeling framework incorporates LFP spectral power to encode and decode a cognitive state. The framework allowed us to capture the temporal evolution of the underlying cognitive processes, which could be key to the development of closed-loop experiments and treatments.

Who May Benefit From On-Demand Control of Deep Brain Stimulation? Noninvasive Evaluation of Parkinson Patients.

BACKGROUND: In closed-loop on-demand control (ODC) of deep brain stimulation (DBS), stimulation is applied only when symptoms appear. Following stimulation of a fixed duration, DBS is switched off until the symptoms reappear. By repeating these demand-driven cycles, the amount of stimulation delivered can be decreased, thereby reducing DBS side-effects and improving battery-life of the pulse-generator. This article introduces Ro metric for quantification of degree of benefit of ODC and explores candidate selection in tremor-dominant Parkinson's disease (PD). METHOD: The study was performed on nine PD patients previously implanted with Medtronic DBS systems. Accelerometer sensor was placed on the tremor-dominant hand to detect onset of tremor. Fixed duration of stimulation (DS) of 20-80 sec was applied. Once the tremor was observed, stimulation was switched on. These trials were repeated during resting, postural, and kinetic conditions. Ro metric was calculated as the ratio of stimulation-off tremor-free period to the DS. Ro calculated at different DS were compared for each patient. RESULTS: We found that for each patient, Ro varied with DS and an optimal DS* gave a higher percentage of stimulation-off time. Average Ro at DS* varied from 0.554 to 4.24 for eight patients giving 35%-80% stimulation-off time. CONCLUSIONS: Ro values can be used for selection of optimal DS* in ODC. Three of nine patients were found to be tremor-free without stimulation for >50% of total time with even up to 80% in one patient. Patients with low Ro may not benefit from ODC in DBS, where the trade-off between having side-effects and using ODC system will need to be assessed.

Anxious/depressed individuals exhibit disrupted frontotemporal synchrony during cognitive conflict encoding.

Anxiety and depressive disorders are associated with cognitive control deficits, yet their underlying neural mechanisms remain poorly understood. Here, we used high-resolution stereotactic EEG (sEEG) to determine how anxiety and/or depression modulates neural and behavioral responses when cognitive control is engaged in individuals with medically refractory epilepsy undergoing sEEG monitoring for surgical evaluation. We analyzed sEEG data recorded from frontotemporal regions of 29 participants (age range: 19-55, mean age: 35.5, female: 16/29) while they performed a Multi-Source Interference Task (MSIT) designed to elicit cognitive conflict. Neurobehavioral interviews, symptom rating scales, and clinical documentation were used to categorize participants as demonstrating anxiety and/or depression symptoms (A/D, n=13) or as epilepsy controls (EC, n=16). Generalized linear mixed-effects (GLME) models were used to analyze behavioral and neural data. Models of oscillatory power were used to identify brain regions within conflict-encoding networks in which coherence and phase locking values (PLV) were examined in A/D and EC. A/D participants demonstrated a greater conflict effect (response time slowing with higher cognitive load), without impairment in response time (RT) or accuracy compared to EC. A/D participants also showed significantly enhanced conflict-evoked theta (4-8Hz) and alpha (8-15Hz) power in the dorsolateral prefrontal cortex (dlPFC) and amygdala as well as widespread broadband activity in the lateral temporal lobe (LTL) compared to EC. Additionally, theta coherence and PLV between dlPFC-LTL and dlPFC-amygdala were reduced by conflict in A/D. Our findings suggest individuals with anxiety/depression symptoms exhibit heightened frontotemporal oscillatory activity and disrupted frontotemporal synchrony during cognitive conflict encoding, which may indicate a greater need for cognitive resources due to ineffective cognitive processing. These results highlight a potential role of frontotemporal circuits in conflict encoding that are altered in anxiety/depression, and may further inform future therapeutic interventions aimed at enhancing cognitive control in these populations.

Modular pipeline for reconstruction and localization of implanted intracranial ECoG and sEEG electrodes.

Implantation of electrodes in the brain has been used as a clinical tool for decades to stimulate and record brain activity. As this method increasingly becomes the standard of care for several disorders and diseases, there is a growing need to quickly and accurately localize the electrodes once they are placed within the brain. We share here a protocol pipeline for localizing electrodes implanted in the brain, which we have applied to more than 260 patients, that is accessible to multiple skill levels and modular in execution. This pipeline uses multiple software packages to prioritize flexibility by permitting multiple different parallel outputs while minimizing the number of steps for each output. These outputs include co-registered imaging, electrode coordinates, 2D and 3D visualizations of the implants, automatic surface and volumetric localizations of the brain regions per electrode, and anonymization and data sharing tools. We demonstrate here some of the pipeline's visualizations and automatic localization algorithms which we have applied to determine appropriate stimulation targets, to conduct seizure dynamics analysis, and to localize neural activity from cognitive tasks in previous studies. Further, the output facilitates the extraction of information such as the probability of grey matter intersection or the nearest anatomic structure per electrode contact across all data sets that go through the pipeline. We expect that this pipeline will be a useful framework for researchers and clinicians alike to localize implanted electrodes in the human brain.

Closed-loop enhancement and neural decoding of cognitive control in humans.

Deficits in cognitive control-that is, in the ability to withhold a default pre-potent response in favour of a more adaptive choice-are common in depression, anxiety, addiction and other mental disorders. Here we report proof-of-concept evidence that, in participants undergoing intracranial epilepsy monitoring, closed-loop direct stimulation of the internal capsule or striatum, especially the dorsal sites, enhances the participants' cognitive control during a conflict task. We also show that closed-loop stimulation upon the detection of lapses in cognitive control produced larger behavioural changes than open-loop stimulation, and that task performance for single trials can be directly decoded from the activity of a small number of electrodes via neural features that are compatible with existing closed-loop brain implants. Closed-loop enhancement of cognitive control might remediate underlying cognitive deficits and aid the treatment of severe mental disorders.

Classification of Seizure Termination Patterns using Deep Learning on intracranial EEG.

Seizure termination has received significantly less attention than initiation and propagation and consequently, remains a poorly understood phase of seizure evolution. Yet, its study may have a significant impact on the development of efficient interventional approaches, i.e., it may be critical for the design of treatments that induce or reproduce termination mechanisms that are triggered in self-terminating seizures. In this work, we aim to study temporal and spectral features of intracranial EEG (iEEG) during epileptic seizures to find time-frequency signatures that can predict the termination patterns. We propose a deep learning model for classification of multi channel iEEG epileptic seizure termination pattern into burst suppression and continuous bursting. We decompose the raw time series seizure data into time-frequency maps using Morlet Wavelet Transform. A Convolution Neural Network (CNN) is then trained on cross-patient time-frequency maps to classify the seizure termination patterns. For evaluation of classification performance, we compared the proposed method with k-Nearest Neighbour (k-NN). The CNN is shown to achieve an accuracy of 90 % and precision of 92 % as compared to 70% and 72% accuracy and precision achieved with the k-NN respectively. The proposed model is thus able to capture the temporal and spatial patterns which results in high performance of the classifier. This method of classification can be used to predict how a particular seizure will end and can potentially inform seizure management and treatment. Clinical relevance- This method establishes a model that can be used to classify seizure termination patterns with an accuracy of 90 % which can assist in better treatment of epilepsy patients.

Neural mass model-based study of frontal-temporal theta oscillations in human subjects during the performance of a cognitive control task.

Cognitive control, the ability to rapidly shift one's attention and behavioral strategy in response to environmental changes, is often compromised across psychiatric disorders. One of the well-validated behavioral paradigms for tapping into the cognitive control circuits is a cognitive interference task, where subjects must suppress a natural response to follow a less intuitive rule. Slower response times on these tasks indicate difficulty exerting control to overcome response conflict. Conflict evokes robust electrophysiological signatures, such as theta (4-8 Hz) oscillations in the prefrontal cortex (PFC). However, the underlying neural mechanisms of conflict-evoked theta oscillations in the PFC are not clear. The objective of this work is to use a neural mass model (NMM) to find feasible cortical networks generating theta oscillations during conflict processing in human subjects. We used intracranial EEG (iEEG) recorded from dorsolateral PFC (dIPFC) and lateral temporal lobe (LTL) of human subjects with intractable epilepsy undergoing invasive monitoring, while they performed a multi-source interference task (MSIT). We used a dynamic causal modeling (DCM) framework to simulate dIPFC-LTL theta using a Jansen-Rit NMM. We found significant evidence for an LTL input into the dlPFC during the initial 500 ms of conflict processing compared to a bidirectional connection between the dlPFC and LTL. We conclude that a neural mass modeling framework can be used to elucidate candidate mechanisms of neural oscillations underlying conflict resolution in human subjects. Clinical Relevance- This can be used to find feasible target mechanisms for designing therapy in patients with compromised cognitive control. This framework can also be expanded to serve as an in-silico test bed for designing and testing neuromodulatory interventions such as electrical stimulation for improving cognitive control in mood/anxiety disorders.

Stochastic modeling of the neuronal activity in the subthalamic nucleus and model parameter identification from Parkinson patient data.

Several stochastic models, with various degrees of complexity, have been proposed to model the neuronal activity from different parts of the human brain. In this article, we use a simple Ornstein-Uhlenbeck process (OUP) to model the spike activity recorded from the subthalamic nucleus of patients suffering from Parkinson's disease at the time of implantation of the electrodes for deep brain stimulation. From the recorded data, which contains information about the spike times of a single neuron, we identify and extract the model parameters of the OUP. We then use these parameters to numerically simulate the inter-spike intervals and the voltage across the neuron membrane. We finally assess how well the proposed mathematical model fits to the measured data and compare it with other commonly adopted stochastic models. We show an excellent agreement between the computer-generated data according to the OUP model and the measured one, as well as the superiority of the OUP model when compared to the Poisson process model and the random walk model; thus, establishing the validity of the OUP as a simple yet biologically plausible model of the neuronal activity recorded from the subthalamic nucleus of Parkinson's disease patients.

Case Report of Dual-Site Neurostimulation and Chronic Recording of Cortico-Striatal Circuitry in a Patient With Treatment Refractory Obsessive Compulsive Disorder.

Psychiatric disorders are increasingly understood as dysfunctions of hyper- or hypoconnectivity in distributed brain circuits. A prototypical example is obsessive compulsive disorder (OCD), which has been repeatedly linked to hyper-connectivity of cortico-striatal-thalamo-cortical (CSTC) loops. Deep brain stimulation (DBS) and lesions of CSTC structures have shown promise for treating both OCD and related disorders involving over-expression of automatic/habitual behaviors. Physiologically, we propose that this CSTC hyper-connectivity may be reflected in high synchrony of neural firing between loop structures, which could be measured as coherent oscillations in the local field potential (LFP). Here we report the results from the pilot patient in an Early Feasibility study (https://clinicaltrials.gov/ct2/show/NCT03184454) in which we use the Medtronic Activa PC+ S device to simultaneously record and stimulate in the supplementary motor area (SMA) and ventral capsule/ventral striatum (VC/VS). We hypothesized that frequency-mismatched stimulation should disrupt coherence and reduce compulsive symptoms. The patient reported subjective improvement in OCD symptoms and showed evidence of improved cognitive control with the addition of cortical stimulation, but these changes were not reflected in primary rating scales specific to OCD and depression, or during blinded cortical stimulation. This subjective improvement was correlated with increased SMA and VC/VS coherence in the alpha, beta, and gamma bands, signals which persisted after correcting for stimulation artifacts. We discuss the implications of this research, and propose future directions for research in network modulation in OCD and more broadly across psychiatric disorders.

Consistent linear and non-linear responses to invasive electrical brain stimulation across individuals and primate species with implanted electrodes.

BACKGROUND: Electrical neuromodulation via implanted electrodes is used in treating numerous neurological disorders, yet our knowledge of how different brain regions respond to varying stimulation parameters is sparse. OBJECTIVE/HYPOTHESIS: We hypothesized that the neural response to electrical stimulation is both region-specific and non-linearly related to amplitude and frequency. METHODS: We examined evoked neural responses following 400 ms trains of 10-400 Hz electrical stimulation ranging from 0.1 to 10 mA. We stimulated electrodes implanted in cingulate cortex (dorsal anterior cingulate and rostral anterior cingulate) and subcortical regions (nucleus accumbens, amygdala) of non-human primates (NHP, N = 4) and patients with intractable epilepsy (N = 15) being monitored via intracranial electrodes. Recordings were performed in prefrontal, subcortical, and temporal lobe locations. RESULTS: In subcortical regions as well as dorsal and rostral anterior cingulate cortex, response waveforms depended non-linearly on frequency (Pearson's linear correlation r < 0.39), but linearly on current (r > 0.58). These relationships between location, and input-output characteristics were similar in homologous brain regions with average Pearson's linear correlation values r > 0.75 between species and linear correlation values between participants r > 0.75 across frequency and current values per brain region. Evoked waveforms could be described by three main principal components (PCs) which allowed us to successfully predict response waveforms across individuals and across frequencies using PC strengths as functions of current and frequency using brain region specific regression models. CONCLUSIONS: These results provide a framework for creation of an atlas of input-output relationships which could be used in the principled selection of stimulation parameters per brain region.

COMPASS: An Open-Source, General-Purpose Software Toolkit for Computational Psychiatry.

Mathematical modeling of behavior during a psychophysical task, referred to as "computational psychiatry," could greatly improve our understanding of mental disorders. One barrier to the broader adoption of computational methods, is that they often require advanced statistical modeling and mathematical skills. Biological and behavioral signals often show skewed or non-Gaussian distributions, and very few toolboxes and analytical platforms are capable of processing such signal categories. We developed the Computational Psychiatry Adaptive State-Space (COMPASS) toolbox, an open-source MATLAB-based software package. This toolbox is easy to use and capable of integrating signals with a variety of distributions. COMPASS has the tools to process signals with continuous-valued and binary measurements, or signals with incomplete-missing or censored-measurements, which makes it well-suited for processing those signals captured during a psychophysical task. After specifying a few parameters in a small set of user-friendly functions, COMPASS allows users to efficiently apply a wide range of computational behavioral models. The model output can be analyzed as an experimental outcome or used as a regressor for neural data and can also be tested using the goodness-of-fit measurement. Here, we demonstrate that COMPASS can replicate two computational behavioral analyses from different groups. COMPASS replicates and can slightly improve on the original modeling results. We also demonstrate the use of COMPASS application in a censored-data problem and compare its performance result with naïve estimation methods. This flexible, general-purpose toolkit should accelerate the use of computational modeling in psychiatric neuroscience.

A neural mass model to predict electrical stimulation evoked responses in human and non-human primate brain.

OBJECTIVE: Deep brain stimulation (DBS) is a valuable tool for ameliorating drug resistant pathologies such as movement disorders and epilepsy. DBS is also being considered for complex neuro-psychiatric disorders, which are characterized by high variability in symptoms and slow responses that hinder DBS setting optimization. The objective of this work was to develop an in silico platform to examine the effects of electrical stimulation in regions neighboring a stimulated brain region. APPROACH: We used the Jansen-Rit neural mass model of single and coupled nodes to simulate the response to a train of electrical current pulses at different frequencies (10-160 Hz) of the local field potential recorded in the amygdala and cortical structures in human subjects and a non-human primate. RESULTS: We found that using a single node model, the evoked responses could be accurately modeled following a narrow range of stimulation frequencies. Including a second coupled node increased the range of stimulation frequencies whose evoked responses could be efficiently modeled. Furthermore, in a chronic recording from a non-human primate, features of the in vivo evoked response remained consistent for several weeks, suggesting that model re-parameterization for chronic stimulation protocols would be infrequent. SIGNIFICANCE: Using a model of neural population activity, we reproduced the evoked response to cortical and subcortical stimulation in human and non-human primate. This modeling framework provides an environment to explore, safely and rapidly, a wide range of stimulation settings not possible in human brain stimulation studies. The model can be trained on a limited dataset of stimulation responses to develop an optimal stimulation strategy for an individual patient.

Treating refractory mental illness with closed-loop brain stimulation: Progress towards a patient-specific transdiagnostic approach.

Mental disorders are a leading cause of disability, morbidity, and mortality among civilian and military populations. Most available treatments have limited efficacy, particularly in disorders where symptoms vary over relatively short time scales. Targeted modulation of neural circuits, particularly through open-loop deep brain stimulation (DBS), showed initial promise but has failed in blinded clinical trials. We propose a new approach, based on targeting neural circuits linked to functional domains that cut across diagnoses. Through that framework, which includes measurement of patients using six psychophysical tasks, we seek to develop a closed-loop DBS system that corrects dysfunctional activity in brain circuits underlying those domains. We present convergent preliminary evidence from functional neuroimaging, invasive human electrophysiology, and human brain stimulation experiments suggesting that this approach is feasible. Using the Emotional Conflict Resolution (ECR) task as an example, we show that emotion-related networks can be identified and modulated in individual patients. Invasive and non-invasive methodologies both identify a network between prefrontal cortex, cingulate cortex, insula, and amygdala. Further, stimulation in cingulate and amygdala changes patients' performance in ways that are linked to the task's emotional content. We present preliminary statistical models that predict this change and allow us to track it at a single-trial level. As these diagnostic and modeling strategies are refined and embodied in an implantable device, they offer the prospect of a new approach to psychiatric treatment and its accompanying neuroscience.

A study of the dynamics of seizure propagation across micro domains in the vicinity of the seizure onset zone.

OBJECTIVE: The use of micro-electrode arrays to measure electrical activity from the surface of the brain is increasingly being investigated as a means to improve seizure onset zone (SOZ) localization. In this work, we used a multivariate autoregressive model to determine the evolution of seizure dynamics in the [Formula: see text] Hz high frequency band across micro-domains sampled by such micro-electrode arrays. We showed that a directed transfer function (DTF) can be used to estimate the flow of seizure activity in a set of simulated micro-electrode data with known propagation pattern. APPROACH: We used seven complex partial seizures recorded from four patients undergoing intracranial monitoring for surgical evaluation to reconstruct the seizure propagation pattern over sliding windows using a DTF measure. MAIN RESULTS: We showed that a DTF can be used to estimate the flow of seizure activity in a set of simulated micro-electrode data with a known propagation pattern. In general, depending on the location of the micro-electrode grid with respect to the clinical SOZ and the time from seizure onset, ictal propagation changed in directional characteristics over a 2-10 s time scale, with gross directionality limited to spatial dimensions of approximately [Formula: see text]. It was also seen that the strongest seizure patterns in the high frequency band and their sources over such micro-domains are more stable over time and across seizures bordering the clinically determined SOZ than inside. SIGNIFICANCE: This type of propagation analysis might in future provide an additional tool to epileptologists for characterizing epileptogenic tissue. This will potentially help narrowing down resection zones without compromising essential brain functions as well as provide important information about targeting anti-epileptic stimulation devices.

Towards closed-loop deep brain stimulation: decision tree-based essential tremor patient's state classifier and tremor reappearance predictor.

Deep Brain Stimulation (DBS) is a surgical procedure to treat some progressive neurological movement disorders, such as Essential Tremor (ET), in an advanced stage. Current FDA-approved DBS systems operate open-loop, i.e., their parameters are unchanged over time. This work develops a Decision Tree (DT) based algorithm that, by using non-invasively measured surface EMG and accelerometer signals as inputs during DBS-OFF periods, classifies the ET patient's state and then predicts when tremor is about to reappear, at which point DBS is turned ON again for a fixed amount of time. The proposed algorithm achieves an overall accuracy of 93.3% and sensitivity of 97.4%, along with 2.9% false alarm rate. Also, the ratio between predicted tremor delay and the actual detected tremor delay is about 0.93, indicating that tremor prediction is very close to the instant where tremor actually reappeared.

Determination of seizure propagation across microdomains using spectral measures of causality.

The use of microelectrode arrays to measure electrical activity from the surface of the brain is increasingly being investigated as a means to improve seizure focus localization. In this work, we determine seizure propagation across microdomains sampled by such microelectrode arrays and compare the results using two widely used frequency domain measures of causality, namely the partial directed coherence and the directed direct transfer function. We show that these two measures produce very similar propagation patterns for simulated microelectrode activity over a relatively smaller number of channels. However as the number of channels increases, partial directed coherence produces better estimates of the actual propagation pattern. Additionally, we apply these two measures to determine seizure propagation over microelectrode arrays measured from a patient undergoing intracranial monitoring for seizure focus localization and find very similar patterns which also agree with a threshold based reconstruction during seizure onset.

On modeling the neuronal activity in movement disorder patients by using the Ornstein Uhlenbeck process.

Mathematical models of the neuronal activity in the affected brain regions of Essential Tremor (ET) and Parkinson's Disease (PD) patients could shed light into the underlying pathophysiology of these diseases, which in turn could help develop personalized treatments including adaptive Deep Brain Stimulation (DBS). In this paper, we use an Ornstein Uhlenbeck Process (OUP) to model the neuronal spiking activity recorded from the brain of ET and PD patients during DBS stereotactic surgery. The parameters of the OUP are estimated based on Inter Spike Interval (ISI) measurements, i.e., the time interval between two consecutive neuronal firings, by means of the Fortet Integral Equation (FIE). The OUP model parameters identified with the FIE method (OUP-FIE) are then used to simulate the ISI distribution resulting from the OUP. Other widely used neuronal activity models, such as the Poisson Process (PP), the Brownian Motion (BM), and the OUP whose parameters are extracted by matching the first two moments of the ISI (OUP-MOM), are also considered. To quantify how close the simulated ISI distribution is to the measured ISI distribution, the Integral Square Error (ISE) criterion is adopted. Amongst all considered stochastic processes, the ISI distribution generated by the OUP-FIE method is shown to produce the least ISE. Finally, a directional Wilcoxon signed rank test is used to show statistically significant reduction in the ISE value obtained from the OUP-FIE compared to the other stochastic processes.