Inferring entire spiking activity from local field potentials

Degradation-aware neural imputation: Advancing decoding stability in brain machine interfaces

Neural signal degradation poses a significant challenge in maintaining stable performance when decoding motor tasks using multiunit activity (MUA) and local field potential (LFP) signals in the implantable brain machine interface (iBMI) applications. Effective methods for imputing degraded or missing signals are essential to restore neural signal integrity, thereby improving decoding accuracy and system robustness over long-term recordings with fluctuating signal quality. This study introduces a confidence-weighted Bayesian linear regression (CW-BLR) approach to impute neural signals affected by degradation, enhancing the robustness and consistency of decoding. The performance of CW-BLR was compared to traditional methods-mean imputation (Mean-imp) and Gaussian-mixture-model-based expectation-maximization (GMM-EM)-using a kernel-sliced inverse regression (kSIR) decoder to evaluate decoding outcomes. Four Wistar rats were trained to perform a forelimb-reaching task while neural activity (MUA and LFPs) was recorded over 27 days. CW-BLR imputed signals degraded during days 8-27. Decoding performance was evaluated using kSIR and compared with Mean-imp and GMM-EM. CW-BLR demonstrated superior performance by effectively preserving both temporal and spatial dependencies within the neural signals. CW-BLR-imputed data significantly improved decoding accuracy over traditional imputation methods, with the kSIR decoder showing consistently higher performance, particularly in maintaining signal quality from the degraded period. CW-BLR offers a robust and effective imputation framework for iBMI applications, addressing signal degradation challenges and maintaining accurate decoding over prolonged recordings. By utilizing confidence-based quality metrics, CW-BLR surpasses traditional methods, providing stable neural decoding across fluctuating signal quality scenarios.

Volitional Regulation and Transferable Patterns of Midbrain Oscillations

Dopaminergic brain areas are crucial for cognition and their dysregulation is linked to neuropsychiatric disorders typically treated with pharmacological interventions. These treatments often have side effects and variable effectiveness, underscoring the need for alternatives. We introduce the first demonstration of neurofeedback using local field potentials (LFP) from the ventral tegmental area (VTA). This approach leverages the real-time temporal resolution of LFP and ability to target deep brain. In our study, two male rhesus macaque monkeys (Macaca mulatta) learned to regulate VTA beta power using a customized normalized metric to stably quantify VTA LFP signal modulation. The subjects demonstrated flexible and specific control with different strategies for specific frequency bands, revealing new insights into the plasticity of VTA neurons contributing to oscillatory activity that is functionally relevant to many aspects of cognition. Excitingly, the subjects showed transferable patterns, a key criterion for clinical applications beyond training settings. This work provides a foundation for neurofeedback-based treatments, which may be a promising alternative to conventional approaches and open new avenues for understanding and managing neuropsychiatric disorders.

Three distinct gamma oscillatory networks within cortical columns in macaque monkeys' area V1

Introduction

A fundamental property of the neocortex is its columnar organization in many species. Generally, neurons of the same column share stimulus preferences and have strong anatomical connections across layers. These features suggest that neurons within a column operate as one unified network. Other features, like the different patterns of input and output connections of neurons located in separate layers and systematic differences in feature tuning, hint at a more segregated and possibly flexible functional organization of neurons within a column.

Methods

To distinguish between these views of columnar processing, we conducted laminar recordings in macaques' area V1 while they performed a demanding attention task. We identified three separate regions with strong gamma oscillatory activity, located in the supragranular, granular, and infragranular laminar domains, based on the current source density (CSD).

Results and Discussion

Their characteristics differed significantly in their dominant gamma frequency and attention-dependent modulation of their gramma power and gamma frequency. In line, spiking activity in the supragranular, infragranular, and upper part of the granular domain exhibited strong phase coherence with the CSD signals of their domain but showed much weaker coherence with the CSD signals of other domains.

Conclusion

These results indicate that columnar processing involves a certain degree of independence between neurons in the three laminar domains, consistent with the assumption of multiple, separate intracolumnar ensembles. Such a functional organization offers various possibilities for dynamic network configuration, indicating that neurons in a column are not restricted to operate as one unified network. Thus, the findings open interesting new possibilities for future concepts and investigations on flexible, dynamic cortical ensemble formation and selective information processing.

Filtered Point Processes Tractably Capture Rhythmic And Broadband Power Spectral Structure in Neural Electrophysiological Recordings

