Local Spiking Activity Research Articles

Flexible modeling of large-scale neural network stimulation: electrical and optical extensions to The Virtual Electrode Recording Tool for EXtracellular Potentials (VERTEX).

Computational models that predict effects of neural stimulation can be used as a preliminary tool to inform in-vivo research, reducing the costs, time, and ethical considerations involved. However, current models do not support the diverse neural stimulation techniques used in-vivo, including the expanding selection of electrodes, stimulation modalities, and stimulation paradigms. To develop a more comprehensive software, we created several extensions to The Virtual Electrode Recording Tool for EXtracellular Potentials (VERTEX), the MATLAB-based neural stimulation tool from Newcastle University. VERTEX simulates input currents in a large population of multi-compartment neurons within a small cortical slice to model electric field stimulation, while recording local field potentials (LFPs) and spiking activity. Our extensions to its existing electric field stimulation framework include multiple pairs of parametrically defined electrodes and biphasic, bipolar stimulation delivered at programmable delays. To support the growing use of optogenetic approaches for targeted neural stimulation, we introduced a feature that models optogenetic stimulation through an additional VERTEX input function that converts irradiance to currents at optogenetically responsive neurons. Finally, we added extensions to allow complex stimulation protocols including paired-pulse, spatiotemporal patterned, and closed-loop stimulation. We demonstrated our novel features using VERTEX's built-in functionalities, illustrating how these extensions can be used to efficiently and systematically test diverse, targeted, and individualized stimulation patterns.

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.

Correlation between desynchrony of hippocampal neural activity and hyperlocomotion in the model mice of schizophrenia and therapeutic effects of aripiprazole.

The hippocampus has been reported to be morphologically and neurochemically altered in schizophrenia (SZ). Hyperlocomotion is a characteristic SZ-associated behavioral phenotype, which is associated with dysregulated dopamine system function induced by hippocampal hyperactivity. However, the neural mechanism of hippocampus underlying hyperlocomotion remains largely unclear. Mouse pups were injected with N-methyl-D-aspartate receptor antagonist (MK-801) or vehicle twice daily on postnatal days (PND) 7-11. In the adulthood phase, one cohort of mice underwent electrode implantation in field CA1 of the hippocampus for the recording local field potentials and spike activity. A separate cohort of mice underwent surgery to allow for calcium imaging of the hippocampus while monitoring the locomotion. Lastly, the effects of atypical antipsychotic (aripiprazole, ARI) were evaluated on hippocampal neural activity. We found that the hippocampal theta oscillations were enhanced in MK-801-treated mice, but the correlation coefficient between the hippocampal spiking activity and theta oscillation was reduced. Consistently, although the rate and amplitude of calcium transients of hippocampal neurons were increased, their synchrony and correlation to locomotion speed were disrupted. ARI ameliorated perturbations produced by the postnatal MK-801 treatment. These results suggest that the disruption of neural coordination may underly the neuropathological mechanism for hyperlocomotion of SZ.

Sequentially activated discrete modules appear as traveling waves in neuronal measurements with limited spatiotemporal sampling.

Numerous studies have identified traveling waves in the cortex and suggested they play important roles in brain processing. These waves are most often measured using macroscopic methods that are unable to assess the local spiking activity underlying wave dynamics. Here, we investigated the possibility that waves may not be traveling at the single neuron scale. We first show that sequentially activating two discrete brain areas can appear as traveling waves in EEG simulations. We next reproduce these results using an analytical model of two sequentially activated regions. Using this model, we were able to generate wave-like activity with variable directions, velocities, and spatial patterns, and to map the discriminability limits between traveling waves and modular sequential activations. Finally, we investigated the link between field potentials and single neuron excitability using large-scale measurements from turtle cortex ex vivo. We found that while field potentials exhibit wave-like dynamics, the underlying spiking activity was better described by consecutively activated spatially adjacent groups of neurons. Taken together, this study suggests caution when interpreting phase delay measurements as continuously propagating wavefronts in two different spatial scales. A careful distinction between modular and wave excitability profiles across scales will be critical for understanding the nature of cortical computations.

