Model-based deconstruction of cortical evoked potentials generated by subthalamic nucleus deep brain stimulation

Neural pathway activation in the subthalamic region depends on stimulation polarity

Deep brain stimulation (DBS) is an effective treatment for Parkinson's disease (PD); however, there is limited understanding of which subthalamic pathways are recruited in response to stimulation. Here, by focusing on the polarity of the stimulus waveform (cathodic vs. anodic), our goal was to elucidate biophysical mechanisms that underlie electrical stimulation in the human brain. In clinical studies, cathodic stimulation more easily triggers behavioral responses, but anodic DBS broadens the therapeutic window. This suggests that neural pathways involved respond preferentially depending on stimulus polarity. To experimentally compare the activation of therapeutically relevant pathways during cathodic and anodic subthalamic nucleus (STN) DBS, pathway activation was quantified by measuring evoked potentials resulting from antidromic or orthodromic activation in 15 PD patients undergoing DBS implantation. Cortical evoked potentials (cEP) were recorded using subdural electrocorticography, DBS local evoked potentials (DLEP) were recorded from non-stimulating contacts and EMG activity was recorded from arm and face muscles. We measured: 1) the amplitude of short-latency cEP, previously demonstrated to reflect activation of the cortico-STN hyperdirect pathway, 2) DLEP amplitude thought to reflect activation of STN-globus pallidus (GP) pathway, and 3) amplitudes of very short-latency cEP and motor evoked potentials (mEP) for activation of cortico-spinal/bulbar tract (CSBT). We constructed recruitment and strength-duration curves for each EP/pathway to compare the excitability for different stimulation polarities. We compared experimental data with the most advanced DBS computational models. Our results provide experimental evidence that subcortical cathodic and anodic stimulation activate the same pathways in the STN region and that cathodic stimulation is in general more efficient. However, relative efficiency varies for different pathways so that anodic stimulation is the least efficient in activating CSBT, more efficient in activating the HDP and as efficient as cathodic in activating STN-GP pathway. Our experiments confirm biophysical model predictions regarding neural activations in the central nervous system and provide evidence that stimulus polarity has differential effects on passing axons, terminal synapses, and local neurons. Comparison of experimental results with clinical DBS studies provides further evidence that the hyperdirect pathway may be involved in the therapeutic mechanisms of DBS.

Personalized depth-specific neuromodulation of the human primary motor cortex via ultrasound

Non-invasive brain stimulation has the potential to boost neuronal plasticity in the primary motor cortex (M1), but it remains unclear whether the stimulation of both superficial and deep layers of the human motor cortex can effectively promote M1 plasticity. Here, we leveraged transcranial ultrasound stimulation (TUS) to precisely target M1 circuits at depths of approximately 5 mm and 16 mm from the cortical surface. Initially, we generated computed tomography images from each participant's individual anatomical magnetic resonance images (MRI), which allowed for the generation of accurate acoustic simulations. This process ensured that personalized TUS was administered exactly to the targeted depths within M1 for each participant. Using long-term depression and long-term potentiation (LTD/LTP) theta-burst stimulation paradigms, we examined whether TUS over distinct depths of M1 could induce LTD/LTP plasticity. Our findings indicated that continuous theta-burst TUS-induced LTD-like plasticity with both superficial and deep M1 stimulation, persisting for at least 30 min. In comparison, sham TUS did not significantly alter M1 excitability. Moreover, intermittent theta-burst TUS did not result in the induction of LTP- or LTD-like plasticity with either superficial or deep M1 stimulation. These findings suggest that the induction of M1 plasticity can be achieved with ultrasound stimulation targeting distinct depths of M1, which is contingent on the characteristics of TUS. KEY POINTS: The study integrated personalized transcranial ultrasound stimulation (TUS) with electrophysiology to determine whether TUS targeting superficial and deep layers of the human motor cortex (M1) could elicit long-term depression (LTD) or long-term potentiation (LTP) plastic changes. Utilizing acoustic simulations derived from individualized pseudo-computed tomography scans, we ensured the precision of TUS delivery to the intended M1 depths for each participant. Continuous theta-burst TUS targeting both the superficial and deep layers of M1 resulted in the emergence of LTD-like plasticity, lasting for at least 30 min. Administering intermittent theta-burst TUS to both the superficial and deep layers of M1 did not lead to the induction of LTP- or LTD-like plastic changes. We suggest that theta-burst TUS targeting distinct depths of M1 can induce plasticity, but this effect is dependent on specific TUS parameters.

Delayed effective connectivity characterizes the epileptogenic zone during stereo-EEG

OBJECTIVE

Single-pulse electrical stimulations (SPES) can elicit normal and abnormal responses that might characterize the epileptogenic zone, including spikes, high-frequency oscillations and cortico-cortical evoked potentials (CCEPs). In this study, we investigate their association with the epileptogenic zone during stereoelectroencephalography (SEEG) in 28 patients with refractory focal epilepsy.

METHODS

Characteristics of CCEPs (distance-corrected or -uncorrected latency, amplitude and the connectivity index) and the occurrence of spikes and ripples were assessed. Responses within the epileptogenic zone and within the non-involved zone were compared using receiver operating characteristics curves and analysis of variance (ANOVA) either in all patients, patients with well-delineated epileptogenic zone, and patients older than 15 years old.

RESULTS

We found an increase in distance-corrected CCEPs latency after stimulation within the epileptogenic zone (area under the curve = 0.71, 0.72, 0.70, ANOVA significant after false discovery rate correction).

CONCLUSIONS

The increased distance-corrected CCEPs latency suggests that neuronal propagation velocity is altered within the epileptogenic network. This association might reflect effective connectivity changes at cortico-cortical or cortico-subcortico-cortical levels. Other responses were not associated with the epileptogenic zone, including the CCEPs amplitude, the connectivity index, the occurrences of induced ripples and spikes. The discrepancy with previous descriptions may be explained by different spatial brain sampling between subdural and depth electrodes.

SIGNIFICANCE

Increased distance-corrected CCEPs latency, indicating delayed effective connectivity, characterizes the epileptogenic zone. This marker could be used to help tailor surgical resection limits after SEEG.

