Average beta burst duration profiles provide a signature of dynamical changes between the ON and OFF medication states in Parkinson's disease

Evoked resonant neural activity long-term dynamics can be reproduced by a computational model with vesicle depletion

Subthalamic deep brain stimulation (DBS) robustly generates high-frequency oscillations known as evoked resonant neural activity (ERNA). Recently the importance of ERNA has been demonstrated through its ability to predict the optimal DBS contact in the subthalamic nucleus in patients with Parkinson's disease. However, the underlying mechanisms of ERNA are not well understood, and previous modelling efforts have not managed to reproduce the wealth of published data describing the dynamics of ERNA. Here, we aim to present a minimal model capable of reproducing the characteristics of the slow ERNA dynamics published to date. We make biophysically-motivated modifications to the Kuramoto model and fit its parameters to the slow dynamics of ERNA obtained from data. Our results demonstrate that it is possible to reproduce the slow dynamics of ERNA (over hundreds of seconds) with a single neuronal population, and, crucially, with vesicle depletion as one of the key mechanisms behind the ERNA frequency decay in our model. We further validate the proposed model against experimental data from Parkinson's disease patients, where it captures the variations in ERNA frequency and amplitude in response to variable stimulation frequency, amplitude, and to stimulation pulse bursting. We provide a series of predictions from the model that could be the subject of future studies for further validation.

How to design optimal brain stimulation to modulate phase-amplitude coupling?

Objective.Phase-amplitude coupling (PAC), the coupling of the amplitude of a faster brain rhythm to the phase of a slower brain rhythm, plays a significant role in brain activity and has been implicated in various neurological disorders. For example, in Parkinson's disease, PAC between the beta (13-30 Hz) and gamma (30-100 Hz) rhythms in the motor cortex is exaggerated, while in Alzheimer's disease, PAC between the theta (4-8 Hz) and gamma rhythms is diminished. Modulating PAC (i.e. reducing or enhancing PAC) using brain stimulation could therefore open new therapeutic avenues. However, while it has been previously reported that phase-locked stimulation can increase PAC, it is unclear what the optimal stimulation strategy to modulate PAC might be. Here, we provide a theoretical framework to narrow down the experimental optimisation of stimulation aimed at modulating PAC, which would otherwise rely on trial and error.Approach.We make analytical predictions using a Stuart-Landau model, and confirm these predictions in a more realistic model of coupled neural populations.Main results.Our framework specifies the critical Fourier coefficients of the stimulation waveform which should be tuned to optimally modulate PAC. Depending on the characteristics of the amplitude response curve of the fast population, these components may include the slow frequency, the fast frequency, combinations of these, as well as their harmonics. We also show that the optimal balance of energy between these Fourier components depends on the relative strength of the endogenous slow and fast rhythms, and that the alignment of fast components with the fast rhythm should change throughout the slow cycle. Furthermore, we identify the conditions requiring to phase-lock stimulation to the fast and/or slow rhythms.Significance.Together, our theoretical framework lays the foundation for guiding the development of innovative and more effective brain stimulation aimed at modulating PAC for therapeutic benefit.

Concurrent inflammation-related brain reorganization in multiple sclerosis and depression

BACKGROUND

Neuroinflammation affects brain tissue integrity in multiple sclerosis (MS) and may have a role in major depressive disorder (MDD). Whether advanced magnetic resonance imaging characteristics of the gray-to-white matter border serve as proxy of neuroinflammatory activity in MDD and MS remain unknown.

METHODS

We included 684 participants (132 MDD patients with recurrent depressive episodes (RDE), 70 MDD patients with a single depressive episode (SDE), 222 MS patients without depressive symptoms (nMS), 58 MS patients with depressive symptoms (dMS), and 202 healthy controls (HC)). 3 T-T1w MRI-derived gray-to-white matter contrast (GWc) was used to reconstruct and characterize connectivity alterations of GWc-covariance networks by means of modularity, clustering coefficient, and degree. A cross-validated support vector machine was used to test the ability of GWc to stratify groups according to their depression symptoms, measured with BDI, at the single-subject level in MS and MDD independently.

FINDINGS

MS and MDD patients showed increased modularity (ANOVA partial-η2 = 0.3) and clustering (partial-η2 = 0.1) compared to HC. In the subgroups, a linear trend analysis attested a gradient of modularity increases in the form: HC, dMS, nMS, SDE, and RDE (ANOVA partial-η2 = 0.28, p < 0.001) while this trend was less evident for clustering coefficient. Reduced morphological integrity (GWc) was seen in patients with increased depressive symptoms (partial-η2 = 0.42, P < 0.001) and was associated with depression scores across patient groups (r = -0.2, P < 0.001). Depressive symptoms in MS were robustly classified (88 %).