Neural electrophysiological recordings arise from interacting rhythmic (oscillatory) and broadband (aperiodic) biological subprocesses. Both rhythmic and broadband processes contribute to the neural power spectrum, which decomposes the variance of a neural recording across frequencies. Although an extensive body of literature has successfully studied rhythms in various diseases and brain states, researchers only recently have systematically studied the characteristics of broadband effects in the power spectrum. Broadband effects can generally be categorized as 1) shifts in power across all frequencies, which correlate with changes in local firing rates and 2) changes in the overall shape of the power spectrum, such as the spectral slope or power law exponent. Shape changes are evident in various conditions and brain states, influenced by factors such as excitation to inhibition balance, age, and various diseases. It is increasingly recognized that broadband and rhythmic effects can interact on a sub-second timescale. For example, broadband power is time-locked to the phase of <1 Hz rhythms in propofol induced unconsciousness. Modeling tools that explicitly deal with both rhythmic and broadband contributors to the power spectrum and that capture their interactions are essential to help improve the interpretability of power spectral effects. Here, we introduce a tractable stochastic forward modeling framework designed to capture both narrowband and broadband spectral effects when prior knowledge or theory about the primary biophysical processes involved is available. Population-level neural recordings are modeled as the sum of filtered point processes (FPPs), each representing the contribution of a different biophysical process such as action potentials or postsynaptic potentials of different types. Our approach builds on prior neuroscience FPP work by allowing multiple interacting processes and time-varying firing rates and by deriving theoretical power spectra and cross-spectra. We demonstrate several properties of the models, including that they divide the power spectrum into frequency ranges dominated by rhythmic and broadband effects, and that they can capture spectral effects across multiple timescales, including sub-second cross-frequency coupling. The framework can be used to interpret empirically observed power spectra and cross-frequency coupling effects in biophysical terms, which bridges the gap between theoretical models and experimental results.

Reducing power requirements for high-accuracy decoding in iBCIs

Objective.Current intracortical brain-computer interfaces (iBCIs) rely predominantly on threshold crossings ('spikes') for decoding neural activity into a control signal for an external device. Spiking data can yield high accuracy online control during complex behaviors; however, its dependence on high-sampling-rate data collection can pose challenges. An alternative signal for iBCI decoding is the local field potential (LFP), a continuous-valued signal that can be acquired simultaneously with spiking activity. However, LFPs are seldom used alone for online iBCI control as their decoding performance has yet to achieve parity with spikes.Approach.Here, we present a strategy to improve the performance of LFP-based decoders by first training a neural dynamics model to use LFPs to reconstruct the firing rates underlying spiking data, and then decoding from the estimated rates. We test these models on previously-collected macaque data during center-out and random-target reaching tasks as well as data collected from a human iBCI participant during attempted speech.Main results.In all cases, training models from LFPs enables firing rate reconstruction with accuracy comparable to spiking-based dynamics models. In addition, LFP-based dynamics models enable decoding performance exceeding that of LFPs alone and approaching that of spiking-based models. In all applications except speech, LFP-based dynamics models also facilitate decoding accuracy exceeding that of direct decoding from spikes.Significance.Because LFP-based dynamics models operate on lower bandwidth and with lower sampling rate than spiking models, our findings indicate that iBCI devices can be designed to operate with lower power requirements than devices dependent on recorded spiking activity, without sacrificing high-accuracy decoding.

Layer-specific dynamics of local field potentials in monkey V1 during electrical stimulation

The mammalian neocortex, organized into six cellular layers or laminae, forms a cortical network within layers. Layer specific computations are crucial for sensory processing of visual stimuli within primary visual cortex. Laminar recordings of local field potentials (LFPs) are a powerful tool to study neural activity within cortical layers. Electric brain stimulation is widely used in basic neuroscience and in a large range of clinical applications. However, the layer-specific effects of electric stimulation on LFPs remain unclear. To address this gap, we conducted laminar LFP recordings of the primary visual cortex in monkeys while presenting a flash visual stimulus. Simultaneously, we applied a low frequency sinusoidal current to the occipital lobe with offset frequency to the flash stimulus repetition rate. We analyzed the modulation of visual-evoked potentials with respect to the applied phase of the electric stimulation. Our results reveal that only the deeper layers, but not the superficial layers, show phase-dependent changes in LFP components with respect to the applied current. Employing a cortical column model, we demonstrate that these in vivo observations can be explained by phase-dependent changes in the driving force within neurons of deeper layers. Our findings offer crucial insight into the selective modulation of cortical layers through electric stimulation, thus advancing approaches for more targeted neuromodulation.