Signatures of Electrical Stimulation Driven Network Interactions in the Human Limbic System

Stimulation-evoked signals are starting to be used as biomarkers to indicate the state and health of brain networks. The human limbic network, often targeted for brain stimulation therapy, is involved in emotion and memory processing. Previous anatomic, neurophysiological, and functional studies suggest distinct subsystems within the limbic network (Rolls, 2015). Studies using intracranial electrical stimulation, however, have emphasized the similarities of the evoked waveforms across the limbic network. We test whether these subsystems have distinct stimulation-driven signatures. In eight patients (four male, four female) with drug-resistant epilepsy, we stimulated the limbic system with single-pulse electrical stimulation. Reliable corticocortical evoked potentials (CCEPs) were measured between hippocampus and the posterior cingulate cortex (PCC) and between the amygdala and the anterior cingulate cortex (ACC). However, the CCEP waveform in the PCC after hippocampal stimulation showed a unique and reliable morphology, which we term the "limbic Hippocampus-Anterior nucleus of the thalamus-Posterior cingulate, HAP-wave." This limbic HAP-wave was visually distinct and separately decoded from the CCEP waveform in ACC after amygdala stimulation. Diffusion MRI data show that the measured end points in the PCC overlap with the end points of the parolfactory cingulum bundle rather than the parahippocampal cingulum, suggesting that the limbic HAP-wave may travel through fornix, mammillary bodies, and the anterior nucleus of the thalamus (ANT). This was further confirmed by stimulating the ANT, which evoked the same limbic HAP-wave but with an earlier latency. Limbic subsystems have unique stimulation-evoked signatures that may be used in the future to help network pathology diagnosis.SIGNIFICANCE STATEMENT The limbic system is often compromised in diverse clinical conditions, such as epilepsy or Alzheimer's disease, and characterizing its typical circuit responses may provide diagnostic insight. Stimulation-evoked waveforms have been used in the motor system to diagnose circuit pathology. We translate this framework to limbic subsystems using human intracranial stereo EEG (sEEG) recordings that measure deeper brain areas. Our sEEG recordings describe a stimulation-evoked waveform characteristic to the memory and spatial subsystem of the limbic network that we term the "limbic HAP-wave." The limbic HAP-wave follows anatomic white matter pathways from hippocampus to thalamus to the posterior cingulum and shows promise as a distinct biomarker of signaling in the human brain memory and spatial limbic network.

Myogenic and cortical evoked potentials vary as a function of stimulus pulse geometry delivered in the subthalamic nucleus of Parkinson's disease patients

Introduction

The therapeutic efficacy of deep brain stimulation (DBS) of the subthalamic nucleus (STN) for Parkinson's disease (PD) may be limited for some patients by the presence of stimulation-related side effects. Such effects are most often attributed to electrical current spread beyond the target region. Prior computational modeling studies have suggested that changing the degree of asymmetry of the individual phases of the biphasic, stimulus pulse may allow for more selective activation of neural elements in the target region. To the extent that different neural elements contribute to the therapeutic vs. side-effect inducing effects of DBS, such improved selectivity may provide a new parameter for optimizing DBS to increase the therapeutic window.

Methods

We investigated the effect of six different pulse geometries on cortical and myogenic evoked potentials in eight patients with PD whose leads were temporarily externalized following STN DBS implant surgery. DBS-cortical evoked potentials were quantified using peak to peak measurements and wavelets and myogenic potentials were quantified using RMS.

Results

We found that the slope of the recruitment curves differed significantly as a function of pulse geometry for both the cortical- and myogenic responses. Notably, this effect was observed most frequently when stimulation was delivered using a monopolar, as opposed to a bipolar, configuration.

Discussion

Manipulating pulse geometry results in differential physiological effects at both the cortical and neuromuscular level. Exploiting these differences may help to expand DBS' therapeutic window and support the potential for incorporating pulse geometry as an additional parameter for optimizing therapeutic benefit.

The Therapeutic Frequency Profile of Subthalamic Nucleus Deep Brain Stimulation in Rats Is Shaped by Antidromic Spike Failure

The therapeutic mechanisms of subthalamic nucleus (STN) deep brain stimulation (DBS) may depend on antidromic activation of cortex via the hyperdirect pathway. However, hyperdirect pathway neurons cannot reliably follow high-stimulation frequencies, and the spike failure rate appears to correlate with symptom relief as a function of stimulation frequency. We hypothesized that antidromic spike failure contributes to the cortical desynchronization caused by DBS. We measured in vivo evoked cortical activity in female Sprague Dawley rats and developed a computational model of cortical activation from STN DBS. We modeled stochastic antidromic spike failure to determine how spike failure affected the desynchronization of pathophysiological oscillatory activity in cortex. We found that high-frequency STN DBS desynchronized pathologic oscillations via the masking of intrinsic spiking through a combination of spike collision, refractoriness, and synaptic depletion. Antidromic spike failure shaped the parabolic relationship between DBS frequency and cortical desynchronization, with maximum desynchronization at ∼130 Hz. These findings reveal that antidromic spike failure plays a critical role in mediating the dependency of symptom relief on stimulation frequency.SIGNIFICANCE STATEMENT Deep brain stimulation (DBS) is a highly effective neuromodulation therapy, yet it remains uncertain why conventionally used stimulation frequencies (e.g., ∼130 Hz) are optimal. In this study, we demonstrate a potential explanation for the stimulation frequency dependency of DBS through a combination of in vivo experimental measurements and computational modeling. We show that high-frequency stimulation can desynchronize pathologic firing patterns in populations of neurons by inducing an informational lesion. However, sporadic spike failure at these high frequencies limits the efficacy of the informational lesion, yielding a parabolic profile with optimal effects at ∼130 Hz. This work provides a potential explanation for the therapeutic mechanism of DBS, and highlights the importance of considering spike failure in mechanistic models of DBS.

Electrophysiologic Evidence That Directional Deep Brain Stimulation Activates Distinct Neural Circuits in Patients With Parkinson Disease

OBJECTIVES

