Real-time brain oscillation detection and phase-locked stimulation using autoregressive spectral estimation and time-series forward prediction

Novel cyclic homogeneous oscillation detection method for high accuracy and specific characterization of neural dynamics

Determining the presence and frequency of neural oscillations is essential to understanding dynamic brain function. Traditional methods that detect peaks over 1/f noise within the power spectrum fail to distinguish between the fundamental frequency and harmonics of often highly non-sinusoidal neural oscillations. To overcome this limitation, we define fundamental criteria that characterize neural oscillations and introduce the cyclic homogeneous oscillation (CHO) detection method. We implemented these criteria based on an autocorrelation approach to determine an oscillation's fundamental frequency. We evaluated CHO by verifying its performance on simulated non-sinusoidal oscillatory bursts and validated its ability to determine the fundamental frequency of neural oscillations in electrocorticographic (ECoG), electroencephalographic (EEG), and stereoelectroencephalographic (SEEG) signals recorded from 27 human subjects. Our results demonstrate that CHO outperforms conventional techniques in accurately detecting oscillations. In summary, CHO demonstrates high precision and specificity in detecting neural oscillations in time and frequency domains. The method's specificity enables the detailed study of non-sinusoidal characteristics of oscillations, such as the degree of asymmetry and waveform of an oscillation. Furthermore, CHO can be applied to identify how neural oscillations govern interactions throughout the brain and to determine oscillatory biomarkers that index abnormal brain function.

Optimizing real-time phase detection in diverse rhythmic biological signals for phase-specific neuromodulation

Closed-loop, phase-specific neurostimulation is a powerful method to modulate ongoing brain activity for clinical and research applications. Phase-specific stimulation relies on estimating the phase of an ongoing oscillation in real time and issuing a control command at a target phase. Phase detection algorithms based on Fast Fourier transform (FFT) are widely used due to their computational efficiency and robustness. However, it is unclear how algorithm performance depends on the spectral properties of the input signal and how algorithm parameters can be optimized. We used offline simulation to evaluate the performance of three algorithms (endpoint-corrected Hilbert Transform, Hilbert Transform and phase mapping) on three rhythmic biological signals with distinct spectral properties (rodent hippocampal theta potential, human EEG alpha and human essential tremor). First, we found that algorithm performance was more strongly influenced by signal amplitude and frequency variation compared with signal to noise ratio. Second, our simulations showed that the size of the data window for phase estimation was critical for the performance of FFT-based algorithms, where the optimal data window corresponds to the period of the oscillation. We validated this prediction with real time phase detection of hippocampal theta oscillations in freely behaving rats performing spatial navigation. Our findings define the relationship between signal properties and algorithm performance and provide a convenient method for optimizing FFT-based phase detection algorithms.

Novel Cyclic Homogeneous Oscillation Detection Method for High Accuracy and Specific Characterization of Neural Dynamics

Detecting temporal and spectral features of neural oscillations is essential to understanding dynamic brain function. Traditionally, the presence and frequency of neural oscillations are determined by identifying peaks over 1/f noise within the power spectrum. However, this approach solely operates within the frequency domain and thus cannot adequately distinguish between the fundamental frequency of a non-sinusoidal oscillation and its harmonics. Non-sinusoidal signals generate harmonics, significantly increasing the false-positive detection rate - a confounding factor in the analysis of neural oscillations. To overcome these limitations, we define the fundamental criteria that characterize a neural oscillation and introduce the Cyclic Homogeneous Oscillation (CHO) detection method that implements these criteria based on an auto-correlation approach that determines the oscillation's periodicity and fundamental frequency. We evaluated CHO by verifying its performance on simulated sinusoidal and non-sinusoidal oscillatory bursts convolved with 1/f noise. Our results demonstrate that CHO outperforms conventional techniques in accurately detecting oscillations. Specifically, we determined the sensitivity and specificity of CHO as a function of signal-to-noise ratio (SNR). We further assessed CHO by testing it on electrocorticographic (ECoG, 8 subjects) and electroencephalographic (EEG, 7 subjects) signals recorded during the pre-stimulus period of an auditory reaction time task and on electrocorticographic signals (6 SEEG subjects and 6 ECoG subjects) collected during resting state. In the reaction time task, the CHO method detected auditory alpha and pre-motor beta oscillations in ECoG signals and occipital alpha and pre-motor beta oscillations in EEG signals. Moreover, CHO determined the fundamental frequency of hippocampal oscillations in the human hippocampus during the resting state (6 SEEG subjects). In summary, CHO demonstrates high precision and specificity in detecting neural oscillations in time and frequency domains. The method's specificity enables the detailed study of non-sinusoidal characteristics of oscillations, such as the degree of asymmetry and waveform of an oscillation. Furthermore, CHO can be applied to identify how neural oscillations govern interactions throughout the brain and to determine oscillatory biomarkers that index abnormal brain function.

Frequency-dependent membrane polarization across neocortical cell types and subcellular elements by transcranial alternating current stimulation

Objective.Transcranial alternating current stimulation (tACS) is a non-invasive brain stimulation technique that directly interacts with ongoing brain oscillations in a frequency-dependent manner. However, it remains largely unclear how the cellular effects of tACS vary between cell types and subcellular elements.Approach.In this study, we use a set of morphologically realistic models of neocortical neurons to simulate the cellular response to uniform oscillating electric fields (EFs). We systematically characterize the membrane polarization in the soma, axons, and dendrites with varying field directions, intensities, and frequencies.Main results.Pyramidal cells are more sensitive to axial EF that is roughly parallel to the cortical column, while interneurons are sensitive to axial EF and transverse EF that is tangent to the cortical surface. Membrane polarization in each subcellular element increases linearly with EF intensity, and its slope, i.e. polarization length, highly depends on the stimulation frequency. At each frequency, pyramidal cells are more polarized than interneurons. Axons usually experience the highest polarization, followed by the dendrites and soma. Moreover, a visible frequency resonance presents in the apical dendrites of pyramidal cells, while the other subcellular elements primarily exhibit low-pass filtering properties. In contrast, each subcellular element of interneurons exhibits complex frequency-dependent polarization. Polarization phase in each subcellular element of cortical neurons lags that of field and exhibits high-pass filtering properties. These results demonstrate that the membrane polarization is not only frequency-dependent, but also cell type- and subcellular element-specific. Through relating effective length and ion mechanism with polarization, we emphasize the crucial role of cell morphology and biophysics in determining the frequency-dependent membrane polarization.Significance.Our findings highlight the diverse polarization patterns across cell types as well as subcellular elements, which provide some insights into the tACS cellular effects and should be considered when understanding the neural spiking activity by tACS.

Targeting motor cortex high-excitability states defined by functional connectivity with real-time EEG-TMS

We tested previous post-hoc findings indicating a relationship between functional connectivity (FC) in the motor network and corticospinal excitability (CsE), in a real-time EEG-TMS experiment in healthy participants. We hypothesized that high FC between left and right motor cortex predicts high CsE. FC was quantified in real-time by single-trial phase-locking value (stPLV), and TMS single pulses were delivered based on the current FC. CsE was indexed by motor-evoked potential (MEP) amplitude in a hand muscle. Possible confounding factors (pre-stimulus μ-power and phase, interstimulus interval) were evaluated post hoc. MEPs were significantly larger during high FC compared to low FC. Post hoc analysis revealed that the FC condition showed a significant interaction with μ-power in the stimulated hemisphere. Further, inter-stimulus interval (ISI) interacted with high vs. low FC conditions. In summary, FC was confirmed to be predictive of CsE, but should not be considered in isolation from μ-power and ISI. Moreover, FC was complementary to μ-phase in predicting CsE. Motor network FC is another marker of real-time accessible CsE beyond previously established markers, in particular phase and power of the μ rhythm, and may help define a more robust composite biomarker of high/low excitability states of human motor cortex.

Repetitive sensorimotor mu-alpha phase-targeted afferent stimulation produces no phase-dependent plasticity related changes in somatosensory evoked potentials or sensory thresholds