CONCLUSIONS

Similar structural network alterations in MDD and MS exist, suggesting possible common inflammatory events like demyelination, neuroinflammation that are caught by GWc analyses. These alterations may vary depending on the severity of symptoms and in the case of MS may elucidate the occurrence of comorbid depression.

Subthalamic stimulation modulates context-dependent effects of beta bursts during fine motor control

Increasing evidence suggests a considerable role of pre-movement beta bursts for motor control and its impairment in Parkinson's disease. However, whether beta bursts occur during precise and prolonged movements and if they affect fine motor control remains unclear. To investigate the role of within-movement beta bursts for fine motor control, we here combine invasive electrophysiological recordings and clinical deep brain stimulation in the subthalamic nucleus in 19 patients with Parkinson's disease performing a context-varying task that comprised template-guided and free spiral drawing. We determined beta bursts in narrow frequency bands around patient-specific peaks and assessed burst amplitude, duration, and their immediate impact on drawing speed. We reveal that beta bursts occur during the execution of drawing movements with reduced duration and amplitude in comparison to rest. Exclusively when drawing freely, they parallel reductions in acceleration. Deep brain stimulation increases the acceleration around beta bursts in addition to a general increase in drawing velocity and improvements of clinical function. These results provide evidence for a diverse and task-specific role of subthalamic beta bursts for fine motor control in Parkinson's disease; suggesting that pathological beta bursts act in a context dependent manner, which can be targeted by clinical deep brain stimulation.

Oscillatory and non-oscillatory features of the magnetoencephalic sensorimotor rhythm in Parkinson's disease

Parkinson's disease (PD) is associated with changes in neural activity in the sensorimotor alpha and beta bands. Using magnetoencephalography (MEG), we investigated the role of spontaneous neuronal activity within the somatosensory cortex in a large cohort of early- to mid-stage PD patients (N = 78) on Parkinsonian medication and age- and sex-matched healthy controls (N = 60) using source reconstructed resting-state MEG. We quantified features of the time series data in terms of oscillatory alpha power and central alpha frequency, beta power and central beta frequency, and 1/f broadband characteristics using power spectral density. Furthermore, we characterised transient oscillatory burst events in the mu-beta band time-domain signals. We examined the relationship between these signal features and the patients' disease state, symptom severity, age, sex, and cortical thickness. PD patients and healthy controls differed on PSD broadband characteristics, with PD patients showing a steeper 1/f exponential slope and higher 1/f offset. PD patients further showed a steeper age-related decrease in the burst rate. Out of all the signal features of the sensorimotor activity, the burst rate was associated with increased severity of bradykinesia, whereas the burst duration was associated with axial symptoms. Our study shows that general non-oscillatory features (broadband 1/f exponent and offset) of the sensorimotor signals are related to disease state and oscillatory burst rate scales with symptom severity in PD.

Abnormally low sensorimotor α band nonlinearity serves as an effective EEG biomarker of Parkinson's disease

Biomarkers obtained from the neurophysiological signals of patients with Parkinson's disease (PD) have objective value in assessing their motor condition for effective diagnosis, monitoring, and clinical intervention. Prominent cortical biomarkers of PD have typically been derived from various β band wave features. This study approached the topic from an alternative perspective and attempted to estimate a recently suggested measure representing α band nonlinear autocorrelative memory from a publicly available EEG dataset that involves 15 patients with earlier-stage PD (dopaminergic medication OFF and ON states) and 16 age-matched healthy controls. The cortical nonlinearity was elevated for the PD ON state compared with the OFF state for bilateral sensorimotor channels C3 and C4 (n = 26; P = 0.003). A similar statistical difference was also identified between PD OFF state and healthy subjects (n = 26; P = 0.049). Analysis over all channels revealed that the α band nonlinearity induced upon medication was constrained to sensorimotor regions. The α nonlinearity measure was compared with a well-accepted cortical biomarker of β-γ phase-amplitude coupling (PAC). They were in moderate negative correlation (r = -0.412; P = 0.036) for only healthy subjects, but not for the patients. The nonlinearity measure was found to be insusceptible to the nonstationary variations within the particular data. Our study provides further evidence that the α band nonlinearity measure can serve as a promising cortical biomarker of PD. The suggested measure can be estimated from a noninvasive low-resolution single scalp EEG channel of patients with relatively early-stage PD, who did not yet need to undergo deep brain stimulation operation.NEW & NOTEWORTHY This study suggests a nonlinearity measure that differentiates Parkinson's disease (PD) dopamine OFF-state scalp EEG data from those of dopamine ON-state patients and healthy subjects. Unlike typical PD cortical biomarkers based on β band activity, this metric shows elevation upon dopaminergic medication in the α band. We provide evidence supporting its potential as an early-stage promising PD biomarker that can be estimated from noninvasive EEG recordings with low resolution and SNR.

