Waveform changes with the evolution of beta bursts in the human subthalamic nucleus

Abnormal beta bursts of depression in the orbitofrontal cortex and its relationship with clinical symptoms

BACKGROUND

Recent researches have reported that frequency-specific patterns of neural activity contain not only rhythmically sustained oscillations but also transient-bursts of isolated events. The aim of this study was to investigated the correlation between beta burst and depression in order to explore depressive disease and the neurological underpinnings of disease-related symptoms.

METHODS

We collected resting-state MEG recordings from 30 depressive patients and a matched 40 healthy controls. A Hidden Markov Model (HMM) was applied on source-space time courses for 78 cortical regions of the AAL atlas and the temporal characteristics of beta burst from the matched HMM states were captured. Group differences were evaluated on these beta burst characteristics after permutation tests and, for the depressive group, associations between burst characteristics and clinical symptom severity were determined using Spearman correlation coefficients.

RESULTS

At a threshold of p=0.05corrected, burst characteristics revealed significant differences between depression patients and controls at the group level, including increased burst amplitude in frontal lobe, decreased burst duration in occipital regions, increased burst rate and decreased burst interval time in some brain regions. Furthermore, burst amplitude in the orbitofrontal cortex (OFC) was positively related to the severity of sleep disturbance and burst rate in the OFC was negatively related to the severity of anxiety in depression patients.

CONCLUSIONS

The findings highlight OFC may be a targeted area responsible for the anxiety and sleep disturbance symptom by abnormal beta burst in depressive patients and beta burst characteristics of OFC might serve as a neuro-marker for the depression.

Beta: bursts of cognition

Beta oscillations are linked to the control of goal-directed processing of sensory information and the timing of motor output. Recent evidence demonstrates they are not sustained but organized into intermittent high-power bursts mediating timely functional inhibition. This implies there is a considerable moment-to-moment variation in the neural dynamics supporting cognition. Beta bursts thus offer new opportunities for studying how sensory inputs are selectively processed, reshaped by inhibitory cognitive operations and ultimately result in motor actions. Recent method advances reveal diversity in beta bursts that provide deeper insights into their function and the underlying neural circuit activity motifs. We propose that brain-wide, spatiotemporal patterns of beta bursting reflect various cognitive operations and that their dynamics reveal nonlinear aspects of cortical processing.

Quantifying local field potential dynamics with amplitude and frequency stability between ON and OFF medication and stimulation in Parkinson's disease

Neural oscillations are critical to understanding the synchronisation of neural activities and their relevance to neurological disorders. For instance, the amplitude of beta oscillations in the subthalamic nucleus has gained extensive attention, as it has been found to correlate with medication status and the therapeutic effects of continuous deep brain stimulation in people with Parkinson's disease. However, the frequency stability of subthalamic nucleus beta oscillations, which has been suggested to be associated with dopaminergic information in brain states, has not been well explored. Moreover, the administration of medicine can have inverse effects on changes in frequency and amplitude. In this study, we proposed a method based on the stationary wavelet transform to quantify the amplitude and frequency stability of subthalamic nucleus beta oscillations and evaluated the method using simulation and real data for Parkinson's disease patients. The results suggest that the amplitude and frequency stability quantification has enhanced sensitivity in distinguishing pathological conditions in Parkinson's disease patients. Our quantification shows the benefit of combining frequency stability information with amplitude and provides a new potential feedback signal for adaptive deep brain stimulation.

Quantifying Sharpness and Nonlinearity in Neonatal Seizure Dynamics

The integration of multiple electrophysiological biomarkers is crucial for monitoring neonatal seizure dynamics. The present study aimed to characterize the temporal dynamics of neonatal seizures by analyzing intrinsic waveforms of epileptic electroencephalogram (EEG) signals. We proposed a complementary set of methods considering envelope power, focal sharpness changes, and nonlinear patterns of EEG signals of 79 neonates with seizures. Features derived from EEG signals were used as input to the machine learning classifier. All three characteristics were significantly elevated during seizure events, as agreed upon by all viewers (P < 0.0001). Envelope power was elevated in the entire seizure period, and the degree of nonlinearity rose at the termination of a seizure event. Epileptic sharpness effectively characterizes an entire seizure event, complementing the role of envelope power in identifying its onset. However, the degree of nonlinearity showed superior discriminability for the termination of a seizure event. The proposed computational methods for intrinsic sharp or nonlinear EEG patterns evolving during neonatal seizure could share some features with envelope power. Current findings may be helpful in developing strategies to improve neonatal seizure monitoring.

