Electrical Stimulation publications

Before selecting a saccadic target, we often acquire partial information about the location of the forthcoming target and preselect a region of visual space even before the target becomes visible. To determine whether the supplementary eye field (SEF) represents information signifying the potential target space, we examined neuronal activity in the SEF of monkeys performing a behavioral task designed to isolate the process of visuospatial preselection under uncertainty from the process of selecting a specified location. Our data showed that the activity of SEF neurons represented information about the potential target space instructed by symbolic cues. Increased activity of visuospatially selective SEF neurons encoded the potential target space, which could be a mechanism facilitating subsequent selection of an appropriate target. Furthermore, electrical stimulation of the SEF during the preselection period disrupted subsequent target selection. These results demonstrate that the SEF contributes to the preselection of potential target spaces based on partial information.

Introduction: Evoked compound action potentials (ECAPs) measured using epidural spinal recordings (ESRs) during epidural spinal cord stimulation (SCS) can help elucidate fundamental mechanisms for the treatment of pain, as well as inform closed-loop control of SCS. Previous studies have used ECAPs to characterize the neural response to various neuromodulation therapies and have demonstrated that ECAPs are highly prone to multiple sources of artifact, including post-stimulus pulse capacitive artifact, electromyography (EMG) bleed-through, and motion artifact resulting from disturbance of the electrode/tissue interface during normal behavior. However, a thorough characterization has yet to be performed for how these sources of artifact may contaminate recordings within the temporal window commonly used to determine activation of A-beta fibers in a large animal model. Methods: We characterized the sources of artifacts that can contaminate the recording of ECAPs in an epidural SCS swine model using the Abbott Octrode™ lead. Muscle paralytics were administered to block muscle activation preventing EMG from contaminating the recorded ECAPs. Concurrent EMG recordings of the longissimus, a long muscle of the back, were used to confirm a 2-4 millisecond (ms) latency source of EMG bleed-through that frequently contaminated the A-beta temporal window. Additionally, we obtained recordings approximately 5-10 minutes post-mortem after clear evoked A-beta and associated EMG responses ceased to characterize the representation of stimulation artifact across the array. Results: Spinal ECAP recordings can be contaminated by capacitive artifact, short latency EMG from nearby long muscles of the back, and motion artifact from multiple sources. In many cases, the capacitive artifact can appear nearly identical in duration and waveshape to evoked A-beta responses. These sources of EMG can have phase shifts across the electrode array, very similar to the phase shift anticipated by propagation of an evoked A-beta fiber response across the array. This short latency EMG is often evident at currents similar to those needed to activate A-beta fibers associated with the treatment of pain. Changes in cerebrospinal fluid between the cord and dura, and motion induced during breathing created a cyclic oscillation in all evoked components of the recorded ECAP signal. Conclusion: Careful controls must be implemented to accurately separate neural signal from the sources of artifact in spinal cord ECAPs. To address this, we suggest experimental procedures and associated reporting requirements necessary to disambiguate the underlying neural response from these confounds. These data are important to better understand the conceptual framework for recorded ESRs, with components such as ECAPs, EMG responses and artifacts, and have important implications for closed-loop control algorithms to account for transient motion such as postural changes and cough.

Objective This study investigated the relationship between DBS evoked potentials (EPs) and diffusion tensor imaging (DTI) in a group of patients with dystonia who underwent DBS treatment. EPs and DTI are both useful methods for studying neural connectivity in the brain but measure different aspects of brain function. EPs provide information on electrical connectivity, while DTI provides information on anatomical pathways connecting regions. Methods This study focused on the pallidum and motor thalamus nuclei, which are common targets for DBS in dystonia. Prior to DBS implantation, DTI images were acquired for each patient, and were processed to obtain DTI coefficients such as length (L), volume (V), and fractional anisotropy (FA) of the fiber tracts. The relationship between the fiber tracts and electrophysiology was examined using a generalized linear model (GLM). Results We showed that the amplitude of EPs correlated with FA and tract volume, while delay correlated with tract length. These findings suggest that DBS signals travel across tracts to affect both local and distant brain regions, and the magnitude of the effect of DBS is determined by the integrity of the white matter tract, while DBS signal delay is affected by the tract length. Our results further suggest that the magnitude and delay of the spread of the DBS signal may be predicted by the DTI connectivity. This provides strong supporting evidence for other studies that have assumed, but have been unable to test, such a relationship. Conclusion Overall, this study suggests that the electrical effects of DBS can be at least partially predicted by noninvasive DTI imaging in patients with dystonia. By combining EPs with DTI, we could investigate the propagation of stimulation pulses through brain regions. While this relationship has been previously hypothesized by the neuroscience community, this is the first study in humans to demonstrate this relationship between DBS EPs and DTI, thereby advancing the field of human brain mapping and enhancing the precision of neurosurgical targeting.