Oscillatory vs. non-oscillatory subthalamic beta activity in Parkinson's disease

Parkinson's disease is characterized by exaggerated beta activity (13-35 Hz) in cortico-basal ganglia motor loops. Beta activity includes both periodic fluctuations (i.e. oscillatory activity) and aperiodic fluctuations reflecting spiking activity and excitation/inhibition balance (i.e. non-oscillatory activity). However, the relative contribution, dopamine dependency and clinical correlations of oscillatory vs. non-oscillatory beta activity remain unclear. We recorded, modelled and analysed subthalamic local field potentials in parkinsonian patients at rest while off or on medication. Autoregressive modelling with additive 1/f noise clarified the relationships between measures of beta activity in the time domain (i.e. amplitude and duration of beta bursts) or in the frequency domain (i.e. power and sharpness of the spectral peak) and oscillatory vs. non-oscillatory activity: burst duration and spectral sharpness are specifically sensitive to oscillatory activity, whereas burst amplitude and spectral power are ambiguously sensitive to both oscillatory and non-oscillatory activity. Our experimental data confirmed the model predictions and assumptions. We subsequently analysed the effect of levodopa, obtaining strong-to-extreme Bayesian evidence that oscillatory beta activity is reduced in patients on vs. off medication, with moderate evidence for absence of modulation of the non-oscillatory component. Finally, specifically the oscillatory component of beta activity correlated with the rate of motor progression of the disease. Methodologically, these results provide an integrative understanding of beta-based biomarkers relevant for adaptive deep brain stimulation. Biologically, they suggest that primarily the oscillatory component of subthalamic beta activity is dopamine dependent and may play a role not only in the pathophysiology but also in the progression of Parkinson's disease. KEY POINTS: Beta activity in Parkinson's disease includes both true periodic fluctuations (i.e. oscillatory activity) and aperiodic fluctuations reflecting spiking activity and synaptic balance (i.e. non-oscillatory activity). The relative contribution, dopamine dependency and clinical correlations of oscillatory vs. non-oscillatory beta activity remain unclear. Burst duration and spectral sharpness are specifically sensitive to oscillatory activity, while burst amplitude and spectral power are ambiguously sensitive to both oscillatory and non-oscillatory activity. Only the oscillatory component of subthalamic beta activity is dopamine-dependent. Stronger beta oscillatory activity correlates with faster motor progression of the disease.

Bursting with Potential: How Sensorimotor Beta Bursts Develop from Infancy to Adulthood

Beta activity is thought to play a critical role in sensorimotor processes. However, little is known about how activity in this frequency band develops. Here, we investigated the developmental trajectory of sensorimotor beta activity from infancy to adulthood. We recorded EEG from 9-month-old, 12-month-old, and adult humans (male and female) while they observed and executed grasping movements. We analyzed "beta burst" activity using a novel method that combines time-frequency decomposition and principal component analysis. We then examined the changes in burst rate and waveform motifs along the selected principal components. Our results reveal systematic changes in beta activity during action execution across development. We found a decrease in beta burst rate during movement execution in all age groups, with the greatest decrease observed in adults. Additionally, we identified three principal components that defined waveform motifs that systematically changed throughout the trial. We found that bursts with waveform shapes closer to the median waveform were not rate-modulated, whereas those with waveform shapes further from the median were differentially rate-modulated. Interestingly, the decrease in the rate of certain burst motifs occurred earlier during movement and was more lateralized in adults than in infants, suggesting that the rate modulation of specific types of beta bursts becomes increasingly refined with age.SIGNIFICANCE STATEMENT We demonstrate that, like in adults, sensorimotor beta activity in infants during reaching and grasping movements occurs in bursts, not oscillations like thought traditionally. Furthermore, different beta waveform shapes were differentially modulated with age, including more lateralization in adults. Aberrant beta activity characterizes various developmental disorders and motor difficulties linked to early brain injury, so looking at burst waveform shape could provide more sensitivity for early identification and treatment of affected individuals before any behavioral symptoms emerge. More generally, comparison of beta burst activity in typical versus atypical motor development will also be instrumental in teasing apart the mechanistic functional roles of different types of beta bursts.

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.