Modeling seizure networks in neuron-glia cultures using microelectrode arrays

Epilepsy is a common neurological disorder, affecting over 65 million people worldwide. Unfortunately, despite resective surgery, over 30 % of patients with drug-resistant epilepsy continue to experience seizures. Retrospective studies considering connectivity using intracranial electrocorticography (ECoG) obtained during neuromonitoring have shown that treatment failure is likely driven by failure to consider critical components of the seizure network, an idea first formally introduced in 2002. However, current studies only capture snapshots in time, precluding the ability to consider seizure network development. Over the past few years, multiwell microelectrode arrays have been increasingly used to study neuronal networks in vitro. As such, we sought to develop a novel in vitro MEA seizure model to allow for study of seizure networks. Specifically, we used 4-aminopyridine (4-AP) to capture hyperexcitable activity, and then show increased network changes after 2 days of chronic treatment. We characterize network changes using functional connectivity measures and a novel technique using dimensionality reduction. We find that 4-AP successfully captures persistently elevated mean firing rate and significant changes in underlying connectivity patterns. We believe this affords a robust in vitro seizure model from which longitudinal network changes can be studied, laying groundwork for future studies exploring seizure network development.

Comparison of alterations in local field potentials and neuronal firing in mouse M1 and CA1 associated with central fatigue induced by high-intensity interval training and moderate-intensity continuous training

Background

The mechanisms underlying central fatigue (CF) induced by high-intensity interval training (HIIT) and moderate-intensity continuous training (MICT) are still not fully understood.

Methods

In order to explore the effects of these exercises on the functioning of cortical and subcortical neural networks, this study investigated the effects of HIIT and MICT on local field potential (LFP) and neuronal firing in the mouse primary motor cortex (M1) and hippocampal CA1 areas. HIIT and MICT were performed on C57BL/6 mice, and simultaneous multichannel recordings were conducted in the M1 motor cortex and CA1 hippocampal region.

Results

A range of responses were elicited, including a decrease in coherence values of LFP rhythms in both areas, and an increase in slow and a decrease in fast power spectral density (PSD, n = 7-9) respectively. HIIT/MICT also decreased the gravity frequency (GF, n = 7-9) in M1 and CA1. Both exercises decreased overall firing rates, increased time lag of firing, declined burst firing rates and the number of spikes in burst, and reduced burst duration (BD) in M1 and CA1 (n = 7-9). While several neuronal firing properties showed a recovery tendency, the alterations of LFP parameters were more sustained during the 10-min post-HIIT/MICT period. MICT appeared to be more effective than HIIT in affecting LFP parameters, neuronal firing rate, and burst firing properties, particularly in CA1. Both exercises significantly affected neural network activities and local neuronal firing in M1 and CA1, with MICT associated with a more substantial and consistent suppression of functional integration between M1 and CA1.

Conclusion

Our study provides valuable insights into the neural mechanisms involved in exercise-induced central fatigue by examining the changes in functional connectivity and coordination between the M1 and CA1 regions. These findings may assist individuals engaged in exercise in optimizing their exercise intensity and timing to enhance performance and prevent excessive fatigue. Additionally, the findings may have clinical implications for the development of interventions aimed at managing conditions related to exercise-induced fatigue.

Multiunit Frontal Eye Field Activity Codes the Visuomotor Transformation, But Not Gaze Prediction or Retrospective Target Memory, in a Delayed Saccade Task

Single-unit (SU) activity-action potentials isolated from one neuron-has traditionally been employed to relate neuronal activity to behavior. However, recent investigations have shown that multiunit (MU) activity-ensemble neural activity recorded within the vicinity of one microelectrode-may also contain accurate estimations of task-related neural population dynamics. Here, using an established model-fitting approach, we compared the spatial codes of SU response fields with corresponding MU response fields recorded from the frontal eye fields (FEFs) in head-unrestrained monkeys (Macaca mulatta) during a memory-guided saccade task. Overall, both SU and MU populations showed a simple visuomotor transformation: the visual response coded target-in-eye coordinates, transitioning progressively during the delay toward a future gaze-in-eye code in the saccade motor response. However, the SU population showed additional secondary codes, including a predictive gaze code in the visual response and retention of a target code in the motor response. Further, when SUs were separated into regular/fast spiking neurons, these cell types showed different spatial code progressions during the late delay period, only converging toward gaze coding during the final saccade motor response. Finally, reconstructing MU populations (by summing SU data within the same sites) failed to replicate either the SU or MU pattern. These results confirm the theoretical and practical potential of MU activity recordings as a biomarker for fundamental sensorimotor transformations (e.g., target-to-gaze coding in the oculomotor system), while also highlighting the importance of SU activity for coding more subtle (e.g., predictive/memory) aspects of sensorimotor behavior.