Hand movements are associated with modulations of neuronal activity across several interconnected cortical areas, including the primary motor cortex (M1) and the dorsal and ventral premotor cortices (PMd and PMv). Local field potentials (LFPs) provide a link between neuronal discharges and synaptic inputs. Our current understanding of how LFPs vary in M1, PMd, and PMv during contralateral and ipsilateral movements is incomplete. To help reveal unique features in the pattern of modulations, we simultaneously recorded LFPs in these areas in two macaque monkeys performing reach and grasp movements with either the right or left hand. The greatest effector-dependent differences were seen in M1, at low (≤13 Hz) and γ frequencies. In premotor areas, differences related to hand use were only present in low frequencies. PMv exhibited the greatest increase in low frequencies during instruction cues and the smallest effector-dependent modulation during movement execution. In PMd, δ oscillations were greater during contralateral reach and grasp, and β activity increased during contralateral grasp. In contrast, β oscillations decreased in M1 and PMv. These results suggest that while M1 primarily exhibits effector-specific LFP activity, premotor areas compute more effector-independent aspects of the task requirements, particularly during movement preparation for PMv and production for PMd. The generation of precise hand movements likely relies on the combination of complementary information contained in the unique pattern of neural modulations contained in each cortical area. Accordingly, integrating LFPs from premotor areas and M1 could enhance the performance and robustness of brain-machine interfaces.

Anxiety is a common symptom across psychiatric disorders, but the neurophysiological underpinnings of these symptoms remain unclear. This knowledge gap has prevented the development of circuit-based treatments that can target the neural substrates underlying anxiety. Here, we conducted an electrophysiological mapping study to identify neurophysiological activity associated with self-reported state anxiety in 17 subjects implanted with intracranial electrodes for seizure localization. Participants had baseline anxiety traits ranging from minimal to severe. Subjects volunteered to participate in an anxiety induction task in which they were temporarily exposed to the threat of unpredictable shock during intracranial recordings. We found that anterior insular beta oscillatory activity was selectively elevated during epochs when unpredictable aversive stimuli were being delivered, and this enhancement in insular beta was correlated with increases in self-reported anxiety. Beta oscillatory activity within the frontoinsular region was also evoked selectively by cues-predictive of threat, but not safety cues. Anterior insular gamma responses were less selective than gamma, strongly evoked by aversive stimuli and had weaker responses to salient threat and safety cues. On longer timescales, this gamma signal also correlated with increased skin conductance, a measure of autonomic state. Lastly, we found that direct electrical stimulation of the anterior insular cortex in a subset of subjects elicited self-reported increases in anxiety that were accompanied by enhanced frontoinsular beta oscillations. Together, these findings suggest that electrophysiologic representations of anxiety- related states and behaviors exist within anterior insular cortex. The findings also suggest the potential of reducing anterior insular beta activity as a therapeutic target for refractory anxiety-spectrum disorders.

Effective neural stimulation requires adequate parametrization. Gaussian-process (GP)-based Bayesian optimization (BO) offers a framework to discover optimal stimulation parameters in real time. Here, we first provide a general protocol to deploy this framework in neurostimulation interventions and follow by exemplifying its use in detail. Specifically, we describe the steps to implant rats with multi-channel electrode arrays in the hindlimb motor cortex. We then detail how to utilize the GP-BO algorithm to maximize evoked target movements, measured as electromyographic responses. For complete details on the use and execution of this protocol, please refer to Bonizzato and colleagues (2023).1.