Sub-harmonic entrainment of cortical gamma oscillations to deep brain stimulation in Parkinson's disease: Model based predictions and validation in three human subjects

OBJECTIVES

The exact mechanisms of deep brain stimulation (DBS) are still an active area of investigation, in spite of its clinical successes. This is due in part to the lack of understanding of the effects of stimulation on neuronal rhythms. Entrainment of brain oscillations has been hypothesised as a potential mechanism of neuromodulation. A better understanding of entrainment might further inform existing methods of continuous DBS, and help refine algorithms for adaptive methods. The purpose of this study is to develop and test a theoretical framework to predict entrainment of cortical rhythms to DBS across a wide range of stimulation parameters.

MATERIALS AND METHODS

We fit a model of interacting neural populations to selected features characterising PD patients' off-stimulation finely-tuned gamma rhythm recorded through electrocorticography. Using the fitted models, we predict basal ganglia DBS parameters that would result in 1:2 entrainment, a special case of sub-harmonic entrainment observed in patients and predicted by theory.

RESULTS

We show that the neural circuit models fitted to patient data exhibit 1:2 entrainment when stimulation is provided across a range of stimulation parameters. Furthermore, we verify key features of the region of 1:2 entrainment in the stimulation frequency/amplitude space with follow-up recordings from the same patients, such as the loss of 1:2 entrainment above certain stimulation amplitudes.

CONCLUSION

Our results reveal that continuous, constant frequency DBS in patients may lead to nonlinear patterns of neuronal entrainment across stimulation parameters, and that these responses can be predicted by modelling. Should entrainment prove to be an important mechanism of therapeutic stimulation, our modelling framework may reduce the parameter space that clinicians must consider when programming devices for optimal benefit.

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.

Translating Pathological Brain Activity Primers in Parkinson's Disease Research

Translational experimental approaches that help us better trace Parkinson's disease (PD) pathophysiological mechanisms leading to new therapeutic targets are urgently needed. In this article, we review recent experimental and clinical studies addressing abnormal neuronal activity and pathological network oscillations, as well as their underlying mechanisms and modulation. Our aim is to enhance our knowledge about the progression of Parkinson's disease pathology and the timing of its symptom's manifestation. Here, we present mechanistic insights relevant for the generation of aberrant oscillatory activity within the cortico-basal ganglia circuits. We summarize recent achievements extrapolated from available PD animal models, discuss their advantages and limitations, debate on their differential applicability, and suggest approaches for transferring knowledge on disease pathology into future research and clinical applications.

Average power and burst analysis revealed complementary information on drug-related changes of motor performance in Parkinson's disease

In patients with Parkinson's disease (PD), suppression of beta and increase in gamma oscillations in the subthalamic nucleus (STN) have been associated with both levodopa treatment and motor functions. Recent results suggest that modulation of the temporal dynamics of theses oscillations (bursting activity) might contain more information about pathological states and behaviour than their average power. Here we directly compared the information provided by power and burst analyses about the drug-related changes in STN activities and their impact on motor performance within PD patients. STN local field potential (LFP) signals were recorded from externalized patients performing self-paced movements ON and OFF levodopa. When normalised across medication states, both power and burst analyses showed an increase in low-beta oscillations in the dopamine-depleted state during rest. When normalised within-medication state, both analyses revealed that levodopa increased movement-related modulation in the alpha and low-gamma bands, with higher gamma activity around movement predicting faster reaches. Finally, burst analyses helped to reveal opposite drug-related changes in low- and high-beta frequency bands, and identified additional within-patient relationships between high-beta bursting and movement performance. Our findings suggest that although power and burst analyses share a lot in common they also provide complementary information on how STN-LFP activity is associated with motor performance, and how levodopa treatment may modify these relationships in a way that helps explain drug-related changes in motor performance. Different ways of normalisation in the power analysis can reveal different information. Similarly, the burst analysis is sensitive to how the threshold is defined - either for separate medication conditions separately, or across pooled conditions. In addition, the burst interpretation has far-reaching implications about the nature of neural oscillations - whether the oscillations happen as isolated burst-events or are they sustained phenomena with dynamic amplitude variations? This can be different for different frequency bands, and different for different medication states even for the same frequency band.

