Inferring Cortical Variability from Local Field Potentials

Intraoperative recordings from the posterior superior insula in awake humans with peripheral neuropathic pain

BACKGROUND AND OBJECTIVES

The activity profile of the posterior insula in neuropathic pain (NeP) remains largely unexplored. To address this and examine its modulation by somatosensory stimulation, we recorded local field potentials (LFP) in awake patients with NeP undergoing deep brain stimulation (DBS) electrode implantation to the posterior-superior insula (PSI) for analgesic purposes.

MATERIALS AND METHODS

Six patients (one woman; 32-45 years), experiencing refractory peripheral NeP and having previously responded to non-invasive stimulation of the PSI underwent stereotactic implantation of DBS electrodes to the PSI as part of a phase II clinical trial. The averaged power of frequencies of LFP and their peaks were calculated during rest and under thermal painful and mechanical non-painful stimulation.

RESULTS

At rest, amplitude peaks within the delta (average min-max.: 2.2 Hz; 1.3-3.7) and theta (6.1 Hz, varying between 5.7 and 6.8 Hz) bands were identified. Compared to rest, both tonic thermal painful, and mechanical non-painful stimulation led to similar mean decreases in gamma power (-24.46 ± 70.56, and -19.56 ± 3.08; respectively). Painful stimuli caused an increase in all the other frequency bands, mainly in alpha and beta ranges, while non-painful stimulation led to decreases in power in all frequencies above 4Hz. Painful tonic stimulation was associated with a significantly greater power variability, both in amplitude and frequency, compared to nonpainful mechanical stimulation.

CONCLUSION

The posterior insula resting state activity in awake patients with chronic NeP was characterized by predominant theta oscillations. Painful and innocuous stimulation led to opposite spectral changes, with a much larger variability across the whole frequency spectrum for painful stimuli, relative to both resting state and non-painful stimulation.

Aperiodic spectral slope tracks the effects of brain state on saliency responses in the human auditory cortex

Alteration of responses to salient stimuli occurs in a wide range of brain disorders and may be rooted in pathophysiological brain state dynamics. Specifically, tonic and phasic modes of activity in the reticular activating system (RAS) influence, and are influenced by, salient stimuli, respectively. The RAS influences the spectral characteristics of activity in the neocortex, shifting the balance between low- and high-frequency fluctuations. Aperiodic '1/f slope' has emerged as a promising composite measure of these brain state dynamics. However, the relationship of 1/f slope to state-dependent processes, such as saliency, is less explored, particularly intracranially in humans. Here, we record pupil diameter as a measure of brain state and intracranial local field potentials in auditory cortical regions of human patients during an auditory oddball stimulus paradigm. We find that phasic high-gamma band responses in auditory cortical regions exhibit an inverted-u shaped relationship to tonic state, as reflected in the 1/f slope. Furthermore, salient stimuli trigger state changes, as indicated by shifts in the 1/f slope. Taken together, these findings suggest that 1/f slope tracks tonic and phasic arousal state dynamics in the human brain, increasing the interpretability of this metric and supporting it as a potential biomarker in brain disorders.

Cortico-striatal beta oscillations as a reward-related signal

The value associated with reward is sensitive to external factors, such as the time between the choice and reward delivery as classically manipulated in temporal discounting tasks. Subjective preference for two reward options is dependent on objective variables of reward magnitude and reward delay. Single neuron correlates of reward value have been observed in regions, including ventral striatum, orbital, and medial prefrontal cortex. Brain imaging studies show cortico-striatal-limbic network activity related to subjective preferences. To explore how oscillatory dynamics represent reward processing across brain regions, we measured local field potentials of rats performing a temporal discounting task. Our goal was to use a data-driven approach to identify an electrophysiological marker that correlates with reward preference. We found that reward-locked oscillations at beta frequencies signaled the magnitude of reward and decayed with longer temporal delays. Electrodes in orbitofrontal/medial prefrontal cortex, anterior insula, ventral striatum, and amygdala individually increased power and were functionally connected at beta frequencies during reward outcome. Beta power during reward outcome correlated with subjective value as defined by a computational model fit to the discounting behavior. These data suggest that cortico-striatal beta oscillations are a reward signal correlated, which may represent subjective value and hold potential to serve as a biomarker and potential therapeutic target.

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.

Flexible low-profile external ventricular drain catheter for real-time brain monitoring

External ventricular drainage is one of the most common neurosurgical procedures in the world for acute hydrocephalus, which must be performed carefully by a neurosurgeon. Although various neuromonitoring external ventricular drain (EVD) catheters have been utilized, they still suffer from rigidity and bulkiness to mitigate post-EVD placement trauma. Here, we introduce a flexible and low-profile smart EVD catheter using a class of technologies with sensitive electrical materials, seamless integration, and flexible mechanics, which serves as a highly soft and minimally invasive device to monitor electrical brain signals. This device reliably captures biopotentials in real time while exhibiting remarkable flexibility and reliability. The seamless integration of its sensory system promises a minimally invasive EVD placement on brain tissue. This work validates the device's distinct characteristics and performances through in vitro experiments and computational analysis. Collectively, this device's exceptional patient- and user-friendly attributes highlight its potential as one of the most practical EVD catheters.

Beyond task response-Pre-stimulus activity modulates contents of consciousness