Deep brain stimulation (DBS) delivered via multicontact leads implanted in the basal ganglia is an established therapy to treat Parkinson disease (PD). However, the different neural circuits that can be modulated through stimulation on different DBS contacts are poorly understood. Evidence shows that electrically stimulating the subthalamic nucleus (STN) causes a therapeutic effect through antidromic activation of the hyperdirect pathway-a monosynaptic connection from the cortex to the STN. Recent studies suggest that stimulating the substantia nigra pars reticulata (SNr) may improve gait. The advent of directional DBS leads now provides a spatially precise means to probe these neural circuits and better understand how DBS affects distinct neural networks.

MATERIALS AND METHODS

We measured cortical evoked potentials (EPs) using electroencephalography (EEG) in response to low-frequency DBS using the different directional DBS contacts in eight patients with PD.

RESULTS

A short-latency EP at 3 milliseconds originating from the primary motor cortex appeared largest in amplitude when stimulating DBS contacts closest to the dorsolateral STN (p < 0.001). A long-latency EP at 10 milliseconds originating from the premotor cortex appeared strongest for DBS contacts closest to the SNr (p < 0.0001).

CONCLUSIONS

Our results show that at the individual patient level, electrical stimulation of different nuclei produces distinct EP signatures. Our approach could be used to identify the functional location of each DBS contact and thus help patient-specific DBS programming.

CLINICAL TRIAL REGISTRATION

The ClinicalTrials.gov registration number for the study is NCT04658641.

Deep Brain Stimulation in the Subthalamic Nucleus Can Improve Skilled Forelimb Movements and Retune Dynamics of Striatal Networks in a Rat Stroke Model

Recovery of upper limb (UL) impairment after stroke is limited in stroke survivors. Since stroke can be considered as a network disorder, neuromodulation may be an approach to improve UL motor dysfunction. Here, we evaluated the effect of high-frequency stimulation (HFS) of the subthalamic nucleus (STN) in rats on forelimb grasping using the single-pellet reaching (SPR) test after stroke and determined costimulated brain regions during STN-HFS using 2-[¹⁸F]Fluoro-2-deoxyglucose-([¹⁸F]FDG)-positron emission tomography (PET). After a 4-week training of SPR, photothrombotic stroke was induced in the sensorimotor cortex of the dominant hemisphere. Thereafter, an electrode was implanted in the STN ipsilateral to the infarction, followed by a continuous STN-HFS or sham stimulation for 7 days. On postinterventional day 2 and 7, an SPR test was performed during STN-HFS. Success rate of grasping was compared between these two time points. [¹⁸F]FDG-PET was conducted on day 2 and 3 after stroke, without and with STN-HFS, respectively. STN-HFS resulted in a significant improvement of SPR compared to sham stimulation. During STN-HFS, a significantly higher [¹⁸F]FDG-uptake was observed in the corticosubthalamic/pallidosubthalamic circuit, particularly ipsilateral to the stimulated side. Additionally, STN-HFS led to an increased glucose metabolism within the brainstem. These data demonstrate that STN-HFS supports rehabilitation of skilled forelimb movements, probably by retuning dysfunctional motor centers within the cerebral network.

The impact of pulse timing on cortical and subthalamic nucleus deep brain stimulation evoked potentials

The impact of pulse timing is an important factor in our understanding of how to effectively modulate the basal ganglia thalamocortical (BGTC) circuit. Single pulse low-frequency DBS-evoked potentials generated through electrical stimulation of the subthalamic nucleus (STN) provide insight into circuit activation, but how the long-latency components change as a function of pulse timing is not well-understood. We investigated how timing between stimulation pulses delivered in the STN region influence the neural activity in the STN and cortex. DBS leads implanted in the STN of five patients with Parkinson's disease were temporarily externalized, allowing for the delivery of paired pulses with inter-pulse intervals (IPIs) ranging from 0.2 to 10 ms. Neural activation was measured through local field potential (LFP) recordings from the DBS lead and scalp EEG. DBS-evoked potentials were computed using contacts positioned in dorsolateral STN as determined through co-registered post-operative imaging. We quantified the degree to which distinct IPIs influenced the amplitude of evoked responses across frequencies and time using the wavelet transform and power spectral density curves. The beta frequency content of the DBS evoked responses in the STN and scalp EEG increased as a function of pulse-interval timing. Pulse intervals <1.0 ms apart were associated with minimal to no change in the evoked response. IPIs from 1.5 to 3.0 ms yielded a significant increase in the evoked response, while those >4 ms produced modest, but non-significant growth. Beta frequency activity in the scalp EEG and STN LFP response was maximal when IPIs were between 1.5 and 4.0 ms. These results demonstrate that long-latency components of DBS-evoked responses are pre-dominantly in the beta frequency range and that pulse interval timing impacts the level of BGTC circuit activation.

Transcranial magnetic stimulation of the brain: What is stimulated? - A consensus and critical position paper

Transcranial (electro)magnetic stimulation (TMS) is currently the method of choice to non-invasively induce neural activity in the human brain. A single transcranial stimulus induces a time-varying electric field in the brain that may evoke action potentials in cortical neurons. The spatial relationship between the locally induced electric field and the stimulated neurons determines axonal depolarization. The induced electric field is influenced by the conductive properties of the tissue compartments and is strongest in the superficial parts of the targeted cortical gyri and underlying white matter. TMS likely targets axons of both excitatory and inhibitory neurons. The propensity of individual axons to fire an action potential in response to TMS depends on their geometry, myelination and spatial relation to the imposed electric field and the physiological state of the neuron. The latter is determined by its transsynaptic dendritic and somatic inputs, intrinsic membrane potential and firing rate. Modeling work suggests that the primary target of TMS is axonal terminals in the crown top and lip regions of cortical gyri. The induced electric field may additionally excite bends of myelinated axons in the juxtacortical white matter below the gyral crown. Neuronal excitation spreads ortho- and antidromically along the stimulated axons and causes secondary excitation of connected neuronal populations within local intracortical microcircuits in the target area. Axonal and transsynaptic spread of excitation also occurs along cortico-cortical and cortico-subcortical connections, impacting on neuronal activity in the targeted network. Both local and remote neural excitation depend critically on the functional state of the stimulated target area and network. TMS also causes substantial direct co-stimulation of the peripheral nervous system. Peripheral co-excitation propagates centrally in auditory and somatosensory networks, but also produces brain responses in other networks subserving multisensory integration, orienting or arousal. The complexity of the response to TMS warrants cautious interpretation of its physiological and behavioural consequences, and a deeper understanding of the mechanistic underpinnings of TMS will be critical for advancing it as a scientific and therapeutic tool.