Introduction Cochlear implants (CIs) are the treatment of choice for severe to profound hearing loss. Variability in CI outcomes remains despite advances in technology and is attributed in part to differences in cortical processing. Studying these differences in CI users is technically challenging. Spectrally degraded stimuli presented to normal-hearing individuals approximate input to the central auditory system in CI users. This study used intracranial electroencephalography (iEEG) to investigate cortical processing of spectrally degraded speech. Methods Participants were adult neurosurgical epilepsy patients. Stimuli were utterances /aba/ and /ada/, spectrally degraded using a noise vocoder (1–4 bands) or presented without vocoding. The stimuli were presented in a two-alternative forced choice task. Cortical activity was recorded using depth and subdural iEEG electrodes. Electrode coverage included auditory core in posteromedial Heschl’s gyrus (HGPM), superior temporal gyrus (STG), ventral and dorsal auditory-related areas, and prefrontal and sensorimotor cortex. Analysis focused on high gamma (70–150 Hz) power augmentation and alpha (8–14 Hz) suppression. Results Chance task performance occurred with 1–2 spectral bands and was near-ceiling for clear stimuli. Performance was variable with 3–4 bands, permitting identification of good and poor performers. There was no relationship between task performance and participants demographic, audiometric, neuropsychological, or clinical profiles. Several response patterns were identified based on magnitude and differences between stimulus conditions. HGPM responded strongly to all stimuli. A preference for clear speech emerged within non-core auditory cortex. Good performers typically had strong responses to all stimuli along the dorsal stream, including posterior STG, supramarginal, and precentral gyrus; a minority of sites in STG and supramarginal gyrus had a preference for vocoded stimuli. In poor performers, responses were typically restricted to clear speech. Alpha suppression was more pronounced in good performers. In contrast, poor performers exhibited a greater involvement of posterior middle temporal gyrus when listening to clear speech. Discussion Responses to noise-vocoded speech provide insights into potential factors underlying CI outcome variability. The results emphasize differences in the balance of neural processing along the dorsal and ventral stream between good and poor performers, identify specific cortical regions that may have diagnostic and prognostic utility, and suggest potential targets for neuromodulation-based CI rehabilitation strategies.

Electronic retinal prostheses for stimulating retinal neurons are promising for vision restoration. However, the rigid electrodes of conventional retinal implants can inflict damage on the soft retina tissue. They also have limited selectivity due to their poor proximity to target cells in the degenerative retina. Here we present a soft artificial retina (thickness, 10 μm) where flexible ultrathin photosensitive transistors are integrated with three-dimensional stimulation electrodes of eutectic gallium–indium alloy. Platinum nanoclusters locally coated only on the tip of these three-dimensional liquid-metal electrodes show advantages in reducing the impedance of the stimulation electrodes. These microelectrodes can enhance the proximity to the target retinal ganglion cells and provide effective charge injections (72.84 mC cm −2 ) to elicit neural responses in the retina. Their low Young’s modulus (234 kPa), owing to their liquid form, can minimize damage to the retina. Furthermore, we used an unsupervised machine learning approach to effectively identify the evoked spikes to grade neural activities within the retinal ganglion cells. Results from in vivo experiments on a retinal degeneration mouse model reveal that the spatiotemporal distribution of neural responses on their retina can be mapped under selective localized illumination areas of light, suggesting the restoration of their vision.

In neurophysiology, achieving precise correlation between physiological responses and anatomic structures is a significant challenge. Therefore, the accuracy of the electrode marking method is crucial. In this study, we describe a tungsten-deposition method, in which tungsten oxide is generated by applying biphasic current pulses to conventional tungsten electrodes. The electrical current used was 40–50 μA, which is similar to that used in electrical microstimulation experiments. The size of the markings ranged from 10 to 100 μm, corresponding to the size of the electrode tip, which is smaller than that of existing marking methods. Despite the small size of the markings, detection is easy as the marking appears in bright red under dark-field observation after Nissl staining. This marking technique resulted in low tissue damage and was maintained in vivo for at least two years. The feasibility of this method was tested in mouse and macaque brains.

Introduction Benzodiazepines (BDZs) are commonly used to treat the symptoms of movement disorders; however, deep brain stimulation (DBS) has become a popular treatment for these disorders. Previous studies have investigated the effects of BDZ on cortical activity, no data are currently available on their effects on deep brain regions, nor on these regions’ responses to DBS. How the BDZ affects the thalamus and basal ganglia in dystonia patients remains unknown. Methods DBS recordings were performed in ventral oralis anterior/posterior (VoaVop), ventral intermediate (VIM) and ventral anterior (VA) thalamic subnuclei, as well as globus pallidus interna (GPi) and subthalamic nucleus (STN). Evoked potentials (EP) and frequency domain analysis were performed to determine the BDZ effect on neural activities compared to the control condition (off-BDZ). Three male pediatric patients with dystonia treated with BDZ and undergoing depth electrode evaluation for clinical targeting were recruited for the study. Stimulation was administered at 25 and 55 Hz frequencies and recordings were simultaneously gathered through pairs of externalized stereoelectroencephalography (sEEG) electrodes. EP amplitude and the effect of stimulation on the frequency spectrum of activity were compared at baseline and following clinical administration of BDZ. Results Frequency analysis showed consistent reductions in activity during BDZ treatment in all studied brain regions for all patients. Evoked potential (EP) analysis showed increased subthalamic nucleus (STN) EP amplitude and decreased ventral intermediate (VIM) and STN EP amplitude during BDZ treatment. Interpretation BDZs reduce thalamic and basal ganglia activity in multiple regions and alter the efficacy of transmission between these regions. While the mechanism is unknown our results confirm the known widespread effects of this class of medications and identify specific areas within the motor system that are directly affected.