The current discussion on the neural correlates of the contents of consciousness (NCCc) focuses mainly on the post-stimulus period of task-related activity. This neglects the substantial impact of the spontaneous or ongoing activity of the brain as manifest in pre-stimulus activity. Does the interaction of pre- and post-stimulus activity shape the contents of consciousness? Addressing this gap in our knowledge, we review and converge two recent lines of findings, that is, pre-stimulus alpha power and pre- and post-stimulus alpha trial-to-trial variability (TTV). The data show that pre-stimulus alpha power modulates post-stimulus activity including specifically the subjective features of conscious contents like confidence and vividness. At the same time, alpha pre-stimulus variability shapes post-stimulus TTV reduction including the associated contents of consciousness. We propose that non-additive rather than merely additive interaction of the internal pre-stimulus activity with the external stimulus in the alpha band is key for contents to become conscious. This is mediated by mechanisms on different levels including neurophysiological, neurocomputational, neurodynamic, neuropsychological and neurophenomenal levels. Overall, considering the interplay of pre-stimulus intrinsic and post-stimulus extrinsic activity across wider timescales, not just evoked responses in the post-stimulus period, is critical for identifying neural correlates of consciousness. This is well in line with both processing and especially the Temporo-spatial theory of consciousness (TTC).

Behavioral Interventions Can Improve Brain Injury-Induced Deficits in Behavioral Flexibility and Impulsivity Linked to Impaired Reward-Feedback Beta Oscillations

Traumatic brain injury (TBI) affects a large population, resulting in severe cognitive impairments. Although cognitive rehabilitation is an accepted treatment for some deficits, studies in patients are limited in ability to probe physiological and behavioral mechanisms. Therefore, animal models are needed to optimize strategies. Frontal TBI in a rat model results in robust and replicable cognitive deficits, making this an ideal candidate for investigating various behavioral interventions. In this study, we report three distinct frontal TBI experiments assessing behavior well into the chronic post-injury period using male Long-Evans rats. First, we evaluated the impact of frontal injury on local field potentials recorded simultaneously from 12 brain regions during a probabilistic reversal learning (PbR) task. Next, a set of rats were tested on a similar PbR task or an impulsivity task (differential reinforcement of low-rate behavior [DRL]) and half received salient cues associated with reinforcement contingencies to encourage engagement in the target behavior. After intervention on the PbR task, brains were stained for markers of activity. On the DRL task, cue relevance was decoupled from outcomes to determine if beneficial effects persisted on impulsive behavior. TBI decreased the ability to detect reinforced outcomes; this was evident in task performance and reward-feedback signals occurring at beta frequencies in lateral orbitofrontal cortex (OFC) and associated frontostriatal regions. The behavioral intervention improved flexibility and increased OFC activity. Intervention also reduced impulsivity, even after cues were decoupled, which was partially mediated by improvements in timing behavior. The current study established a platform to begin investigating cognitive rehabilitation in rats and identified a strong role for dysfunctional OFC signaling in probabilistic learning after frontal TBI.

Hierarchical VAEs provide a normative account of motion processing in the primate brain

The relationship between perception and inference, as postulated by Helmholtz in the 19th century, is paralleled in modern machine learning by generative models like Variational Autoencoders (VAEs) and their hierarchical variants. Here, we evaluate the role of hierarchical inference and its alignment with brain function in the domain of motion perception. We first introduce a novel synthetic data framework, Retinal Optic Flow Learning (ROFL), which enables control over motion statistics and their causes. We then present a new hierarchical VAE and test it against alternative models on two downstream tasks: (i) predicting ground truth causes of retinal optic flow (e.g., self-motion); and (ii) predicting the responses of neurons in the motion processing pathway of primates. We manipulate the model architectures (hierarchical versus non-hierarchical), loss functions, and the causal structure of the motion stimuli. We find that hierarchical latent structure in the model leads to several improvements. First, it improves the linear decodability of ground truth factors and does so in a sparse and disentangled manner. Second, our hierarchical VAE outperforms previous state-of-the-art models in predicting neuronal responses and exhibits sparse latent-to-neuron relationships. These results depend on the causal structure of the world, indicating that alignment between brains and artificial neural networks depends not only on architecture but also on matching ecologically relevant stimulus statistics. Taken together, our results suggest that hierarchical Bayesian inference underlines the brain's understanding of the world, and hierarchical VAEs can effectively model this understanding.

Brain signatures indexing variation in internal processing during perceptual decision-making

Brain activity is highly variable during a task. Discovering, characterizing, and linking variability in brain activity to internal processes has primarily relied on experimental manipulations. However, changes in internal processing could arise from many factors independent of experimental conditions. Here we utilize a data-driven clustering method based on modularity-maximation to identify consistent spatial-temporal EEG activity patterns across individual trials. Subjects (N = 25) performed a motion discrimination task with six interleaved levels of coherence. Clustering identified two discrete subtypes of trials with different patterns of activity. Surprisingly, Subtype 1 occurred more frequently in trials with lower motion coherence but was associated with faster response times. Computational modeling suggests that Subtype 1 was characterized by a lower threshold for reaching a decision. These results highlight across-trial variability in decision processes traditionally hidden to experimenters and provide a method for identifying endogenous brain state variability relevant to cognition and behavior.

Resting state electroencephalographic rhythms are affected by immediately preceding memory demands in cognitively unimpaired elderly and patients with mild cognitive impairment