The cortical evoked potential corresponds with deep brain stimulation efficacy in rats

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) antidromically activates the motor cortex (M1), and this cortical activation appears to play a role in the treatment of hypokinetic motor behaviors (Gradinaru V, Mogri M, Thompson KR, Henderson JM, Deisseroth K. Science 324: 354-359, 2009; Yu C, Cassar IR, Sambangi J, Grill WM. J Neurosci 40: 4323-4334, 2020). The synchronous antidromic activation takes the form of a short-latency cortical evoked potential (cEP) in electrocorticography (ECoG) recordings of M1. We assessed the utility of the cEP as a biomarker for STN DBS in unilateral 6-hydroxydopamine-lesioned female Sprague Dawley rats, with stimulating electrodes implanted in the STN and the ECoG recorded above M1. We quantified the correlations of the cEP magnitude and latency with changes in motor behavior from DBS and compared them to the correlation between motor behaviors and several commonly used spectral-based biomarkers. The cEP features correlated strongly with motor behaviors and were highly consistent across animals, whereas the spectral biomarkers correlated weakly with motor behaviors and were highly variable across animals. The cEP may thus be a useful biomarker for assessing the therapeutic efficacy of DBS parameters, as its features strongly correlate with motor behavior, it is consistent across time and subjects, it can be recorded under anesthesia, and it is simple to quantify with a large signal-to-noise ratio, enabling rapid, real-time evaluation. Additionally, our work provides further evidence that antidromic cortical activation mediates changes in motor behavior from STN DBS and that the dependence of DBS efficacy on stimulation frequency may be related to antidromic spike failure.NEW & NOTEWORTHY We characterize a new potential biomarker for deep brain stimulation (DBS), the cortical evoked potential (cEP), and demonstrate that it exhibits a robust correlation with motor behaviors as a function of stimulation frequency. The cEP may thus be a useful clinical biomarker for changes in motor behavior. This work also provides insight into the cortical mechanisms of DBS, suggesting that motor behaviors are strongly affected by the rate of antidromic spike failure during DBS.

Subthalamic deep brain stimulation of an anatomically detailed model of the human hyperdirect pathway

The motor hyperdirect pathway (HDP) is considered a key target in the treatment of Parkinson's disease with subthalamic deep brain stimulation (DBS). This hypothesis is partially derived from the association of HDP activation with evoked potentials (EPs) generated in the motor cortex and subthalamic nucleus (STN) after a DBS pulse. However, the biophysical details of how and when DBS-induced action potentials (APs) in HDP neurons reach their terminations in the cortex or STN remain unclear. Therefore, we used an anatomically detailed representation of the motor HDP, as well as the internal capsule (IC), in a model of human subthalamic DBS to explore AP activation and transmission in the HDP and IC. Our results show that small diameter HDP axons exhibited AP initiation in their subthalamic terminal arbor, which resulted in relatively long transmission latencies to cortex (∼3.5-8 ms). Alternatively, large diameter HDP axons were most likely to be directly activated in the capsular region, which resulted in short transmission times to the cortex (∼1-3 ms). However, those large diameter HDP antidromic APs would be indistinguishable from any other IC axons that were also activated by the stimulus. Conversely, DBS-induced APs in both small and large diameter HDP axons reached their synaptic boutons in the STN with similar timings, but both spanned a wide temporal range (∼0.5-5 ms). We also found that using anodic or bipolar stimulation helped to bias activation of the HDP over the IC. These computational results provide useful information for linking HDP activation with EP recordings in clinical experiments.NEW & NOTEWORTHY We used biophysical models to study pathway recruitment and conduction latencies of the hyperdirect pathway (HDP) in response to subthalamic deep brain stimulation (DBS). The model system allowed us to assess the influence of increased anatomical realism on pathway activity and the possibility of identifying HDP activity in evoked potentials (EPs) recorded in either the subthalamic nucleus (STN) or cortex. The model predicts that HDP activation is accentuated by complex axonal branching in the STN.

Hyperdirect connectivity of opercular speech network to the subthalamic nucleus

How the basal ganglia participate in the uniquely human behavior of speech is poorly understood, despite their known role in modulating critical aspects of cognitive and motor behavior. The subthalamic nucleus (STN) is well positioned to facilitate basal ganglia functions critical for speech. Using electrocorticography in patients undergoing awake deep brain stimulation (DBS) surgery, evidence is reported for a left opercular hyperdirect pathway in humans via stimulating the STN and examining antidromic-evoked activity in the left temporal, parietal, and frontal opercular cortex. These high-resolution cortical and subcortical mapping data provide evidence for hyperdirect connectivity between the inferior frontal gyrus and the STN. In addition, evoked potential data are consistent with the presence of monosynaptic projections from areas of the opercular ections may be unique to humans, evolving alongside the ability for speech.

Using electrical stimulation of the subthalamic nucleus and simultaneous cortical recordings in individuals undergoing deep brain stimulation, Jorge et al. provide electrophysiological evidence for a hyperdirect pathway to the basal ganglia from cortical areas that control sensory and motor-planning aspects of speech.

Electrophysiological characterization of the hyperdirect pathway and its functional relevance for subthalamic deep brain stimulation

The subthalamic nucleus (STN) receives input from various cortical areas via hyperdirect pathway (HDP) which bypasses the basal-ganglia loop. Recently, the HDP has gained increasing interest, because of its relevance for STN deep brain stimulation (DBS). To understand the HDP's role cortical responses evoked by STN-DBS have been investigated. These responses have short (<2 ms), medium (2-15 ms), and long (20-70 ms) latencies. Medium-latency responses are supposed to represent antidromic cortical activations via HDP. Together with long-latency responses the medium responses can potentially be used as biomarker of DBS efficacy as well as side effects. We here propose that the activation sequence of the cortical evoked responses can be conceptualized as high frequency oscillations (HFO) for signal analysis. HFO might therefore serve as marker for antidromic activation. Using existing knowledge on HFO recordings, this approach allows data analyses and physiological modeling to advance the pathophysiological understanding of cortical DBS-evoked high-frequency activity.

