Sleep restores an optimal computational regime in cortical networks.

Sleep is assumed to subserve homeostatic processes in the brain; however, the set point around which sleep tunes circuit computations is unknown. Slow-wave activity (SWA) is commonly used to reflect the homeostatic aspect of sleep; although it can indicate sleep pressure, it does not explain why animals need sleep. This study aimed to assess whether criticality may be the computational set point of sleep. By recording cortical neuron activity continuously for 10-14 d in freely behaving rats, we show that normal waking experience progressively disrupts criticality and that sleep functions to restore critical dynamics. Criticality is perturbed in a context-dependent manner, and waking experience is causal in driving these effects. The degree of deviation from criticality predicts future sleep/wake behavior more accurately than SWA, behavioral history or other neural measures. Our results demonstrate that perturbation and recovery of criticality is a network homeostatic mechanism consistent with the core, restorative function of sleep.

Tauopathy severely disrupts homeostatic set-points in emergent neural dynamics but not in the activity of individual neurons.

The homeostatic regulation of neuronal activity is essential for robust computation; key set-points, such as firing rate, are actively stabilized to compensate for perturbations. From this perspective, the disruption of brain function central to neurodegenerative disease should reflect impairments of computationally essential set-points. Despite connecting neurodegeneration to functional outcomes, the impact of disease on set-points in neuronal activity is unknown. Here we present a comprehensive, theory-driven investigation of the effects of tau-mediated neurodegeneration on homeostatic set-points in neuronal activity. In a mouse model of tauopathy, we examine 27,000 hours of hippocampal recordings during free behavior throughout disease progression. Contrary to our initial hypothesis that tauopathy would impact set-points in spike rate and variance, we found that cell-level set-points are resilient to even the latest stages of disease. Instead, we find that tauopathy disrupts neuronal activity at the network-level, which we quantify using both pairwise measures of neuron interactions as well as measurement of the network's nearness to criticality, an ideal computational regime that is known to be a homeostatic set-point. We find that shifts in network criticality 1) track with symptoms, 2) predict underlying anatomical and molecular pathology, 3) occur in a sleep/wake dependent manner, and 4) can be used to reliably classify an animal's genotype. Our data suggest that the critical set-point is intact, but that homeostatic machinery is progressively incapable of stabilizing hippocampal networks, particularly during waking. This work illustrates how neurodegenerative processes can impact the computational capacity of neurobiological systems, and suggest an important connection between molecular pathology, circuit function, and animal behavior.

Stable representation of a naturalistic movie emerges from episodic activity with gain variability.

Visual cortical responses are known to be highly variable across trials within an experimental session. However, the long-term stability of visual cortical responses is poorly understood. Here using chronic imaging of V1 in mice we show that neural responses to repeated natural movie clips are unstable across weeks. Individual neuronal responses consist of sparse episodic activity which are stable in time but unstable in gain across weeks. Further, we find that the individual episode, instead of neuron, serves as the basic unit of the week-to-week fluctuation. To investigate how population activity encodes the stimulus, we extract a stable one-dimensional representation of the time in the natural movie, using an unsupervised method. Most week-to-week fluctuation is perpendicular to the stimulus encoding direction, thus leaving the stimulus representation largely unaffected. We propose that precise episodic activity with coordinated gain changes are keys to maintain a stable stimulus representation in V1.

Precision multidimensional neural population code recovered from single intracellular recordings.

Neurons in sensory cortices are more naturally and deeply integrated than any current neural population recording tools (e.g. electrode arrays, fluorescence imaging). Two concepts facilitate efforts to observe population neural code with single-cell recordings. First, even the highest quality single-cell recording studies find a fraction of the stimulus information in high-dimensional population recordings. Finding any of this missing information provides proof of principle. Second, neurons and neural populations are understood as coupled nonlinear differential equations. Therefore, fitted ordinary differential equations provide a basis for single-trial single-cell stimulus decoding. We obtained intracellular recordings of fluctuating transmembrane current and potential in mouse visual cortex during stimulation with drifting gratings. We use mean deflection from baseline when comparing to prior single-cell studies because action potentials are too sparse and the deflection response to drifting grating stimuli (e.g. tuning curves) are well studied. Equation-based decoders allowed more precise single-trial stimulus discrimination than tuning-curve-base decoders. Performance varied across recorded signal types in a manner consistent with population recording studies and both classification bases evinced distinct stimulus-evoked phases of population dynamics, providing further corroboration. Naturally and deeply integrated observations of population dynamics would be invaluable. We offer proof of principle and a versatile framework.