The rate model of basal ganglia function predicts that muscle activity in dystonia is due to disinhibition of thalamus resulting from decreased inhibitory input from pallidum. We seek to test this hypothesis in children with dyskinetic cerebral palsy undergoing evaluation for deep brain stimulation (DBS) to analyze movement-related activity in different brain regions. The results revealed prominent beta-band frequency peaks in the globus pallidus interna (GPi), ventral oralis anterior/posterior (VoaVop) subnuclei of the thalamus, and subthalamic nucleus (STN) during movement but not at rest. Connectivity analysis indicated stronger coupling between STN-VoaVop and STN-GPi compared to GPi-STN. These findings contradict the hypothesis of decreased thalamic inhibition in dystonia, suggesting that abnormal patterns of inhibition and disinhibition, rather than reduced GPi activity, contribute to the disorder. Additionally, the study implies that correcting abnormalities in GPi function may explain the effectiveness of DBS targeting the STN and GPi in treating dystonia.

We report on the power spectral analysis of multiunit activity (MUA) of lumbar spinal interneurons during a locomotor task. In line with prior studies, we found evidence of greater power at the frequency of locomotion in high lumbar segments and during the flexion phase. Our results also confirm prior observations from our laboratory that the rhythmically active MUA behaves as a longitudinal standing wave of neural activation that is flexor dominant.

Application of closed-loop approaches in systems neuroscience and brain-computer interfaces holds great promise for revolutionizing our understanding of the brain and for developing novel neuromodulation strategies to restore lost function. The anterior forebrain mesocircuit (AFM) of the mammalian brain is hypothesized to underlie arousal regulation of the cortex and striatum, and support cognitive functions during wakefulness. Dysfunction of arousal regulation is hypothesized to contribute to cognitive dysfunctions in various neurological disorders, and most prominently in patients following traumatic brain injury (TBI). Several clinical studies have explored the use of daily central thalamic deep brain stimulation (CT-DBS) within the AFM to restore consciousness and executive attention in TBI patients. In this study, we explored the use of closed-loop CT-DBS in order to episodically regulate arousal of the AFM of a healthy non-human primate (NHP) with the goal of restoring behavioral performance. We used pupillometry and near real-time analysis of ECoG signals to episodically initiate closed-loop CT-DBS and here we report on our ability to enhance arousal and restore the animal's performance. The initial computer based approach was then experimentally validated using a customized clinical-grade DBS device, the DyNeuMo-X, a bi-directional research platform used for rapidly testing closed-loop DBS. The successful implementation of the DyNeuMo-X in a healthy NHP supports ongoing clinical trials employing the internal DyNeuMo system (NCT05437393, NCT05197816) and our goal of developing and accelerating the deployment of novel neuromodulation approaches to treat cognitive dysfunction in patients with structural brain injuries and other etiologies.

Intracortical microstimulation (ICMS) is commonly used in many experimental and clinical paradigms; however, its effects on the activation of neurons are still not completely understood. To document the responses of cortical neurons in awake nonhuman primates to stimulation, we recorded single-unit activity while delivering single-pulse stimulation via Utah arrays implanted in primary motor cortex (M1) of three macaque monkeys. Stimuli between 5 and 50 μA delivered to single channels reliably evoked spikes in neurons recorded throughout the array with delays of up to 12 ms. ICMS pulses also induced a period of inhibition lasting up to 150 ms that typically followed the initial excitatory response. Higher current amplitudes led to a greater probability of evoking a spike and extended the duration of inhibition. The likelihood of evoking a spike in a neuron was dependent on the spontaneous firing rate as well as the delay between its most recent spike time and stimulus onset. Tonic repetitive stimulation between 2 and 20 Hz often modulated both the probability of evoking spikes and the duration of inhibition; high-frequency stimulation was more likely to change both responses. On a trial-by-trial basis, whether a stimulus evoked a spike did not affect the subsequent inhibitory response; however, their changes over time were often positively or negatively correlated. Our results document the complex dynamics of cortical neural responses to electrical stimulation that need to be considered when using ICMS for scientific and clinical applications.

