Features for detection of Parkinson's disease tremor from local field potentials of the subthalamic nucleus

Deep Brain Stimulation (DBS) is a treatment routinely used to alleviate the symptoms of Parkinson's disease (PD). In this type of treatment, electrical pulses are applied through electrodes implanted into the basal ganglia of the patient. As the symptoms are not permanent in most patients, it is desirable to develop an on-demand stimulator, applying pulses only when onset of the symptoms is detected. This study evaluates a feature set created for the detection of tremor - a cardinal symptom of PD. The designed feature set was based on standard signal features and researched properties of the electrical signals recorded from subthalamic nucleus (STN) within the basal ganglia, which together included temporal, spectral, statistical, autocorrelation and fractal properties. The most characterized tremor related features were selected using statistical testing and backward algorithms then used for classification on unseen patient signals. The spectral features were among the most efficient at detecting tremor, notably spectral bands 3.5-5.5 Hz and 0-1 Hz proved to be highly significant. The classification results for determination of tremor achieved 94% sensitivity with specificity equaling one.

Balanced excitation and inhibition: Model based analysis of local field potentials

We have developed a model of the local field potential (LFP) based on the conservation of charge, the independence principle of ionic flows and the classical Hodgkin–Huxley (HH) type intracellular model of synaptic activity. Insights were gained through the simulation of the HH intracellular model on the nonlinear relationship between the balance of synaptic conductances and that of post-synaptic currents. The latter is dependent not only on the former, but also on the temporal lag between the excitatory and inhibitory conductances, as well as the strength of the afferent signal. The proposed LFP model provides a method for decomposing the LFP recordings near the soma of layer IV pyramidal neurons in the barrel cortex of anaesthetised rats into two highly correlated components with opposite polarity. The temporal dynamics and the proportional balance of the two components are comparable to the excitatory and inhibitory post-synaptic currents computed from the HH model. This suggests that the two components of the LFP reflect the underlying excitatory and inhibitory post-synaptic currents of the local neural population. We further used the model to decompose a sequence of evoked LFP responses under repetitive electrical stimulation (5 Hz) of the whisker pad. We found that as neural responses adapted, the excitatory and inhibitory components also adapted proportionately, while the temporal lag between the onsets of the two components increased during frequency adaptation. Our results demonstrated that the balance between neural excitation and inhibition can be investigated using extracellular recordings. Extension of the model to incorporate multiple compartments should allow more quantitative interpretations of surface Electroencephalography (EEG) recordings into components reflecting the excitatory, inhibitory and passive ionic current flows generated by local neural populations.

Parkinsonian tremor identification with multiple local field potential feature classification

This paper explores the development of multi-feature classification techniques used to identify tremor-related characteristics in the Parkinsonian patient. Local field potentials were recorded from the subthalamic nucleus and the globus pallidus internus of eight Parkinsonian patients through the implanted electrodes of a Deep brain stimulation (DBS) device prior to device internalization. A range of signal processing techniques were evaluated with respect to their tremor detection capability and used as inputs in a multi-feature neural network classifier to identify the activity of Parkinsonian tremor. The results of this study show that a trained multi-feature neural network is able, under certain conditions, to achieve excellent detection accuracy on patients unseen during training. Overall the tremor detection accuracy was mixed, although an accuracy of over 86% was achieved in four out of the eight patients.

Cannabidivarin is anticonvulsant in mouse and rat in vitro and in seizure models

Summary Background and purpose: Phytocannabinoids in Cannabis sativa have diverse pharmacological targets extending beyond cannabinoid receptors and several exert notable anticonvulsant effects. For the first time, we investigated the anticonvulsant profile of the phytocannabinoid cannabidivarin (CBDV) in vitro and in in vivo seizure models. Experimental approach: The effect of CBDV (1-100μM) on epileptiform local field potentials (LFPs) induced in rat hippocampal brain slices by 4-AP application or Mg2+-free conditions was assessed by in vitro multi-electrode array recordings. Additionally, the anticonvulsant profile of CBDV (50-200 mg kg-1) in vivo was investigated in four rodent seizure models: maximal electroshock (mES) and audiogenic seizures in mice, and pentylenetetrazole (PTZ) and pilocarpine-induced seizures in rat. CBDV effects in combination with commonly-used antiepileptic drugs were investigated in rat seizures. Finally, the motor side effect profile of CBDV was investigated using static beam and gripstrength assays. Key results: CDBV significantly attenuated status epilepticus-like epileptiform LFPs induced by 4-AP and Mg2+-free conditions. CBDV had significant anticonvulsant effects in mES (≥100 mg kg-1), audiogenic (≥50 mg kg-1) and PTZ-induced seizures (≥100 mg kg-1). CBDV alone had no effect against pilocarpine-induced seizures, but significantly attenuated these seizures when administered with valproate or phenobarbital at 200 mg kg-1 CBDV. CBDV had no effect on motor function. Conclusions and Implications: These results indicate that CBDV is an effective anticonvulsant across a broad range of seizure models, does not significantly affect normal motor function and therefore merits further investigation in chronic epilepsy models to justify human trials.