Stability of motor cortex network states during learning-associated neural reorganizations.

A substantial reorganization of neural activity and neuron-to-movement relationship in motor cortical circuits accompanies the emergence of reproducible movement patterns during motor learning. Little is known about how this tempest of neural activity restructuring impacts the stability of network states in recurrent cortical circuits. To investigate this issue, we reanalyzed data in which we recorded for 14 days via population calcium imaging the activity of the same neural populations of pyramidal neurons in layer 2/3 and layer 5 of forelimb motor and premotor cortex in mice during the daily learning of a lever-press task. We found that motor cortex network states remained stable with respect to the critical network state during the extensive reorganization of both neural population activity and its relation to lever movement throughout learning. Specifically, layer 2/3 cortical circuits unceasingly displayed robust evidence for operating at the critical network state, a regime that maximizes information capacity and transmission and provides a balance between network robustness and flexibility. In contrast, layer 5 circuits operated away from the critical network state for all 14 days of recording and learning. In conclusion, this result indicates that the wide-ranging malleability of synapses, neurons, and neural connectivity during learning operates within the constraint of a stable and layer-specific network state regarding dynamic criticality, and suggests that different cortical layers operate under distinct constraints because of their specialized goals.NEW & NOTEWORTHY The neural activity reorganizes throughout motor learning, but how this reorganization impacts the stability of network states is unclear. We used two-photon calcium imaging to investigate how the network states in layer 2/3 and layer 5 of forelimb motor and premotor cortex are modulated by motor learning. We show that motor cortex network states are layer-specific and constant regarding criticality during neural activity reorganization, and suggests that layer-specific constraints could be motivated by different functions.

Diverse coactive neurons encode stimulus-driven and stimulus-independent variables.

Both experimenter-controlled stimuli and stimulus-independent variables impact cortical neural activity. A major hurdle to understanding neural representation is distinguishing between qualitatively different causes of the fluctuating population activity. We applied an unsupervised low-rank tensor decomposition analysis to the recorded population activity in the visual cortex of awake mice in response to repeated presentations of naturalistic visual stimuli. We found that neurons covaried largely independently of individual neuron stimulus response reliability and thus encoded both stimulus-driven and stimulus-independent variables. Importantly, a neuron's response reliability and the neuronal coactivation patterns substantially reorganized for different external visual inputs. Analysis of recurrent balanced neural network models revealed that both the stimulus specificity and the mixed encoding of qualitatively different variables can arise from clustered external inputs. These results establish that coactive neurons with diverse response reliability mediate a mixed representation of stimulus-driven and stimulus-independent variables in the visual cortex.NEW & NOTEWORTHY V1 neurons covary largely independently of individual neuron's response reliability. A single neuron's response reliability imposes only a weak constraint on its encoding capabilities. Visual stimulus instructs a neuron's reliability and coactivation pattern. Network models revealed using clustered external inputs.

Pre-Synaptic Pool Modification (PSPM): A supervised learning procedure for recurrent spiking neural networks.