Sequentially activated discrete modules appear as traveling waves in neuronal measurements with limited spatiotemporal sampling

Numerous studies have identified traveling waves in the cortex and suggested they play important roles in brain processing. These waves are most often measured using macroscopic methods that are unable to assess the local spiking activity underlying wave dynamics. Here, we investigated the possibility that waves may not be traveling at the single neuron scale. We first show that sequentially activating two discrete brain areas can appear as traveling waves in EEG simulations. We next reproduce these results using an analytical model of two sequentially activated regions. Using this model, we were able to generate wave-like activity with variable directions, velocities, and spatial patterns, and to map the discriminability limits between traveling waves and modular sequential activations. Finally, we investigated the link between field potentials and single neuron excitability using large-scale measurements from turtle cortex ex vivo. We found that while field potentials exhibit wave-like dynamics, the underlying spiking activity was better described by consecutively activated spatially adjacent groups of neurons. Taken together, this study suggests caution when interpreting phase delay measurements as continuously propagating wavefronts in two different spatial scales. A careful distinction between modular and wave excitability profiles across scales will be critical for understanding the nature of cortical computations.

Local field potentials, spiking activity, and receptive fields in human visual cortex.

The concept of receptive field (RF) is central to sensory neuroscience. Neuronal RF properties have been substantially studied in animals, while those in humans remain nearly unexplored. Here, we measured neuronal RFs with intracranial local field potentials (LFPs) and spiking activity in human visual cortex (V1/V2/V3). We recorded LFPs via macro-contacts and discovered that RF sizes estimated from low-frequency activity (LFA, 0.5-30 Hz) were larger than those estimated from low-gamma activity (LGA, 30-60 Hz) and high-gamma activity (HGA, 60-150 Hz). We then took a rare opportunity to record LFPs and spiking activity via microwires in V1 simultaneously. We found that RF sizes and temporal profiles measured from LGA and HGA closely matched those from spiking activity. In sum, this study reveals that spiking activity of neurons in human visual cortex could be well approximated by LGA and HGA in RF estimation and temporal profile measurement, implying the pivotal functions of LGA and HGA in early visual information processing.

Rapid Input-Output Transformation between Local Field Potential and Spiking Activity during Sensation but not Action in the Superior Colliculus.

Sensorimotor transformation is the sequential process of registering a sensory signal in the environment and then responding with the relevant movement at an appropriate time. For visually guided eye movements, neural signatures in the form of spiking activity of neurons have been extensively studied along the dorsoventral axis of the superior colliculus (SC). In contrast, the local field potential (LFP), which represents the putative input to a region, remains largely unexplored in the SC. We therefore compared amplitude levels and onset times of both spike bursts and LFP modulations recorded simultaneously with a laminar probe along the dorsoventral axis of SC in 3 male monkeys performing the visually guided delayed saccade task. Both signals displayed a gradual transition from sensory activity in the superficial layers to a predominantly motor response in the deeper layers, although the transition from principally sensory to mostly motor response occurred ∼500 μm deeper for the LFP. For the sensory response, LFP modulation preceded spike burst onset by <5 ms in the superficial and intermediate layers and only when data were analyzed on a trial-by-trial basis. The motor burst in the spiking activity led LFP modulation by >25 ms in the deeper layers. The results reveal a fast and efficient input-output transformation between LFP modulation and spike burst in the visually responsive layers activity during sensation but not during action. The spiking pattern observed during the movement phase is likely dominated by intracollicular processing that is not captured in the LFP.SIGNIFICANCE STATEMENT What is the sequence of events between local field potential (LFP) modulation and spiking activity during sensorimotor transformation? A trial-by-trial analysis reveals that the LFP activity leads the spike burst in the superficial and intermediate layers of the superior colliculus during visual processing, while both trial-by-trial and the average analyses show that the spike burst leads the LFP modulation during movement generation. These results suggest an almost instantaneous LFP input, spike burst output transformation in the visually responsive layers of the superior colliculus when registering the stimulus. In contrast, substantial intracollicular processing likely results in a saccade-related spike burst that leads LFP modulation.