In humans, the level of consciousness is assessed by quantifying the spatiotemporal complexity of cortical responses using Perturbational Complexity Index (PCI) and related PCIst (st, state transitions). Here we validate PCIst in freely moving rats and mice by showing that it is lower in NREM sleep and slow wave anesthesia than in wake or REM sleep, as in humans. We then show that (1) low PCIst is associated with the occurrence of an OFF period of neuronal silence; (2) stimulation of deep, but not superficial, cortical layers leads to reliable PCIst changes across sleep/wake and anesthesia; (3) consistent PCIst changes are independent of which single area is being stimulated or recorded, except for recordings in mouse prefrontal cortex. These experiments show that PCIst can reliably measure vigilance states in unresponsive animals and support the hypothesis that it is low when an OFF period disrupts causal interactions in cortical networks.

The contribution of different brain regions to movement abnormalities in children with dystonia is unknown. Three awake subjects undergoing depth electrode implantation for assessments of potential deep brain recording targets performed a rhythmic figure-8 drawing task. Two subjects had dystonia, one was undergoing testing for treatment of Tourette Syndrome and had neither dystonia nor abnormal movements during testing. Movement-related signals were evaluated by determining the magnitude of task-related frequency components. Brain signals were recorded in globus pallidus internus (GPi), the ventral oralis anterior/posterior (VoaVop) and the ventral intermediate (Vim) nuclei of the thalamus. In comparison to the subject without dystonia, both children with dystonia showed increased task-related activity in GPi and Vim. This finding is consistent with a role of both basal ganglia and cerebellar outputs in the pathogenesis of dystonia. Our results further suggest that frequency analysis of brain recordings during cyclic movements may be a useful tool for analysis of the presence of movement-related signals in various brain regions.

The primate amygdala is a complex consisting of over a dozen nuclei that have been implicated in a host of cognitive functions, individual differences, and psychiatric illnesses. These functions are implemented through distinct connectivity profiles, which have been documented in animals but remain largely unknown in humans. Here we present results from 25 neurosurgical patients who had concurrent electrical stimulation of the amygdala with intracranial electroencephalography (electrical stimulation tract-tracing; es-TT), or fMRI (electrical stimulation fMRI; es-fMRI), methods providing strong inferences about effective connectivity of amygdala subdivisions with the rest of the brain. We quantified functional connectivity with medial and lateral amygdala, the temporal order of these connections on the timescale of milliseconds, and also detail second-order effective connectivity among the key nodes. These findings provide a uniquely detailed characterization of human amygdala functional connectivity that will inform functional neuroimaging studies in healthy and clinical populations.

Animals can capitalize on invariance in the environment by learning and automating highly consistent actions; however, they must also remain flexible and adapt to environmental changes. It remains unclear how primary motor cortex (M1) can drive precise movements, yet also support behavioral exploration when faced with consistent errors. Using a reach-to-grasp task in rats, along with simultaneous electrophysiological monitoring in M1 and dorsolateral striatum (DLS), we find that behavioral exploration to overcome consistent task errors is closely associated with tandem increases in M1 and DLS neural variability; subsequently, consistent ensemble patterning returns with convergence to a new successful strategy. We also show that compared to reliably patterned intracranial microstimulation in M1, variable stimulation patterns result in significantly greater movement variability. Our results thus indicate that motor and striatal areas can flexibly transition between two modes, reliable neural pattern generation for automatic and precise movements versus variable neural patterning for behavioral exploration.

Spinal interneurons play a critical role in motor output. A given interneuron may receive convergent input from several different sensory modalities and descending centers and relay this information to just as many targets. Therefore, there is a critical need to quantify populations of spinal interneurons simultaneously. Here, we quantify the functional connectivity of spinal neurons through the concurrent recording of populations of lumbar interneurons and hindlimb motor units in the in vivo cat model during activation of either the ipsilateral sural nerve or contralateral tibial nerve. Two microelectrode arrays were placed into lamina VII, one at L3 and a second at L6/7, while an electrode array was placed on the surface of the exposed muscle. Stimulation of tibial and sural nerves elicited similar changes in the discharge rate of both interneurons and motor units. However, these same neurons showed highly significant differences in prevalence and magnitude of correlated activity underlying these two forms of afferent drive. Activation of the ipsilateral sural nerve resulted in highly correlated activity, particularly at the caudal array. In contrast, the contralateral tibial nerve resulted in less, but more widespread correlated activity at both arrays. These data suggest that the ipsilateral sural nerve has dense projections onto caudal lumbar spinal neurons, while contralateral tibial nerve has a sparse pattern of projections.

We report on the ensemble organization of interneuronal activity in the spinal cord during locomotor movements and show that lumbar intermediate zone interneurons organize in two groups related to the two major phases of walking: stance and swing. Ensemble organization is also shown to better correlate with muscular output than single-cell activity, although ensemble membership does not appear to be somatotopically organized within the spinal cord.