Experiments on event-related electroencephalographic oscillations in aged people typically include blocks of cognitive tasks with a few minutes of interval between them. The present exploratory study tested the effect of being engaged on cognitive tasks over the resting state cortical arousal after task completion, and whether it differs according to the level of the participant’s cognitive decline. To investigate this issue, we used a local database including data in 30 healthy cognitively unimpaired (CU) persons and 40 matched patients with amnestic mild cognitive impairment (aMCI). They had been involved in 2 memory tasks for about 40 min and underwent resting-state electroencephalographic (rsEEG) recording after 5 min from the task end. eLORETA freeware estimated rsEEG alpha source activity as an index of general cortical arousal. In the CU but not aMCI group, there was a negative correlation between memory tasks performance and posterior rsEEG alpha source activity. The better the memory tasks performance, the lower the posterior alpha activity (i.e., higher cortical arousal). There was also a negative correlation between neuropsychological test scores of global cognitive status and alpha source activity. These results suggest that engagement in memory tasks may perturb background brain arousal for more than 5 min after the tasks end, and that this effect are dependent on participants global cognitive status. Future studies in CU and aMCI groups may cross-validate and extend these results with experiments including (1) rsEEG recordings before memory tasks and (2) post-tasks rsEEG recordings after 5, 15, and 30 min.

Binocular rivalry reveals an out-of-equilibrium neural dynamics suited for decision-making

In ambiguous or conflicting sensory situations, perception is often 'multistable' in that it perpetually changes at irregular intervals, shifting abruptly between distinct alternatives. The interval statistics of these alternations exhibits quasi-universal characteristics, suggesting a general mechanism. Using binocular rivalry, we show that many aspects of this perceptual dynamics are reproduced by a hierarchical model operating out of equilibrium. The constitutive elements of this model idealize the metastability of cortical networks. Independent elements accumulate visual evidence at one level, while groups of coupled elements compete for dominance at another level. As soon as one group dominates perception, feedback inhibition suppresses supporting evidence. Previously unreported features in the serial dependencies of perceptual alternations compellingly corroborate this mechanism. Moreover, the proposed out-of-equilibrium dynamics satisfies normative constraints of continuous decision-making. Thus, multistable perception may reflect decision-making in a volatile world: integrating evidence over space and time, choosing categorically between hypotheses, while concurrently evaluating alternatives.

Model-based signal processing enables bidirectional inferring between local field potential and spikes evoked by noxious stimulation

BACKGROUND

Recording spontaneous and evoked activities by means of unitary extracellular recordings and local field potential (LFP) are key understanding the mechanisms of neural coding. The LFP is one of the most popular and easy methods to measure the activity of a population of neurons. LFP is also a composite signal known to be difficult to interpret and model. There is a growing need to highlight the relationship between spiking activity and LFP. Here, we hypothesized that LFP could be inferred from spikes under evoked noxious conditions.

METHOD

Recording was performed from the medullary dorsal horn (MDH) in deeply anesthetized rats. We detail a process to highlight the C-fiber (nociceptive) evoked activity, by removing the A-fiber evoked activity using a model-based approach. Then, we applied the convolution kernel theory and optimization algorithms to infer the C-fiber LFP from the single cell spikes. Finally, we used a probability density function and an optimization algorithm to infer the spikes distribution from the LFP.

RESULTS

We successfully extracted C-fiber LFP in all data recordings. We observed that C-fibers spikes preceded the C-fiber LFP and were rather correlated to the LFP derivative. Finally, we inferred LFP from spikes with excellent correlation coefficient (r = 0.9) and reverse generated the spikes distribution from LFP with good correlation coefficients (r = 0.7) on spikes number.

CONCLUSION

We introduced the kernel convolution theory to successfully infer the LFP from spikes, and we demonstrated that we could generate the spikes distribution from the LFP.

Behavior needs neural variability

Human and non-human animal behavior is highly malleable and adapts successfully to internal and external demands. Such behavioral success stands in striking contrast to the apparent instability in neural activity (i.e., variability) from which it arises. Here, we summon the considerable evidence across scales, species, and imaging modalities that neural variability represents a key, undervalued dimension for understanding brain-behavior relationships at inter- and intra-individual levels. We believe that only by incorporating a specific focus on variability will the neural foundation of behavior be comprehensively understood.

Alpha/beta power decreases during episodic memory formation predict the magnitude of alpha/beta power decreases during subsequent retrieval

Episodic memory retrieval is characterised by the vivid reinstatement of information about a personally-experienced event. Growing evidence suggests that this reinstatement is supported by reductions in the spectral power of alpha/beta activity. Given that the amount of information that can be recalled depends on the amount of information that was originally encoded, information-based accounts of alpha/beta activity would suggest that retrieval-related alpha/beta power decreases similarly depend upon decreases in alpha/beta power during encoding. To test this hypothesis, seventeen human participants completed a sequence-learning task while undergoing concurrent MEG recordings. Regression-based analyses were then used to estimate how alpha/beta power decreases during encoding predicted alpha/beta power decreases during retrieval on a trial-by-trial basis. When subjecting these parameter estimates to group-level analysis, we find evidence to suggest that retrieval-related alpha/beta (7-15Hz) power decreases fluctuate as a function of encoding-related alpha/beta power decreases. These results suggest that retrieval-related alpha/beta power decreases are contingent on the decrease in alpha/beta power that arose during encoding. Subsequent analysis uncovered no evidence to suggest that these alpha/beta power decreases reflect stimulus identity, indicating that the contingency between encoding- and retrieval-related alpha/beta power reflects the reinstatement of a neurophysiological operation, rather than neural representation, during episodic memory retrieval.