Learning synaptic weights of spiking neural network (SNN) models that can reproduce target spike trains from provided neural firing data is a central problem in computational neuroscience and spike-based computing. The discovery of the optimal weight values can be posed as a supervised learning task wherein the weights of the model network are chosen to maximize the similarity between the target spike trains and the model outputs. It is still largely unknown whether optimizing spike train similarity of highly recurrent SNNs produces weight matrices similar to those of the ground truth model. To this end, we propose flexible heuristic supervised learning rules, termed Pre-Synaptic Pool Modification (PSPM), that rely on stochastic weight updates in order to produce spikes within a short window of the desired times and eliminate spikes outside of this window. PSPM improves spike train similarity for all-to-all SNNs and makes no assumption about the post-synaptic potential of the neurons or the structure of the network since no gradients are required. We test whether optimizing for spike train similarity entails the discovery of accurate weights and explore the relative contributions of local and homeostatic weight updates. Although PSPM improves similarity between spike trains, the learned weights often differ from the weights of the ground truth model, implying that connectome inference from spike data may require additional constraints on connectivity statistics. We also find that spike train similarity is sensitive to local updates, but other measures of network activity such as avalanche distributions, can be learned through synaptic homeostasis.

Dynamics and sources of response variability and its coordination in visual cortex.

The trial-to-trial response variability in sensory cortices and the extent to which this variability can be coordinated among cortical units have strong implications for cortical signal processing. Yet, little is known about the relative contributions and dynamics of defined sources to the cortical response variability and their correlations across cortical units. To fill this knowledge gap, here we obtained and analyzed multisite local field potential (LFP) recordings from visual cortex of turtles in response to repeated naturalistic movie clips and decomposed cortical across-trial LFP response variability into three defined sources, namely, input, network, and local fluctuations. We found that input fluctuations dominate cortical response variability immediately following stimulus onset, whereas network fluctuations dominate the response variability in the steady state during continued visual stimulation. Concurrently, we found that the network fluctuations dominate the correlations of the variability during the ongoing and steady-state epochs, but not immediately following stimulus onset. Furthermore, simulations of various model networks indicated that (i) synaptic time constants, leading to oscillatory activity, and (ii) synaptic clustering and synaptic depression, leading to spatially constrained pockets of coherent activity, are both essential features of cortical circuits to mediate the observed relative contributions and dynamics of input, network, and local fluctuations to the cortical LFP response variability and their correlations across recording sites. In conclusion, these results show how a mélange of multiscale thalamocortical circuit features mediate a complex stimulus-modulated cortical activity that, when naively related to the visual stimulus alone, appears disguised as high and coordinated across-trial response variability.

Cortical Circuit Dynamics Are Homeostatically Tuned to Criticality In Vivo.

Homeostatic mechanisms stabilize neuronal activity in vivo, but whether this process gives rise to balanced network dynamics is unknown. Here, we continuously monitored the statistics of network spiking in visual cortical circuits in freely behaving rats for 9 days. Under control conditions in light and dark, networks were robustly organized around criticality, a regime that maximizes information capacity and transmission. When input was perturbed by visual deprivation, network criticality was severely disrupted and subsequently restored to criticality over 48 h. Unexpectedly, the recovery of excitatory dynamics preceded homeostatic plasticity of firing rates by >30 h. We utilized model investigations to manipulate firing rate homeostasis in a cell-type-specific manner at the onset of visual deprivation. Our results suggest that criticality in excitatory networks is established by inhibitory plasticity and architecture. These data establish that criticality is consistent with a homeostatic set point for visual cortical dynamics and suggest a key role for homeostatic regulation of inhibition.

Single-Cell Membrane Potential Fluctuations Evince Network Scale-Freeness and Quasicriticality.