Diverse beta burst waveform motifs characterize movement-related cortical dynamics

Classical analyses of induced, frequency-specific neural activity typically average band-limited power over trials. More recently, it has become widely appreciated that in individual trials, beta band activity occurs as transient bursts rather than amplitude-modulated oscillations. Most studies of beta bursts treat them as unitary, and having a stereotyped waveform. However, we show there is a wide diversity of burst shapes. Using a biophysical model of burst generation, we demonstrate that waveform variability is predicted by variability in the synaptic drives that generate beta bursts. We then use a novel, adaptive burst detection algorithm to identify bursts from human MEG sensor data recorded during a joystick-based reaching task, and apply principal component analysis to burst waveforms to define a set of dimensions, or motifs, that best explain waveform variance. Finally, we show that bursts with a particular range of waveform motifs, ones not fully accounted for by the biophysical model, differentially contribute to movement-related beta dynamics. Sensorimotor beta bursts are therefore not homogeneous events and likely reflect distinct computational processes.

Clinical neurophysiological interrogation of motor slowing: A critical step towards tuning adaptive deep brain stimulation

OBJECTIVE

Subthalamic nucleus (STN) beta activity (13-30 Hz) is the most accepted biomarker for adaptive deep brain stimulation (aDBS) for Parkinson's disease (PD). We hypothesize that different frequencies within the beta range may exhibit distinct temporal dynamics and, as a consequence, different relationships to motor slowing and adaptive stimulation patterns. We aim to highlight the need for an objective method to determine the aDBS feedback signal.

METHODS

STN LFPs were recorded in 15 PD patients at rest and while performing a cued motor task. The impact of beta bursts on motor performance was assessed for different beta candidate frequencies: the individual frequency strongest associated with motor slowing, the individual beta peak frequency, the frequency most modulated by movement execution, as well as the entire-, low- and high beta band. How these candidate frequencies differed in their bursting dynamics and theoretical aDBS stimulation patterns was further investigated.

RESULTS

The individual motor slowing frequency often differs from the individual beta peak or beta-related movement-modulation frequency. Minimal deviations from a selected target frequency as feedback signal for aDBS leads to a substantial drop in the burst overlapping and in the alignment of the theoretical onset of stimulation triggers (to ∼ 75% for 1 Hz, to ∼ 40% for 3 Hz deviation).

CONCLUSIONS

Clinical-temporal dynamics within the beta frequency range are highly diverse and deviating from a reference biomarker frequency can result in altered adaptive stimulation patterns.

SIGNIFICANCE

A clinical-neurophysiological interrogation could be helpful to determine the patient-specific feedback signal for aDBS.

The potential of convolutional neural networks for identifying neural states based on electrophysiological signals: experiments on synthetic and real patient data

Processing incoming neural oscillatory signals in real-time and decoding from them relevant behavioral or pathological states is often required for adaptive Deep Brain Stimulation (aDBS) and other brain-computer interface (BCI) applications. Most current approaches rely on first extracting a set of predefined features, such as the power in canonical frequency bands or various time-domain features, and then training machine learning systems that use those predefined features as inputs and infer what the underlying brain state is at each given time point. However, whether this algorithmic approach is best suited to extract all available information contained within the neural waveforms remains an open question. Here, we aim to explore different algorithmic approaches in terms of their potential to yield improvements in decoding performance based on neural activity such as measured through local field potentials (LFPs) recordings or electroencephalography (EEG). In particular, we aim to explore the potential of end-to-end convolutional neural networks, and compare this approach with other machine learning methods that are based on extracting predefined feature sets. To this end, we implement and train a number of machine learning models, based either on manually constructed features or, in the case of deep learning-based models, on features directly learnt from the data. We benchmark these models on the task of identifying neural states using simulated data, which incorporates waveform features previously linked to physiological and pathological functions. We then assess the performance of these models in decoding movements based on local field potentials recorded from the motor thalamus of patients with essential tremor. Our findings, derived from both simulated and real patient data, suggest that end-to-end deep learning-based methods may surpass feature-based approaches, particularly when the relevant patterns within the waveform data are either unknown, difficult to quantify, or when there may be, from the point of view of the predefined feature extraction pipeline, unidentified features that could contribute to decoding performance. The methodologies proposed in this study might hold potential for application in adaptive deep brain stimulation (aDBS) and other brain-computer interface systems.