Phase-dependent plasticity has been proposed as a neurobiological mechanism by which oscillatory phase-amplitude cross-frequency coupling mediates memory process in the brain. Mimicking this mechanism, real-time EEG oscillatory phase-triggered transcranial magnetic stimulation (TMS) has successfully induced LTP-like changes in corticospinal excitability in the human motor cortex. Here we asked whether EEG phase-triggered afferent stimulation alone, if repetitively applied to the peaks, troughs, or random phases of the sensorimotor mu-alpha rhythm, would be sufficient to modulate the strength of thalamocortical synapses as assessed by changes in somatosensory evoked potential (SEP) N20 and P25 amplitudes and sensory thresholds (ST). Specifically, we applied 100 Hz triplets of peripheral electrical stimulation (PES) to the thumb, middle, and little finger of the right hand in pseudorandomized trials, with the afferent input from each finger repetitively and consistently arriving either during the cortical mu-alpha trough or peak or at random phases. No significant changes in SEP amplitudes or ST were observed across the phase-dependent PES intervention. We discuss potential limitations of the study and argue that suboptimal stimulation parameter choices rather than a general lack of phase-dependent plasticity in thalamocortical synapses are responsible for this null finding. Future studies should further explore the possibility of phase-dependent sensory stimulation.

Multifunctional fibers enable modulation of cortical and deep brain activity during cognitive behavior in macaques

Recording and modulating neural activity in vivo enables investigations of the neurophysiology underlying behavior and disease. However, there is a dearth of translational tools for simultaneous recording and localized receptor-specific modulation. We address this limitation by translating multifunctional fiber neurotechnology previously only available for rodent studies to enable cortical and subcortical neural recording and modulation in macaques. We record single-neuron and broader oscillatory activity during intracranial GABA infusions in the premotor cortex and putamen. By applying state-space models to characterize changes in electrophysiology, we uncover that neural activity evoked by a working memory task is reshaped by even a modest local inhibition. The recordings provide detailed insight into the electrophysiological effect of neurotransmitter receptor modulation in both cortical and subcortical structures in an awake macaque. Our results demonstrate a first-time application of multifunctional fibers for causal studies of neuronal activity in behaving nonhuman primates and pave the way for clinical translation of fiber-based neurotechnology.

A wearable EEG system for closed-loop neuromodulation of sleep-related oscillations

Objective.Healthy sleep plays a critical role in general well-being. Enhancement of slow-wave sleep by targeting acoustic stimuli to particular phases of delta (0.5-2 Hz) waves has shown promise as a non-invasive approach to improve sleep quality. Closed-loop stimulation during other sleep phases targeting oscillations at higher frequencies such as theta (4-7 Hz) or alpha (8-12 Hz) could be another approach to realize additional health benefits. However, systems to track and deliver stimulation relative to the instantaneous phase of electroencephalogram (EEG) signals at these higher frequencies have yet to be demonstrated outside of controlled laboratory settings.Approach.Here we examine the feasibility of using an endpoint-corrected version of the Hilbert transform (ecHT) algorithm implemented on a headband wearable device to measure alpha phase and deliver phase-locked auditory stimulation during the transition from wakefulness to sleep, during which alpha power is greatest. First, the ecHT algorithm is implementedin silicoto evaluate the performance characteristics of this algorithm across a range of sleep-related oscillatory frequencies. Secondly, a pilot sleep study tests feasibility to use the wearable device by users in the home setting for measurement of EEG activity during sleep and delivery of real-time phase-locked stimulation.Main results.The ecHT is capable of computing the instantaneous phase of oscillating signals with high precision, allowing auditory stimulation to be delivered at the intended phases of neural oscillations with low phase error. The wearable system was capable of measuring sleep-related neural activity with sufficient fidelity for sleep stage scoring during the at-home study, and phase-tracking performance matched simulated results. Users were able to successfully operate the system independently using the companion smartphone app to collect data and administer stimulation, and presentation of auditory stimuli during sleep initiation did not negatively impact sleep onset.Significance.This study demonstrates the feasibility of closed-loop real-time tracking and neuromodulation of a range of sleep-related oscillations using a wearable EEG device. Preliminary results suggest that this approach could be used to deliver non-invasive neuromodulation across all phases of sleep.

Real-time low latency estimation of brain rhythms with deep neural networks

Objective.Neurofeedback and brain-computer interfacing technology open the exciting opportunity for establishing interactive closed-loop real-time communication with the human brain. This requires interpreting brain's rhythmic activity and generating timely feedback to the brain. Lower delay between neuronal events and the appropriate feedback increases the efficacy of such interaction. Novel more efficient approaches capable of tracking brain rhythm's phase and envelope are needed for scenarios that entail instantaneous interaction with the brain circuits.Approach.Isolating narrow-band signals incurs fundamental delays. To some extent they can be compensated using forecasting models. Given the high quality of modern time series forecasting neural networks we explored their utility for low-latency extraction of brain rhythm parameters. We tested five neural networks with conceptually distinct architectures in forecasting synthetic EEG rhythms. The strongest architecture was then trained to simultaneously filter and forecast EEG data. We compared it against the state-of-the-art techniques using synthetic and real data from 25 subjects.Main results.The temporal convolutional network (TCN) remained the strongest forecasting model that achieved in the majority of testing scenarios>90% rhythm's envelope correlation with<10 ms effective delay and<20∘circular standard deviation of phase estimates. It also remained stable enough to noise level perturbations. Trained to filter and predict the TCN outperformed the cFIR, the Kalman filter based state-space estimation technique and remained on par with the larger Conv-TasNet architecture.Significance.Here we have for the first time demonstrated the utility of the neural network approach for low-latency narrow-band filtering of brain activity signals. Our proposed approach coupled with efficient implementation enhances the effectiveness of brain-state dependent paradigms across various applications. Moreover, our framework for forecasting EEG signals holds promise for investigating the predictability of brain activity, providing valuable insights into the fundamental questions surrounding the functional organization and hierarchical information processing properties of the brain.

Closing the loop between brain and electrical stimulation: towards precision neuromodulation treatments

One of the most critical challenges in using noninvasive brain stimulation (NIBS) techniques for the treatment of psychiatric and neurologic disorders is inter- and intra-individual variability in response to NIBS. Response variations in previous findings suggest that the one-size-fits-all approach does not seem the most appropriate option for enhancing stimulation outcomes. While there is a growing body of evidence for the feasibility and effectiveness of individualized NIBS approaches, the optimal way to achieve this is yet to be determined. Transcranial electrical stimulation (tES) is one of the NIBS techniques showing promising results in modulating treatment outcomes in several psychiatric and neurologic disorders, but it faces the same challenge for individual optimization. With new computational and methodological advances, tES can be integrated with real-time functional magnetic resonance imaging (rtfMRI) to establish closed-loop tES-fMRI for individually optimized neuromodulation. Closed-loop tES-fMRI systems aim to optimize stimulation parameters based on minimizing differences between the model of the current brain state and the desired value to maximize the expected clinical outcome. The methodological space to optimize closed-loop tES fMRI for clinical applications includes (1) stimulation vs. data acquisition timing, (2) fMRI context (task-based or resting-state), (3) inherent brain oscillations, (4) dose-response function, (5) brain target trait and state and (6) optimization algorithm. Closed-loop tES-fMRI technology has several advantages over non-individualized or open-loop systems to reshape the future of neuromodulation with objective optimization in a clinically relevant context such as drug cue reactivity for substance use disorder considering both inter and intra-individual variations. Using multi-level brain and behavior measures as input and desired outcomes to individualize stimulation parameters provides a framework for designing personalized tES protocols in precision psychiatry.

The Effect of the Peristimulus α Phase on Visual Perception through Real-Time Phase-Locked Stimulus Presentation

The α phase has been theorized to reflect fluctuations in cortical excitability and thereby impose a cyclic influence on visual perception. Despite its appeal, this notion is not fully substantiated, as both supporting and opposing evidence has been recently reported. In contrast to previous research, this study examined the effect of the peristimulus instead of prestimulus phase on visual detection through a real-time phase-locked stimulus presentation (PLSP) approach. Specifically, we monitored phase data from magnetoencephalography (MEG) recordings over time, with a newly developed algorithm based on adaptive Kalman filtering (AKF). This information guided online presentations of masked stimuli that were phased-locked to different stages of the α cycle while healthy humans concurrently performed detection tasks. Behavioral evidence showed that the overall detection rate did not significantly vary according to the four predetermined peristimulus α phases. Nevertheless, the follow-up analyses highlighted that the phase at 90° relative to 180° likely enhanced detection. Corroborating neural parietal activity showed that early interaction between α phases and incoming stimuli orchestrated the neural representation of the hits and misses of the stimuli. This neural representation varied according to the phase and in turn shaped the behavioral outcomes. In addition to directly investigating to what extent fluctuations in perception can be ascribed to the α phases, this study suggests that phase-dependent perception is not as robust as previously presumed, and might also depend on how the stimuli are differentially processed as a result of a stimulus-phase interaction, in addition to reflecting alternations of the perceptual states between phases.