Spiking Laguerre Volterra networks-predicting neuronal activity from local field potentials

Objective.Understanding the generative mechanism between local field potentials (LFP) and neuronal spiking activity is a crucial step for understanding information processing in the brain. Up to now, most approaches have relied on simply quantifying the coupling between LFP and spikes. However, very few have managed to predict the exact timing of spike occurrence based on LFP variations.Approach.Here, we fill this gap by proposing novel spiking Laguerre-Volterra network (sLVN) models to describe the dynamic LFP-spike relationship. Compared to conventional artificial neural networks, the sLVNs are interpretable models that provide explainable features of the underlying dynamics.Main results.The proposed networks were applied on extracellular microelectrode recordings of Parkinson's Disease patients during deep brain stimulation (DBS) surgery. Based on the predictability of the LFP-spike pairs, we detected three neuronal populations with unique signal characteristics and sLVN model features.Significance.These clusters were indirectly associated with motor score improvement following DBS surgery, warranting further investigation into the potential of spiking activity predictability as an intraoperative biomarker for optimal DBS lead placement.

Modeling and Analysis of Environmental Electromagnetic Interference in Multiple-Channel Neural Recording Systems for High Common-Mode Interference Rejection Performance

Environmental electromagnetic interference (EMI) has always been a major interference source for multiple-channel neural recording systems, and little theoretical work has been attempted to address it. In this paper, equivalent circuit models are proposed to model both electromagnetic interference sources and neural signals in such systems, and analysis has been performed to generate the design guidelines for neural probes and the subsequent recording circuit towards higher common-mode interference (CMI) rejection performance while maintaining the recorded neural action potential (AP) signal quality. In vivo animal experiments with a configurable 32-channel neural recording system are carried out to validate the proposed models and design guidelines. The results show the power spectral density (PSD) of environmental 50 Hz EMI interference is reduced by three orders from 4.43 × 10-3 V2/Hz to 4.04 × 10-6 V2/Hz without affecting the recorded AP signal quality in an unshielded experiment environment.

Large-scale neurophysiology and single-cell profiling in human neuroscience

Advances in large-scale single-unit human neurophysiology, single-cell RNA sequencing, spatial transcriptomics and long-term ex vivo tissue culture of surgically resected human brain tissue have provided an unprecedented opportunity to study human neuroscience. In this Perspective, we describe the development of these paradigms, including Neuropixels and recent brain-cell atlas efforts, and discuss how their convergence will further investigations into the cellular underpinnings of network-level activity in the human brain. Specifically, we introduce a workflow in which functionally mapped samples of human brain tissue resected during awake brain surgery can be cultured ex vivo for multi-modal cellular and functional profiling. We then explore how advances in human neuroscience will affect clinical practice, and conclude by discussing societal and ethical implications to consider. Potential findings from the field of human neuroscience will be vast, ranging from insights into human neurodiversity and evolution to providing cell-type-specific access to study and manipulate diseased circuits in pathology. This Perspective aims to provide a unifying framework for the field of human neuroscience as we welcome an exciting era for understanding the functional cytoarchitecture of the human brain.

Granule cells perform frequency-dependent pattern separation in a computational model of the dentate gyrus

Mnemonic discrimination (MD) may be dependent on oscillatory perforant path input frequencies to the hippocampus in a "U"-shaped fashion, where some studies show that slow and fast input frequencies support MD, while other studies show that intermediate frequencies disrupt MD. We hypothesize that pattern separation (PS) underlies frequency-dependent MD performance. We aim to study, in a computational model of the hippocampal dentate gyrus (DG), the network and cellular mechanisms governing this putative "U"-shaped PS relationship. We implemented a biophysical model of the DG that produces the hypothesized "U"-shaped input frequency-PS relationship, and its associated oscillatory electrophysiological signatures. We subsequently evaluated the network's PS ability using an adapted spatiotemporal task. We undertook systematic lesion studies to identify the network-level mechanisms driving the "U"-shaped input frequency-PS relationship. A minimal circuit of a single granule cell (GC) stimulated with oscillatory inputs was also used to study potential cellular-level mechanisms. Lesioning synapses onto GCs did not impact the "U"-shaped input frequency-PS relationship. Furthermore, GC inhibition limits PS performance for fast frequency inputs, while enhancing PS for slow frequency inputs. GC interspike interval was found to be input frequency dependent in a "U"-shaped fashion, paralleling frequency-dependent PS observed at the network level. Additionally, GCs showed an attenuated firing response for fast frequency inputs. We conclude that independent of network-level inhibition, GCs may intrinsically be capable of producing a "U"-shaped input frequency-PS relationship. GCs may preferentially decorrelate slow and fast inputs via spike timing reorganization and high frequency filtering.