Deep Brain Stimulation Can Differentiate Subregions of the Human Subthalamic Nucleus Area by EEG Biomarkers

Introduction: Precise lead localization is crucial for an optimal clinical outcome of subthalamic nucleus (STN) deep brain stimulation (DBS) treatment in patients with Parkinson's disease (PD). Currently, anatomical measures, as well as invasive intraoperative electrophysiological recordings, are used to locate DBS electrodes. The objective of this study was to find an alternative electrophysiology tool for STN DBS lead localization.

Methods: Sixty-one postoperative electrophysiology recording sessions were obtained from 17 DBS-treated patients with PD. An intraoperative physiological method automatically detected STN borders and subregions. Postoperative EEG cortical activity was measured, while STN low frequency stimulation (LFS) was applied to different areas inside and outside the STN. Machine learning models were used to differentiate stimulation locations, based on EEG analysis of engineered features.

Results: A machine learning algorithm identified the top 25 evoked response potentials (ERPs), engineered features that can differentiate inside and outside STN stimulation locations as well as within STN stimulation locations. Evoked responses in the medial and ipsilateral fronto-central areas were found to be most significant for predicting the location of STN stimulation. Two-class linear support vector machine (SVM) predicted the inside (dorso-lateral region, DLR, and ventro-medial region, VMR) vs. outside [zona incerta, ZI, STN stimulation classification with an accuracy of 0.98 and 0.82 for ZI vs. VMR and ZI vs. DLR, respectively, and an accuracy of 0.77 for the within STN (DLR vs. VMR)]. Multiclass linear SVM predicted all areas with an accuracy of 0.82 for the outside and within STN stimulation locations (ZI vs. DLR vs. VMR).

Conclusions: Electroencephalogram biomarkers can use low-frequency STN stimulation to localize STN DBS electrodes to ZI, DLR, and VMR STN subregions. These models can be used for both intraoperative electrode localization and postoperative stimulation programming sessions, and have a potential to improve STN DBS clinical outcomes.

Three-Dimensional Heterogeneity of Cerebellar Interposed Nucleus-Recipient Zones in the Thalamic Nuclei

The cerebellum is conceptualized as a processor of complex movements and is also endowed with roles in cognitive and emotional behaviors. Although the axons of deep cerebellar nuclei are known to project to primary thalamic nuclei, macroscopic investigation of the characteristics of these projections, such as the spatial distribution of recipient zones, is lacking. Here, we studied the output of the cerebellar interposed nucleus (IpN) to the ventrolateral (VL) and centrolateral (CL) thalamic nuclei using electrophysiological recording in vivo and trans-synaptic viral tracing. We found that IpN stimulation induced mono-synaptic evoked potentials (EPs) in the VL but not the CL region. Furthermore, both the EPs induced by the IpN and the innervation of IpN projections displayed substantial heterogeneity across the VL region in three-dimensional space. These findings indicate that the recipient zones of IpN inputs vary between and within thalamic nuclei and may differentially control thalamo-cortical networks.

Histology-driven model of the macaque motor hyperdirect pathway

Emerging appreciation for the hyperdirect pathway (HDP) as an important cortical glutamatergic input to the subthalamic nucleus (STN) has motivated a wide range of recent investigations on its role in motor control, as well as the mechanisms of subthalamic deep brain stimulation (DBS). However, the pathway anatomy and terminal arbor morphometry by which the HDP links cortical and subthalamic activity are incompletely understood. One critical hindrance to advancing understanding is the lack of anatomically detailed population models which can help explain how HDP pathway anatomy and neuronal biophysics give rise to spatiotemporal patterns of stimulus-response activity observed in vivo. Therefore, the goal of this study was to establish a population model of motor HDP axons through application of generative algorithms constrained by recent histology and imaging data. The products of this effort include a de novo macaque brain atlas, detailed statistical analysis of histological reconstructions of macaque motor HDP axons, and the generation of 10,000 morphometrically constrained synthetic motor HDP axons. The synthetic HDP axons exhibited a 3.8% mean error with respect to parametric distributions of the fiber target volume, total length, number of bifurcations, bifurcation angles, meander angles, and segment lengths measured in BDA-labeled HDP axon reconstructions. As such, this large population of synthetic motor HDP axons represents an anatomically based foundation for biophysical simulations that can be coupled to electrophysiological and/or behavioral measurements, with the goal of better understanding the role of the HDP in motor system activity.

Axono-cortical evoked potentials as a new method of IONM for preserving the motor control network: a first study in three cases

BACKGROUND

White matter stimulation in an awake patient is currently the gold standard for identification of functional pathways. Despite the robustness and reproducibility of this method, very little is known about the electrophysiological mechanisms underlying the functional disruption. Axono-cortical evoked potentials (ACEPs) provide a reliable technique to explore these mechanisms.

OBJECTIVE

To describe the shape and spatial patterns of ACEPs recorded when stimulating the white matter of the caudal part of the right superior frontal gyrus while recording in the precentral gyrus.

METHODS

We report on three patients operated on under awake condition for a right superior frontal diffuse low-grade glioma. Functional sites were identified in the posterior wall of the cavity, whose 2-3-mA stimulation generated an arrest of movement. Once the resection was done, axono-cortical potentials were evoked: recording electrodes were put over the precentral gyrus, while stimulating at 1 Hz the white matter functional sites during 30-60 s. Unitary evoked potentials were averaged off-line. Waveform was visually analyzed, defining peaks and troughs, with quantitative measurements of their amplitudes and latencies. Spatial patterns of ACEPs were compared with patients' own and HCP-derived structural connectomics.

RESULTS

Axono-cortical evoked potentials (ACEPs) were obtained and exhibited complex shapes and spatial patterns that correlated only partially with structural connectivity patterns.