What information single neurons receive about general neural circuit activity is a fundamental question for neuroscience. Somatic membrane potential (Vm) fluctuations are driven by the convergence of synaptic inputs from a diverse cross-section of upstream neurons. Furthermore, neural activity is often scale-free, implying that some measurements should be the same, whether taken at large or small scales. Together, convergence and scale-freeness support the hypothesis that single Vm recordings carry useful information about high-dimensional cortical activity. Conveniently, the theory of "critical branching networks" (one purported explanation for scale-freeness) provides testable predictions about scale-free measurements that are readily applied to Vm fluctuations. To investigate, we obtained whole-cell current-clamp recordings of pyramidal neurons in visual cortex of turtles with unknown genders. We isolated fluctuations in Vm below the firing threshold and analyzed them by adapting the definition of "neuronal avalanches" (i.e., spurts of population spiking). The Vm fluctuations which we analyzed were scale-free and consistent with critical branching. These findings recapitulated results from large-scale cortical population data obtained separately in complementary experiments using microelectrode arrays described previously (Shew et al., 2015). Simultaneously recorded single-unit local field potential did not provide a good match, demonstrating the specific utility of Vm Modeling shows that estimation of dynamical network properties from neuronal inputs is most accurate when networks are structured as critical branching networks. In conclusion, these findings extend evidence of critical phenomena while also establishing subthreshold pyramidal neuron Vm fluctuations as an informative gauge of high-dimensional cortical population activity.SIGNIFICANCE STATEMENT The relationship between membrane potential (Vm) dynamics of single neurons and population dynamics is indispensable to understanding cortical circuits. Just as important to the biophysics of computation are emergent properties such as scale-freeness, where critical branching networks offer insight. This report makes progress on both fronts by comparing statistics from single-neuron whole-cell recordings with population statistics obtained with microelectrode arrays. Not only are fluctuations of somatic Vm scale-free, they match fluctuations of population activity. Thus, our results demonstrate appropriation of the brain's own subsampling method (convergence of synaptic inputs) while extending the range of fundamental evidence for critical phenomena in neural systems from the previously observed mesoscale (fMRI, LFP, population spiking) to the microscale, namely, Vm fluctuations.

The turtle visual system mediates a complex spatiotemporal transformation of visual stimuli into cortical activity.

The three-layered visual cortex of turtle is characterized by extensive intracortical axonal projections and receives non-retinotopic axonal projections from lateral geniculate nucleus. What spatiotemporal transformation of visual stimuli into cortical activity arises from such tangle of malleable cortical inputs and intracortical connections? To address this question, we obtained band-pass filtered extracellular recordings of neural activity in turtle dorsal cortex during visual stimulation of the retina. We discovered important spatial and temporal features of stimulus-modulated cortical local field potential (LFP) recordings. Spatial receptive fields span large areas of the visual field, have an intricate internal structure, and lack directional tuning. The receptive field structure varies across recording sites in a distant-dependent manner. Such composite spatial organization of stimulus-modulated cortical activity is accompanied by an equally multifaceted temporal organization. Cortical visual responses are delayed, persistent, and oscillatory. Further, prior cortical activity contributes globally to adaptation in turtle visual cortex. In conclusion, these results demonstrate convoluted spatiotemporal transformations of visual stimuli into stimulus-modulated cortical activity that, at present, largely evade computational frameworks.

Coupling of synaptic inputs to local cortical activity differs among neurons and adapts after stimulus onset.

Cortical activity contributes significantly to the high variability of sensory responses of interconnected pyramidal neurons, which has crucial implications for sensory coding. Yet, largely because of technical limitations of in vivo intracellular recordings, the coupling of a pyramidal neuron's synaptic inputs to the local cortical activity has evaded full understanding. Here we obtained excitatory synaptic conductance ( g) measurements from putative pyramidal neurons and local field potential (LFP) recordings from adjacent cortical circuits during visual processing in the turtle whole brain ex vivo preparation. We found a range of g-LFP coupling across neurons. Importantly, for a given neuron, g-LFP coupling increased at stimulus onset and then relaxed toward intermediate values during continued visual stimulation. A model network with clustered connectivity and synaptic depression reproduced both the diversity and the dynamics of g-LFP coupling. In conclusion, these results establish a rich dependence of single-neuron responses on anatomical, synaptic, and emergent network properties. NEW & NOTEWORTHY Cortical neurons are strongly influenced by the networks in which they are embedded. To understand sensory processing, we must identify the nature of this influence and its underlying mechanisms. Here we investigate synaptic inputs to cortical neurons, and the nearby local field potential, during visual processing. We find a range of neuron-to-network coupling across cortical neurons. This coupling is dynamically modulated during visual processing via biophysical and emergent network properties.