Inferring synaptic inputs from spikes with a conductance-based neural encoding model

Descriptive statistical models of neural responses generally aim to characterize the mapping from stimuli to spike responses while ignoring biophysical details of the encoding process. Here, we introduce an alternative approach, the conductance-based encoding model (CBEM), which describes a mapping from stimuli to excitatory and inhibitory synaptic conductances governing the dynamics of sub-threshold membrane potential. Remarkably, we show that the CBEM can be fit to extracellular spike train data and then used to predict excitatory and inhibitory synaptic currents. We validate these predictions with intracellular recordings from macaque retinal ganglion cells. Moreover, we offer a novel quasi-biophysical interpretation of the Poisson generalized linear model (GLM) as a special case of the CBEM in which excitation and inhibition are perfectly balanced. This work forges a new link between statistical and biophysical models of neural encoding and sheds new light on the biophysical variables that underlie spiking in the early visual pathway.

Alpha/beta power decreases track the fidelity of stimulus-specific information

Massed synchronised neuronal firing is detrimental to information processing. When networks of task-irrelevant neurons fire in unison, they mask the signal generated by task-critical neurons. On a macroscopic level, such synchronisation can contribute to alpha/beta (8-30 Hz) oscillations. Reducing the amplitude of these oscillations, therefore, may enhance information processing. Here, we test this hypothesis. Twenty-one participants completed an associative memory task while undergoing simultaneous EEG-fMRI recordings. Using representational similarity analysis, we quantified the amount of stimulus-specific information represented within the BOLD signal on every trial. When correlating this metric with concurrently-recorded alpha/beta power, we found a significant negative correlation which indicated that as post-stimulus alpha/beta power decreased, stimulus-specific information increased. Critically, we found this effect in three unique tasks: visual perception, auditory perception, and visual memory retrieval, indicating that this phenomenon transcends both stimulus modality and cognitive task. These results indicate that alpha/beta power decreases parametrically track the fidelity of both externally-presented and internally-generated stimulus-specific information represented within the cortex.

Integrated open-source software for multiscale electrophysiology

The methods for electrophysiology in neuroscience have evolved tremendously over the recent years with a growing emphasis on dense-array signal recordings. Such increased complexity and augmented wealth in the volume of data recorded, have not been accompanied by efforts to streamline and facilitate access to processing methods, which too are susceptible to grow in sophistication. Moreover, unsuccessful attempts to reproduce peer-reviewed publications indicate a problem of transparency in science. This growing problem could be tackled by unrestricted access to methods that promote research transparency and data sharing, ensuring the reproducibility of published results. Here, we provide a free, extensive, open-source software that provides data-analysis, data-management and multi-modality integration solutions for invasive neurophysiology. Users can perform their entire analysis through a user-friendly environment without the need of programming skills, in a tractable (logged) way. This work contributes to open-science, analysis standardization, transparency and reproducibility in invasive neurophysiology.

Data-Driven Approaches to Understanding Visual Neuron Activity

With modern neurophysiological methods able to record neural activity throughout the visual pathway in the context of arbitrarily complex visual stimulation, our understanding of visual system function is becoming limited by the available models of visual neurons that can be directly related to such data. Different forms of statistical models are now being used to probe the cellular and circuit mechanisms shaping neural activity, understand how neural selectivity to complex visual features is computed, and derive the ways in which neurons contribute to systems-level visual processing. However, models that are able to more accurately reproduce observed neural activity often defy simple interpretations. As a result, rather than being used solely to connect with existing theories of visual processing, statistical modeling will increasingly drive the evolution of more sophisticated theories.

Characterizing and dissociating multiple time-varying modulatory computations influencing neuronal activity

In many brain areas, sensory responses are heavily modulated by factors including attentional state, context, reward history, motor preparation, learned associations, and other cognitive variables. Modelling the effect of these modulatory factors on sensory responses has proven challenging, mostly due to the time-varying and nonlinear nature of the underlying computations. Here we present a computational model capable of capturing and dissociating multiple time-varying modulatory effects on neuronal responses on the order of milliseconds. The model’s performance is tested on extrastriate perisaccadic visual responses in nonhuman primates. Visual neurons respond to stimuli presented around the time of saccades differently than during fixation. These perisaccadic changes include sensitivity to the stimuli presented at locations outside the neuron’s receptive field, which suggests a contribution of multiple sources to perisaccadic response generation. Current computational approaches cannot quantitatively characterize the contribution of each modulatory source in response generation, mainly due to the very short timescale on which the saccade takes place. In this study, we use a high spatiotemporal resolution experimental paradigm along with a novel extension of the generalized linear model framework (GLM), termed the sparse-variable GLM, to allow for time-varying model parameters representing the temporal evolution of the system with a resolution on the order of milliseconds. We used this model framework to precisely map the temporal evolution of the spatiotemporal receptive field of visual neurons in the middle temporal area during the execution of a saccade. Moreover, an extended model based on a factorization of the sparse-variable GLM allowed us to disassociate and quantify the contribution of individual sources to the perisaccadic response. Our results show that our novel framework can precisely capture the changes in sensitivity of neurons around the time of saccades, and provide a general framework to quantitatively track the role of multiple modulatory sources over time.