When long-range zero-lag synchronization is feasible in cortical networks

Many studies have reported long-range synchronization of neuronal activity between brain areas, in particular in the beta and gamma bands with frequencies in the range of 14–30 and 40–80 Hz, respectively. Several studies have reported synchrony with zero phase lag, which is remarkable considering the synaptic and conduction delays inherent in the connections between distant brain areas. This result has led to many speculations about the possible functional role of zero-lag synchrony, such as for neuronal communication, attention, memory, and feature binding. However, recent studies using recordings of single-unit activity and local field potentials report that neuronal synchronization may occur with non-zero phase lags. This raises the questions whether zero-lag synchrony can occur in the brain and, if so, under which conditions. We used analytical methods and computer simulations to investigate which connectivity between neuronal populations allows or prohibits zero-lag synchrony. We did so for a model where two oscillators interact via a relay oscillator. Analytical results and computer simulations were obtained for both type I Mirollo–Strogatz neurons and type II Hodgkin–Huxley neurons. We have investigated the dynamics of the model for various types of synaptic coupling and importantly considered the potential impact of Spike-Timing Dependent Plasticity (STDP) and its learning window. We confirm previous results that zero-lag synchrony can be achieved in this configuration. This is much easier to achieve with Hodgkin–Huxley neurons, which have a biphasic phase response curve, than for type I neurons. STDP facilitates zero-lag synchrony as it adjusts the synaptic strengths such that zero-lag synchrony is feasible for a much larger range of parameters than without STDP.

Classification of cortical microcircuits based on micro-electrode-array data from slices of rat barrel cortex

The bewildering complexity of cortical microcircuits at the single cell level gives rise to surprisingly robust emergent activity patterns at the level of laminar and columnar local field potentials (LFPs) in response to targeted local stimuli. Here we report the results of our multivariate data-analytic approach based on simultaneous multi-site recordings using micro-electrode-array chips for investigation of the microcircuitary of rat somatosensory (barrel) cortex. We find high repeatability of stimulus-induced responses, and typical spatial distributions of LFP responses to stimuli in supragranular, granular, and infragranular layers, where the last form a particularly distinct class. Population spikes appear to travel with about 33 cm/s from granular to infragranular layers. Responses within barrel related columns have different profiles than those in neighbouring columns to the left or interchangeably to the right. Variations between slices occur, but can be minimized by strictly obeying controlled experimental protocols. Cluster analysis on normalized recordings indicates specific spatial distributions of time series reflecting the location of sources and sinks independent of the stimulus layer. Although the precise correspondences between single cell activity and LFPs are still far from clear, a sophisticated neuroinformatics approach in combination with multi-site LFP recordings in the standardized slice preparation is suitable for comparing normal conditions to genetically or pharmacologically altered situations based on real cortical microcircuitry.

A dynamic causal model of the coupling between pulse stimulation and neural activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

The resting-state neurovascular coupling relationship: rapid changes in spontaneous neural activity in the somatosensory cortex are associated with haemodynamic fluctuations that resemble stimulus-evoked haemodynamics

Although promise exists for patterns of resting-state blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) brain connectivity to be used as biomarkers of early brain pathology, a full understanding of the nature of the relationship between neural activity and spontaneous fMRI BOLD fluctuations is required before such data can be correctly interpreted. To investigate this issue, we combined electrophysiological recordings of rapid changes in multi-laminar local field potentials from the somatosensory cortex of anaesthetized rats with concurrent two-dimensional optical imaging spectroscopy measurements of resting-state haemodynamics that underlie fluctuations in the BOLD fMRI signal. After neural ‘events’ were identified, their time points served to indicate the start of an epoch in the accompanying haemodynamic fluctuations. Multiple epochs for both neural ‘events’ and the accompanying haemodynamic fluctuations were averaged. We found that the averaged epochs of resting-state haemodynamic fluctuations taken after neural ‘events’ closely resembled the temporal profile of stimulus-evoked cortical haemodynamics. Furthermore, we were able to demonstrate that averaged epochs of resting-state haemodynamic fluctuations resembling the temporal profile of stimulus-evoked haemodynamics could also be found after peaks in neural activity filtered into specific electroencephalographic frequency bands (theta, alpha, beta, and gamma). This technique allows investigation of resting-state neurovascular coupling using methodologies that are directly comparable to that developed for investigating stimulus-evoked neurovascular responses.