Hippocampal 4-Hz oscillations emerge during stationary running in a wheel and are resistant to medial septum inactivation.

Recent studies described 2-4 Hz oscillations in the hippocampus of rats performing stationary locomotion on treadmills and other apparatus. Since the 2-4 Hz rhythm shares common features with theta (5-12 Hz) oscillations-such as a positive amplitude-running speed relationship and modulation of spiking activity-many have questioned whether these rhythms are related or independently generated. Here, we analyzed local field potentials and spiking activity from the dorsal CA1 of rats executing a spatial alternation task and running for ~15 s in a wheel during the intertrial intervals both before and after muscimol injection into the medial septum. We observed remarkable 4-Hz oscillations during wheel runs, which presented amplitude positively correlated with running speed. Surprisingly, the amplitude of 4-Hz and theta oscillations were inversely related. Medial septum inactivation abolished hippocampal theta but preserved 4-Hz oscillations. It also affected the entrainment of pyramidal cells and interneurons by 4-Hz rhythmic activity. In all, these results dissociate the underlying mechanism of 4-Hz and theta oscillations in the rat hippocampus.

Intrinsic timescales in the visual cortex change with selective attention and reflect spatial connectivity

Intrinsic timescales characterize dynamics of endogenous fluctuations in neural activity. Variation of intrinsic timescales across the neocortex reflects functional specialization of cortical areas, but less is known about how intrinsic timescales change during cognitive tasks. We measured intrinsic timescales of local spiking activity within columns of area V4 in male monkeys performing spatial attention tasks. The ongoing spiking activity unfolded across at least two distinct timescales, fast and slow. The slow timescale increased when monkeys attended to the receptive fields location and correlated with reaction times. By evaluating predictions of several network models, we found that spatiotemporal correlations in V4 activity were best explained by the model in which multiple timescales arise from recurrent interactions shaped by spatially arranged connectivity, and attentional modulation of timescales results from an increase in the efficacy of recurrent interactions. Our results suggest that multiple timescales may arise from the spatial connectivity in the visual cortex and flexibly change with the cognitive state due to dynamic effective interactions between neurons.

Mechanisms and functions of respiration-driven gamma oscillations in the primary olfactory cortex.

Gamma oscillations are believed to underlie cognitive processes by shaping the formation of transient neuronal partnerships on a millisecond scale. These oscillations are coupled to the phase of breathing cycles in several brain areas, possibly reflecting local computations driven by sensory inputs sampled at each breath. Here, we investigated the mechanisms and functions of gamma oscillations in the piriform (olfactory) cortex of awake mice to understand their dependence on breathing and how they relate to local spiking activity. Mechanistically, we find that respiration drives gamma oscillations in the piriform cortex, which correlate with local feedback inhibition and result from recurrent connections between local excitatory and inhibitory neuronal populations. Moreover, respiration-driven gamma oscillations are triggered by the activation of mitral/tufted cells in the olfactory bulb and are abolished during ketamine/xylazine anesthesia. Functionally, we demonstrate that they locally segregate neuronal assemblies through a winner-take-all computation leading to sparse odor coding during each breathing cycle. Our results shed new light on the mechanisms of gamma oscillations, bridging computation, cognition, and physiology.

Emergent Low-Frequency Activity in Cortico-Cerebellar Networks with Motor Skill Learning.

The motor cortex controls skilled arm movement by recruiting a variety of targets in the nervous system, and it is important to understand the emergent activity in these regions as refinement of a motor skill occurs. One fundamental projection of the motor cortex (M1) is to the cerebellum. However, the emergent activity in the motor cortex and the cerebellum that appears as a dexterous motor skill is consolidated is incompletely understood. Here, we report on low-frequency oscillatory (LFO) activity that emerges in cortico-cerebellar networks with learning the reach-to-grasp motor skill. We chronically recorded the motor and the cerebellar cortices in rats, which revealed the emergence of coordinated movement-related activity in the local-field potentials as the reaching skill consolidated. Interestingly, we found this emergent activity only in the rats that gained expertise in the task. We found that the local and cross-area spiking activity was coordinated with LFOs in proficient rats. Finally, we also found that these neural dynamics were more prominently expressed during accurate behavior in the M1. This work furthers our understanding on emergent dynamics in the cortico-cerebellar loop that underlie learning and execution of precise skilled movement.