Delay Analysis in Closed-Loop EEG Phase-Triggered Transcranial Magnetic Stimulation

The recent development of closed-loop EEG phase-triggered transcranial magnetic stimulation (TMS) has advanced potential applications of adaptive neuromodulation based on the current brain state. Closed-loop TMS involves instantaneous acquisition of the EEG rhythm, timing prediction of the target phase, and triggering of TMS. However, the accuracy of EEG phase prediction algorithms is largely influenced by the system's transport delay, and their relationship is rarely considered in related work. This paper proposes a delay analysis that considers the delay of the closed-loop EEG phase-triggered TMS system as a primary factor in the validation of phase prediction algorithms. An in-silico validation using real EEG data was performed to compare the performance of commonly used algorithms. The experimental results indicate a significant influence of the total delay on the algorithm performance, and the performance ranking among algorithms varies at different levels of delay. We conclude that the delay analysis framework should be widely adopted in the design and validation of phase prediction algorithms for closed-loop EEG phase-triggered TMS systems.

Manipulating single-unit theta phase-locking with PhaSER: An open-source tool for real-time phase estimation and manipulation

The precise timing of neuronal spiking relative to the brain's endogenous oscillations (i.e., phase-locking or spike-phase coupling) has long been hypothesized to coordinate cognitive processes and maintain excitatory-inhibitory homeostasis. Indeed, disruptions in theta phase-locking have been described in models of neurological diseases with associated cognitive deficits and seizures, such as Alzheimer's disease, temporal lobe epilepsy, and autism spectrum disorders. However, due to technical limitations, determining if phase-locking causally contributes to these disease phenotypes has not been possible until recently. To fill this gap and allow for the flexible manipulation of single-unit phase-locking to on-going endogenous oscillations, we developed PhaSER, an open-source tool that allows for phase-specific manipulations. PhaSER can deliver optogenetic stimulation at defined phases of theta in order to shift the preferred firing phase of neurons relative to theta in real-time. Here, we describe and validate this tool in a subpopulation of inhibitory neurons that express somatostatin (SOM) in the CA1 and dentate gyrus (DG) regions of the dorsal hippocampus. We show that PhaSER is able to accurately deliver a photo-manipulation that activates opsin+ SOM neurons at specified phases of theta in real-time in awake, behaving mice. Further, we show that this manipulation is sufficient to alter the preferred firing phase of opsin+ SOM neurons without altering the referenced theta power or phase. All software and hardware requirements to implement real-time phase manipulations during behavior are available online (https://github.com/ShumanLab/PhaSER).

Optimization of phase prediction for brain-state dependent stimulation: a grid-search approach

Objective.Sources of heterogeneity in non-invasive brain stimulation literature can be numerous, with underlying brain states and protocol differences at the top of the list. Yet, incoherent results from brain-state-dependent stimulation experiments suggest that there are further factors adding to the variance. Hypothesizing that different signal processing pipelines might be partly responsible for heterogeneity; we investigated their effects on brain-state forecasting approaches.Approach.A grid-search was used to determine the fastest and most-accurate combination of preprocessing parameters and phase-forecasting algorithms. The grid-search was applied on a synthetic dataset and validated on electroencephalographic (EEG) data from a healthy (n= 18) and stroke (n= 31) cohort.Main results.Differences in processing pipelines led to different results; the grid-search chosen pipelines significantly increased the accuracy of published forecasting methods. The accuracy achieved in healthy was comparably high in stroke patients.Significance.This systematic offline analysis highlights the importance of the specific EEG processing and forecasting pipelines used for online state-dependent setups where precision in phase prediction is critical. Moreover, successful results in the stroke cohort pave the way to test state-dependent interventional treatment approaches.

Modulation effect of mouse hippocampal neural oscillations by closed-loop transcranial ultrasound stimulation

Objective. Closed-loop transcranial ultrasound stimulation (TUS) can be applied at a specific time according to the state of neural activity to achieve timely and precise neuromodulation and improve the modulation effect. In a previous study, we found that closed-loop TUS at the peaks and troughs of the theta rhythm in the mouse hippocampus was able to increase the absolute power and decrease the relative power of the theta rhythm of local field potentials (LFPs) independent of the peaks and troughs of the stimulus. However, it remained unclear whether the modulation effect of this closed-loop TUS-induced mouse hippocampal neural oscillation depended on the peaks and troughs of the theta rhythm.Approach. In this study, we used ultrasound with different stimulation modes and durations to stimulate the peaks (peak stimulation) and troughs (trough stimulation) of the hippocampal theta rhythm. The LFPs in the area of ultrasound stimulation were recorded and the amplitudes and power spectra of the theta rhythm before and after ultrasound stimulation were analyzed.Main results. The results showed that (a) the relative change in amplitude of theta rhythm decreases as the number of stimulation trials under peak stimulation increases; (b) the relative change in the absolute power of the theta rhythm decreases as the number of stimulation trials under peak stimulation increases; (c) the relative change in amplitude of the theta rhythm increases nonlinearly with the stimulation duration (SD) under peak stimulation, and; (d) the relative change in absolute power exhibits a nonlinear increase with SD under peak stimulation.Significance. These results suggest that the modulation effect of closed-loop TUS on theta rhythm is dependent on the stimulation mode and duration under peak stimulation. TUS has the potential to precisely modulate theta rhythm-related neural activity.

The Portiloop: A deep learning-based open science tool for closed-loop brain stimulation

Closed-loop brain stimulation refers to capturing neurophysiological measures such as electroencephalography (EEG), quickly identifying neural events of interest, and producing auditory, magnetic or electrical stimulation so as to interact with brain processes precisely. It is a promising new method for fundamental neuroscience and perhaps for clinical applications such as restoring degraded memory function; however, existing tools are expensive, cumbersome, and offer limited experimental flexibility. In this article, we propose the Portiloop, a deep learning-based, portable and low-cost closed-loop stimulation system able to target specific brain oscillations. We first document open-hardware implementations that can be constructed from commercially available components. We also provide a fast, lightweight neural network model and an exploration algorithm that automatically optimizes the model hyperparameters to the desired brain oscillation. Finally, we validate the technology on a challenging test case of real-time sleep spindle detection, with results comparable to off-line expert performance on the Massive Online Data Annotation spindle dataset (MODA; group consensus). Software and plans are available to the community as an open science initiative to encourage further development and advance closed-loop neuroscience research [https://github.com/Portiloop].

Seizure Forecasting by High-Frequency Activity (80-170 Hz) in Long-term Continuous Intracranial EEG Recordings

BACKGROUND AND OBJECTIVES

Reliable seizure forecasting has important implications in epilepsy treatment and improving the quality of lives for people with epilepsy. High-frequency activity (HFA) is a biomarker that has received significant attention over the past 2 decades, but its predictive value in seizure forecasting remains uncertain. This work aimed to determine the utility of HFA in seizure forecasting.

METHODS

We used seizure data and HFA (80-170 Hz) data obtained from long-term, continuous intracranial EEG recordings of patients with drug-resistant epilepsy. Instantaneous rates and phases of HFA cycles were used as features for seizure forecasting. Seizure forecasts based on each individual HFA feature, and with the use of a combined approach, were generated pseudo-prospectively (causally). To compute the instantaneous phases for pseudo-prospective forecasting, real-time phase estimation based on an autoregressive model was used. Features were combined with a weighted average approach. The performance of seizure forecasting was primarily evaluated by the area under the curve (AUC).

RESULTS

Of 15 studied patients (median recording duration 557 days, median seizures 151), 12 patients with >10 seizures after 100 recording days were included in the pseudo-prospective analysis. The presented real-time phase estimation is feasible and can causally estimate the instantaneous phases of HFA cycles with high accuracy. Pseudo-prospective seizure forecasting based on HFA rates and phases performed significantly better than chance in 11 of 12 patients, although there were patient-specific differences. Combining rate and phase information improved forecasting performance compared to using either feature alone. The combined forecast using the best-performing channel yielded a median AUC of 0.70, a median sensitivity of 0.57, and a median specificity of 0.77.

DISCUSSION

These findings show that HFA could be useful for seizure forecasting and represent proof of concept for using prior information of patient-specific relationships between HFA and seizures in pseudo-prospective forecasting. Future seizure forecasting algorithms might benefit from the inclusion of HFA, and the real-time phase estimation approach can be extended to other biomarkers.

CLASSIFICATION OF EVIDENCE

This study provides Class IV evidence that HFA (80-170 Hz) in long-term continuous intracranial EEG can be useful to forecast seizures in patients with refractory epilepsy.

Temporally Targeted Interactions With Pathologic Oscillations as Therapeutical Targets in Epilepsy and Beyond

Self-organized neuronal oscillations rely on precisely orchestrated ensemble activity in reverberating neuronal networks. Chronic, non-malignant disorders of the brain are often coupled to pathological neuronal activity patterns. In addition to the characteristic behavioral symptoms, these disturbances are giving rise to both transient and persistent changes of various brain rhythms. Increasing evidence support the causal role of these "oscillopathies" in the phenotypic emergence of the disease symptoms, identifying neuronal network oscillations as potential therapeutic targets. While the kinetics of pharmacological therapy is not suitable to compensate the disease related fine-scale disturbances of network oscillations, external biophysical modalities (e.g., electrical stimulation) can alter spike timing in a temporally precise manner. These perturbations can warp rhythmic oscillatory patterns via resonance or entrainment. Properly timed phasic stimuli can even switch between the stable states of networks acting as multistable oscillators, substantially changing the emergent oscillatory patterns. Novel transcranial electric stimulation (TES) approaches offer more reliable neuronal control by allowing higher intensities with tolerable side-effect profiles. This precise temporal steerability combined with the non- or minimally invasive nature of these novel TES interventions make them promising therapeutic candidates for functional disorders of the brain. Here we review the key experimental findings and theoretical background concerning various pathological aspects of neuronal network activity leading to the generation of epileptic seizures. The conceptual and practical state of the art of temporally targeted brain stimulation is discussed focusing on the prevention and early termination of epileptic seizures.

Phase- targeted stimulation modulates phase-amplitude coupling in the motor cortex of the human brain

BACKGROUND

Phase-amplitude coupling (PAC) in which the amplitude of a faster field potential oscillation is coupled to the phase of a slower rhythm, is one of the most well-studied interactions between oscillations at different frequency bands. In a healthy brain, PAC accompanies cognitive functions such as learning and memory, and changes in PAC have been associated with neurological diseases including Parkinson's disease (PD), schizophrenia, obsessive-compulsive disorder, Alzheimer's disease, and epilepsy.

OBJECTIVE

/Hypothesis: In PD, normalization of PAC in the motor cortex has been reported in the context of effective treatments such as dopamine replacement therapy and deep brain stimulation (DBS), but the possibility of normalizing PAC through intervention at the cortex has not been shown in humans. Phase-targeted stimulation (PDS) has a strong potential to modulate PAC levels and potentially normalize it.

METHODS

We applied stimulation pulses triggered by specific phases of the beta oscillations, the low frequency oscillations that define phase of gamma amplitude in beta-gamma PAC, to the motor cortex of seven PD patients at rest during DBS lead placement surgery We measured the effect on PAC modulation in the motor cortex relative to stimulation-free periods.

RESULTS

We describe a system for phase-targeted stimulation locked to specific phases of a continuously updated slow local field potential oscillation (in this case, beta band oscillations) prediction. Stimulation locked to the phase of the peak of beta oscillations increased beta-gamma coupling both during and after stimulation in the motor cortex, and the opposite phase (trough) stimulation reduced the magnitude of coupling after stimulation.

CONCLUSION

These results demonstrate the capacity of cortical phase-targeted stimulation to modulate PAC without evoking motor activation, which could allow applications in the treatment of neurological disorders associated with abnormal PAC, such as PD.

A state-informed stimulation approach with real-time estimation of the instantaneous phase of neural oscillations by a Kalman filter

Objective.We propose a novel method to estimate the instantaneous oscillatory phase to implement a real-time system for state-informed sensory stimulation in electroencephalography (EEG) experiments.Approach.The method uses Kalman filter-based prediction to estimate current and future EEG signals. We tested the performance of our method in a real-time situation.Main results.Our method showed higher accuracy in predicting the EEG phase than the conventional autoregressive (AR) model-based method.Significance.A Kalman filter allows us to easily estimate the instantaneous phase of EEG oscillations based on the automatically estimated AR model implemented in a real-time signal processing machine. The proposed method has a potential for versatile applications targeting the modulation of EEG phase dynamics and the plasticity of brain networks in relation to perceptual or cognitive functions.

Causal roles of prefrontal cortex during spontaneous perceptual switching are determined by brain state dynamics

The prefrontal cortex (PFC) is thought to orchestrate cognitive dynamics. However, in tests of bistable visual perception, no direct evidence supporting such presumable causal roles of the PFC has been reported except for a recent work. Here, using a novel brain-state-dependent neural stimulation system, we identified causal effects on percept dynamics in three PFC activities—right frontal eye fields, dorsolateral PFC (DLPFC), and inferior frontal cortex (IFC). The causality is behaviourally detectable only when we track brain state dynamics and modulate the PFC activity in brain-state-/state-history-dependent manners. The behavioural effects are underpinned by transient neural changes in the brain state dynamics, and such neural effects are quantitatively explainable by structural transformations of the hypothetical energy landscapes. Moreover, these findings indicate distinct functions of the three PFC areas: in particular, the DLPFC enhances the integration of two PFC-active brain states, whereas IFC promotes the functional segregation between them. This work resolves the controversy over the PFC roles in spontaneous perceptual switching and underlines brain state dynamics in fine investigations of brain-behaviour causality.

A cube that seems to shift its spatial arrangement as you keep looking; the elegant silhouette of a pirouetting dancer, which starts to spin in the opposite direction the more you stare at it; an illustration that shows two profiles – or is it a vase? These optical illusions are examples of bistable visual perception. Beyond their entertaining aspect, they provide a way for scientists to explore the dynamics of human consciousness, and the neural regions involved in this process.

Some studies show that bistable visual perception is associated with the activation of the prefrontal cortex, a brain area involved in complex cognitive processes. However, it is unclear whether this region is required for the illusions to emerge. Some research has showed that even if sections of the prefrontal cortex are temporally deactivated, participants can still experience the illusions.

Instead, Takamitsu Watanabe proposes that bistable visual perception is a process tied to dynamic brain states – that is, that distinct regions of the prefontal cortex are required for this fluctuating visual awareness, depending on the state of the whole brain. Such causal link cannot be observed if brain activity is not tracked closely.

To investigate this, the brain states of 65 participants were recorded as individuals were experiencing the optical illusions; the activity of their various brain regions could therefore be mapped, and then areas of the prefrontal cortex could precisely be inhibited at the right time using transcranial magnetic stimulation. This revealed that, indeed, prefrontal cortex regions were necessary for bistable visual perception, but not in a simple way. Instead, which ones were required and when depended on activity dynamics taking place in the whole brain. Overall, these results indicate that monitoring brain states is necessary to better understand – and ultimately, control – the neural pathways underlying perception and behaviour.

A state space modeling approach to real-time phase estimation

Brain rhythms have been proposed to facilitate brain function, with an especially important role attributed to the phase of low-frequency rhythms. Understanding the role of phase in neural function requires interventions that perturb neural activity at a target phase, necessitating estimation of phase in real-time. Current methods for real-time phase estimation rely on bandpass filtering, which assumes narrowband signals and couples the signal and noise in the phase estimate, adding noise to the phase and impairing detections of relationships between phase and behavior. To address this, we propose a state space phase estimator for real-time tracking of phase. By tracking the analytic signal as a latent state, this framework avoids the requirement of bandpass filtering, separately models the signal and the noise, accounts for rhythmic confounds, and provides credible intervals for the phase estimate. We demonstrate in simulations that the state space phase estimator outperforms current state-of-the-art real-time methods in the contexts of common confounds such as broadband rhythms, phase resets, and co-occurring rhythms. Finally, we show applications of this approach to in vivo data. The method is available as a ready-to-use plug-in for the Open Ephys acquisition system, making it widely available for use in experiments.

Real-time estimation of phase and amplitude with application to neural data

Computation of the instantaneous phase and amplitude via the Hilbert Transform is a powerful tool of data analysis. This approach finds many applications in various science and engineering branches but is not proper for causal estimation because it requires knowledge of the signal's past and future. However, several problems require real-time estimation of phase and amplitude; an illustrative example is phase-locked or amplitude-dependent stimulation in neuroscience. In this paper, we discuss and compare three causal algorithms that do not rely on the Hilbert Transform but exploit well-known physical phenomena, the synchronization and the resonance. After testing the algorithms on a synthetic data set, we illustrate their performance computing phase and amplitude for the accelerometer tremor measurements and a Parkinsonian patient's beta-band brain activity.

Prediction of cortical theta oscillations in humans for phase-locked visual stimulation

BACKGROUND

The timing of an event within an oscillatory phase is considered to be one of the key strategies used by the brain to code and process neural information. Whereas existing methods of studying this phenomenon are chiefly based on retrospective analysis of electroencephalography (EEG) data, we now present a method to study it prospectively. New method: We present a system that allows for the delivery of visual stimuli at a specific phase of the cortical theta oscillation by fitting a sine to raw surface EEG data to estimate and predict the phase. One noteworthy feature of the method is that it can minimize potentially confounding effects of previous trials by using only a short sequence of past data.

RESULTS

In a trial with 10 human participants we achieved a significant phase locking with an inter-trial phase coherence of 0.39. We demonstrated successful phase locking on synthetic signals with a signal-to-noise ratio of less than - 20 dB. Comparison with existing method(s): We compared the new method to an autoregressive method published in the literature and found the new method was superior in mean phase offset, circular standard deviation, and prediction latency.

CONCLUSIONS

By fitting sine waves to raw EEG traces, we locked visual stimuli to arbitrary phases within the theta oscillatory cycle of healthy humans.

Phase-dependent offline enhancement of human motor memory

BACKGROUND

Skill learning engages offline activity in the primary motor cortex (M1). Sensorimotor cortical activity oscillates between excitatory trough and inhibitory peak phases of the mu (8-12 Hz) rhythm. We recently showed that these mu phases influence the magnitude and direction of neuroplasticity induction within M1. However, the contribution of M1 activity during mu peak and trough phases to human skill learning has not been investigated.

OBJECTIVE

To evaluate the effects of phase-dependent TMS during mu peak and trough phases on offline learning of a newly-acquired motor skill.

METHODS

On Day 1, three groups of healthy adults practiced an explicit motor sequence learning task with their non-dominant left hand. After practice, phase-dependent TMS was applied to the right M1 during either mu peak or mu trough phases. The third group received sham TMS during random mu phases. On Day 2, all subjects were re-tested on the same task to evaluate offline learning.

RESULTS

Subjects who received phase-dependent TMS during mu trough phases showed increased offline skill learning compared to those who received phase-dependent TMS during mu peak phases or sham TMS during random mu phases. Additionally, phase-dependent TMS during mu trough phases elicited stronger whole-brain broadband oscillatory power responses than phase-dependent TMS during mu peak phases.

CONCLUSIONS

We conclude that sensorimotor mu trough phases reflect brief windows of opportunity during which TMS can strengthen newly-acquired skill memories.

EEG, MEG and neuromodulatory approaches to explore cognition: Current status and future directions

Neural oscillations and their association with brain states and cognitive functions have been object of extensive investigation over the last decades. Several electroencephalography (EEG) and magnetoencephalography (MEG) analysis approaches have been explored and oscillatory properties have been identified, in parallel with the technical and computational advancement. This review provides an up-to-date account of how EEG/MEG oscillations have contributed to the understanding of cognition. Methodological challenges, recent developments and translational potential, along with future research avenues, are discussed.

Non-invasive suppression of essential tremor via phase-locked disruption of its temporal coherence

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.

Real-Time Implementation of EEG Oscillatory Phase-Informed Visual Stimulation Using a Least Mean Square-Based AR Model

It is a technically challenging problem to assess the instantaneous brain state using electroencephalography (EEG) in a real-time closed-loop setup because the prediction of future signals is required to define the current state, such as the instantaneous phase and amplitude. To accomplish this in real-time, a conventional Yule–Walker (YW)-based autoregressive (AR) model has been used. However, the brain state-dependent real-time implementation of a closed-loop system employing an adaptive method has not yet been explored. Our primary purpose was to investigate whether time-series forward prediction using an adaptive least mean square (LMS)-based AR model would be implementable in a real-time closed-loop system or not. EEG state-dependent triggers synchronized with the EEG peaks and troughs of alpha oscillations in both an open-eyes resting state and a visual task. For the resting and visual conditions, statistical results showed that the proposed method succeeded in giving triggers at a specific phase of EEG oscillations for all participants. These individual results showed that the LMS-based AR model was successfully implemented in a real-time closed-loop system targeting specific phases of alpha oscillations and can be used as an adaptive alternative to the conventional and machine-learning approaches with a low computational load.

Phase-dependent Stimulation for Modulating Phase-amplitude Coupling: A Computational Modeling Approach

Phase-amplitude coupling (PAC), in which the amplitude of a faster neural oscillation couples to the phase of a slower rhythm, is one of the most common representations of complex neuronal rhythmic activities. In a healthy brain, PAC accompanies cognitive function, and abnormal patterns of PAC have been linked to several neurological disorders. Among the various brain neuromodulation techniques, phase-dependent stimulation has a strong potential to modulate PAC levels. In this study, we utilize a computational model in the NEURON environment based on a detailed mathematical model of neuronal populations, consisting of networks with both excitatory and inhibitory neurons, to simulate PAC generation. The model was then used to investigate the modulatory effects of phase-dependent stimulation on the generated PAC. Simulated data from the model shows that stimulation locked to the phase of slower rhythms increased PAC level during stimulation. These results demonstrate the capacity of phase-dependent stimulation to modulate PAC, which could allow for applications in the treatment of neurological disorders associated with abnormal PAC, such as Parkinson's disease.Clinical Relevance- Analyzing the origins of neuronal PAC and developing a brain stimulation technique for modulating the level of PAC can facilitate the development of novel treatment methods for neurological disorders associated with abnormal cross-frequency coupling.

A user-friendly algorithm for adaptive closed-loop phase-locked stimulation

BACKGROUND

Closed-loop phase-locked stimulation experiments are rare due to the unavailability of user-friendly algorithms and devices. Our goal is to provide an algorithm for the detection of oscillatory activity in local field potentials (LFPs) and phase prediction, which is user-friendly and robust to non-stationarities in LFPs of behaving animals.

NEW METHOD

We propose an algorithm that only requires specification of the frequency range within which oscillatory episodes are tracked. Frequency-specific detection thresholds and filter parameters are adjusted automatically based on the short-time LFP power spectrum. Estimates of instantaneous frequency and instantaneous phase are used for phase extrapolation, taking advantage of Bayesian estimation. We used real LFP signals, recorded from a variety of different species and different brain areas, as well as artificial LFP signals with known properties to assess the detection and prediction performance of our algorithm and three previously published reference algorithms under various conditions.

RESULTS AND COMPARISON WITH EXISTING METHODS

Our algorithm, while significantly more user-friendly than previous approaches, provides a solid detection and prediction performance over a wide range of realistic conditions and, in many cases, has a longer prediction horizon than the reference algorithms. Due to its ability to adjust to changes in the signal, the algorithm is well-prepared to deal with non-stationarities in oscillation frequency, even in the presence of multiple oscillation components.

CONCLUSIONS

We have created a universal algorithm for oscillation detection and phase prediction, which performs well and is user-friendly at the same time, making closed-loop phase-locked stimulation experiments easier to accomplish.

Perspective: Phase Amplitude Coupling-Based Phase-Dependent Neuromodulation in Parkinson's Disease

Deep brain stimulation (DBS) is an effective surgical therapy for Parkinson's disease (PD). However, limitations of the DBS systems have led to great interest in adaptive neuromodulation systems that can dynamically adjust stimulation parameters to meet concurrent therapeutic demand. Constant high-frequency motor cortex stimulation has not been remarkably efficacious, which has led to greater focus on modulation of subcortical targets. Understanding of the importance of timing in both cortical and subcortical stimulation has generated an interest in developing more refined, parsimonious stimulation techniques based on critical oscillatory activities of the brain. Concurrently, much effort has been put into identifying biomarkers of both parkinsonian and physiological patterns of neuronal activities to drive next generation of adaptive brain stimulation systems. One such biomarker is beta-gamma phase amplitude coupling (PAC) that is detected in the motor cortex. PAC is strongly correlated with parkinsonian specific motor signs and symptoms and respond to therapies in a dose-dependent manner. PAC may represent the overall state of the parkinsonian motor network and have less instantaneously dynamic fluctuation during movement. These findings raise the possibility of novel neuromodulation paradigms that are potentially less invasiveness than DBS. Successful application of PAC in neuromodulation may necessitate phase-dependent stimulation technique, which aims to deliver precisely timed stimulation pulses to a specific phase to predictably modulate to selectively modulate pathological network activities and behavior in real time. Overcoming current technical challenges can lead to deeper understanding of the parkinsonian pathophysiology and development of novel neuromodulatory therapies with potentially less side-effects and higher therapeutic efficacy.

Digital filters for low-latency quantification of brain rhythms in real time

OBJECTIVE

The rapidly developing paradigm of closed-loop neuroscience has extensively employed brain rhythms as the signal forming real-time neurofeedback, triggering brain stimulation, or governing stimulus selection. However, the efficacy of brain rhythm contingent paradigms suffers from significant delays related to the process of extraction of oscillatory parameters from broad-band neural signals with conventional methods. To this end, real-time algorithms are needed that would shorten the delay while maintaining an acceptable speed-accuracy trade-off.

APPROACH

Here we evaluated a family of techniques based on the application of the least-squares complex-valued filter (LSCF) design to real-time quantification of brain rhythms. These techniques allow for explicit optimization of the speed-accuracy trade-off when quantifying oscillatory patterns. We used EEG data collected from 10 human participants to systematically compare LSCF approach to the other commonly used algorithms. Each method being evaluated was optimized by scanning through the grid of its hyperparameters using independent data samples.

MAIN RESULTS

When applied to the task of estimating oscillatory envelope and phase, the LSCF techniques outperformed in speed and accuracy both conventional Fourier transform and rectification based methods as well as more advanced techniques such as those that exploit autoregressive extrapolation of narrow-band filtered signals. When operating at zero latency, the weighted LSCF approach yielded 75% accuracy when detecting alpha-activity episodes, as defined by the amplitude crossing of the 95th-percentile threshold.

SIGNIFICANCE

The LSCF approaches are easily applicable to low-delay quantification of brain rhythms. As such, these methods are useful in a variety of neurofeedback, brain-computer-interface and other experimental paradigms that require rapid monitoring of brain rhythms.

Experimental evaluation of methods for real-time EEG phase-specific transcranial magnetic stimulation

OBJECTIVE

Real-time approaches for transcranial magnetic stimulation (TMS) based on a specific EEG phase are a promising avenue for more precise neuromodulation interventions. However, optimal approaches to reliably extract the EEG phase in a frequency band of interest to inform TMS are still to be identified. Here, we implement a new real-time phase detection method for closed-loop EEG-TMS for robust phase extraction. We compare this algorithm with state-of-the-art methods and evaluate its performance both in silico and experimentally.

APPROACH

We propose a new robust algorithm (Educated Temporal Prediction) for delivering real-time EEG phase-specific stimulation based on short prerecorded EEG training data. This method estimates the interpeak period from a training period and applies a bias correction to predict future peaks. We compare the accuracy and computation speed of the ETP algorithm with two existing methods (Fourier based, Autoregressive Prediction) using prerecorded resting EEG data and real-time experiments.

MAIN RESULTS

We found that Educated Temporal Prediction performs with higher accuracy than Fourier-based or Autoregressive methods both in silico and in vivo while being computationally more efficient. Further, we document the dependency of the EEG signal-to-noise ratio (SNR) on algorithm accuracy across all algorithms.

SIGNIFICANCE

Our results give important insights for real-time EEG-TMS technical development as well as experimental design. Due to its robustness and computational efficiency, our method can find broad use in experimental research or clinical applications. Through open sharing of code for all three methods, we enable broad access of TMS-EEG real-time algorithms to the community.

The shaky ground truth of real-time phase estimation

Instantaneous phase of brain oscillations in electroencephalography (EEG) is a measure of brain state that is relevant to neuronal processing and modulates evoked responses. However, determining phase at the time of a stimulus with standard signal processing methods is not possible due to the stimulus artifact masking the future part of the signal. Here, we quantify the degree to which signal-to-noise ratio and instantaneous amplitude of the signal affect the variance of phase estimation error and the precision with which "ground truth" phase is even defined, using both the variance of equivalent estimators and realistic simulated EEG data with known synthetic phase. Necessary experimental conditions are specified in which pre-stimulus phase estimation is meaningfully possible based on instantaneous amplitude and signal-to-noise ratio of the oscillation of interest. An open source toolbox is made available for causal (using pre-stimulus signal only) phase estimation along with a EEG dataset consisting of recordings from 140 participants and a best practices workflow for algorithm optimization and benchmarking. As an illustration, post-hoc sorting of open-loop transcranial magnetic stimulation (TMS) trials according to pre-stimulus sensorimotor μ-rhythm phase is performed to demonstrate modulation of corticospinal excitability, as indexed by the amplitude of motor evoked potentials.

Time-Series Prediction of the Oscillatory Phase of EEG Signals Using the Least Mean Square Algorithm-Based AR Model

Neural oscillations are vital for the functioning of a central nervous system because they assist in brain communication across a huge network of neurons. Alpha frequency oscillations are believed to depict idling or inhibition of task-irrelevant cortical activities. However, recent studies on alpha oscillations (particularly alpha phase) hypothesize that they have an active and direct role in the mechanisms of attention and working memory. To understand the role of alpha oscillations in several cognitive processes, accurate estimations of phase, amplitude, and frequency are required. Herein, we propose an approach for time-series forward prediction by comparing an autoregressive (AR) model and an adaptive method (least mean square (LMS)-based AR model). This study tested both methods for two prediction lengths of data. Our results indicate that for shorter data segments (prediction of 128 ms), the AR model outperforms the LMS-based AR model, while for longer prediction lengths (256 ms), the LMS- based AR model surpasses the AR model. LMS with low computational cost can aid in electroencephalography (EEG) phase prediction (alpha oscillations) in basic research to reveal the functional role of the oscillatory phase as well as for applications for brain-computer interfaces.

Phase-dependence of response curves to deep brain stimulation and their relationship: from essential tremor patient data to a Wilson-Cowan model

Essential tremor manifests predominantly as a tremor of the upper limbs. One therapy option is high-frequency deep brain stimulation, which continuously delivers electrical stimulation to the ventral intermediate nucleus of the thalamus at about 130 Hz. Constant stimulation can lead to side effects, it is therefore desirable to find ways to stimulate less while maintaining clinical efficacy. One strategy, phase-locked deep brain stimulation, consists of stimulating according to the phase of the tremor. To advance methods to optimise deep brain stimulation while providing insights into tremor circuits, we ask the question: can the effects of phase-locked stimulation be accounted for by a canonical Wilson-Cowan model? We first analyse patient data, and identify in half of the datasets significant dependence of the effects of stimulation on the phase at which stimulation is provided. The full nonlinear Wilson-Cowan model is fitted to datasets identified as statistically significant, and we show that in each case the model can fit to the dynamics of patient tremor as well as to the phase response curve. The vast majority of top fits are stable foci. The model provides satisfactory prediction of how patient tremor will react to phase-locked stimulation by predicting patient amplitude response curves although they were not explicitly fitted. We also approximate response curves of the significant datasets by providing analytical results for the linearisation of a stable focus model, a simplification of the Wilson-Cowan model in the stable focus regime. We report that the nonlinear Wilson-Cowan model is able to describe response to stimulation more precisely than the linearisation.

Large-scale cortical travelling waves predict localized future cortical signals

Predicting future brain signal is highly sought-after, yet difficult to achieve. To predict the future phase of cortical activity at localized ECoG and MEG recording sites, we exploit its predominant, large-scale, spatiotemporal dynamics. The dynamics are extracted from the brain signal through Fourier analysis and principal components analysis (PCA) only, and cast in a data model that predicts future signal at each site and frequency of interest. The dominant eigenvectors of the PCA that map the large-scale patterns of past cortical phase to future ones take the form of smoothly propagating waves over the entire measurement array. In ECoG data from 3 subjects and MEG data from 20 subjects collected during a self-initiated motor task, mean phase prediction errors were as low as 0.5 radians at local sites, surpassing state-of-the-art methods of within-time-series or event-related models. Prediction accuracy was highest in delta to beta bands, depending on the subject, was more accurate during episodes of high global power, but was not strongly dependent on the time-course of the task. Prediction results did not require past data from the to-be-predicted site. Rather, best accuracy depended on the availability in the model of long wavelength information. The utility of large-scale, low spatial frequency traveling waves in predicting future phase activity at local sites allows estimation of the error introduced by failing to account for irreducible trajectories in the activity dynamics.

Prediction is an important step in scientific progress, often leading to real-world applications. Prediction of future brain activity could lead to improvements in detecting driver and pilot error or real-time brain testing using transcranial magnetic stimulation. Previous studies have either supposed that the ‘noise’ level in the cortex is high, setting the prediction bar rather low; or used localized measurements to predict future activity, with modest success. A long-held but controversial hypothesis is that the cortex is best characterized as a multi-scale dynamic structure, in which the flow of activity at one scale, say, the area responsible for motor control, is inextricably tied to activity at smaller and larger scales, for example within a cortical column and the whole cortex. We test this hypothesis by analyzing large-scale traveling waves of cortical activity. Like waves arriving on a beach, the ongoing wave motion allows better prediction of future activity compared to monitoring the local rise and fall; in the best cases the future wave cycle is predicted with as low as 20° average error angle. The prediction techniques developed for the present research rely on mathematics related to quantifying large-scale weather patterns or analysis of fluid dynamics.

Continuous Phase Estimation for Phase-Locked Neural Stimulation Using an Autoregressive Model for Signal Prediction

Neural oscillations enable communication between brain regions. Closed-loop brain stimulation attempts to modify this activity by stimulation locked to the phase of concurrent neural oscillations. If successful, this may be a major step forward for clinical brain stimulation therapies. The challenge for effective phase-locked systems is accurately calculating the phase of a source oscillation in real time. The basic operations of filtering the source signal to a frequency band of interest and extracting its phase cannot be performed in real time without distortion. We present a method for continuously estimating phase that reduces this distortion by using an autoregressive model to predict the future of a filtered signal before passing it though the Hilbert transform. This method outperforms published approaches on real data and is available as a reusable open-source module. We also examine the challenge of compensating for the filter phase response and outline promising directions of future study.

Pulsed Facilitation of Corticospinal Excitability by the Sensorimotor μ-Alpha Rhythm

Alpha oscillations (8-14 Hz) are assumed to gate information flow in the brain by means of pulsed inhibition; that is, the phasic suppression of cortical excitability and information processing once per alpha cycle, resulting in stronger net suppression for larger alpha amplitudes due to the assumed amplitude asymmetry of the oscillation. While there is evidence for this hypothesis regarding occipital alpha oscillations, it is less clear for the central sensorimotor μ-alpha rhythm. Probing corticospinal excitability via transcranial magnetic stimulation (TMS) of the primary motor cortex and the measurement of motor evoked potentials (MEPs), we have previously demonstrated that corticospinal excitability is modulated by both amplitude and phase of the sensorimotor μ-alpha rhythm. However, the direction of this modulation, its proposed asymmetry, and its underlying mechanisms remained unclear. We therefore used real-time EEG-triggered single- and paired-pulse TMS in healthy humans of both sexes to assess corticospinal excitability and GABA-A-receptor mediated short-latency intracortical inhibition (SICI) at rest during spontaneous high amplitude μ-alpha waves at different phase angles (peaks, troughs, rising and falling flanks) and compared them to periods of low amplitude (desynchronized) μ-alpha. MEP amplitude was facilitated during troughs and rising flanks, but no phasic suppression was observed at any time, nor any modulation of SICI. These results are best compatible with sensorimotor μ-alpha reflecting asymmetric pulsed facilitation but not pulsed inhibition of motor cortical excitability. The asymmetric excitability with respect to rising and falling flanks of the μ-alpha cycle further reveals that voltage differences alone cannot explain the impact of phase.SIGNIFICANCE STATEMENT The pulsed inhibition hypothesis, which assumes that alpha oscillations actively inhibit neuronal processing in a phasic manner, is highly influential and has substantially shaped our understanding of these oscillations. However, some of its basic assumptions, in particular its asymmetry and inhibitory nature, have rarely been tested directly. Here, we explicitly investigated the asymmetry of modulation and its direction for the human sensorimotor μ-alpha rhythm. We found clear evidence of pulsed facilitation, but not inhibition, in the human motor cortex, challenging the generalizability of the pulsed inhibition hypothesis and advising caution when interpreting sensorimotor μ-alpha changes in the sensorimotor system. This study also demonstrates how specific assumptions about the neurophysiological underpinnings of cortical oscillations can be experimentally tested noninvasively in humans.

Brain oscillation-synchronized stimulation of the left dorsolateral prefrontal cortex in depression using real-time EEG-triggered TMS

BACKGROUND

Repetitive transcranial magnetic stimulation (rTMS) of the left dorsolateral prefrontal cortex (DLPFC) is an effective treatment for major depressive disorder (MDD), but response rates are low and effect sizes small. Synchronizing TMS pulses with instantaneous brain oscillations can reduce variability and increase efficacy of TMS-induced plasticity.

OBJECTIVE

To study whether brain oscillation-synchronized rTMS is feasible, safe and has neuromodulatory effects when targeting the DLPFC of patients with MDD.

METHODS

Using real-time EEG-triggered TMS we conducted a pseudo-randomized controlled single-session crossover trial of brain oscillation-synchronized rTMS of left DLPFC in 17 adult patients with antidepressant-resistant MDD. Stimulation conditions in separate sessions were: (1) rTMS triggered at the negative EEG peak of instantaneous alpha oscillations (alpha-synchronized rTMS), (2) a variation of intermittent theta-burst stimulation (modified iTBS), and (3) a random alpha phase control condition.

RESULTS

Triggering TMS at the negative peak of instantaneous alpha oscillations by real-time analysis of the electrode F5 EEG signal was successful in 15 subjects. Two subjects reported mild transient discomfort at the site of stimulation during stimulation; no serious adverse events were reported. Alpha-synchronized rTMS, but not modified iTBS or the random alpha phase control condition, reduced resting-state alpha activity in left DLPFC and increased TMS-induced beta oscillations over frontocentral channels.

CONCLUSIONS

Alpha-synchronized rTMS of left DLPFC is feasible, safe and has specific single-session neuromodulatory effects in patients with antidepressant-resistant MDD. Future studies need to further elucidate the mechanisms, optimize the parameters and investigate the therapeutic potential and efficacy of brain oscillation-synchronized rTMS in MDD.

Causally Investigating Cortical Dynamics and Signal Processing by Targeting Natural System Attractors With Precisely Timed (Electrical) Stimulation

Electrical stimulation is a promising tool for interacting with neuronal dynamics to identify neural mechanisms that underlie cognitive function. Since effects of a single short stimulation pulse typically vary greatly and depend on the current network state, many experimental paradigms have rather resorted to continuous or periodic stimulation in order to establish and maintain a desired effect. However, such an approach explicitly leads to forced and “unnatural” brain activity. Further, continuous stimulation can make it hard to parse the recorded activity and separate neural signal from stimulation artifacts. In this study we propose an alternate strategy: by monitoring a system in realtime, we use the existing preferred states or attractors of the network and apply short and precise pulses in order to switch between those states. When pushed into one of its attractors, one can use the natural tendency of the system to remain in such a state to prolong the effect of a stimulation pulse, opening a larger window of opportunity to observe the consequences on cognitive processing. To elaborate on this idea, we consider flexible information routing in the visual cortex as a prototypical example. When processing a stimulus, neural populations in the visual cortex have been found to engage in synchronized gamma activity. In this context, selective signal routing is achieved by changing the relative phase between oscillatory activity in sending and receiving populations (communication through coherence, CTC). In order to explore how perturbations interact with CTC, we investigate a network of interneuronal gamma (ING) oscillators composed of integrate-and-fire neurons exhibiting similar synchronization and signal routing phenomena. We develop a closed-loop stimulation paradigm based on the phase-response characteristics of the network and demonstrate its ability to establish desired synchronization states. By measuring information content throughout the model, we evaluate the effect of signal contamination caused by the stimulation in relation to the magnitude of the injected pulses and intrinsic noise in the system. Finally, we demonstrate that, up to a critical noise level, precisely timed perturbations can be used to artificially induce the effect of attention by selectively routing visual signals to higher cortical areas.

Cross-Frequency Coupling Based Neuromodulation for Treating Neurological Disorders

Synchronous, rhythmic changes in the membrane polarization of neurons form oscillations in local field potentials. It is hypothesized that high-frequency brain oscillations reflect local cortical information processing, and low-frequency brain oscillations project information flow across larger cortical networks. This provides complex forms of information transmission due to interactions between oscillations at different frequency bands, which can be rendered with cross-frequency coupling (CFC) metrics. Phase-amplitude coupling (PAC) is one of the most common representations of the CFC. PAC reflects the coupling of the phase of oscillations in a specific frequency band to the amplitude of oscillations in another frequency band. In a normal brain, PAC accompanies multi-item working memory in the hippocampus, and changes in PAC have been associated with diseases such as schizophrenia, obsessive-compulsive disorder (OCD), Alzheimer disease (AD), epilepsy, and Parkinson's disease (PD). The purpose of this article is to explore CFC across the central nervous system and demonstrate its correlation to neurological disorders. Results from previously published studies are reviewed to explore the significant role of CFC in large neuronal network communication and its abnormal behavior in neurological disease. Specifically, the association of effective treatment in PD such as dopaminergic medication and deep brain stimulation with PAC changes is described. Lastly, CFC analysis of the electrocorticographic (ECoG) signals recorded from the motor cortex of a Parkinson's disease patient and the parahippocampal gyrus of an epilepsy patient are demonstrated. This information taken together illuminates possible roles of CFC in the nervous system and its potential as a therapeutic target in disease states. This will require new neural interface technologies such as phase-dependent stimulation triggered by PAC changes, for the accurate recording, monitoring, and modulation of the CFC signal.

Phase Synchronicity of μ-Rhythm Determines Efficacy of Interhemispheric Communication Between Human Motor Cortices

The theory of communication through coherence predicts that effective connectivity between nodes in a distributed oscillating neuronal network depends on their instantaneous excitability state and phase synchronicity (Fries, 2005). Here, we tested this prediction by using state-dependent millisecond-resolved real-time electroencephalography-triggered dual-coil transcranial magnetic stimulation (EEG-TMS) (Zrenner et al., 2018) to target the EEG-negative (high-excitability state) versus EEG-positive peak (low-excitability state) of the sensorimotor μ-rhythm in the left (conditioning) and right (test) motor cortex (M1) of 16 healthy human subjects (9 female, 7 male). Effective connectivity was tested by short-interval interhemispheric inhibition (SIHI); that is, the inhibitory effect of the conditioning TMS pulse given 10-12 ms before the test pulse on the test motor-evoked potential. We compared the four possible combinations of excitability states (negative peak, positive peak) and phase relations (in-phase, out-of-phase) of the μ-rhythm in the conditioning and test M1 and a random phase condition. Strongest SIHI was found when the two M1 were in phase for the high-excitability state (negative peak of the μ-rhythm), whereas the weakest SIHI occurred when they were out of phase and the conditioning M1 was in the low-excitability state (positive peak). Phase synchronicity contributed significantly to SIHI variation, with stronger SIHI in the in-phase than out-of-phase conditions. These findings are in exact accord with the predictions of the theory of communication through coherence. They open a translational route for highly effective modification of brain connections by repetitive stimulation at instants in time when nodes in the network are phase synchronized and excitable.SIGNIFICANCE STATEMENT The theory of communication through coherence predicts that effective connectivity between nodes in distributed oscillating brain networks depends on their instantaneous excitability and phase relation. We tested this hypothesis in healthy human subjects by real-time analysis of brain states by electroencephalography in combination with transcranial magnetic stimulation of left and right motor cortex. We found that short-interval interhemispheric inhibition, a marker of interhemispheric effective connectivity, was maximally expressed when the two motor cortices were in phase for a high-excitability state (the trough of the sensorimotor μ-rhythm). We conclude that findings are consistent with the theory of communication through coherence. They open a translational route to highly effectively modify brain connections by repetitive stimulation at instants in time of phase-synchronized high-excitability states.

A Fast EEG Forecasting Algorithm for Phase-Locked Transcranial Electrical Stimulation of the Human Brain

A growing body of research suggests that non-invasive electrical brain stimulation can more effectively modulate neural activity when phase-locked to the underlying brain rhythms. Transcranial alternating current stimulation (tACS) can potentially stimulate the brain in-phase to its natural oscillations as recorded by electroencephalography (EEG), but matching these oscillations is a challenging problem due to the complex and time-varying nature of the EEG signals. Here we address this challenge by developing and testing a novel approach intended to deliver tACS phase-locked to the activity of the underlying brain region in real-time. This novel approach extracts phase and frequency from a segment of EEG, then forecasts the signal to control the stimulation. A careful tuning of the EEG segment length and prediction horizon is required and has been investigated here for different EEG frequency bands. The algorithm was tested on EEG data from 5 healthy volunteers. Algorithm performance was quantified in terms of phase-locking values across a variety of EEG frequency bands. Phase-locking performance was found to be consistent across individuals and recording locations. With current parameters, the algorithm performs best when tracking oscillations in the alpha band (8–13 Hz), with a phase-locking value of 0.77 ± 0.08. Performance was maximized when the frequency band of interest had a dominant frequency that was stable over time. The algorithm performs faster, and provides better phase-locked stimulation, compared to other recently published algorithms devised for this purpose. The algorithm is suitable for use in future studies of phase-locked tACS in preclinical and clinical applications.

Gamma-Rhythmic Gain Modulation

Cognition requires the dynamic modulation of effective connectivity, i.e., the modulation of the postsynaptic neuronal response to a given input. If postsynaptic neurons are rhythmically active, this might entail rhythmic gain modulation, such that inputs synchronized to phases of high gain benefit from enhanced effective connectivity. We show that visually induced gamma-band activity in awake macaque area V4 rhythmically modulates responses to unpredictable stimulus events. This modulation exceeded a simple additive superposition of a constant response onto ongoing gamma-rhythmic firing, demonstrating the modulation of multiplicative gain. Gamma phases leading to strongest neuronal responses also led to shortest behavioral reaction times, suggesting functional relevance of the effect. Furthermore, we find that constant optogenetic stimulation of anesthetized cat area 21a produces gamma-band activity entailing a similar gain modulation. As the gamma rhythm in area 21a did not spread backward to area 17, this suggests that postsynaptic gamma is sufficient for gain modulation.

High-definition transcranial direct current stimulation induces both acute and persistent changes in broadband cortical synchronization: a simultaneous tDCS-EEG study

The goal of this study was to develop methods for simultaneously acquiring electrophysiological data during high-definition transcranial direct current stimulation (tDCS) using high-resolution electroencephalography (EEG). Previous studies have pointed to the after-effects of tDCS on both motor and cognitive performance, and there appears to be potential for using tDCS in a variety of clinical applications. However, little is known about the real-time effects of tDCS on rhythmic cortical activity in humans due to the technical challenges of simultaneously obtaining electrophysiological data during ongoing stimulation. Furthermore, the mechanisms of action of tDCS in humans are not well understood. We have conducted a simultaneous tDCS-EEG study in a group of healthy human subjects. Significant acute and persistent changes in spontaneous neural activity and event-related synchronization (ERS) were observed during and after the application of high-definition tDCS over the left sensorimotor cortex. Both anodal and cathodal stimulation resulted in acute global changes in broadband cortical activity which were significantly different than the changes observed in response to sham stimulation. For the group of eight subjects studied, broadband individual changes in spontaneous activity during stimulation were apparent both locally and globally. In addition, we found that high-definition tDCS of the left sensorimotor cortex can induce significant ipsilateral and contralateral changes in event-related desynchronization and ERS during motor imagination following the end of the stimulation period. Overall, our results demonstrate the feasibility of acquiring high-resolution EEG during high-definition tDCS and provide evidence that tDCS in humans directly modulates rhythmic cortical synchronization during and after its administration.