Oscillatory beta dynamics inform biomarker-driven treatment optimization for Parkinson's disease

Parkinson's disease (PD) is a neurodegenerative disorder characterized by loss of dopaminergic neurons and dysregulation of the basal ganglia. Cardinal motor symptoms include bradykinesia, rigidity, and tremor. Deep brain stimulation (DBS) of select subcortical nuclei is standard of care for medication-refractory PD. Conventional open-loop DBS delivers continuous stimulation with fixed parameters that do not account for a patient's dynamic activity state or medication cycle. In comparison, closed-loop DBS, or adaptive DBS (aDBS), adjusts stimulation based on biomarker feedback that correlates with clinical state. Recent work has identified several neurophysiological biomarkers in local field potential recordings from PD patients, the most promising of which are 1) elevated beta (∼13-30 Hz) power in the subthalamic nucleus (STN), 2) increased beta synchrony throughout basal ganglia-thalamocortical circuits, notably observed as coupling between the STN beta phase and cortical broadband gamma (∼50-200 Hz) amplitude, and 3) prolonged beta bursts in the STN and cortex. In this review, we highlight relevant frequency and time domain features of STN beta measured in PD patients and summarize how spectral beta power, oscillatory beta synchrony, phase-amplitude coupling, and temporal beta bursting inform PD pathology, neurosurgical targeting, and DBS therapy. We then review how STN beta dynamics inform predictive, biomarker-driven aDBS approaches for optimizing PD treatment. We therefore provide clinically useful and actionable insight that can be applied toward aDBS implementation for PD.

Evoked resonant neural activity in subthalamic local field potentials reflects basal ganglia network dynamics

Evoked resonant neural activity (ERNA) is induced by subthalamic deep brain stimulation (DBS) and was recently suggested as a marker of lead placement and contact selection in Parkinson's disease. Yet, its underlying mechanisms and how it is modulated by stimulation parameters are unclear. Here, we recorded local field potentials from 27 Parkinson's disease patients, while leads were externalised to scrutinise the ERNA. First, we show that ERNA in the time series waveform and spectrogram likely represent the same activity, which was contested before. Second, our results show that the ERNA has fast and slow dynamics during stimulation, consistent with the synaptic failure hypothesis. Third, we show that ERNA parameters are modulated by different DBS frequencies, intensities, medication states and stimulation modes (continuous DBS vs. adaptive DBS). These results suggest the ERNA might prove useful as a predictor of the best DBS frequency and lowest effective intensity in addition to contact selection. Changes with levodopa and DBS mode suggest that the ERNA may indicate the state of the cortico-basal ganglia circuit making it a putative biomarker to track clinical state in adaptive DBS.

The aperiodic exponent of subthalamic field potentials reflects excitation/inhibition balance in Parkinsonism

Periodic features of neural time-series data, such as local field potentials (LFPs), are often quantified using power spectra. While the aperiodic exponent of spectra is typically disregarded, it is nevertheless modulated in a physiologically relevant manner and was recently hypothesised to reflect excitation/inhibition (E/I) balance in neuronal populations. Here, we used a cross-species in vivo electrophysiological approach to test the E/I hypothesis in the context of experimental and idiopathic Parkinsonism. We demonstrate in dopamine-depleted rats that aperiodic exponents and power at 30-100 Hz in subthalamic nucleus (STN) LFPs reflect defined changes in basal ganglia network activity; higher aperiodic exponents tally with lower levels of STN neuron firing and a balance tipped towards inhibition. Using STN-LFPs recorded from awake Parkinson's patients, we show that higher exponents accompany dopaminergic medication and deep brain stimulation (DBS) of STN, consistent with untreated Parkinson's manifesting as reduced inhibition and hyperactivity of STN. These results suggest that the aperiodic exponent of STN-LFPs in Parkinsonism reflects E/I balance and might be a candidate biomarker for adaptive DBS.

When do bursts matter in the primary motor cortex? Investigating changes in the intermittencies of beta rhythms associated with movement states

Brain activity exhibits significant temporal structure that is not well captured in the power spectrum. Recently, attention has shifted to characterising the properties of intermittencies in rhythmic neural activity (i.e. bursts), yet the mechanisms that regulate them are unknown. Here, we present evidence from electrocorticography recordings made over the motor cortex to show that the statistics of bursts, such as duration or amplitude, in the beta frequency (14-30 Hz) band, significantly aid the classification of motor states such as rest, movement preparation, execution, and imagery. These features reflect nonlinearities not detectable in the power spectrum, with states increasing in nonlinearity from movement execution to preparation to rest. Further, we show using a computational model of the cortical microcircuit, constrained to account for burst features, that modulations of laminar specific inhibitory interneurons are responsible for the temporal organisation of activity. Finally, we show that the temporal characteristics of spontaneous activity can be used to infer the balance of cortical integration between incoming sensory information and endogenous activity. Critically, we contribute to the understanding of how transient brain rhythms may underwrite cortical processing, which in turn, could inform novel approaches for brain state classification, and modulation with novel brain-computer interfaces.