CONCLUSION

ACEPs is a new IONM methodology that could both contribute to elucidate the propagation of neuronal activity within a distributed network when stimulating white matter and provide a new technique for preserving motor control abilities during brain tumor resections.

Image-based biophysical modeling predicts cortical potentials evoked with subthalamic deep brain stimulation

BACKGROUND

Subthalamic deep brain stimulation (DBS) is an effective surgical treatment for Parkinson's disease and continues to advance technologically with an enormous parameter space. As such, in-silico DBS modeling systems have become common tools for research and development, but their underlying methods have yet to be standardized and validated.

OBJECTIVE

Evaluate the accuracy of patient-specific estimates of neural pathway activations in the subthalamic region against intracranial, cortical evoked potential (EP) recordings.

METHODS

Pathway activations were modeled in eleven patients using the latest advances in connectomic modeling of subthalamic DBS, focusing on the hyperdirect pathway (HDP) and corticospinal/bulbar tract (CSBT) for their relevance in human research studies. Correlations between pathway activations and respective EP amplitudes were quantified.

RESULTS

Good model performance required accurate lead localization and image fusions, as well as appropriate selection of fiber diameter in the biophysical model. While optimal model parameters varied across patients, good performance could be achieved using a global set of parameters that explained 60% and 73% of electrophysiologic activations of CSBT and HDP, respectively. Moreover, restricted models fit to only EP amplitudes of eight standard (monopolar and bipolar) electrode configurations were able to extrapolate variation in EP amplitudes across other directional electrode configurations and stimulation parameters, with no significant reduction in model performance across the cohort.

CONCLUSIONS

Our findings demonstrate that connectomic models of DBS with sufficient anatomical and electrical details can predict recruitment dynamics of white matter. These results will help to define connectomic modeling standards for preoperative surgical targeting and postoperative patient programming applications.

Single-pulse electrical stimulation methodology in freely moving rat

BACKGROUND

Cortico-cortical evoked potentials (CCEP) are becoming popular to infer brain connectivity and cortical excitability in implanted refractory epilepsy patients. Our goal was to transfer this methodology to the freely moving rodent.

NEW METHOD

CCEP were recorded on freely moving Sprague-Dawley rats, from cortical and subcortical areas using depth electrodes. Electrical stimulation was applied using 1 ms biphasic current pulse, cathodic first, at a frequency of 0.5 Hz, with intensities ranging from 0.2 to 0.8 mA. Data were then processed in a similar fashion to human clinical studies, which included epoch selection, artefact correction and smart averaging.

RESULTS

For a large range of tested intensities, we recorded CCEPs with very good signal to noise ratio and reproducibility between animals, without any behavioral modification. The CCEP were composed of different components according to recorded and stimulated sites, similarly to human recordings.

COMPARISON WITH EXISTING METHODS

We minimally adapted a clinically-motivated methodology to a freely moving rodent model to achieve high translational relevance of future preclinical studies.

CONCLUSIONS

Our results indicate that the CCEP methodology can be applied to freely moving rodents and transferred to preclinical research. This will be of interest to address various neuroscientific questions, in physiological and pathological conditions.

Real-time suppression and amplification of frequency-specific neural activity using stimulation evoked oscillations

BACKGROUND

Approaches to predictably control neural oscillations are needed to understand their causal role in brain function in healthy or diseased states and to advance the development of neuromodulation therapies.

OBJECTIVE

We present a closed-loop neural control and optimization framework to actively suppress or amplify low-frequency neural oscillations observed in local field potentials in real-time by using electrical stimulation. The rationale behind this control approach and our working hypothesis is that neural oscillatory activity evoked by electrical pulses can suppress or amplify spontaneous oscillations via destructive or constructive interference when the pulses are continuously delivered with appropriate amplitudes and at precise phases of the modulated oscillations in a closed-loop scheme.

METHODS

We tested our hypothesis in two nonhuman primates that exhibited a robust increase in low-frequency (8-30 Hz) oscillatory power in the subthalamic nucleus (STN) following administration of the neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP). To test our neural control approach, we targeted 8-17 Hz oscillations and used electrode arrays and electrical stimulation waveforms similar to those used in humans chronically implanted with brain stimulation systems. Stimulation parameters that maximize the suppression or amplification of neural oscillations were predicted using mathematical models of the stimulation evoked oscillations.

RESULTS

Our neural control and optimization approach was capable of actively and robustly suppressing or amplifying oscillations in the targeted frequency band (8-17 Hz) in real-time in the studied subjects.

CONCLUSIONS

The results from this study support our hypothesis and suggest that the proposed neural control framework allows one to characterize in controlled experiments the functional role of frequency-specific neural oscillations by using electrodes and stimulation waveforms currently being employed in humans.

A comprehensive model-based framework for optimal design of biomimetic patterns of electrical stimulation for prosthetic sensation

OBJECTIVE

Touch and proprioception are essential to motor function as shown by the movement deficits that result from the loss of these senses, e.g. due to neuropathy of sensory nerves. To achieve a high-performance brain-controlled prosthetic arm/hand thus requires the restoration of somatosensation, perhaps through intracortical microstimulation (ICMS) of somatosensory cortex (S1). The challenge is to generate patterns of neuronal activation that evoke interpretable percepts. We present a framework to design optimal spatiotemporal patterns of ICMS (STIM) that evoke naturalistic patterns of neuronal activity and demonstrate performance superior to four previous approaches.

APPROACH

We recorded multiunit activity from S1 during a center-out reach task (from proprioceptive neurons in Brodmann's area 2) and during application of skin indentations (from cutaneous neurons in Brodmann's area 1). We implemented a computational model of a cortical hypercolumn and used a genetic algorithm to design STIM that evoked patterns of model neuron activity that mimicked their experimentally-measured counterparts. Finally, from the ICMS patterns, the evoked neuronal activity, and the stimulus parameters that gave rise to it, we trained a recurrent neural network (RNN) to learn the mapping function between the physical stimulus and the biomimetic stimulation pattern, i.e. the sensory encoder to be integrated into a neuroprosthetic device.

MAIN RESULTS