The sensory responses of neurons in many brain areas, particularly those in higher prefrontal or parietal areas, are strongly influenced by factors including task rules, attentional state, context, reward history, motor preparation, learned associations, and other cognitive variables. These modulations often occur in combination, or on fast timescales which present a challenge for both experimental and modelling approaches aiming to describe the underlying mechanisms or computations. Here we present a computational model capable of capturing and dissociating multiple time-varying modulatory effects on spiking responses on the order of milliseconds. The model’s performance is evaluated by testing its ability to reproduce and dissociate multiple changes in visual sensitivity occurring in extrastriate visual cortex around the time of rapid eye movements. No previous model is capable of capturing these changes with as fine a resolution as that presented here. Our model both provides specific insight into the nature and time course of changes in visual sensitivity around the time of eye movements, and offers a general framework applicable to a wide variety of contexts in which sensory processing is modulated dynamically by multiple time-varying cognitive or behavioral factors, to understand the neuronal computations underpinning these modulations and make predictions about the underlying mechanisms.

The quest for interpretable models of neural population activity

Many aspects of brain function arise from the coordinated activity of large populations of neurons. Recent developments in neural recording technologies are providing unprecedented access to the activity of such populations during increasingly complex experimental contexts; however, extracting scientific insights from such recordings requires the concurrent development of analytical tools that relate this population activity to system-level function. This is a primary motivation for latent variable models, which seek to provide a low-dimensional description of population activity that can be related to experimentally controlled variables, as well as uncontrolled variables such as internal states (e.g. attention and arousal) and elements of behavior. While deriving an understanding of function from traditional latent variable methods relies on low-dimensional visualizations, new approaches are targeting more interpretable descriptions of the components underlying system-level function.

Single-Trial Decoding from Local Field Potential Using Bag of Word Representation

Neural decoding allows us to study the brain functions by investigating the relationship between a stimulus and the corresponding response. Recently, the local field potential (LFP) has been targeted as a hallmark of brain activity for neural decoding. Despite several decoding methods, there is still a lack of a comprehensive framework to decode cognitive functions in an integrated structure. Here, we addressed this issue by developing a dictionary-based method to represent the LFP signals via a bag-of-words (BOW) approach. First, we defined a general dictionary consisting of various Gabor wavelets as the words which enabled us to represent LFPs in word domain. For each trial, the LFP signal was convolved with the dictionary words. The integral of the absolute value and the mean phase of the complex output were considered as histogram weights. In the next step, using cross-validation leave-one-out method, the trials were split into the training and test sets. The weights of each individual word were swapped across trials within a certain category of the training set while the sequential order was maintained. Finally, the test trial was classified using label voting in the k-nearest training trials. We conducted the proposed method on two independent LFP data sets, recorded from the rat primary auditory cortex (A1) and monkey middle temporal area in order to evaluate its efficiency. In addition to the chance level, the proposed method was compared with a standard BOW approach that has been extended recently for biomedical signals classification. Results show a high efficiency (~ 15% improvement in decoding accuracy) of the proposed method. Together, the aforementioned method provides a comprehensive framework for single-trial decoding from short-length electrophysiological signals.

Phase shifts in high-beta- and low-gamma-band local field potentials predict the focus of visual spatial attention

The local field potential (LFP) contains rich information about activity in local neuronal populations. However, it has been challenging to establish direct links between LFP modulations and task-relevant behavior or cognitive processes, such as attention. We sought to determine whether LFP amplitude or phase modulations are predictive of the allocation of visual spatial attention. LFPs were recorded simultaneously in multiple early visual brain structures of alert macaque monkeys performing attention-demanding detection and discrimination tasks. Attention directed toward the receptive field of recorded neurons generated systematically larger phase shifts in high-beta- and low-gamma-frequency LFPs compared with LFP phase shifts on trials in which attention was directed away from the receptive field. This attention-mediated temporal advance corresponded to ~10 ms. LFP phase shifts also correlated with reaction times when monkeys were engaged in the tasks. Importantly, attentional modulation of LFP phase was consistent across monkeys, tasks, visual brain structures, and cortical layers. In contrast, attentional modulation of LFP amplitude varied across frequency bands, visual structures/layers, and tasks. Because LFP phase shifts were robust, consistent, and predictive of spatial attention, they could serve as a reliable marker for attention signals in the brain. NEW & NOTEWORTHY Local field potentials (LFPs) reflect the activity of spatially localized populations of neurons. Whether alterations in LFP activity are indicative of cognitive processes, such as attention, is unclear. We found that shifts in the phase of LFPs measured in multiple visual brain areas reliably predicted the focus of spatial attention. LFP phase shifts could therefore serve as a marker for behaviorally relevant attention signals in the brain.

Modeling the Correlated Activity of Neural Populations: A Review

The principles of neural encoding and computations are inherently collective and usually involve large populations of interacting neurons with highly correlated activities. While theories of neural function have long recognized the importance of collective effects in populations of neurons, only in the past two decades has it become possible to record from many cells simultaneously using advanced experimental techniques with single-spike resolution and to relate these correlations to function and behavior. This review focuses on the modeling and inference approaches that have been recently developed to describe the correlated spiking activity of populations of neurons. We cover a variety of models describing correlations between pairs of neurons, as well as between larger groups, synchronous or delayed in time, with or without the explicit influence of the stimulus, and including or not latent variables. We discuss the advantages and drawbacks or each method, as well as the computational challenges related to their application to recordings of ever larger populations.

Reducing the Effect of Spurious Phase Variations in Neural Oscillatory Signals