Cross-Frequency Coupling and Intelligent Neuromodulation

Cross-frequency coupling (CFC) reflects (nonlinear) interactions between signals of different frequencies. Evidence from both patient and healthy participant studies suggests that CFC plays an essential role in neuronal computation, interregional interaction, and disease pathophysiology. The present review discusses methodological advances and challenges in the computation of CFC with particular emphasis on potential solutions to spurious coupling, inferring intrinsic rhythms in a targeted frequency band, and causal interferences. We specifically focus on the literature exploring CFC in the context of cognition/memory tasks, sleep, and neurological disorders, such as Alzheimer's disease, epilepsy, and Parkinson's disease. Furthermore, we highlight the implication of CFC in the context and for the optimization of invasive and noninvasive neuromodulation and rehabilitation. Mainly, CFC could support advancing the understanding of the neurophysiology of cognition and motor control, serve as a biomarker for disease symptoms, and leverage the optimization of therapeutic interventions, e.g., closed-loop brain stimulation. Despite the evident advantages of CFC as an investigative and translational tool in neuroscience, further methodological improvements are required to facilitate practical and correct use in cyborg and bionic systems in the field.

A Novel Framework in Quantifying Oscillatory Coupling to Gait Disturbance in Parkinson's Disease

Phase-amplitude coupling (PAC) based on the uniform phase empirical mode decomposition (UPEMD) is proposed to improve the accuracy of PAC assessment. The framework is applied to investigate the mechanism and improvement measure of gait disturbance for Parkinson's disease (PD). Hβ modulation is suppressed at the time of contralateral heel strikes and rebounds when the contralateral foot rests on the ground and the ipsilateral foot is raised. Prominent PACs exist between δ and Lβ/Hβ activities. Auditory cue improves the gait; meanwhile, it enhances the Hβ modulation, and suppresses the δ-Lβ/Hβ PACs, which may rebound toward the before-cue stage afterward. Our findings suggest the proposed UPEMD-PAC is a useful framework in quantifying PAC with pre-determined frequencies, whereas the δ-Lβ/Hβ PACs in the subthalamic nucleus serve as potential biomarkers for gait disturbance in PD. Clinical Relevance- This manifests the efficacy of auditory cues on gait disturbance. The proposed framework may be useful in diagnosing the severity of motor impairment.

Subthalamic low beta bursts differ in Parkinson's disease phenotypes

OBJECTIVE

Parkinson's disease (PD) patients may be categorized into tremor-dominant (TD) and postural-instability and gait disorder (PIGD) motor phenotypes, but the dynamical aspects of subthalamic nucleus local field potentials (STN-LFP) and the neural correlates of this phenotypical classification remain unclear.

METHODS

35 STN-LFP (20 PIGD and 15 TD) were investigated through continuous wavelet transform and machine-learning-based methods. The beta oscillation - the main band associated with motor impairment in PD - dynamics was characterized through beta burst parameters across phenotypes and burst intervals under specific proposed criteria for optimal burst threshold definition.

RESULTS

Low-frequency (13-22 Hz) beta burst probability was the best predictor for PD phenotypes (75% accuracy). PIGD patients presented higher average burst duration (p = 0.018), while TD patients exhibited higher burst probability (p = 0.014). Categorization into shorter and longer than 400 ms bursts led to significant interaction between burst length categories and the phenotypes (p < 0.050) as revealed by mixed-effects models. Long burst durations and short bursts probability positively correlated, respectively, with rigidity-bradykinesia (p = 0.029) and tremor (p = 0.038) scores.

CONCLUSIONS

Subthalamic low-frequency beta bursts differed between TD and PIGD phenotypes and correlated with motor symptoms.

SIGNIFICANCE

These findings improve the PD phenotypes' electrophysiological characterization and may define new criteria for adaptive deep brain stimulation.

Balance between competing spectral states in subthalamic nucleus is linked to motor impairment in Parkinson's disease