Visually induced γ band rhythm in spatial summation beyond the receptive field in mouse primary visual cortex.

Recent studies in many kinds of mammals have established the existence of multiple γ rhythms in the cerebral cortex subserving different functions. In the primary visual cortex (V1), visually induced γ rhythms are dependent on stimulus features. However, experimental findings of γ power induced by varying the size of the drifting grating are inconsistent. Here, we reinvestigated the spatial summation properties of visually induced spike and γ rhythm activities in mouse V1. Our results show that drifting sinusoidal grating stimuli mainly induce 2 γ band rhythms, including a low-frequency band (25-45Hz) and a high-frequency band (55-75Hz). Unlike previous findings, we discovered that visually induced γ power could also exhibit extrareceptive field (ERF) modulatory properties. The modulation by ERF stimulation could be either suppressive, countersuppressive, or nonsuppressive, mostly similar to the local spike activity. Moreover, further analysis of the neuron group exhibiting surround suppression in both spike and γ activity revealed that the strength of the surround suppression and the receptive field size showed moderate correlations between measurements by spike and γ rhythm activity. Our findings improve the understanding of the characteristics and neural mechanisms of induced γ rhythms in visual spatial summation.

Psilocybin reduces low frequency oscillatory power and neuronal phase-locking in the anterior cingulate cortex of awake rodents

Psilocybin is a hallucinogenic compound that is showing promise in the ability to treat neurological conditions such as depression and post-traumatic stress disorder. There have been several investigations into the neural correlates of psilocybin administration using non-invasive methods, however, there has yet to be an invasive study of the mechanism of action in awake rodents. Using multi-unit extracellular recordings, we recorded local field potential and spiking activity from populations of neurons in the anterior cingulate cortex of awake mice during the administration of psilocybin (2 mg/kg). The power of low frequency bands in the local field potential was found to significantly decrease in response to psilocybin administration, whilst gamma band activity trended towards an increase. The population firing rate was found to increase overall, with just under half of individual neurons showing a significant increase. Psilocybin significantly decreased the level of phase modulation of cells with each neural frequency band except high-gamma oscillations, consistent with a desynchronization of cortical populations. Furthermore, bursting behavior was altered in a subset of cells, with both positive and negative changes in the rate of bursting. Neurons that increased their burst firing following psilocybin administration were highly likely to transition from a phase-modulated to a phase unmodulated state. Taken together, psilocybin reduces low frequency oscillatory power, increases overall firing rates and desynchronizes local neural activity. These findings are consistent with dissolution of the default mode network under psilocybin, and may be indicative of disruption of top-down processing in the acute psychedelic state.

Hippocampal–prefrontal interactions during spatial decision‐making

AbstractThe hippocampus has been linked to memory encoding and spatial navigation, while the prefrontal cortex is associated with cognitive functions such as decision‐making. These regions are hypothesized to communicate in tasks that demand both spatial navigation and decision‐making processes. However, the electrophysiological signatures underlying this communication remain to be better elucidated. To investigate the dynamics of the hippocampal–prefrontal interactions, we have analyzed their local field potentials and spiking activity recorded from rats performing a spatial alternation task on a figure eight‐shaped maze. We found that the phase coherence of theta peaked around the choice point area of the maze. Moreover, Granger causality revealed a hippocampus → prefrontal cortex directionality of information flow at theta frequency, peaking at starting areas of the maze, and on the reverse direction at delta frequency, peaking near the turn onset. Additionally, the patterns of phase‐amplitude cross‐frequency coupling within and between the regions also showed spatial selectivity, and hippocampal theta and prefrontal delta modulated not only gamma amplitude but also inter‐regional gamma synchrony. Finally, we found that the theta rhythm dynamically modulated neurons in both regions, with the highest modulation at the choice area; interestingly, prefrontal cortex neurons were more strongly modulated by the hippocampal theta rhythm than by their local field rhythm. In all, our results reveal maximum electrophysiological interactions between the hippocampus and the prefrontal cortex near the decision‐making period of the spatial alternation task, corroborating the hypothesis that a dynamic interplay between these regions takes place during spatial decisions.