The phase-reset model of oscillatory EEG activity has received a lot of attention in the last decades for decoding different cognitive processes. Based on this model, the ERPs are assumed to be generated as a result of phase reorganization in ongoing EEG. Alignment of the phase of neuronal activities can be observed within or between different assemblies of neurons across the brain. Phase synchronization has been used to explore and understand perception, attentional binding and considering it in the domain of neuronal correlates of consciousness. The importance of the topic and its vast exploration in different domains of the neuroscience presses the need for appropriate tools and methods for measuring the level of phase synchronization of neuronal activities. Measuring the level of instantaneous phase (IP) synchronization has been used extensively in numerous studies of ERPs as well as oscillatory activity for a better understanding of the underlying cognitive binding with regard to different set of stimulations such as auditory and visual. However, the reliability of results can be challenged as a result of noise artifact in IP. Phase distortion due to environmental noise artifacts as well as different pre-processing steps on signals can lead to generation of artificial phase jumps. One of such effects presented recently is the effect of low envelope on the IP of signal. It has been shown that as the instantaneous envelope of the analytic signal approaches zero, the variations in the phase increase, effectively leading to abrupt transitions in the phase. These abrupt transitions can distort the phase synchronization results as they are not related to any neurophysiological effect. These transitions are called spurious phase variation. In this study, we present a model to remove generated artificial phase variations due to the effect of low envelope. The proposed method is based on a simplified form of a Kalman smoother, that is able to model the IP behavior in narrow-bandpassed oscillatory signals. In this work we first explain the details of the proposed Kalman smoother for modeling the dynamics of the phase variations in narrow-bandpassed signals and then evaluate it on a set of synthetic signals. Finally, we apply the model on ongoing-EEG signals to assess the removal of spurious phase variations.

A video-driven model of response statistics in the primate middle temporal area

Neurons in the primate middle temporal area (MT) encode information about visual motion and binocular disparity. MT has been studied intensively for decades, so there is a great deal of information in the literature about MT neuron tuning. In this study, our goal is to consolidate some of this information into a statistical model of the MT population response. The model accepts arbitrary stereo video as input. It uses computer-vision methods to calculate known correlates of the responses (such as motion velocity), and then predicts activity using a combination of tuning functions that have previously been used to describe data in various experiments. To construct the population response, we also estimate the distributions of many model parameters from data in the electrophysiology literature. We show that the model accounts well for a separate dataset of MT speed tuning that was not used in developing the model. The model may be useful for studying relationships between MT activity and behavior in ethologically relevant tasks. As an example, we show that the model can provide regression targets for internal activity in a deep convolutional network that performs a visual odometry task, so that its representations become more physiologically realistic.

Relating the Structure of Noise Correlations in Macaque Primary Visual Cortex to Decoder Performance

Noise correlations in neuronal responses can have a strong influence on the information available in large populations. In addition, the structure of noise correlations may have a great impact on the utility of different algorithms to extract this information that may depend on the specific algorithm, and hence may affect our understanding of population codes in the brain. Thus, a better understanding of the structure of noise correlations and their interplay with different readout algorithms is required. Here we use eigendecomposition to investigate the structure of noise correlations in populations of about 50–100 simultaneously recorded neurons in the primary visual cortex of anesthetized monkeys, and we relate this structure to the performance of two common decoders: the population vector and the optimal linear estimator. Our analysis reveals a non-trivial correlation structure, in which the eigenvalue spectrum is composed of several distinct large eigenvalues that represent different shared modes of fluctuation extending over most of the population, and a semi-continuous tail. The largest eigenvalue represents a uniform collective mode of fluctuation. The second and third eigenvalues typically show either a clear functional (i.e., dependent on the preferred orientation of the neurons) or spatial structure (i.e., dependent on the physical position of the neurons). We find that the number of shared modes increases with the population size, being roughly 10% of that size. Furthermore, we find that the noise in each of these collective modes grows linearly with the population. This linear growth of correlated noise power can have limiting effects on the utility of averaging neuronal responses across large populations, depending on the readout. Specifically, the collective modes of fluctuation limit the accuracy of the population vector but not of the optimal linear estimator.

A Streaming PCA VLSI Chip for Neural Data Compression

Neural recording system miniaturization and integration with low-power wireless technologies require compressing neural data before transmission. Feature extraction is a procedure to represent data in a low-dimensional space; its integration into a recording chip can be an efficient approach to compress neural data. In this paper, we propose a streaming principal component analysis algorithm and its microchip implementation to compress multichannel local field potential (LFP) and spike data. The circuits have been designed in a 65-nm CMOS technology and occupy a silicon area of 0.06 mm. Throughout the experiments, the chip compresses LFPs by 10 at the expense of as low as 1% reconstruction errors and 144-nW/channel power consumption; for spikes, the achieved compression ratio is 25 with 8% reconstruction errors and 3.05-W/channel power consumption. In addition, the algorithm and its hardware architecture can swiftly adapt to nonstationary spiking activities, which enables efficient hardware sharing among multiple channels to support a high-channel count recorder.

The HTM Spatial Pooler-A Neocortical Algorithm for Online Sparse Distributed Coding