Exaggerated local field potential bursts of activity at frequencies in the low beta band are a well-established phenomenon in the subthalamic nucleus of patients with Parkinson's disease. However, such activity is only moderately correlated with motor impairment. Here we test the hypothesis that beta bursts are just one of several dynamic states in the subthalamic nucleus local field potential in Parkinson's disease, and that together these different states predict motor impairment with high fidelity. Local field potentials were recorded in 32 patients (64 hemispheres) undergoing deep brain stimulation surgery targeting the subthalamic nucleus. Recordings were performed following overnight withdrawal of anti-parkinsonian medication, and after administration of levodopa. Local field potentials were analysed using hidden Markov modelling to identify transient spectral states with frequencies under 40 Hz. Findings in the low beta frequency band were similar to those previously reported; levodopa reduced occurrence rate and duration of low beta states, and the greater the reductions, the greater the improvement in motor impairment. However, additional local field potential states were distinguished in the theta, alpha and high beta bands, and these behaved in an opposite manner. They were increased in occurrence rate and duration by levodopa, and the greater the increases, the greater the improvement in motor impairment. In addition, levodopa favoured the transition of low beta states to other spectral states. When all local field potential states and corresponding features were considered in a multivariate model it was possible to predict 50% of the variance in patients' hemibody impairment OFF medication, and in the change in hemibody impairment following levodopa. This only improved slightly if signal amplitude or gamma band features were also included in the multivariate model. In addition, it compares with a prediction of only 16% of the variance when using beta bursts alone. We conclude that multiple spectral states in the subthalamic nucleus local field potential have a bearing on motor impairment, and that levodopa-induced shifts in the balance between these states can predict clinical change with high fidelity. This is important in suggesting that some states might be upregulated to improve parkinsonism and in suggesting how local field potential feedback can be made more informative in closed-loop deep brain stimulation systems.

An Impending Paradigm Shift in Motor Imagery Based Brain-Computer Interfaces

The development of reliable assistive devices for patients that suffer from motor impairments following central nervous system lesions remains a major challenge in the field of non-invasive Brain-Computer Interfaces (BCIs). These approaches are predominated by electroencephalography and rely on advanced signal processing and machine learning methods to extract neural correlates of motor activity. However, despite tremendous and still ongoing efforts, their value as effective clinical tools remains limited. We advocate that a rather overlooked research avenue lies in efforts to question neurophysiological markers traditionally targeted in non-invasive motor BCIs. We propose an alternative approach grounded by recent fundamental advances in non-invasive neurophysiology, specifically subject-specific feature extraction of sensorimotor bursts of activity recorded via (possibly magnetoencephalography-optimized) electroencephalography. This path holds promise in overcoming a significant proportion of existing limitations, and could foster the wider adoption of online BCIs in rehabilitation protocols.

Within-cycle instantaneous frequency profiles report oscillatory waveform dynamics

The nonsinusoidal waveform is emerging as an important feature of neuronal oscillations. However, the role of single-cycle shape dynamics in rapidly unfolding brain activity remains unclear. Here, we develop an analytical framework that isolates oscillatory signals from time series using masked empirical mode decomposition to quantify dynamical changes in the shape of individual cycles (along with amplitude, frequency, and phase) with instantaneous frequency. We show how phase-alignment, a process of projecting cycles into a regularly sampled phase grid space, makes it possible to compare cycles of different durations and shapes. "Normalized shapes" can then be constructed with high temporal detail while accounting for differences in both duration and amplitude. We find that the instantaneous frequency tracks nonsinusoidal shapes in both simulated and real data. Notably, in local field potential recordings of mouse hippocampal CA1, we find that theta oscillations have a stereotyped slow-descending slope in the cycle-wise average yet exhibit high variability on a cycle-by-cycle basis. We show how principal component analysis allows identification of motifs of theta cycle waveform that have distinct associations to cycle amplitude, cycle duration, and animal movement speed. By allowing investigation into oscillation shape at high temporal resolution, this analytical framework will open new lines of inquiry into how neuronal oscillations support moment-by-moment information processing and integration in brain networks.NEW & NOTEWORTHY We propose a novel analysis approach quantifying nonsinusoidal waveform shape. The approach isolates oscillations with empirical mode decomposition before waveform shape is quantified using phase-aligned instantaneous frequency. This characterizes the full shape profile of individual cycles while accounting for between-cycle differences in duration, amplitude, and timing. We validated in simulations before applying to identify a range of data-driven nonsinusoidal shape motifs in hippocampal theta oscillations.