Novel approaches for quantifying beta synchrony in Parkinson's disease

Despite the clinical and financial burden of Parkinson's disease (PD), there is no standardized, reliable biomarker to diagnose and track PD progression. Instead, PD is primarily assessed using subjective clinical rating scales and patient self-report. Such approaches can be imprecise, hindering diagnosis and disease monitoring. An objective biomarker would be beneficial for clinical care, refining diagnosis, and treatment. Due to widespread electrophysiological abnormalities both within and between brain structures in PD, development of electrophysiologic biomarkers may be feasible. Basal ganglia recordings acquired with neurosurgical approaches have revealed elevated power in the beta frequency range (13-30 Hz) in PD, suggesting that beta power could be a putative PD biomarker. However, there are limitations to the use of beta power as a biomarker. Recent advances in analytic approaches have led to novel methods to quantify oscillatory synchrony in the beta frequency range. Here we describe some of these novel approaches in the context of PD and explore how they may serve as electrophysiological biomarkers. These novel signatures include (1) interactions between beta phase and broadband (> 50 Hz, "gamma") amplitude (i.e., phase amplitude coupling, PAC), (2) asymmetries in waveform shape, (3) beta coherence, and (4) beta "bursts." Development of a robust, reliable, and readily accessible electrophysiologic biomarker would represent a major step towards more precise and personalized care in PD.

Controlling Clinical States Governed by Different Temporal Dynamics With Closed-Loop Deep Brain Stimulation: A Principled Framework

Closed-loop strategies for deep brain stimulation (DBS) are paving the way for improving the efficacy of existing neuromodulation therapies across neurological disorders. Unlike continuous DBS, closed-loop DBS approaches (cl-DBS) optimize the delivery of stimulation in the temporal domain. However, clinical and neurophysiological manifestations exhibit highly diverse temporal properties and evolve over multiple time-constants. Moreover, throughout the day, patients are engaged in different activities such as walking, talking, or sleeping that may require specific therapeutic adjustments. This broad range of temporal properties, along with inter-dependencies affecting parallel manifestations, need to be integrated in the development of therapies to achieve a sustained, optimized control of multiple symptoms over time. This requires an extended view on future cl-DBS design. Here we propose a conceptual framework to guide the development of multi-objective therapies embedding parallel control loops. Its modular organization allows to optimize the personalization of cl-DBS therapies to heterogeneous patient profiles. We provide an overview of clinical states and symptoms, as well as putative electrophysiological biomarkers that may be integrated within this structure. This integrative framework may guide future developments and become an integral part of next-generation precision medicine instruments.

Inference of brain networks with approximate Bayesian computation - assessing face validity with an example application in Parkinsonism

This paper describes and validates a novel framework using the Approximate Bayesian Computation (ABC) algorithm for parameter estimation and model selection in models of mesoscale brain network activity. We provide a proof of principle, first pass validation of this framework using a set of neural mass models of the cortico-basal ganglia thalamic circuit inverted upon spectral features from experimental, in vivo recordings. This optimization scheme relaxes an assumption of fixed-form posteriors (i.e. the Laplace approximation) taken in previous approaches to inverse modelling of spectral features. This enables the exploration of model dynamics beyond that approximated from local linearity assumptions and so fit to explicit, numerical solutions of the underlying non-linear system of equations. In this first paper, we establish a face validation of the optimization procedures in terms of: (i) the ability to approximate posterior densities over parameters that are plausible given the known causes of the data; (ii) the ability of the model comparison procedures to yield posterior model probabilities that can identify the model structure known to generate the data; and (iii) the robustness of these procedures to local minima in the face of different starting conditions. Finally, as an illustrative application we show (iv) that model comparison can yield plausible conclusions given the known neurobiology of the cortico-basal ganglia-thalamic circuit in Parkinsonism. These results lay the groundwork for future studies utilizing highly nonlinear or brittle models that can explain time dependant dynamics, such as oscillatory bursts, in terms of the underlying neural circuits.