We identified ICMS patterns that evoked simulated responses that closely approximated the measured responses for neurons within 50 µm of the electrode tip. The RNN-based sensory encoder generalized well to untrained limb movements or skin indentations. STIM designed using the model-based optimization approach outperformed STIM designed using existing linear and nonlinear mappings.

SIGNIFICANCE

The proposed framework produces an encoder that converts limb state or patterns of pressure exerted onto the prosthetic hand into STIM that evoke naturalistic patterns of neuronal activation.

Concurrent Deep Brain Stimulation Reduces the Direct Cortical Stimulation Necessary for Motor Output

BACKGROUND

Converging literatures suggest that deep brain stimulation (DBS) in Parkinson's disease affects multiple circuit mechanisms. One proposed mechanism is the normalization of primary motor cortex (M1) pathophysiology via effects on the hyperdirect pathway.

OBJECTIVES

We hypothesized that DBS would reduce the current intensity necessary to modulate motor-evoked potentials from focally applied direct cortical stimulation (DCS).

METHODS

Intraoperative subthalamic DBS, DCS, and preoperative diffusion tensor imaging data were acquired in 8 patients with Parkinson's disease.

RESULTS

In 7 of 8 patients, DBS significantly reduced the M1 DCS current intensity required to elicit motor-evoked potentials. This neuromodulation was specific to select DBS bipolar configurations. In addition, the volume of activated tissue models of these configurations were significantly associated with overlap of the hyperdirect pathway.

CONCLUSIONS

DBS reduces the current necessary to elicit a motor-evoked potential using DCS. This supports a circuit mechanism of DBS effectiveness, potentially involving the hyperdirect pathway that speculatively may underlie reductions in hypokinetic abnormalities in Parkinson's disease. © 2020 International Parkinson and Movement Disorder Society.

Proceedings of the Seventh Annual Deep Brain Stimulation Think Tank: Advances in Neurophysiology, Adaptive DBS, Virtual Reality, Neuroethics and Technology

The Seventh Annual Deep Brain Stimulation (DBS) Think Tank held on September 8th of 2019 addressed the most current: (1) use and utility of complex neurophysiological signals for development of adaptive neurostimulation to improve clinical outcomes; (2) Advancements in recent neuromodulation techniques to treat neuropsychiatric disorders; (3) New developments in optogenetics and DBS; (4) The use of augmented Virtual reality (VR) and neuromodulation; (5) commercially available technologies; and (6) ethical issues arising in and from research and use of DBS. These advances serve as both "markers of progress" and challenges and opportunities for ongoing address, engagement, and deliberation as we move to improve the functional capabilities and translational value of DBS. It is in this light that these proceedings are presented to inform the field and initiate ongoing discourse. As consistent with the intent, and spirit of this, and prior DBS Think Tanks, the overarching goal is to continue to develop multidisciplinary collaborations to rapidly advance the field and ultimately improve patient outcomes.

Latency of subthalamic nucleus deep brain stimulation-evoked cortical activity as a potential biomarker for postoperative motor side effects

OBJECTIVE

Here, we investigate whether cortical activation predicts motor side effects of deep brain stimulation (DBS) and whether these potential biomarkers have utility under general anesthesia.

METHODS

We recorded scalp potentials elicited by DBS during surgery (n = 11), both awake and under general anesthesia, and in an independent ambulatory cohort (n = 8). Across a range of stimulus configurations, we measured the amplitude and timing of short- and long-latency response components and linked them to motor side effects.

RESULTS

Regardless of anesthesia state, in both cohorts, DBS settings with capsular side effects elicited early responses with peak latencies clustering at <1 ms. This early response was preserved under anesthesia in all participants (11/11). In contrast, the long-latency components were suppressed completely in 6/11 participants. Finally, the latency of the earliest response could predict the presence of postoperative motor side effects both awake and under general anesthesia (84.8% and 75.8% accuracy, awake and under anesthesia, respectively).

CONCLUSION

DBS elicits short-latency cortical activation, both awake and under general anesthesia, which appears to reveal interactions between the stimulus and the corticospinal tract.

SIGNIFICANCE

Short-latency evoked cortical activity can potentially be used to aid both DBS lead placement and post-operative programming.

Biophysical reconstruction of the signal conduction underlying short-latency cortical evoked potentials generated by subthalamic deep brain stimulation

OBJECTIVE

Direct activation of the hyperdirect (HD) pathway has been linked to therapeutic benefit from subthalamic deep brain stimulation (DBS) for the treatment of Parkinson's disease (PD). We sought to quantify the axonal conduction biophysics of corticofugal axons directly stimulated by subthalamic DBS and reconcile those findings with short-latency cortical evoked potential (EP) results.

METHODS

We used a detailed computational model of human subthalamic DBS to quantify axonal activation and conduction. Signal propagation to cortex was evaluated for medium (5.7 µm), large (10.0 µm), and exceptionally large (15.0 µm) diameter corticofugal axons associated with either internal capsule (IC) fibers of passage or the HD pathway. We then compared the modeling results to human cortical EP measurements that have described an exceptionally fast component (EP0) occurring ~1 ms after the stimulus pulse, a fast component (EP1) at ~3 ms, and a slower component (EP2) at ~5 ms.

RESULTS

Subthalamic stimulation of the HD pathway with large and medium diameter axons propagated action potentials to cortex with timings that coincide with the EP1 and EP2 signals, respectively. Only direct activation of exceptionally large diameter fibers in the IC generated signals that could approach the EP0 timing. However, the action potential biophysics do not generally support the existence of a cortical EP less than 1.5 ms after DBS onset.

CONCLUSIONS

The EP1 and EP2 signals can be biophysically linked to antidromic activation of the HD pathway.

SIGNIFICANCE

Theoretical reconstruction of cortical EPs from subthalamic DBS demonstrate a convergence of anatomical, biophysical, and electrophysiological results.

Phase-Dependent Suppression of Beta Oscillations in Parkinson's Disease Patients

