skip to main content


Title: Non-invasive suppression of essential tremor via phase-locked disruption of its temporal coherence
Abstract Aberrant neural oscillations hallmark numerous brain disorders. Here, we first report a method to track the phase of neural oscillations in real-time via endpoint-corrected Hilbert transform (ecHT) that mitigates the characteristic Gibbs distortion. We then used ecHT to show that the aberrant neural oscillation that hallmarks essential tremor (ET) syndrome, the most common adult movement disorder, can be transiently suppressed via transcranial electrical stimulation of the cerebellum phase-locked to the tremor. The tremor suppression is sustained shortly after the end of the stimulation and can be phenomenologically predicted. Finally, we use feature-based statistical-learning and neurophysiological-modelling to show that the suppression of ET is mechanistically attributed to a disruption of the temporal coherence of the aberrant oscillations in the olivocerebellar loop, thus establishing its causal role. The suppression of aberrant neural oscillation via phase-locked driven disruption of temporal coherence may in the future represent a powerful neuromodulatory strategy to treat brain disorders.  more » « less
Award ID(s):
1845348
NSF-PAR ID:
10215413
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Nature Communications
Volume:
12
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Deep brain stimulation (DBS) is a commonly used treatment for medication resistant Parkinson’s disease and is an emerging treatment for other neurological disorders. More recently, phase-specific adaptive DBS (aDBS), whereby the application of stimulation is locked to a particular phase of tremor, has been proposed as a strategy to improve therapeutic efficacy and decrease side effects. In this work, in the context of these phase-specific aDBS strategies, we investigate the dynamical behavior of large populations of coupled neurons in response to near-periodic stimulation, namely, stimulation that is periodic except for a slowly changing amplitude and phase offset that can be used to coordinate the timing of applied input with a specified phase of model oscillations. Using an adaptive phase-amplitude reduction strategy, we illustrate that for a large population of oscillatory neurons, the temporal evolution of the associated phase distribution in response to near-periodic forcing can be captured using a reduced order model with four state variables. Subsequently, we devise and validate a closed-loop control strategy to disrupt synchronization caused by coupling. Additionally, we identify strategies for implementing the proposed control strategy in situations where underlying model equations are unavailable by estimating the necessary terms of the reduced order equations in real-time from observables.

     
    more » « less
  2. Essential tremor (ET) is among the most prevalent movement disorders, but its origins are elusive. The inferior olivary nucleus (ION) has been hypothesized as the prime generator of tremor because of the pacemaker properties of ION neurons, but structural and functional changes in ION are unlikely under ET. Abnormalities have instead been reported in the cerebello-thalamo-cortical network, including dysfunctions of the GABAergic projections from the cerebellar cortex to the dentate nucleus. It remains unclear, though, how tremor would relate to a dysfunction of cerebellar connectivity. To address this question, we built a computational model of the cortico-cerebello-thalamo-cortical loop. We simulated the effects of a progressive loss of GABA A α 1 -receptor subunits and up-regulation of α 2/3 -receptor subunits in the dentate nucleus, and correspondingly, we studied the evolution of the firing patterns along the loop. The model closely reproduced experimental evidence for each structure in the loop. It showed that an alteration of amplitudes and decay times of the GABAergic currents to the dentate nucleus can facilitate sustained oscillatory activity at tremor frequency throughout the network as well as a robust bursting activity in the thalamus, which is consistent with observations of thalamic tremor cells in ET patients. Tremor-related oscillations initiated in small neural populations and spread to a larger network as the synaptic dysfunction increased, while thalamic high-frequency stimulation suppressed tremor-related activity in thalamus but increased the oscillation frequency in the olivocerebellar loop. These results suggest a mechanism for tremor generation under cerebellar dysfunction, which may explain the origin of ET. 
    more » « less
  3. Introduction:Current brain-computer interfaces (BCIs) primarily rely on visual feedback. However, visual feedback may not be sufficient for applications such as movement restoration, where somatosensory feedback plays a crucial role. For electrocorticography (ECoG)-based BCIs, somatosensory feedback can be elicited by cortical surface electro-stimulation [1]. However, simultaneous cortical stimulation and recording is challenging due to stimulation artifacts. Depending on the orientation of stimulating electrodes, their distance to the recording site, and the stimulation intensity, these artifacts may overwhelm the neural signals of interest and saturate the recording bioamplifiers, making it impossible to recover the underlying information [2]. To understand how these factors affect artifact propagation, we performed a preliminary characterization of ECoG signals during cortical stimulation.Materials/Methods/ResultsECoG electrodes were implanted in a 39-year old epilepsy patient as shown in Fig. 1. Pairs of adjacent electrodes were stimulated as a part of language cortical mapping. For each stimulating pair, a charge-balanced biphasic square pulse train of current at 50 Hz was delivered for five seconds at 2, 4, 6, 8 and 10 mA. ECoG signals were recorded at 512 Hz. The signals were then high-pass filtered (≥1.5 Hz, zero phase), and the 5-second stimulation epochs were segmented. Within each epoch, artifact-induced peaks were detected for each electrode, except the stimulating pair, where signals were clipped due to amplifier saturation. These peaks were phase-locked across electrodes and were 20 ms apart, thus matching the pulse train frequency. The response was characterized by calculating the median peak within the 5-second epochs. Fig. 1 shows a representative response of the right temporal grid (RTG), with the stimulation channel at RTG electrodes 14 and 15. It also shows a hypothetical amplifier saturation contour of an implantable, bi-directional, ECoG-based BCI prototype [2], assuming the supply voltage of 2.2 V and a gain of 66 dB. Finally, we quantify the worstcase scenario by calculating the largest distance between the saturation contour and the midpoint of each stimulating channel.Discussion:Our results indicate that artifact propagation follows a dipole potential distribution with the extent of the saturation region (the interior of the white contour) proportional to the stimulation amplitude. In general, the artifacts propagated farthest when a 10 mA current was applied with the saturation regions extending from 17 to 32 mm away from the midpoint of the dipole. Consistent with the electric dipole model, this maximum spread happened along the direction of the dipole moment. An exception occurred at stimulation channel RTG11-16, for which an additional saturation contour emerged away from the dipole contour (not shown), extending the saturation region to 41 mm. Also, the worst-case scenario was observed at 6 mA stimulation amplitude. This departure could be a sign of a nonlinear, switch-like behavior, wherein additional conduction pathways could become engaged in response to sufficiently high stimulation.Significance:While ECoG stimulation is routinely performed in the clinical setting, quantitative studies of the resulting signals are lacking. Our preliminary study demonstrates that stimulation artifacts largely obey dipole distributions, suggesting that the dipole model could be used to predict artifact propagation. Further studies are necessary to ascertain whether these results hold across other subjects and combinations of stimulation/recording grids. Once completed, these studies will reveal practical design constraints for future implantable bi-directional ECoG-based BCIs. These include parameters such as the distances between and relative orientations of the stimulating and recording electrodes, the choice of the stimulating electrodes, the optimal placement of the reference electrode, and the maximum stimulation amplitude. These findings would also have important implications for the design of custom, low-power bioamplifiers for implantable bi-directional ECoG-based BCIs.References:[1] Hiremath, S. V., et al. "Human perception of electrical stimulation on the surface of somatosensory cortex." PloS one 12.5 (2017): e0176020.[2] Rouse, A. G., et al. "A chronic generalized bi-directional brain-machine interface." Journal of Neural Engineering 8.3 (2011): 036018 
    more » « less
  4. SUMMARY

    Protracted episodes of 0.5–7 Hz pre-eruptive volcanic tremor (PVT) are common at active stratovolcanoes. Reliable links to processes related to magma movement consequently enable a potential to use properties of PVT as diagnostic eruptive precursors. A challenging feature of PVT is that generic spectral and amplitude properties of the signal evolve similarly, independent of widely varying volcano structures and conduit geometries on which most physical models rely. The ‘magma wagging’ model introduced in Jellinek & Bercovici (2011) and extended by Bercovici et al. (2013), Liao et al. and Liao & Bercovici (2018) makes progress because it depends on magma dynamics that are only weakly sensitive to volcano architecture: The flow of gas through a permeable foamy annulus of gas bubbles excites, modulates and maintains a wagging oscillation of a central magma column rising in an erupting conduit. ‘Magma wagging’ and resulting PVT are driven through an energy transfer from a ‘Bernoulli mode’ related to azimuthal variations in annular gas flow speeds. Consistent with observations, spectral and amplitude properties of PVT are predicted to evolve before an eruption as the width of the annulus decreases with increased gas fluxes. To confirm this critical Bernoulli-to-wagging energy transfer we use extensive experiments and restricted numerical simulations on wagging oscillations excited on analogue viscoelastic columns by annular air flows. We also explore sensitivities of the spatial and temporal characters of wagging to asymmetric annular air flows that are intractable in the existing magma wagging model and expected to occur in nature with spatial variations in annulus permeability. From high-resolution time-series of linear and orbital displacements of analogue column tops and time-series of axial deflections and accelerations of the column centre line, we characterize the excitation, evolution, and steady-state oscillations in unprecedented detail over a broad range of conditions. We show that the Bernoulli mode corresponds to the timescale for the buildup of axial elastic bending stresses in response to pressure variations related to air flows over the heights of columns. We identify three distinct wagging modes: (i) rotational (cf. Liao et al. 2018); (ii) mixed-mode and (iii) chaotic. Rotational modes are favoured for symmetric, high intensity forcing and a maximal delivery of mechanical energy to the fundamental magma wagging mode. Mixed-mode oscillations regimes are favoured for a symmetric, intermediate intensity forcing. Chaotic modes, involving the least efficient delivery of energy to the fundamental mode, occur for asymmetric forcing and where the intensity of imposed airflow is low. Numerical simulations also show that where forcing frequencies are comparable to a natural mode of free oscillation, power delivered by peripheral air flows is concentrated at the lowest frequency fundamental mode generally and spread among higher frequency natural modes where air pressure and column elastic forces are comparable. Our combined experimental and numerical results make qualitative predictions for the evolution of the character of volcanic tremor and its expression in seismic or infrasound arrays during natural events that is testable in field-based studies of PVT and syn-eruptive volcanic tremor.

     
    more » « less
  5. Cortical stimulationviaelectrocorticography (ECoG) may be an effective method for inducing artificial sensation in bi-directional brain-computer interfaces (BD-BCIs). However, strong electrical artifacts caused by electrostimulation may significantly degrade or obscure neural information. A detailed understanding of stimulation artifact propagation through relevant tissues may improve existing artifact suppression techniques or inspire the development of novel artifact mitigation strategies. Our work thus seeks to comprehensively characterize and model the propagation of artifacts in subdural ECoG stimulation. To this end, we collected and analyzed data from eloquent cortex mapping procedures of four subjects with epilepsy who were implanted with subdural ECoG electrodes. From this data, we observed that artifacts exhibited phase-locking and ratcheting characteristics in the time domain across all subjects. In the frequency domain, stimulation caused broadband power increases, as well as power bursts at the fundamental stimulation frequency and its super-harmonics. The spatial distribution of artifacts followed the potential distribution of an electric dipole with a median goodness-of-fit ofR2= 0.80 across all subjects and stimulation channels. Artifacts as large as ±1,100 μV appeared anywhere from 4.43 to 38.34 mm from the stimulation channel. These temporal, spectral and spatial characteristics can be utilized to improve existing artifact suppression techniques, inspire new strategies for artifact mitigation, and aid in the development of novel cortical stimulation protocols. Taken together, these findings deepen our understanding of cortical electrostimulation and provide critical design specifications for future BD-BCI systems.

     
    more » « less