Induced cortical oscillations in turtle cortex are coherent at the mesoscale of population activity, but not at the microscale of the membrane potential of neurons.

Bursts of oscillatory neural activity have been hypothesized to be a core mechanism by which remote brain regions can communicate. Such a hypothesis raises the question to what extent oscillations are coherent across spatially distant neural populations. To address this question, we obtained local field potential (LFP) and membrane potential recordings from the visual cortex of turtle in response to visual stimulation of the retina. The time-frequency analysis of these recordings revealed pronounced bursts of oscillatory neural activity and a large trial-to-trial variability in the spectral and temporal properties of the observed oscillations. First, local bursts of oscillations varied from trial to trial in both burst duration and peak frequency. Second, oscillations of a given recording site were not autocoherent; i.e., the phase did not progress linearly in time. Third, LFP oscillations at spatially separate locations within the visual cortex were more phase coherent in the presence of visual stimulation than during ongoing activity. In contrast, the membrane potential oscillations from pairs of simultaneously recorded pyramidal neurons showed smaller phase coherence, which did not change when switching from black screen to visual stimulation. In conclusion, neuronal oscillations at distant locations in visual cortex are coherent at the mesoscale of population activity, but coherence is largely absent at the microscale of the membrane potential of neurons.NEW & NOTEWORTHY Coherent oscillatory neural activity has long been hypothesized as a potential mechanism for communication across locations in the brain. In this study we confirm the existence of coherent oscillations at the mesoscale of integrated cortical population activity. However, at the microscopic level of neurons, we find no evidence for coherence among oscillatory membrane potential fluctuations. These results raise questions about the applicability of the communication through coherence hypothesis to the level of the membrane potential.

Criticality predicts maximum irregularity in recurrent networks of excitatory nodes.

A rigorous understanding of brain dynamics and function requires a conceptual bridge between multiple levels of organization, including neural spiking and network-level population activity. Mounting evidence suggests that neural networks of cerebral cortex operate at a critical regime, which is defined as a transition point between two phases of short lasting and chaotic activity. However, despite the fact that criticality brings about certain functional advantages for information processing, its supporting evidence is still far from conclusive, as it has been mostly based on power law scaling of size and durations of cascades of activity. Moreover, to what degree such hypothesis could explain some fundamental features of neural activity is still largely unknown. One of the most prevalent features of cortical activity in vivo is known to be spike irregularity of spike trains, which is measured in terms of the coefficient of variation (CV) larger than one. Here, using a minimal computational model of excitatory nodes, we show that irregular spiking (CV > 1) naturally emerges in a recurrent network operating at criticality. More importantly, we show that even at the presence of other sources of spike irregularity, being at criticality maximizes the mean coefficient of variation of neurons, thereby maximizing their spike irregularity. Furthermore, we also show that such a maximized irregularity results in maximum correlation between neuronal firing rates and their corresponding spike irregularity (measured in terms of CV). On the one hand, using a model in the universality class of directed percolation, we propose new hallmarks of criticality at single-unit level, which could be applicable to any network of excitable nodes. On the other hand, given the controversy of the neural criticality hypothesis, we discuss the limitation of this approach to neural systems and to what degree they support the criticality hypothesis in real neural networks. Finally, we discuss the limitations of applying our results to real networks and to what degree they support the criticality hypothesis.

Network activity influences the subthreshold and spiking visual responses of pyramidal neurons in the three-layer turtle cortex.