Resting tremor classification and detection in Parkinson's disease patients

Parkinson is a neurodegenerative disease, in which tremor is the main symptom. This paper investigates the use of different classification methods to identify tremors experienced by Parkinsonian patients.Some previous research has focussed tremor analysis on external body signals (e.g., electromyography, accelerometer signals, etc.). Our advantage is that we have access to sub-cortical data, which facilitates the applicability of the obtained results into real medical devices since we are dealing with brain signals directly. Local field potentials (LFP) were recorded in the subthalamic nucleus of 7 Parkinsonian patients through the implanted electrodes of a deep brain stimulation (DBS) device prior to its internalization. Measured LFP signals were preprocessed by means of splinting, down sampling, filtering, normalization and rec-tification. Then, feature extraction was conducted through a multi-level decomposition via a wavelettrans form. Finally, artificial intelligence techniques were applied to feature selection, clustering of tremor types, and tremor detection.The key contribution of this paper is to present initial results which indicate, to a high degree of certainty, that there appear to be two distinct subgroups of patients within the group-1 of patients according to the Consensus Statement of the Movement Disorder Society on Tremor. Such results may well lead to different resultant treatments for the patients involved, depending on how their tremor has been classified. Moreover, we propose a new approach for demand driven stimulation, in which tremor detection is also based on the subtype of tremor the patient has. Applying this knowledge to the tremor detection problem, it can be concluded that the results improve when patient clustering is applied prior to detection.

Field line resonances as a trigger and a tracer for substorm onset

In this paper, we show that periodic auroral arc structures are seen at the location of one particular auroral substorm onset for the 15 min preceding onset, suggesting that field line resonances should be considered a strong candidate for triggering substorm onset. Irrespective of whether this field line resonance is coincidentally or causally linked to this substorm onset, the characteristics of the field line resonance can be used to remote sense the characteristics of the geomagnetic field line that supports substorm onset. In this instance, the eigenfrequency of this resonance is around 12 mHz. Interestingly, however, there is no evidence of this field line resonance in a seven satellite major Time History of Events and Macroscale Interactions during Substorms (THEMIS)-GOES conjunction, ranging from geosynchronous orbit to ~30 RE. However, using space-based cross-phase measurements of the local field line eigenfrequency at the inner THEMIS locations, we find that the local field line eigenfrequency is 6–10 mHz. Hence, we can reliably say that this 12 mHz Field Line Resonance (FLR) must lie inside of THEMIS locations. Our conclusion is that a high-m field line resonance can both represent a strong candidate for a trigger for substorm onset, as first proposed by Samson et al. (1992), and that its characteristics can provide invaluable information as to where substorm onset occurs in the magnetosphere.

Cannabidiol displays antiepileptiform and antiseizure properties in vitro and in vivo

Plant-derived cannabinoids (phytocannabinoids) are compounds with emerging therapeutic potential. Early studies suggested that cannabidiol (CBD) has anticonvulsant properties in animal models and reduced seizure frequency in limited human trials. Here, we examine the anti-epileptiform and anti-seizure potential of CBD using in vitro electrophysiology and an in vivo animal seizure model, respectively. CBD (0.01-100 muM) effects were assessed in vitro using the Mg(2+)-free and 4-aminopyridine (4-AP) models of status epilepticus-like epileptiform activity in hippocampal brain slices via multi-electrode array (MEA) recordings. In the Mg(2+)-free model, CBD decreased epileptiform local field potential (LFP) burst amplitude (in CA1 and dentate gyrus (DG) regions) and burst duration (in all regions) and increased burst frequency (in all regions). In the 4-AP model, CBD decreased LFP burst amplitude (in CA1, only at 100 muM CBD), burst duration (in CA3 and DG), and burst frequency (in all regions). CBD (1, 10 and 100 mg/kg) effects were also examined in vivo using the pentylenetetrazole (PTZ) model of generalised seizures. CBD (100 mg/kg) exerted clear anticonvulsant effects with significant decreases in incidence of severe seizures and mortality in comparison to vehicle-treated animals. Finally, CBD acted with only low affinity at cannabinoid CB(1) receptors and displayed no agonist activity in [(35)S]GTPgammaS assays in cortical membranes. These findings suggest that CBD acts to inhibit epileptiform activity in vitro and seizure severity in vivo. Thus, we demonstrate the potential of CBD as a novel anti-epileptic drug (AED) in the unmet clinical need associated with generalised seizures.

Cholinergic modulation of intrinsic fibre-evoked excitatory transmission contains a nicotinic component in immature but not adult rat piriform cortex, in vitro