Enhancing Prediction of Forelimb Movement Trajectory through a Calibrating-Feedback Paradigm Incorporating RAT Primary Motor and Agranular Cortical Ensemble Activity in the Goal-Directed Reaching Task

Complete reaching movements involve target sensing, motor planning, and arm movement execution, and this process requires the integration and communication of various brain regions. Previously, reaching movements have been decoded successfully from the motor cortex (M1) and applied to prosthetic control. However, most studies attempted to decode neural activities from a single brain region, resulting in reduced decoding accuracy during visually guided reaching motions. To enhance the decoding accuracy of visually guided forelimb reaching movements, we propose a parallel computing neural network using both M1 and medial agranular cortex (AGm) neural activities of rats to predict forelimb-reaching movements. The proposed network decodes M1 neural activities into the primary components of the forelimb movement and decodes AGm neural activities into internal feedforward information to calibrate the forelimb movement in a goal-reaching movement. We demonstrate that using AGm neural activity to calibrate M1 predicted forelimb movement can improve decoding performance significantly compared to neural decoders without calibration. We also show that the M1 and AGm neural activities contribute to controlling forelimb movement during goal-reaching movements, and we report an increase in the power of the local field potential (LFP) in beta and gamma bands over AGm in response to a change in the target distance, which may involve sensorimotor transformation and communication between the visual cortex and AGm when preparing for an upcoming reaching movement. The proposed parallel computing neural network with the internal feedback model improves prediction accuracy for goal-reaching movements.

Hippocampal memory reactivation during sleep is correlated with specific cortical states of the retrosplenial and prefrontal cortices

Episodic memories are thought to be stabilized through the coordination of cortico-hippocampal activity during sleep. However, the timing and mechanism of this coordination remain unknown. To investigate this, we studied the relationship between hippocampal reactivation and slow-wave sleep up and down states of the retrosplenial cortex (RTC) and prefrontal cortex (PFC). We found that hippocampal reactivations are strongly correlated with specific cortical states. Reactivation occurred during sustained cortical Up states or during the transition from up to down state. Interestingly, the most prevalent interaction with memory reactivation in the hippocampus occurred during sustained up states of the PFC and RTC, while hippocampal reactivation and cortical up-to-down state transition in the RTC showed the strongest coordination. Reactivation usually occurred within 150-200 msec of a cortical Up state onset, indicating that a buildup of excitation during cortical Up state activity influences the probability of memory reactivation in CA1. Conversely, CA1 reactivation occurred 30-50 msec before the onset of a cortical down state, suggesting that memory reactivation affects down state initiation in the RTC and PFC, but the effect in the RTC was more robust. Our findings provide evidence that supports and highlights the complexity of bidirectional communication between cortical regions and the hippocampus during sleep.

Neural Decoding for Intracortical Brain-Computer Interfaces

Brain-computer interfaces have revolutionized the field of neuroscience by providing a solution for paralyzed patients to control external devices and improve the quality of daily life. To accurately and stably control effectors, it is important for decoders to recognize an individual's motor intention from neural activity either by noninvasive or intracortical neural recording. Intracortical recording is an invasive way of measuring neural electrical activity with high temporal and spatial resolution. Herein, we review recent developments in neural signal decoding methods for intracortical brain-computer interfaces. These methods have achieved good performance in analyzing neural activity and controlling robots and prostheses in nonhuman primates and humans. For more complex paradigms in motor rehabilitation or other clinical applications, there remains more space for further improvements of decoders.

Neuronal selectivity for stimulus information determines prefrontal LFP gamma power regardless of task execution

Local field potential (LFP) power in the gamma frequency is modulated by cognitive variables during task execution. We sought to examine whether such modulations only emerge when task rules are established. We therefore analyzed neuronal firing and LFPs in different prefrontal subdivisions before and after the same monkeys were trained to perform cognitive tasks. Prior to task rule learning, sites containing neurons selective for stimuli already exhibited increased gamma power during and after the passive viewing of stimuli compared to the baseline period. Unexpectedly, when the same monkeys learned to maintain these stimuli in working memory, the elevation of gamma power above the baseline was diminished, despite an overall increase in firing rate. Learning and executing the task further decoupled LFP power from single neuron firing. Gamma power decreased at the time when subjects needed to make a judgment about whether two stimuli were the same or not, and differential gamma power was observed for matching and nonmatching stimuli. Our results indicate that prefrontal gamma power emerges spontaneously, not necessarily tied to a cognitive task being executed.