A primary goal of systems neuroscience is to understand cortical function, typically by studying spontaneous and stimulus-modulated cortical activity. Mounting evidence suggests a strong and complex relationship exists between the ongoing and stimulus-modulated cortical state. To date, most work in this area has been based on spiking in populations of neurons. While advantageous in many respects, this approach is limited in scope: it records the activity of a minority of neurons and gives no direct indication of the underlying subthreshold dynamics. Membrane potential recordings can fill these gaps in our understanding, but stable recordings are difficult to obtain in vivo. Here, we recorded subthreshold cortical visual responses in the ex vivo turtle eye-attached whole brain preparation, which is ideally suited for such a study. We found that, in the absence of visual stimulation, the network was "synchronous"; neurons displayed network-mediated transitions between hyperpolarized (Down) and depolarized (Up) membrane potential states. The prevalence of these slow-wave transitions varied across turtles and recording sessions. Visual stimulation evoked similar Up states, which were on average larger and less reliable when the ongoing state was more synchronous. Responses were muted when immediately preceded by large, spontaneous Up states. Evoked spiking was sparse, highly variable across trials, and mediated by concerted synaptic inputs that were, in general, only very weakly correlated with inputs to nearby neurons. Together, these results highlight the multiplexed influence of the cortical network on the spontaneous and sensory-evoked activity of individual cortical neurons.NEW & NOTEWORTHY Most studies of cortical activity focus on spikes. Subthreshold membrane potential recordings can provide complementary insight, but stable recordings are difficult to obtain in vivo. Here, we recorded the membrane potentials of cortical neurons during ongoing and visually evoked activity. We observed a strong relationship between network and single-neuron evoked activity spanning multiple temporal scales. The membrane potential perspective of cortical dynamics thus highlights the influence of intrinsic network properties on visual processing.

Adaptation modulates correlated subthreshold response variability in visual cortex.

Cortical sensory responses are highly variable across stimulus presentations. This variability can be correlated across neurons (due to some combination of dense intracortical connectivity, cortical activity level, and cortical state), with fundamental implications for population coding. Yet the interpretation of correlated response variability (or "noise correlation") has remained fraught with difficulty, in part because of the restriction to extracellular neuronal spike recordings. Here, we measured response variability and its correlation at the most microscopic level of electrical neural activity, the membrane potential, by obtaining dual whole cell recordings from pairs of cortical pyramidal neurons during visual processing in the turtle whole brain ex vivo preparation. We found that during visual stimulation, correlated variability adapts toward an intermediate level and that this correlation dynamic is likely mediated by intracortical mechanisms. A model network with external inputs, synaptic depression, and structure reproduced the observed dynamics of correlated variability. These results suggest that intracortical adaptation self-organizes cortical circuits toward a balanced regime at which correlated variability is maintained at an intermediate level.NEW & NOTEWORTHY Correlated response variability has profound implications for stimulus encoding, yet our understanding of this phenomenon is based largely on spike data. Here, we investigate the dynamics and mechanisms of membrane potential-correlated variability (CC) in visual cortex with a combined experimental and computational approach. We observe a visually evoked increase in CC, followed by a fast return to baseline. Our results further suggest a link between this observation and the adaptation-mediated dynamics of emergent network phenomena.

Neocortical activity is stimulus- and scale-invariant.

Mounting evidence supports the hypothesis that the cortex operates near a critical state, defined as the transition point between order (large-scale activity) and disorder (small-scale activity). This criticality is manifested by power law distribution of the size and duration of spontaneous cascades of activity, which are referred as neuronal avalanches. The existence of such neuronal avalanches has been confirmed by several studies both in vitro and in vivo, among different species and across multiple spatial scales. However, despite the prevalence of scale free activity, still very little is known concerning whether and how the scale-free nature of cortical activity is altered during external stimulation. To address this question, we performed in vivo two-photon population calcium imaging of layer 2/3 neurons in primary visual cortex of behaving mice during visual stimulation and conducted statistical analyses on the inferred spike trains. Our investigation for each mouse and condition revealed power law distributed neuronal avalanches, and irregular spiking individual neurons. Importantly, both the avalanche and the spike train properties remained largely unchanged for different stimuli, while the cross-correlation structure varied with stimuli. Our results establish that microcircuits in the visual cortex operate near the critical regime, while rearranging functional connectivity in response to varying sensory inputs.