The piriform cortex (PC) is highly prone to epileptogenesis, particularly in immature animals, where decreased muscarinic modulation of PC intrinsic fibre excitatory neurotransmission is implicated as a likely cause. However, whether higher levels of acetylcholine (ACh) release occur in immature vs. adult PC remains unclear. We investigated this using in vitro extracellular electrophysiological recording techniques. Intrinsic fibre-evoked extracellular field potentials (EFPs) were recorded from layers II to III in PC brain slices prepared from immature (P14-18) and adult (P>40) rats. Adult and immature PC EFPs were suppressed by eserine (1muM) or neostigmine (1muM) application, with a greater suppression in immature ( approximately 40%) than adult ( approximately 30%) slices. Subsequent application of atropine (1muM) reversed EFP suppression, producing supranormal ( approximately 12%) recovery in adult slices, suggesting that suppression was solely muscarinic ACh receptor-mediated and that some 'basal' cholinergic 'tone' was present. Conversely, atropine only partially reversed anticholinesterase effects in immature slices, suggesting the presence of additional non-muscarinic modulation. Accordingly, nicotine (50muM) caused immature field suppression ( approximately 30%) that was further enhanced by neostigmine, whereas it had no effect on adult EFPs. Unlike atropine, nicotinic antagonists, mecamylamine and methyllycaconitine, induced immature supranormal field recovery ( approximately 20%) following anticholinesterase-induced suppression (with no effect on adult slices), confirming that basal cholinergic 'tone' was also present. We suggest that nicotinic inhibitory cholinergic modulation occurs in the immature rat PC intrinsic excitatory fibre system, possibly to complement the existing, weak muscarinic modulation, and could be another important developmentally regulated system governing immature PC susceptibility towards epileptogenesis.

Prediction of Parkinson’s Disease tremor onset using a radial basis function neural network based on particle swarm optimization

Deep Brain Stimulation (DBS) has been successfully used throughout the world for the treatment of Parkinson's disease symptoms. To control abnormal spontaneous electrical activity in target brain areas DBS utilizes a continuous stimulation signal. This continuous power draw means that its implanted battery power source needs to be replaced every 18–24 months. To prolong the life span of the battery, a technique to accurately recognize and predict the onset of the Parkinson's disease tremors in human subjects and thus implement an on-demand stimulator is discussed here. The approach is to use a radial basis function neural network (RBFNN) based on particle swarm optimization (PSO) and principal component analysis (PCA) with Local Field Potential (LFP) data recorded via the stimulation electrodes to predict activity related to tremor onset. To test this approach, LFPs from the subthalamic nucleus (STN) obtained through deep brain electrodes implanted in a Parkinson patient are used to train the network. To validate the network's performance, electromyographic (EMG) signals from the patient's forearm are recorded in parallel with the LFPs to accurately determine occurrences of tremor, and these are compared to the performance of the network. It has been found that detection accuracies of up to 89% are possible. Performance comparisons have also been made between a conventional RBFNN and an RBFNN based on PSO which show a marginal decrease in performance but with notable reduction in computational overhead.

Parkinson’s Disease tremor classification – a comparison between Support Vector Machines and neural networks

Deep Brain Stimulation has been used in the study of and for treating Parkinson’s Disease (PD) tremor symptoms since the 1980s. In the research reported here we have carried out a comparative analysis to classify tremor onset based on intraoperative microelectrode recordings of a PD patient’s brain Local Field Potential (LFP) signals. In particular, we compared the performance of a Support Vector Machine (SVM) with two well known artificial neural network classifiers, namely a Multiple Layer Perceptron (MLP) and a Radial Basis Function Network (RBN). The results show that in this study, using specifically PD data, the SVM provided an overall better classification rate achieving an accuracy of 81% recognition.

A dynamic model of neurovascular coupling: Implications for blood vessel dilation and constriction

Neurovascular coupling in response to stimulation of the rat barrel cortex was investigated using concurrent multichannel electrophysiology and laser Doppler flowmetry. The data were used to build a linear dynamic model relating neural activity to blood flow. Local field potential time series were subject to current source density analysis, and the time series of a layer IV sink of the barrel cortex was used as the input to the model. The model output was the time series of the changes in regional cerebral blood flow (CBF). We show that this model can provide excellent fit of the CBF responses for stimulus durations of up to 16 s. The structure of the model consisted of two coupled components representing vascular dilation and constriction. The complex temporal characteristics of the CBF time series were reproduced by the relatively simple balance of these two components. We show that the impulse response obtained under the 16-s duration stimulation condition generalised to provide a good prediction to the data from the shorter duration stimulation conditions. Furthermore, by optimising three out of the total of nine model parameters, the variability in the data can be well accounted for over a wide range of stimulus conditions. By establishing linearity, classic system analysis methods can be used to generate and explore a range of equivalent model structures (e.g., feed-forward or feedback) to guide the experimental investigation of the control of vascular dilation and constriction following stimulation.