Hierarchical temporal memory (HTM) provides a theoretical framework that models several key computational principles of the neocortex. In this paper, we analyze an important component of HTM, the HTM spatial pooler (SP). The SP models how neurons learn feedforward connections and form efficient representations of the input. It converts arbitrary binary input patterns into sparse distributed representations (SDRs) using a combination of competitive Hebbian learning rules and homeostatic excitability control. We describe a number of key properties of the SP, including fast adaptation to changing input statistics, improved noise robustness through learning, efficient use of cells, and robustness to cell death. In order to quantify these properties we develop a set of metrics that can be directly computed from the SP outputs. We show how the properties are met using these metrics and targeted artificial simulations. We then demonstrate the value of the SP in a complete end-to-end real-world HTM system. We discuss the relationship with neuroscience and previous studies of sparse coding. The HTM spatial pooler represents a neurally inspired algorithm for learning sparse representations from noisy data streams in an online fashion.

Sensory coding accuracy and perceptual performance are improved during the desynchronized cortical state

Cortical activity changes continuously during the course of the day. At a global scale, population activity varies between the ‘synchronized’ state during sleep and ‘desynchronized’ state during waking. However, whether local fluctuations in population synchrony during wakefulness modulate the accuracy of sensory encoding and behavioral performance is poorly understood. Here, we show that populations of cells in monkey visual cortex exhibit rapid fluctuations in synchrony ranging from desynchronized responses, indicative of high alertness, to highly synchronized responses. These fluctuations are local and control the trial variability in population coding accuracy and behavioral performance in a discrimination task. When local population activity is desynchronized, the correlated variability between neurons is reduced, and network and behavioral performance are enhanced. These findings demonstrate that the structure of variability in local cortical populations is not noise but rather controls how sensory information is optimally integrated with ongoing processes to guide network coding and behavior.

Changes in synchrony of cortical populations are observed across the sleep-wake cycle, however the effect of fluctuations in synchrony during wakefulness is not understood. Here the authors show that visual cortical neurons have improved sensory encoding accuracy as well as improved perceptual performance during periods of local population desynchrony.

Transcranial Direct Current Stimulation Facilitates Associative Learning and Alters Functional Connectivity in the Primate Brain

There has been growing interest in transcranial direct current stimulation (tDCS), a non-invasive technique purported to modulate neural activity via weak, externally applied electric fields. Although some promising preliminary data have been reported for applications ranging from stroke rehabilitation to cognitive enhancement, little is known about how tDCS affects the human brain, and some studies have concluded that it may have no effect at all. Here, we describe a macaque model of tDCS that allows us to simultaneously examine the effects of tDCS on brain activity and behavior. We find that applying tDCS to right prefrontal cortex improves monkeys' performance on an associative learning task. While firing rates do not change within the targeted area, tDCS does induce large low-frequency oscillations in the underlying tissue. These oscillations alter functional connectivity, both locally and between distant brain areas, and these long-range changes correlate with tDCS's effects on behavior. Together, these results are consistent with the idea that tDCS leads to widespread changes in brain activity and suggest that it may be a valuable method for cheaply and non-invasively altering functional connectivity in humans.

Dynamic Structure of Neural Variability in the Cortical Representation of Speech Sounds

Accurate sensory discrimination is commonly believed to require precise representations in the nervous system; however, neural stimulus responses can be highly variable, even to identical stimuli. Recent studies suggest that cortical response variability decreases during stimulus processing, but the implications of such effects on stimulus discrimination are unclear. To address this, we examined electrocorticographic cortical field potential recordings from the human nonprimary auditory cortex (superior temporal gyrus) while subjects listened to speech syllables. Compared with a prestimulus baseline, activation variability decreased upon stimulus onset, similar to findings from microelectrode recordings in animal studies. We found that this decrease was simultaneous with encoding and spatially specific for those electrodes that most strongly discriminated speech sounds. We also found that variability was predominantly reduced in a correlated subspace across electrodes. We then compared signal and variability (noise) correlations and found that noise correlations reduce more for electrodes with strong signal correlations. Furthermore, we found that this decrease in variability is strongest in the high gamma band, which correlates with firing rate response. Together, these findings indicate that the structure of single-trial response variability is shaped to enhance discriminability despite non-stimulus-related noise.

SIGNIFICANCE STATEMENT

Cortical responses can be highly variable to auditory speech sounds. Despite this, sensory perception can be remarkably stable. Here, we recorded from the human superior temporal gyrus, a high-order auditory cortex, and studied the changes in the cortical representation of speech stimuli across multiple repetitions. We found that neural variability is reduced upon stimulus onset across electrodes that encode speech sounds.

Variability and Correlations in Primary Visual Cortical Neurons Driven by Fixational Eye Movements

The ability to distinguish between elements of a sensory neuron's activity that are stimulus independent versus driven by the stimulus is critical for addressing many questions in systems neuroscience. This is typically accomplished by measuring neural responses to repeated presentations of identical stimuli and identifying the trial-variable components of the response as noise. In awake primates, however, small "fixational" eye movements (FEMs) introduce uncontrolled trial-to-trial differences in the visual stimulus itself, potentially confounding this distinction. Here, we describe novel analytical methods that directly quantify the stimulus-driven and stimulus-independent components of visual neuron responses in the presence of FEMs. We apply this approach, combined with precise model-based eye tracking, to recordings from primary visual cortex (V1), finding that standard approaches that ignore FEMs typically miss more than half of the stimulus-driven neural response variance, creating substantial biases in measures of response reliability. We show that these effects are likely not isolated to the particular experimental conditions used here, such as the choice of visual stimulus or spike measurement time window, and thus will be a more general problem for V1 recordings in awake primates. We also demonstrate that measurements of the stimulus-driven and stimulus-independent correlations among pairs of V1 neurons can be greatly biased by FEMs. These results thus illustrate the potentially dramatic impact of FEMs on measures of signal and noise in visual neuron activity and also demonstrate a novel approach for controlling for these eye-movement-induced effects.