Brain signal predictions from multi-scale networks using a linearized framework

Simulations of neural activity at different levels of detail are ubiquitous in modern neurosciences, aiding the interpretation of experimental data and underlying neural mechanisms at the level of cells and circuits. Extracellular measurements of brain signals reflecting transmembrane currents throughout the neural tissue remain commonplace. The lower frequencies (≲ 300Hz) of measured signals generally stem from synaptic activity driven by recurrent interactions among neural populations and computational models should also incorporate accurate predictions of such signals. Due to limited computational resources, large-scale neuronal network models (≳ 10⁶ neurons or so) often require reducing the level of biophysical detail and account mainly for times of action potentials (‘spikes’) or spike rates. Corresponding extracellular signal predictions have thus poorly accounted for their biophysical origin.

Here we propose a computational framework for predicting spatiotemporal filter kernels for such extracellular signals stemming from synaptic activity, accounting for the biophysics of neurons, populations, and recurrent connections. Signals are obtained by convolving population spike rates by appropriate kernels for each connection pathway and summing the contributions. Our main results are that kernels derived via linearized synapse and membrane dynamics, distributions of cells, conduction delay, and volume conductor model allow for accurately capturing the spatiotemporal dynamics of ground truth extracellular signals from conductance-based multicompartment neuron networks. One particular observation is that changes in the effective membrane time constants caused by persistent synapse activation must be accounted for.

The work also constitutes a major advance in computational efficiency of accurate, biophysics-based signal predictions from large-scale spike and rate-based neuron network models drastically reducing signal prediction times compared to biophysically detailed network models. This work also provides insight into how experimentally recorded low-frequency extracellular signals of neuronal activity may be approximately linearly dependent on spiking activity. A new software tool LFPykernels serves as a reference implementation of the framework.

Understanding the brain’s function and activity in healthy and pathological states across spatial scales and times spanning entire lives is one of humanity’s great undertakings. In experimental and clinical work probing the brain’s activity, a variety of electric and magnetic measurement techniques are routinely applied. However interpreting the extracellularly measured signals remains arduous due to multiple factors, mainly the large number of neurons contributing to the signals and complex interactions occurring in recurrently connected neuronal circuits. To understand how neurons give rise to such signals, mechanistic modeling combined with forward models derived using volume conductor theory has proven to be successful, but this approach currently does not scale to the systems level (encompassing millions of neurons or more) where simplified or abstract neuron representations typically are used. Motivated by experimental findings implying approximately linear relationships between times of neuronal action potentials and extracellular population signals, we provide a biophysics-based method for computing causal filters relating spikes and extracellular signals that can be applied with spike times or rates of large-scale neuronal network models for predictions of population signals without relying on ad hoc approximations.

A machine learning approach to characterize sequential movement-related states in premotor and motor cortices

Non-human primate (NHP) movement kinematics have been decoded from spikes and local field potentials (LFPs) recorded during motor tasks. However, the potential of LFPs to provide network-like characterizations of neural dynamics during planning and execution of sequential movements requires further exploration. Is the aggregate nature of LFPs suitable to construct informative brain state descriptors of movement preparation and execution? To investigate this, we developed a framework to process LFPs based on machine learning classifiers and analyzed LFP from a primate, implanted with several microelectrode arrays covering the premotor cortex in both hemispheres and the primary motor cortex on one side. The primate performed a reach and grasp task, consisting of five consecutive states, starting from rest until a rewarding target (food) was attained. We use this five-state task to characterize brain activity and connectivity within eight frequency bands, using spectral amplitude and pair-wise correlations across electrodes as features. Our results show that we could best distinguish all five movement-related states using the highest frequency band (200-500Hz), yielding an 87% accuracy with spectral amplitude, and 60% with pair-wise electrode correlation. Further analyses characterized each movement-related state, showing differential neuronal population activity at above-gamma frequencies during the various stages of movement. In summary, our results show that the concerted use of novel machine-learning techniques with coarse grained queue broad signals such as LFPs may be successfully used to track and decode fine reach and grasp movement aspects across several brain regions.