Synchronized oscillations within and between brain areas facilitate normal processing, but are often amplified in disease. A prominent example is the abnormally sustained beta-frequency (∼20 Hz) oscillations recorded from the cortex and subthalamic nucleus of Parkinson's disease patients. Computational modeling suggests that the amplitude of such oscillations could be modulated by applying stimulation at a specific phase. Such a strategy would allow selective targeting of the oscillation, with relatively little effect on other activity parameters. Here, activity was recorded from 10 awake, parkinsonian patients (6 male, 4 female human subjects) undergoing functional neurosurgery. We demonstrate that stimulation arriving on a particular patient-specific phase of the beta oscillation over consecutive cycles could suppress the amplitude of this pathophysiological activity by up to 40%, while amplification effects were relatively weak. Suppressive effects were accompanied by a reduction in the rhythmic output of subthalamic nucleus (STN) neurons and synchronization with the mesial cortex. While stimulation could alter the spiking pattern of STN neurons, there was no net effect on firing rate, suggesting that reduced beta synchrony was a result of alterations to the relative timing of spiking activity, rather than an overall change in excitability. Together, these results identify a novel intrinsic property of cortico-basal ganglia synchrony that suggests the phase of ongoing neural oscillations could be a viable and effective control signal for the treatment of Parkinson's disease. This work has potential implications for other brain diseases with exaggerated neuronal synchronization and for probing the function of rhythmic activity in the healthy brain.SIGNIFICANCE STATEMENT In Parkinson's disease (PD), movement impairment is correlated with exaggerated beta frequency oscillations in the cerebral cortex and subthalamic nucleus (STN). Using a novel method of stimulation in PD patients undergoing neurosurgery, we demonstrate that STN beta oscillations can be suppressed when consecutive electrical pulses arrive at a specific phase of the oscillation. This effect is likely because of interrupting the timing of neuronal activity rather than excitability, as stimulation altered the firing pattern of STN spiking without changing overall rate. These findings show the potential of oscillation phase as an input for "closed-loop" stimulation, which could provide a valuable neuromodulation strategy for the treatment of brain disorders and for elucidating the role of neuronal oscillations in the healthy brain.

Cortical Potentials Evoked by Subthalamic Stimulation Demonstrate a Short Latency Hyperdirect Pathway in Humans

A monosynaptic projection from the cortex to the subthalamic nucleus is thought to have an important role in basal ganglia function and in the mechanism of therapeutic subthalamic deep-brain stimulation, but in humans the evidence for its existence is limited. We sought physiological confirmation of the cortico-subthalamic hyperdirect pathway using invasive recording techniques in patients with Parkinson's disease (9 men, 1 woman). We measured sensorimotor cortical evoked potentials using a temporary subdural strip electrode in response to low-frequency deep-brain stimulation in patients undergoing awake subthalamic or pallidal lead implantations. Evoked potentials were grouped into very short latency (<2 ms), short latency (2-10 ms), and long latency (10-100 ms) from the onset of the stimulus pulse. Subthalamic and pallidal stimulation resulted in very short-latency evoked potentials at 1.5 ms in the primary motor cortex accompanied by EMG-evoked potentials consistent with corticospinal tract activation. Subthalamic, but not pallidal stimulation, resulted in three short-latency evoked potentials at 2.8, 5.8, and 7.7 ms in a widespread cortical distribution, consistent with antidromic activation of the hyperdirect pathway. Long-latency potentials were evoked by both targets, with subthalamic responses lagging pallidal responses by 10-20 ms, consistent with orthodromic activation of the thalamocortical pathway. The amplitude of the first short-latency evoked potential was predictive of the chronic therapeutic stimulation contact.SIGNIFICANCE STATEMENT This is the first physiological demonstration of the corticosubthalamic hyperdirect pathway and its topography at high spatial resolution in humans. We studied cortical potentials evoked by deep-brain stimulation in patients with Parkinson's disease undergoing awake lead implantation surgery. Subthalamic stimulation resulted in multiple short-latency responses consistent with activation of hyperdirect pathway, whereas no such response was present during pallidal stimulation. We contrast these findings with very short latency, direct corticospinal tract activations, and long-latency responses evoked through polysynaptic orthodromic projections. These findings underscore the importance of incorporating the hyperdirect pathway into models of human basal ganglia function.

Patterned low-frequency deep brain stimulation induces motor deficits and modulates cortex-basal ganglia neural activity in healthy rats

Deep brain stimulation (DBS) is an effective therapy for movement disorders, including Parkinson's disease (PD), although the mechanisms of action remain unclear. Abnormal oscillatory neural activity is correlated with motor symptoms, and pharmacological or DBS treatment that alleviates motor symptoms appears to suppress abnormal oscillations. However, whether such oscillatory activity is causal of motor deficits such as tremor remains unclear. Our goal was to generate abnormal oscillatory activity in the cortex-basal ganglia loop using patterned subthalamic nucleus DBS and to quantify motor behavior in awake healthy rats. Stimulation patterns were designed via model-based optimization to increase power in the low-frequency (7-11 Hz) band because these oscillations are associated with the emergence of motor symptoms in the 6-hydroxydopamine lesioned rat model of parkinsonism. We measured motor activity using a head-mounted accelerometer, as well as quantified neural activity in cortex and globus pallidus (GP), in response to 5 stimulation patterns that generated a range of 7- to 11-Hz spectral power. Stimulation patterns induced oscillatory activity in the low-frequency band in the cortex and GP and caused tremor, whereas control patterns and regular 50-Hz DBS did not generate any such effects. Neural and motor-evoked responses observed during stimulation were synchronous and time-locked to stimulation bursts within the patterns. These results identified elements of irregular patterns of stimulation that were correlated with tremor and tremor-related neural activity in the cortex and basal ganglia and may lead to the identification of the oscillatory activity and structures associated with the generation of tremor activity. NEW & NOTEWORTHY Subthalamic nucleus deep brain stimulation is a promising therapy for movement disorders such as Parkinson's disease. Several groups reported correlation between suppression of abnormal oscillatory activity in the cortex-basal ganglia and motor symptoms, but it remains unclear whether such oscillations play a causal role in the emergence of motor symptoms. We demonstrate generation of tremor and pathological oscillatory activity in otherwise healthy rats by stimulation with patterns that produced increases in low-frequency oscillatory activity.