Adaptation towards scale-free dynamics improves cortical stimulus discrimination at the cost of reduced detection.

Fundamental to the function of nervous systems is the ability to reorganize to cope with changing sensory input. Although well-studied in single neurons, how such adaptive versatility manifests in the collective population dynamics and function of cerebral cortex remains unknown. Here we measured population neural activity with microelectrode arrays in turtle visual cortex while visually stimulating the retina. First, we found that, following the onset of stimulation, adaptation tunes the collective population dynamics towards a special regime with scale-free spatiotemporal activity, after an initial large-scale transient response. Concurrently, we observed an adaptive tradeoff between two important aspects of population coding-sensory detection and discrimination. As adaptation tuned the cortex toward scale-free dynamics, stimulus discrimination was enhanced, while stimulus detection was reduced. Finally, we used a network-level computational model to show that short-term synaptic depression was sufficient to mechanistically explain our experimental results. In the model, scale-free dynamics emerge only when the model operates near a special regime called criticality. Together our model and experimental results suggest unanticipated functional benefits and costs of adaptation near criticality in visual cortex.

Coherent and intermittent ensemble oscillations emerge from networks of irregular spiking neurons.

Local field potential (LFP) recordings from spatially distant cortical circuits reveal episodes of coherent gamma oscillations that are intermittent, and of variable peak frequency and duration. Concurrently, single neuron spiking remains largely irregular and of low rate. The underlying potential mechanisms of this emergent network activity have long been debated. Here we reproduce such intermittent ensemble oscillations in a model network, consisting of excitatory and inhibitory model neurons with the characteristics of regular-spiking (RS) pyramidal neurons, and fast-spiking (FS) and low-threshold spiking (LTS) interneurons. We find that fluctuations in the external inputs trigger reciprocally connected and irregularly spiking RS and FS neurons in episodes of ensemble oscillations, which are terminated by the recruitment of the LTS population with concurrent accumulation of inhibitory conductance in both RS and FS neurons. The model qualitatively reproduces experimentally observed phase drift, oscillation episode duration distributions, variation in the peak frequency, and the concurrent irregular single-neuron spiking at low rate. Furthermore, consistent with previous experimental studies using optogenetic manipulation, periodic activation of FS, but not RS, model neurons causes enhancement of gamma oscillations. In addition, increasing the coupling between two model networks from low to high reveals a transition from independent intermittent oscillations to coherent intermittent oscillations. In conclusion, the model network suggests biologically plausible mechanisms for the generation of episodes of coherent intermittent ensemble oscillations with irregular spiking neurons in cortical circuits.

Inferring presynaptic population spiking from single-trial membrane potential recordings.

BACKGROUND: The time-varying membrane potential of a cortical neuron contains important information about the network activity. Extracting this information requires separating excitatory and inhibitory synaptic inputs from single-trial membrane potential recordings without averaging across trials. NEW METHOD: We propose a method to extract the time course of excitatory and inhibitory synaptic inputs to a neuron from a single-trial membrane potential recording. The method takes advantage of the differences in the time constants and the reversal potentials of the excitatory and inhibitory synaptic currents, which allows the untangling of the two conductance types. RESULTS: We evaluate the applicability of the method on a leaky integrate-and-fire model neuron and find high quality of estimation of excitatory synaptic conductance changes and presynaptic population spikes. Application of the method to a real cortical neuron with known synaptic inputs in a brain slice returns high-quality estimation of the time course of the excitatory synaptic conductance. Application of the method to membrane potential recordings from a cortical pyramidal neuron of an intact brain reveals complex network activity. COMPARISON WITH EXISTING METHODS: Existing methods are based on repeated trials and thus are limited to estimating the statistical features of synaptic conductance changes, or, when based on single trials, are limited to special cases, have low temporal resolution, or are impractically complicated. CONCLUSIONS: We propose and test an efficient method for estimating the full time course of excitatory and inhibitory synaptic conductances from single-trial membrane potential recordings. The method is sufficiently simple to ensure widespread use in neuroscience.