SIGNIFICANCE STATEMENT

Distinguishing between the signal and noise in a sensory neuron's activity is typically accomplished by measuring neural responses to repeated presentations of an identical stimulus. For recordings from the visual cortex of awake animals, small "fixational" eye movements (FEMs) inevitably introduce trial-to-trial variability in the visual stimulus, potentially confounding such measures. Here, we show that FEMs often have a dramatic impact on several important measures of response variability for neurons in primary visual cortex. We also present an analytical approach for quantifying signal and noise in visual neuron activity in the presence of FEMs. These results thus highlight the importance of controlling for FEMs in studies of visual neuron function, and demonstrate novel methods for doing so.

Revealing unobserved factors underlying cortical activity with a rectified latent variable model applied to neural population recordings

The activity of sensory cortical neurons is not only driven by external stimuli but also shaped by other sources of input to the cortex. Unlike external stimuli, these other sources of input are challenging to experimentally control, or even observe, and as a result contribute to variability of neural responses to sensory stimuli. However, such sources of input are likely not "noise" and may play an integral role in sensory cortex function. Here we introduce the rectified latent variable model (RLVM) in order to identify these sources of input using simultaneously recorded cortical neuron populations. The RLVM is novel in that it employs nonnegative (rectified) latent variables and is much less restrictive in the mathematical constraints on solutions because of the use of an autoencoder neural network to initialize model parameters. We show that the RLVM outperforms principal component analysis, factor analysis, and independent component analysis, using simulated data across a range of conditions. We then apply this model to two-photon imaging of hundreds of simultaneously recorded neurons in mouse primary somatosensory cortex during a tactile discrimination task. Across many experiments, the RLVM identifies latent variables related to both the tactile stimulation as well as nonstimulus aspects of the behavioral task, with a majority of activity explained by the latter. These results suggest that properly identifying such latent variables is necessary for a full understanding of sensory cortical function and demonstrate novel methods for leveraging large population recordings to this end.NEW & NOTEWORTHY The rapid development of neural recording technologies presents new opportunities for understanding patterns of activity across neural populations. Here we show how a latent variable model with appropriate nonlinear form can be used to identify sources of input to a neural population and infer their time courses. Furthermore, we demonstrate how these sources are related to behavioral contexts outside of direct experimental control.

Neurons in Primate Visual Cortex Alternate between Responses to Multiple Stimuli in Their Receptive Field

A fundamental question concerning representation of the visual world in our brain is how a cortical cell responds when presented with more than a single stimulus. We find supportive evidence that most cells presented with a pair of stimuli respond predominantly to one stimulus at a time, rather than a weighted average response. Traditionally, the firing rate is assumed to be a weighted average of the firing rates to the individual stimuli (response-averaging model) (Bundesen et al., 2005). Here, we also evaluate a probability-mixing model (Bundesen et al., 2005), where neurons temporally multiplex the responses to the individual stimuli. This provides a mechanism by which the representational identity of multiple stimuli in complex visual scenes can be maintained despite the large receptive fields in higher extrastriate visual cortex in primates. We compare the two models through analysis of data from single cells in the middle temporal visual area (MT) of rhesus monkeys when presented with two separate stimuli inside their receptive field with attention directed to one of the two stimuli or outside the receptive field. The spike trains were modeled by stochastic point processes, including memory effects of past spikes and attentional effects, and statistical model selection between the two models was performed by information theoretic measures as well as the predictive accuracy of the models. As an auxiliary measure, we also tested for uni- or multimodality in interspike interval distributions, and performed a correlation analysis of simultaneously recorded pairs of neurons, to evaluate population behavior.

The Population Tracking Model: A Simple, Scalable Statistical Model for Neural Population Data

Our understanding of neural population coding has been limited by a lack of analysis methods to characterize spiking data from large populations. The biggest challenge comes from the fact that the number of possible network activity patterns scales exponentially with the number of neurons recorded ([Formula: see text]). Here we introduce a new statistical method for characterizing neural population activity that requires semi-independent fitting of only as many parameters as the square of the number of neurons, requiring drastically smaller data sets and minimal computation time. The model works by matching the population rate (the number of neurons synchronously active) and the probability that each individual neuron fires given the population rate. We found that this model can accurately fit synthetic data from up to 1000 neurons. We also found that the model could rapidly decode visual stimuli from neural population data from macaque primary visual cortex about 65 ms after stimulus onset. Finally, we used the model to estimate the entropy of neural population activity in developing mouse somatosensory cortex and, surprisingly, found that it first increases, and then decreases during development. This statistical model opens new options for interrogating neural population data and can bolster the use of modern large-scale in vivo Ca[Formula: see text] and voltage imaging tools.

Artefactual origin of biphasic cortical spike-LFP correlation

Electrophysiological data acquisition systems introduce various distortions into the signals they record. While such distortions were discussed previously, their effects are often not appreciated. Here I show that the biphasic shape of cortical spike-triggered LFP average (stLFP), reported in multiple studies, is likely an artefact introduced by high-pass filter of the neural data acquisition system when the actual stLFP has a single trough around the zero lag.