Self-Organized Criticality in the Brain

Self-organized criticality (SOC) refers to the ability of complex systems to evolve toward a second-order phase transition at which interactions between system components lead to scale-invariant events that are beneficial for system performance. For the last two decades, considerable experimental evidence has accumulated that the mammalian cortex with its diversity in cell types, interconnectivity, and plasticity might exhibit SOC. Here, we review the experimental findings of isolated, layered cortex preparations to self-organize toward four dynamical motifs presently identified in the intact cortex in vivo: up-states, oscillations, neuronal avalanches, and coherence potentials. During up-states, the synchronization observed for nested theta/gamma oscillations embeds scale-invariant neuronal avalanches, which can be identified by robust power law scaling in avalanche sizes with a slope of −3/2 and a critical branching parameter of 1. This precise dynamical coordination, tracked in the negative transients of the local field potential (nLFP) and spiking activity of pyramidal neurons using two-photon imaging, emerges autonomously in superficial layers of organotypic cortex cultures and acute cortex slices, is homeostatically regulated, exhibits separation of time scales, and reveals unique size vs. quiet time dependencies. A subclass of avalanches, the coherence potentials, exhibits precise maintenance of the time course in propagated local synchrony. Avalanches emerge in superficial layers of the cortex under conditions of strong external drive. The balance of excitation and inhibition (E/I), as well as neuromodulators such as dopamine, establishes powerful control parameters for avalanche dynamics. This rich dynamical repertoire is not observed in dissociated cortex cultures, which lack the differentiation into cortical layers and exhibit a dynamical phenotype expected for a first-order phase transition. The precise interactions between up-states, nested oscillations, and avalanches in superficial layers of the cortex provide compelling evidence for SOC in the brain.

Frequency of theta rhythm is controlled by acceleration, but not speed, in running rats

SummaryThe theta rhythm organizes neural activity across hippocampus and entorhinal cortex. A role for theta oscillations in spatial navigation is supported by half a century of research reporting that theta frequency encodes running speed linearly so that displacement can be estimated through theta frequency integration. We show that this relationship is an artifact caused by the fact that the speed of freely moving animals could not be systematically disentangled from acceleration. Using an experimental procedure that clamps running speed at pre-set values, we find that the theta frequency of local field potentials and spike activity is linearly related to positive acceleration, but not negative acceleration or speed. The modulation by positive-only acceleration makes rhythmic activity at theta frequency unfit as a code to compute displacement or any other kinematic variable. Temporally precise variations in theta frequency may instead serve as a mechanism for speeding up entorhinal-hippocampal computations during accelerated movement.

Theta Synchrony Is Increased near Neural Populations That Are Active When Initiating Instructed Movement.

Theta oscillations (3–8 Hz) in the human brain have been linked to perception, cognitive control, and spatial memory, but their relation to the motor system is less clear. We tested the hypothesis that theta oscillations coordinate distributed behaviorally relevant neural representations during movement using intracranial electroencephalography (iEEG) recordings from nine patients (n = 490 electrodes) as they performed a simple instructed movement task. Using high frequency activity (HFA; 70–200 Hz) as a marker of local spiking activity, we identified electrodes that were positioned near neural populations that showed increased activity during instruction and movement. We found that theta synchrony was widespread throughout the brain but was increased near regions that showed movement-related increases in neural activity. These results support the view that theta oscillations represent a general property of brain activity that may also play a specific role in coordinating widespread neural activity when initiating voluntary movement.

Blockade of Retinal Oscillations by Benzodiazepines Improves Efficiency of Electrical Stimulation in the Mouse Model of RP, rd10.

PurposeIn RP, photoreceptors degenerate. Retinal prostheses are considered a suitable strategy to restore vision. In animal models of RP, a pathologic rhythmic activity seems to compromise the efficiency of retinal ganglion cell stimulation by an electrical prosthesis. We, therefore, strove to eliminate this pathologic activity.MethodsElectrophysiologic recordings of local field potentials and spike activity of retinal ganglion cells were obtained in vitro from retinae of wild-type and rd10 mice using multielectrode arrays. Retinae were stimulated electrically.ResultsThe efficiency of electrical stimulation was lower in rd10 retina than in wild-type retina and this was highly correlated with the presence of oscillations in retinal activity. Glycine and GABA, as well as the benzodiazepines diazepam, lorazepam, and flunitrazepam, abolished retinal oscillations and, most important, increased the efficiency of electrical stimulation to values similar to those in wild-type retina.ConclusionsTreatment of patients with these benzodiazepines may offer a way to improve the performance of retinal implants in cases with poor implant proficiency. This study may open the way to a therapy that supports electrical stimulation by prostheses with pharmacologic treatment.

Adaptation Modulates Spike-Phase Coupling Tuning Curve in the Rat Primary Auditory Cortex.

Adaptation is an important mechanism that causes a decrease in the neural response both in terms of local field potentials (LFP) and spiking activity. We previously showed this reduction effect in the tuning curve of the primary auditory cortex. Moreover, we revealed that a repeated stimulus reduces the neural response in terms of spike-phase coupling (SPC). In the current study, we examined the effect of adaptation on the SPC tuning curve. To this end, employing the phase-locking value (PLV) method, we estimated the spike-LFP coupling. The data was obtained by a simultaneous recording from four single-electrodes in the primary auditory cortex of 15 rats. We first investigated whether the neural system may use spike-LFP phase coupling in the primary auditory cortex to encode sensory information. Secondly, we investigated the effect of adaptation on this potential SPC tuning. Our data showed that the coupling between spikes’ times and the LFP phase in beta oscillations represents sensory information (different stimulus frequencies), with an inverted bell-shaped tuning curve. Furthermore, we showed that adaptation to a specific frequency modulates SPC tuning curve of the adapter and its neighboring frequencies. These findings could be useful for interpretation of feature representation in terms of SPC and the underlying neural mechanism of adaptation.

Aberrant Auditory Steady-State Response of Awake Mice After Single Application of the NMDA Receptor Antagonist MK-801 Into the Medial Geniculate Body

BackgroundSystemic administration of noncompetitive N-methyl-D-aspartate receptor (NMDAR) antagonists such as MK-801 is widely used to model psychosis of schizophrenia (SZ). Acute systemic MK-801 in rodents caused an increase of the auditory steady-state responses (ASSRs), the oscillatory neural responses to periodic auditory stimulation, while most studies in patients with SZ reported a decrease of ASSRs. This inconsistency may be attributable to the comprehensive effects of systemic administration of MK-801. Here, we examined how the ASSR is affected by selectively blocking NMDAR in the thalamus.MethodsWe implanted multiple electrodes in the auditory cortex (AC) and prefrontal cortex to simultaneously record the local field potential and spike activity (SA) of multiple sites from awake mice. Click-trains at a 40-Hz repetition rate were used to evoke the ASSR. We compared the mean trial power and phase-locking factor and the firing rate of SA before and after microinjection of MK-801 (1.5 µg) into the medial geniculate body (MGB).ResultsWe found that both the AC and prefrontal cortex showed a transient local field potential response at the onset of click-train stimulus, which was less affected by the application of MK-801 in the MGB. Following the onset response, the AC also showed a response continuing throughout the stimulus period, corresponding to the ASSR, which was suppressed by the application of MK-801.ConclusionOur data suggest that the MGB is one of the generators of ASSR, and NMDAR hypofunction in the thalamocortical projection may account for the ASSR deficits in SZ.