Inhibitory control of correlated intrinsic variability in cortical networks

Sensory Encoding Alternates With Hippocampal Ripples across Cycles of Forebrain Spiking Cascades

The brain's response to external events depends on its internal arousal states, which are dynamically governed by neuromodulatory systems and have recently been linked to coordinated spike timing cascades in widespread brain networks. At rest, both arousal fluctuations and spiking cascades are evident throughout the forebrain and play out over multisecond time scales. Here, by analyzing large-scale neural recording data collected by the Allen Institute, it is demonstrated that these intrinsic processes persist across the mouse brain even during periods of continuous visual stimulation. In the stationary animal, each quasi-periodic cascade cycle systematically influenced 1) the efficacy of encoding in visually responsive brain areas and 2) the incidence of memory-related hippocampal ripples. During this cycle, the phase of high arousal is marked by high efficiency in visual encoding whereas the phase of low arousal is marked by the occurrence of hippocampal ripples. However, during bouts of active locomotion, this cycle is abolished and brain maintained a state of elevated visual coding efficiency, with ripples being nearly absent. It is hypothesized that the infra-slow cascade dynamics reflect an adaptive cycle of alternating exteroceptive sensory sampling and internal mnemonic function that persistently pervades the forebrain, only stopping during active exploration of the environment.

Subthreshold moment analysis of neuronal populations driven by synchronous synaptic inputs

Even when driven by the same stimulus, neuronal responses are well-known to exhibit a striking level of spiking variability. In-vivo electrophysiological recordings also reveal a surprisingly large degree of variability at the subthreshold level. In prior work, we considered biophysically relevant neuronal models to account for the observed magnitude of membrane voltage fluctuations. We found that accounting for these fluctuations requires weak but nonzero synchrony in the spiking activity, in amount that are consistent with experimentally measured spiking correlations. Here we investigate whether such synchrony can explain additional statistical features of the measured neural activity, including neuronal voltage covariability and voltage skewness. Addressing this question involves conducting a generalized moment analysis of conductance-based neurons in response to input drives modeled as correlated jump processes. Technically, we perform such an analysis using fixed-point techniques from queuing theory that are applicable in the stationary regime of activity. We found that weak but nonzero synchrony can consistently explain the experimentally reported voltage covariance and skewness. This confirms the role of synchrony as a primary driver of cortical variability and supports that physiological neural activity emerges as a population-level phenomenon, especially in the spontaneous regime.

Subthreshold variability of neuronal populations driven by synchronous synaptic inputs

Even when driven by the same stimulus, neuronal responses are well-known to exhibit a striking level of spiking variability. In-vivo electrophysiological recordings also reveal a surprisingly large degree of variability at the subthreshold level. In prior work, we considered biophysically relevant neuronal models to account for the observed magnitude of membrane voltage fluctuations. We found that accounting for these fluctuations requires weak but nonzero synchrony in the spiking activity, in amount that are consistent with experimentally measured spiking correlations. Here we investigate whether such synchrony can explain additional statistical features of the measured neural activity, including neuronal voltage covariability and voltage skewness. Addressing this question involves conducting a generalized moment analysis of conductance-based neurons in response to input drives modeled as correlated jump processes. Technically, we perform such an analysis using fixed-point techniques from queuing theory that are applicable in the stationary regime of activity. We found that weak but nonzero synchrony can consistently explain the experimentally reported voltage covariance and skewness. This confirms the role of synchrony as a primary driver of cortical variability and supports that physiological neural activity emerges as a population-level phenomenon, especially in the spontaneous regime.

Role of short-term plasticity and slow temporal dynamics in enhancing time series prediction with a brain-inspired recurrent neural network

Typical reservoir networks are based on random connectivity patterns that differ from brain circuits in two important ways. First, traditional reservoir networks lack synaptic plasticity among recurrent units, whereas cortical networks exhibit plasticity across all neuronal types and cortical layers. Second, reservoir networks utilize random Gaussian connectivity, while cortical networks feature a heavy-tailed distribution of synaptic strengths. It is unclear what are the computational advantages of these features for predicting complex time series. In this study, we integrated short-term plasticity (STP) and lognormal connectivity into a novel recurrent neural network (RNN) framework. The model exhibited rich patterns of population activity characterized by slow coordinated fluctuations. Using graph spectral decomposition, we show that weighted networks with lognormal connectivity and STP yield higher complexity than several graph types. When tested on various tasks involving the prediction of complex time series data, the RNN model outperformed a baseline model with random connectivity as well as several other network architectures. Overall, our results underscore the potential of incorporating brain-inspired features such as STP and heavy-tailed connectivity to enhance the robustness and performance of artificial neural networks in complex data prediction and signal processing tasks.

Touch-evoked traveling waves establish a translaminar spacetime code

Linking sensory-evoked traveling waves to underlying circuit patterns is critical to understanding the neural basis of sensory perception. To form this link, we performed simultaneous electrophysiology and two-photon calcium imaging through transparent NeuroGrids and mapped touch-evoked traveling waves and underlying microcircuit dynamics. In awake mice, both passive and active whisker touch elicited traveling waves within and across barrels, with a fast early component followed by a late wave that lasted hundreds of milliseconds poststimulus. Notably, late waves were modulated by perceived value and predicted behavioral choice in a two-whisker discrimination task. We found that the late wave feature was (i) modulated by motor feedback, (ii) differentially engaged a sparse ensemble reactivation pattern across layer 2/3, which a balanced-state network model reconciled via feedback-induced inhibitory stabilization, and (iii) aligned to regenerative layer 5 apical dendritic Ca2+ events. Our results reveal that translaminar spacetime patterns organized by cortical feedback support sparse touch-evoked traveling waves.

Automated customization of large-scale spiking network models to neuronal population activity

Understanding brain function is facilitated by constructing computational models that accurately reproduce aspects of brain activity. Networks of spiking neurons capture the underlying biophysics of neuronal circuits, yet their activity's dependence on model parameters is notoriously complex. As a result, heuristic methods have been used to configure spiking network models, which can lead to an inability to discover activity regimes complex enough to match large-scale neuronal recordings. Here we propose an automatic procedure, Spiking Network Optimization using Population Statistics (SNOPS), to customize spiking network models that reproduce the population-wide covariability of large-scale neuronal recordings. We first confirmed that SNOPS accurately recovers simulated neural activity statistics. Then, we applied SNOPS to recordings in macaque visual and prefrontal cortices and discovered previously unknown limitations of spiking network models. Taken together, SNOPS can guide the development of network models, thereby enabling deeper insight into how networks of neurons give rise to brain function.

Pathological cell assembly dynamics in a striatal MSN network model

Under normal conditions the principal cells of the striatum, medium spiny neurons (MSNs), show structured cell assembly activity patterns which alternate sequentially over exceedingly long timescales of many minutes. It is important to understand this activity since it is characteristically disrupted in multiple pathologies, such as Parkinson's disease and dyskinesia, and thought to be caused by alterations in the MSN to MSN lateral inhibitory connections and in the strength and distribution of cortical excitation to MSNs. To understand how these long timescales arise we extended a previous network model of MSN cells to include synapses with short-term plasticity, with parameters taken from a recent detailed striatal connectome study. We first confirmed the presence of sequentially switching cell clusters using the non-linear dimensionality reduction technique, Uniform Manifold Approximation and Projection (UMAP). We found that the network could generate non-stationary activity patterns varying extremely slowly on the order of minutes under biologically realistic conditions. Next we used Simulation Based Inference (SBI) to train a deep net to map features of the MSN network generated cell assembly activity to MSN network parameters. We used the trained SBI model to estimate MSN network parameters from ex-vivo brain slice calcium imaging data. We found that best fit network parameters were very close to their physiologically observed values. On the other hand network parameters estimated from Parkinsonian, decorticated and dyskinetic ex-vivo slice preparations were different. Our work may provide a pipeline for diagnosis of basal ganglia pathology from spiking data as well as for the design pharmacological treatments.

A Translaminar Spacetime Code Supports Touch-Evoked Traveling Waves

Linking sensory-evoked traveling waves to underlying circuit patterns is critical to understanding the neural basis of sensory perception. To form this link, we performed simultaneous electrophysiology and two-photon calcium imaging through transparent NeuroGrids and mapped touch-evoked cortical traveling waves and their underlying microcircuit dynamics. In awake mice, both passive and active whisker touch elicited traveling waves within and across barrels, with a fast early component followed by a variable late wave that lasted hundreds of milliseconds post-stimulus. Strikingly, late-wave dynamics were modulated by stimulus value and correlated with task performance. Mechanistically, the late wave component was i) modulated by motor feedback, ii) complemented by a sparse ensemble pattern across layer 2/3, which a balanced-state network model reconciled via inhibitory stabilization, and iii) aligned to regenerative Layer-5 apical dendritic Ca 2+ events. Our results reveal a translaminar spacetime pattern organized by cortical feedback in the sensory cortex that supports touch-evoked traveling waves.

GRAPHICAL ABSTRACT AND HIGHLIGHTS

Whisker touch evokes both early- and late-traveling waves in the barrel cortex over 100's of millisecondsReward reinforcement modulates wave dynamics Late wave emergence coincides with network sparsity in L23 and time-locked L5 dendritic Ca 2+ spikes Experimental and computational results link motor feedback to distinct translaminar spacetime patterns.

Rapid compensatory plasticity revealed by dynamic correlated activity in monkeys in vivo

To produce adaptive behavior, neural networks must balance between plasticity and stability. Computational work has demonstrated that network stability requires plasticity mechanisms to be counterbalanced by rapid compensatory processes. However, such processes have yet to be experimentally observed. Here we demonstrate that repeated optogenetic activation of excitatory neurons in monkey visual cortex (area V1) induces a population-wide dynamic reduction in the strength of neuronal interactions over the timescale of minutes during the awake state, but not during rest. This new form of rapid plasticity was observed only in the correlation structure, with firing rates remaining stable across trials. A computational network model operating in the balanced regime confirmed experimental findings and revealed that inhibitory plasticity is responsible for the decrease in correlated activity in response to repeated light stimulation. These results provide the first experimental evidence for rapid homeostatic plasticity that primarily operates during wakefulness, which stabilizes neuronal interactions during strong network co-activation.

VIP interneurons regulate cortical size tuning and visual perception

Cortical circuit function is regulated by extensively interconnected, diverse populations of GABAergic interneurons that may play key roles in shaping circuit operation according to behavioral context. A specialized population of interneurons that co-express vasoactive intestinal peptides (VIP-INs) are activated during arousal and innervate other INs and pyramidal neurons (PNs). Although state-dependent modulation of VIP-INs has been extensively studied, their role in regulating sensory processing is less well understood. We examined the impact of VIP-INs in the primary visual cortex of awake behaving mice. Loss of VIP-IN activity alters the behavioral state-dependent modulation of somatostatin-expressing INs (SST-INs) but not PNs. In contrast, reduced VIP-IN activity globally disrupts visual feature selectivity for stimulus size. Moreover, the impact of VIP-INs on perceptual behavior varies with context and is more acute for small than large visual cues. VIP-INs thus contribute to both state-dependent modulation of cortical activity and sensory context-dependent perceptual performance.

Automated customization of large-scale spiking network models to neuronal population activity

Understanding brain function is facilitated by constructing computational models that accurately reproduce aspects of brain activity. Networks of spiking neurons capture the underlying biophysics of neuronal circuits, yet the dependence of their activity on model parameters is notoriously complex. As a result, heuristic methods have been used to configure spiking network models, which can lead to an inability to discover activity regimes complex enough to match large-scale neuronal recordings. Here we propose an automatic procedure, Spiking Network Optimization using Population Statistics (SNOPS), to customize spiking network models that reproduce the population-wide covariability of large-scale neuronal recordings. We first confirmed that SNOPS accurately recovers simulated neural activity statistics. Then, we applied SNOPS to recordings in macaque visual and prefrontal cortices and discovered previously unknown limitations of spiking network models. Taken together, SNOPS can guide the development of network models and thereby enable deeper insight into how networks of neurons give rise to brain function.

Interneuron-Targeted Disruption of SYNGAP1 Alters Sensory Representations in the Neocortex and Impairs Sensory Learning

SYNGAP1 haploinsufficiency in humans leads to severe neurodevelopmental disorders characterized by intellectual disability, autism, epilepsy, and sensory processing deficits. However, the circuit mechanisms underlying these disorders are not well understood. In mice, a decrease of SynGAP levels results in cognitive deficits by interfering with the development of excitatory glutamatergic connections. Recent evidence suggests that SynGAP also plays a crucial role in the development and function of GABAergic inhibitory interneurons. Nevertheless, it remains uncertain whether and to what extent the expression of SYNGAP1 in inhibitory interneurons contributes to cortical circuit function and related behaviors. The activity of cortical neurons has not been measured simultaneously with behavior. To address these gaps, we recorded from layer 2/3 neurons in the primary whisker somatosensory cortex (wS1) of mice while they learned to perform a whisker tactile detection task. Our results demonstrate that mice with interneuron-specific SYNGAP1 haploinsufficiency exhibit learning deficits characterized by heightened behavioral responses in the absence of relevant sensory input and premature responses to unrelated sensory stimuli not associated with reward acquisition. These behavioral deficits are accompanied by specific circuit abnormalities within wS1. Interneuron-specific SYNGAP1 haploinsufficiency increases detrimental neuronal correlations directly related to task performance and enhances responses to irrelevant sensory stimuli unrelated to the reward acquisition. In summary, our findings indicate that a reduction of SynGAP in inhibitory interneurons impairs sensory representation in the primary sensory cortex by disrupting neuronal correlations, which likely contributes to the observed cognitive deficits in mice with pan-neuronal SYNGAP1 haploinsufficiency.SIGNIFICANCE STATEMENT SYNGAP1 haploinsufficiency leads to severe neurodevelopmental disorders. The exact nature of neural circuit dysfunction caused by SYNGAP1 haploinsufficiency remains poorly understood. SynGAP plays a critical role in the function of GABAergic inhibitory interneurons as well as glutamatergic pyramidal neurons in the neocortex. Whether and how decreasing SYNGAP1 level in inhibitory interneurons disrupts a behaviorally relevant circuit remains unclear. We measure neural activity and behavior in mice learning a perceptual task. Mice with interneuron-targeted disruption of SYNGAP1 display increased detrimental neuronal correlations and elevated responses to irrelevant sensory inputs, which are related to impaired task performance. These results show that cortical interneuron dysfunction contributes to sensory deficits in SYNGAP1 haploinsufficiency with important implications for identifying therapeutic targets.

VIP interneurons regulate cortical size tuning and visual perception

Local cortical circuit function is regulated by diverse populations of GABAergic interneurons with distinct properties and extensive interconnectivity. Inhibitory-to-inhibitory interactions between interneuron populations may play key roles in shaping circuit operation according to behavioral context. A specialized population of GABAergic interneurons that co-express vasoactive intestinal peptide (VIP-INs) are activated during arousal and locomotion and innervate other local interneurons and pyramidal neurons. Although modulation of VIP-IN activity by behavioral state has been extensively studied, their role in regulating information processing and selectivity is less well understood. Using a combination of cellular imaging, short and long-term manipulation, and perceptual behavior, we examined the impact of VIP-INs on their synaptic target populations in the primary visual cortex of awake behaving mice. We find that loss of VIP-IN activity alters the behavioral state-dependent modulation of somatostatin-expressing interneurons (SST-INs) but not pyramidal neurons (PNs). In contrast, reduced VIP-IN activity disrupts visual feature selectivity for stimulus size in both populations. Inhibitory-to inhibitory interactions thus directly shape the selectivity of GABAergic interneurons for sensory stimuli. Moreover, the impact of VIP-IN activity on perceptual behavior varies with visual context and is more acute for small than large visual cues. VIP-INs thus contribute to both state-dependent modulation of cortical circuit activity and sensory context-dependent perceptual performance.

A deep generative adversarial network capturing complex spiral waves in disinhibited circuits of the cerebral cortex

BACKGROUND

In the cerebral cortex, disinhibited activity is characterized by propagating waves that spread across neural tissue. In this pathological state, a widely reported form of activity are spiral waves that travel in a circular pattern around a fixed spatial locus termed the center of mass. Spiral waves exhibit stereotypical activity and involve broad patterns of co-fluctuations, suggesting that they may be of lower complexity than healthy activity.

RESULTS

To evaluate this hypothesis, we performed dense multi-electrode recordings of cortical networks where disinhibition was induced by perfusing a pro-epileptiform solution containing 4-Aminopyridine as well as increased potassium and decreased magnesium. Spiral waves were identified based on a spatially delimited center of mass and a broad distribution of instantaneous phases across electrodes. Individual waves were decomposed into "snapshots" that captured instantaneous neural activation across the entire network. The complexity of these snapshots was examined using a measure termed the participation ratio. Contrary to our expectations, an eigenspectrum analysis of these snapshots revealed a broad distribution of eigenvalues and an increase in complexity compared to baseline networks. A deep generative adversarial network was trained to generate novel exemplars of snapshots that closely captured cortical spiral waves. These synthetic waves replicated key features of experimental data including a tight center of mass, a broad eigenvalue distribution, spatially-dependent correlations, and a high complexity. By adjusting the input to the model, new samples were generated that deviated in systematic ways from the experimental data, thus allowing the exploration of a broad range of states from healthy to pathologically disinhibited neural networks.

CONCLUSIONS

Together, results show that the complexity of population activity serves as a marker along a continuum from healthy to disinhibited brain states. The proposed generative adversarial network opens avenues for replicating the dynamics of cortical seizures and accelerating the design of optimal neurostimulation aimed at suppressing pathological brain activity.

Size tuning of neural response variability in laminar circuits of macaque primary visual cortex

A defining feature of the cortex is its laminar organization, which is likely critical for cortical information processing. For example, visual stimuli of different size evoke distinct patterns of laminar activity. Visual information processing is also influenced by the response variability of individual neurons and the degree to which this variability is correlated among neurons. To elucidate laminar processing, we studied how neural response variability across the layers of macaque primary visual cortex is modulated by visual stimulus size. Our laminar recordings revealed that single neuron response variability and the shared variability among neurons are tuned for stimulus size, and this size-tuning is layer-dependent. In all layers, stimulation of the receptive field (RF) reduced single neuron variability, and the shared variability among neurons, relative to their pre-stimulus values. As the stimulus was enlarged beyond the RF, both single neuron and shared variability increased in supragranular layers, but either did not change or decreased in other layers. Surprisingly, we also found that small visual stimuli could increase variability relative to baseline values. Our results suggest multiple circuits and mechanisms as the source of variability in different layers and call for the development of new models of neural response variability.

It's about time: Linking dynamical systems with human neuroimaging to understand the brain

Most human neuroscience research to date has focused on statistical approaches that describe stationary patterns of localized neural activity or blood flow. While these patterns are often interpreted in light of dynamic, information-processing concepts, the static, local, and inferential nature of the statistical approach makes it challenging to directly link neuroimaging results to plausible underlying neural mechanisms. Here, we argue that dynamical systems theory provides the crucial mechanistic framework for characterizing both the brain's time-varying quality and its partial stability in the face of perturbations, and hence, that this perspective can have a profound impact on the interpretation of human neuroimaging results and their relationship with behavior. After briefly reviewing some key terminology, we identify three key ways in which neuroimaging analyses can embrace a dynamical systems perspective: by shifting from a local to a more global perspective, by focusing on dynamics instead of static snapshots of neural activity, and by embracing modeling approaches that map neural dynamics using "forward" models. Through this approach, we envisage ample opportunities for neuroimaging researchers to enrich their understanding of the dynamic neural mechanisms that support a wide array of brain functions, both in health and in the setting of psychopathology.

[Auditory response patterns of mouse primary auditory cortex to sound stimuli]

OBJECTIVE

To investigate the auditory response patterns of mouse primary auditory cortex (A1) neurons.

METHODS

In vivo cell-attached recordings and neural network modeling were performed to detect the changes in response patterns of A1 neurons of awake C57BL/6J mice to sound stimulation with varying lengths. A1 neuron signals were recorded for 216 neurons in 20 awake mice using a target sound stimulation sequence, and the classification and response characteristics of A1 neuron response patterns were examined using post-stimulus spike time histograms. To simulate the diversity of the A1 neuron response patterns, an A1 neuron model was established based on the Wilson-Cowan model and integral-firing model. The neuron connection weight parameters in the model were calculated by examining the micro loop structure of the pyramidal neurons, parvalbumin neurons, and somatostatin neurons in the A1 region, and the A1 neural network information coding model was constructed.

RESULTS

The Onset response neurons only had fast spike response within 10 to 40 ms after the beginning of noise stimulation (122 neurons). The Sustained response neurons had spike response continuously during the noise stimulation (26 neurons). The On-off response neurons had fast spike response after the beginning and the end of noise stimulation (40 neurons). The Offset response neurons only had fast spike response within 10 to 40 ms after the end of noise stimulation (22 neurons). In the neural network model, the Onset peak neural activities of A1 pyramidal neurons, parvalbumin neurons, and somatostatin neurons were 0.7483, 0.5236 and 0.9427, respectively, and their response half peak widths were 18.5 ms, 12 ms and 31 ms during the 100 ms noise stimulation, respectively. By changing the feedforward excitation and synaptic inhibition time constants in the model, the neurons generated numerous different types of spike train.

CONCLUSION

The auditory response of mouse A1 neurons to sound stimuli shows mainly the Onset, Sustained, On-off, and Offset response patterns.

Not so spontaneous: Multi-dimensional representations of behaviors and context in sensory areas

Sensory areas are spontaneously active in the absence of sensory stimuli. This spontaneous activity has long been studied; however, its functional role remains largely unknown. Recent advances in technology, allowing large-scale neural recordings in the awake and behaving animal, have transformed our understanding of spontaneous activity. Studies using these recordings have discovered high-dimensional spontaneous activity patterns, correlation between spontaneous activity and behavior, and dissimilarity between spontaneous and sensory-driven activity patterns. These findings are supported by evidence from developing animals, where a transition toward these characteristics is observed as the circuit matures, as well as by evidence from mature animals across species. These newly revealed characteristics call for the formulation of a new role for spontaneous activity in neural sensory computation.

Secondary auditory cortex mediates a sensorimotor mechanism for action timing

The ability to accurately determine when to perform an action is a fundamental brain function and vital to adaptive behavior. The behavioral mechanism and neural circuit for action timing, however, remain largely unknown. Using a new, self-paced action timing task in mice, we found that deprivation of auditory, but not somatosensory or visual input, disrupts learned action timing. The hearing effect was dependent on the auditory feedback derived from the animal's own actions, rather than passive environmental cues. Neuronal activity in the secondary auditory cortex was found to be both correlated with and necessary for the proper execution of learned action timing. Closed-loop, action-dependent optogenetic stimulation of the specific task-related neuronal population within the secondary auditory cortex rescued the key features of learned action timing under auditory deprivation. These results unveil a previously underappreciated sensorimotor mechanism in which the secondary auditory cortex transduces self-generated audiomotor feedback to control action timing.

A flexible Bayesian framework for unbiased estimation of timescales

Timescales characterize the pace of change for many dynamic processes in nature. Timescales are usually estimated by fitting the exponential decay of data autocorrelation in the time or frequency domain. Here we show that this standard procedure often fails to recover the correct timescales due to a statistical bias arising from the finite sample size. We develop an alternative approach which estimates timescales by fitting the sample autocorrelation or power spectrum with a generative model based on a mixture of Ornstein-Uhlenbeck (OU) processes using adaptive approximate Bayesian computations (aABC). Our method accounts for finite sample size and noise in data and returns a posterior distribution of timescales that quantifies the estimation uncertainty and can be used for model selection. We demonstrate the accuracy of our method on synthetic data and illustrate its application to recordings from primate cortex. We provide a customizable Python package implementing our framework with different generative models suitable for diverse applications.

Cortical state dynamics and selective attention define the spatial pattern of correlated variability in neocortex

Correlated activity fluctuations in the neocortex influence sensory responses and behavior. Neural correlations reflect anatomical connectivity but also change dynamically with cognitive states such as attention. Yet, the network mechanisms defining the population structure of correlations remain unknown. We measured correlations within columns in the visual cortex. We show that the magnitude of correlations, their attentional modulation, and dependence on lateral distance are explained by columnar On-Off dynamics, which are synchronous activity fluctuations reflecting cortical state. We developed a network model in which the On-Off dynamics propagate across nearby columns generating spatial correlations with the extent controlled by attentional inputs. This mechanism, unlike previous proposals, predicts spatially non-uniform changes in correlations during attention. We confirm this prediction in our columnar recordings by showing that in superficial layers the largest changes in correlations occur at intermediate lateral distances. Our results reveal how spatially structured patterns of correlated variability emerge through interactions of cortical state dynamics, anatomical connectivity, and attention.

Structure and function of axo-axonic inhibition

Inhibitory neurons in mammalian cortex exhibit diverse physiological, morphological, molecular, and connectivity signatures. While considerable work has measured the average connectivity of several interneuron classes, there remains a fundamental lack of understanding of the connectivity distribution of distinct inhibitory cell types with synaptic resolution, how it relates to properties of target cells, and how it affects function. Here, we used large-scale electron microscopy and functional imaging to address these questions for chandelier cells in layer 2/3 of the mouse visual cortex. With dense reconstructions from electron microscopy, we mapped the complete chandelier input onto 153 pyramidal neurons. We found that synapse number is highly variable across the population and is correlated with several structural features of the target neuron. This variability in the number of axo-axonic ChC synapses is higher than the variability seen in perisomatic inhibition. Biophysical simulations show that the observed pattern of axo-axonic inhibition is particularly effective in controlling excitatory output when excitation and inhibition are co-active. Finally, we measured chandelier cell activity in awake animals using a cell-type-specific calcium imaging approach and saw highly correlated activity across chandelier cells. In the same experiments, in vivo chandelier population activity correlated with pupil dilation, a proxy for arousal. Together, these results suggest that chandelier cells provide a circuit-wide signal whose strength is adjusted relative to the properties of target neurons.

Sensory coding and contrast invariance emerge from the control of plastic inhibition over emergent selectivity

Visual stimuli are represented by a highly efficient code in the primary visual cortex, but the development of this code is still unclear. Two distinct factors control coding efficiency: Representational efficiency, which is determined by neuronal tuning diversity, and metabolic efficiency, which is influenced by neuronal gain. How these determinants of coding efficiency are shaped during development, supported by excitatory and inhibitory plasticity, is only partially understood. We investigate a fully plastic spiking network of the primary visual cortex, building on phenomenological plasticity rules. Our results suggest that inhibitory plasticity is key to the emergence of tuning diversity and accurate input encoding. We show that inhibitory feedback (random and specific) increases the metabolic efficiency by implementing a gain control mechanism. Interestingly, this led to the spontaneous emergence of contrast-invariant tuning curves. Our findings highlight that (1) interneuron plasticity is key to the development of tuning diversity and (2) that efficient sensory representations are an emergent property of the resulting network.

Synaptic plasticity is crucial for the development of efficient input representation in the different sensory cortices, such as the primary visual cortex. Efficient visual representation is determined by two factors: representational efficiency, i.e. how many different input features can be represented, and metabolic efficiency, i.e. how many spikes are required to represent a specific feature. Previous research has pointed out the importance of plasticity at excitatory synapses to achieve high representational efficiency and feedback inhibition as a gain control mechanism for controlling metabolic efficiency. However, it is only partially understood how the influence of inhibitory plasticity on excitatory plasticity can lead to an efficient representation. Using a spiking neural network, we show that plasticity at feed-forward and feedback inhibitory synapses is necessary for the emergence of well-distributed neuronal selectivity to improve representational efficiency. Further, the emergent balance between excitatory and inhibitory currents improves the metabolic efficiency, and leads to contrast-invariant tuning as an inherent network property. Extending previous work, our simulation results highlight the importance of plasticity at inhibitory synapses.

Collapse of complexity of brain and body activity due to excessive inhibition and MeCP2 disruption

Complex body movements require complex dynamics and coordination among neurons in motor cortex. Conversely, a long-standing theoretical notion supposes that if many neurons in motor cortex become excessively synchronized, they may lack the necessary complexity for healthy motor coding. However, direct experimental support for this idea is rare and underlying mechanisms are unclear. Here we recorded three-dimensional body movements and spiking activity of many single neurons in motor cortex of rats with enhanced synaptic inhibition and a transgenic rat model of Rett syndrome (RTT). For both cases, we found a collapse of complexity in the motor system. Reduced complexity was apparent in lower-dimensional, stereotyped brain-body interactions, neural synchrony, and simpler behavior. Our results demonstrate how imbalanced inhibition can cause excessive synchrony among movement-related neurons and, consequently, a stereotyped motor code. Excessive inhibition and synchrony may underlie abnormal motor function in RTT.

How Cortical Circuits Implement Cortical Computations: Mouse Visual Cortex as a Model

The mouse, as a model organism to study the brain, gives us unprecedented experimental access to the mammalian cerebral cortex. By determining the cortex's cellular composition, revealing the interaction between its different components, and systematically perturbing these components, we are obtaining mechanistic insight into some of the most basic properties of cortical function. In this review, we describe recent advances in our understanding of how circuits of cortical neurons implement computations, as revealed by the study of mouse primary visual cortex. Further, we discuss how studying the mouse has broadened our understanding of the range of computations performed by visual cortex. Finally, we address how future approaches will fulfill the promise of the mouse in elucidating fundamental operations of cortex.

Somatostatin interneurons activated by 5-HT2A receptor suppress slow oscillations in medial entorhinal cortex

Serotonin (5-HT) is one of the major neuromodulators present in the mammalian brain and has been shown to play a role in multiple physiological processes. The mechanisms by which 5-HT modulates cortical network activity, however, are not yet fully understood. We investigated the effects of 5-HT on slow oscillations (SOs), a synchronized cortical network activity universally present across species. SOs are observed during anesthesia and are considered to be the default cortical activity pattern. We discovered that (±)3,4-methylenedioxymethamphetamine (MDMA) and fenfluramine, two potent 5-HT releasers, inhibit SOs within the entorhinal cortex (EC) in anesthetized mice. Combining opto- and pharmacogenetic manipulations with in vitro electrophysiological recordings, we uncovered that somatostatin-expressing (Sst) interneurons activated by the 5-HT_2A receptor (5-HT_2AR) play an important role in the suppression of SOs. Since 5-HT_2AR signaling is involved in the etiology of different psychiatric disorders and mediates the psychological effects of many psychoactive serotonergic drugs, we propose that the newly discovered link between Sst interneurons and 5-HT will contribute to our understanding of these complex topics.

Hippocampal hyperactivity in a rat model of Alzheimer's disease

Neuronal network dysfunction is a hallmark of Alzheimer's disease (AD). However, the underlying pathomechanisms remain unknown. We analyzed the hippocampal micronetwork in transgenic McGill-R-Thy1-APP rats (APPtg) at the beginning of extracellular amyloid beta (Aβ) deposition. We established two-photon Ca²⁺ -imaging in vivo in the hippocampus of rats and found hyperactivity of CA1 neurons. Patch-clamp recordings in brain slices in vitro revealed increased neuronal input resistance and prolonged action potential width in CA1 pyramidal neurons. We did neither observe changes in synaptic inhibition, nor in excitation. Our data support the view that increased intrinsic excitability of CA1 neurons may precede inhibitory dysfunction at an early stage of Aβ-deposition and disease progression.

Low rank mechanisms underlying flexible visual representations

Neuronal population responses to sensory stimuli are remarkably flexible. The responses of neurons in visual cortex have heterogeneous dependence on stimulus properties (e.g., contrast), processes that affect all stages of visual processing (e.g., adaptation), and cognitive processes (e.g., attention or task switching). Understanding whether these processes affect similar neuronal populations and whether they have similar effects on entire populations can provide insight into whether they utilize analogous mechanisms. In particular, it has recently been demonstrated that attention has low rank effects on the covariability of populations of visual neurons, which impacts perception and strongly constrains mechanistic models. We hypothesized that measuring changes in population covariability associated with other sensory and cognitive processes could clarify whether they utilize similar mechanisms or computations. Our experimental design included measurements in multiple visual areas using four distinct sensory and cognitive processes. We found that contrast, adaptation, attention, and task switching affect the variability of responses of populations of neurons in primate visual cortex in a similarly low rank way. These results suggest that a given circuit may use similar mechanisms to perform many forms of modulation and likely reflects a general principle that applies to a wide range of brain areas and sensory, cognitive, and motor processes.

Bayesian inference for biophysical neuron models enables stimulus optimization for retinal neuroprosthetics

While multicompartment models have long been used to study the biophysics of neurons, it is still challenging to infer the parameters of such models from data including uncertainty estimates. Here, we performed Bayesian inference for the parameters of detailed neuron models of a photoreceptor and an OFF- and an ON-cone bipolar cell from the mouse retina based on two-photon imaging data. We obtained multivariate posterior distributions specifying plausible parameter ranges consistent with the data and allowing to identify parameters poorly constrained by the data. To demonstrate the potential of such mechanistic data-driven neuron models, we created a simulation environment for external electrical stimulation of the retina and optimized stimulus waveforms to target OFF- and ON-cone bipolar cells, a current major problem of retinal neuroprosthetics.

Tuning network dynamics from criticality to an asynchronous state

According to many experimental observations, neurons in cerebral cortex tend to operate in an asynchronous regime, firing independently of each other. In contrast, many other experimental observations reveal cortical population firing dynamics that are relatively coordinated and occasionally synchronous. These discrepant observations have naturally led to competing hypotheses. A commonly hypothesized explanation of asynchronous firing is that excitatory and inhibitory synaptic inputs are precisely correlated, nearly canceling each other, sometimes referred to as 'balanced' excitation and inhibition. On the other hand, the 'criticality' hypothesis posits an explanation of the more coordinated state that also requires a certain balance of excitatory and inhibitory interactions. Both hypotheses claim the same qualitative mechanism-properly balanced excitation and inhibition. Thus, a natural question arises: how are asynchronous population dynamics and critical dynamics related, how do they differ? Here we propose an answer to this question based on investigation of a simple, network-level computational model. We show that the strength of inhibitory synapses relative to excitatory synapses can be tuned from weak to strong to generate a family of models that spans a continuum from critical dynamics to asynchronous dynamics. Our results demonstrate that the coordinated dynamics of criticality and asynchronous dynamics can be generated by the same neural system if excitatory and inhibitory synapses are tuned appropriately.

Training deep neural density estimators to identify mechanistic models of neural dynamics

Mechanistic modeling in neuroscience aims to explain observed phenomena in terms of underlying causes. However, determining which model parameters agree with complex and stochastic neural data presents a significant challenge. We address this challenge with a machine learning tool which uses deep neural density estimators-trained using model simulations-to carry out Bayesian inference and retrieve the full space of parameters compatible with raw data or selected data features. Our method is scalable in parameters and data features and can rapidly analyze new data after initial training. We demonstrate the power and flexibility of our approach on receptive fields, ion channels, and Hodgkin-Huxley models. We also characterize the space of circuit configurations giving rise to rhythmic activity in the crustacean stomatogastric ganglion, and use these results to derive hypotheses for underlying compensation mechanisms. Our approach will help close the gap between data-driven and theory-driven models of neural dynamics.

Dynamic representations in networked neural systems

A group of neurons can generate patterns of activity that represent information about stimuli; subsequently, the group can transform and transmit activity patterns across synapses to spatially distributed areas. Recent studies in neuroscience have begun to independently address the two components of information processing: the representation of stimuli in neural activity and the transmission of information in networks that model neural interactions. Yet only recently are studies seeking to link these two types of approaches. Here we briefly review the two separate bodies of literature; we then review the recent strides made to address this gap. We continue with a discussion of how patterns of activity evolve from one representation to another, forming dynamic representations that unfold on the underlying network. Our goal is to offer a holistic framework for understanding and describing neural information representation and transmission while revealing exciting frontiers for future research.

Converging intracortical signatures of two separated processing timescales in human early auditory cortex

Neural oscillations in auditory cortex are argued to support parsing and representing speech constituents at their corresponding temporal scales. Yet, how incoming sensory information interacts with ongoing spontaneous brain activity, what features of the neuronal microcircuitry underlie spontaneous and stimulus-evoked spectral fingerprints, and what these fingerprints entail for stimulus encoding, remain largely open questions. We used a combination of human invasive electrophysiology, computational modeling and decoding techniques to assess the information encoding properties of brain activity and to relate them to a plausible underlying neuronal microarchitecture. We analyzed intracortical auditory EEG activity from 10 patients while they were listening to short sentences. Pre-stimulus neural activity in early auditory cortical regions often exhibited power spectra with a shoulder in the delta range and a small bump in the beta range. Speech decreased power in the beta range, and increased power in the delta-theta and gamma ranges. Using multivariate machine learning techniques, we assessed the spectral profile of information content for two aspects of speech processing: detection and discrimination. We obtained better phase than power information decoding, and a bimodal spectral profile of information content with better decoding at low (delta-theta) and high (gamma) frequencies than at intermediate (beta) frequencies. These experimental data were reproduced by a simple rate model made of two subnetworks with different timescales, each composed of coupled excitatory and inhibitory units, and connected via a negative feedback loop. Modeling and experimental results were similar in terms of pre-stimulus spectral profile (except for the iEEG beta bump), spectral modulations with speech, and spectral profile of information content. Altogether, we provide converging evidence from both univariate spectral analysis and decoding approaches for a dual timescale processing infrastructure in human auditory cortex, and show that it is consistent with the dynamics of a simple rate model.

GABA-mediated tonic inhibition differentially modulates gain in functional subtypes of cortical interneurons

The binding of GABA (γ-aminobutyric acid) to extrasynaptic GABA_A receptors generates tonic inhibition that acts as a powerful modulator of cortical network activity. Despite GABA being present throughout the extracellular space of the brain, previous work has shown that GABA may differentially modulate the excitability of neuron subtypes according to variation in chloride gradient. Here, using biophysically detailed neuron models, we predict that tonic inhibition can differentially modulate the excitability of neuron subtypes according to variation in electrophysiological properties. Surprisingly, tonic inhibition increased the responsiveness (or gain) in models with features typical for somatostatin interneurons but decreased gain in models with features typical for parvalbumin interneurons. Patch-clamp recordings from cortical interneurons supported these predictions, and further in silico analysis was then performed to seek a putative mechanism underlying gain modulation. We found that gain modulation in models was dependent upon the magnitude of tonic current generated at depolarized membrane potential-a property associated with outward rectifying GABA_A receptors. Furthermore, tonic inhibition produced two biophysical changes in models of relevance to neuronal excitability: 1) enhanced action potential repolarization via increased current flow into the dendritic compartment, and 2) reduced activation of voltage-dependent potassium channels. Finally, we show theoretically that reduced potassium channel activation selectively increases gain in models possessing action potential dynamics typical for somatostatin interneurons. Potassium channels in parvalbumin-type models deactivate rapidly and are unavailable for further modulation. These findings show that GABA can differentially modulate interneuron excitability and suggest a mechanism through which this occurs in silico via differences of intrinsic electrophysiological properties.

Inferring and validating mechanistic models of neural microcircuits based on spike-train data

The interpretation of neuronal spike train recordings often relies on abstract statistical models that allow for principled parameter estimation and model selection but provide only limited insights into underlying microcircuits. In contrast, mechanistic models are useful to interpret microcircuit dynamics, but are rarely quantitatively matched to experimental data due to methodological challenges. Here we present analytical methods to efficiently fit spiking circuit models to single-trial spike trains. Using derived likelihood functions, we statistically infer the mean and variance of hidden inputs, neuronal adaptation properties and connectivity for coupled integrate-and-fire neurons. Comprehensive evaluations on synthetic data, validations using ground truth in-vitro and in-vivo recordings, and comparisons with existing techniques demonstrate that parameter estimation is very accurate and efficient, even for highly subsampled networks. Our methods bridge statistical, data-driven and theoretical, model-based neurosciences at the level of spiking circuits, for the purpose of a quantitative, mechanistic interpretation of recorded neuronal population activity.

It is difficult to fit mechanistic, biophysically constrained circuit models to spike train data from in vivo extracellular recordings. Here the authors present analytical methods that enable efficient parameter estimation for integrate-and-fire circuit models and inference of the underlying connectivity structure in subsampled networks.

The Low-Dimensional Neural Architecture of Cognitive Complexity Is Related to Activity in Medial Thalamic Nuclei

Cognitive activity emerges from large-scale neuronal dynamics that are constrained to a low-dimensional manifold. How this low-dimensional manifold scales with cognitive complexity, and which brain regions regulate this process, are not well understood. We addressed this issue by analyzing sub-second high-field fMRI data acquired during performance of a task that systematically varied the complexity of cognitive reasoning. We show that task performance reconfigures the low-dimensional manifold and that deviations from these patterns relate to performance errors. We further demonstrate that individual differences in thalamic activity relate to reconfigurations of the low-dimensional architecture during task engagement.

The Dynamical Regime of Sensory Cortex: Stable Dynamics around a Single Stimulus-Tuned Attractor Account for Patterns of Noise Variability

Correlated variability in cortical activity is ubiquitously quenched following stimulus onset, in a stimulus-dependent manner. These modulations have been attributed to circuit dynamics involving either multiple stable states (“attractors”) or chaotic activity. Here we show that a qualitatively different dynamical regime, involving fluctuations about a single, stimulus-driven attractor in a loosely balanced excitatory-inhibitory network (the stochastic “stabilized supralinear network”), best explains these modulations. Given the supralinear input/output functions of cortical neurons, increased stimulus drive strengthens effective network connectivity. This shifts the balance from interactions that amplify variability to suppressive inhibitory feedback, quenching correlated variability around more strongly driven steady states. Comparing to previously published and original data analyses, we show that this mechanism, unlike previous proposals, uniquely accounts for the spatial patterns and fast temporal dynamics of variability suppression. Specifying the cortical operating regime is key to understanding the computations underlying perception.

•

A simple network model explains stimulus-tuning of cortical variability suppression

•

Inhibition stabilizes recurrently interacting neurons with supralinear I/O functions

•

Stimuli strengthen inhibitory stabilization around a stable state, quenching variability

•

Single-trial V1 data are compatible with this model and rules out competing proposals

Probing spatial inhomogeneity of cholinergic changes in cortical state in rat

Acetylcholine (ACh) plays an essential role in cortical information processing. Cholinergic changes in cortical state can fundamentally change how the neurons encode sensory input and motor output. Traditionally, ACh distribution in cortex and associated changes in cortical state have been assumed to be spatially diffuse. However, recent studies demonstrate a more spatially inhomogeneous structure of cholinergic projections to cortex. Moreover, many experimental manipulations of ACh have been done at a single spatial location, which inevitably results in spatially non-uniform ACh distribution. Such non-uniform application of ACh across the spatial extent of a cortical microcircuit could have important impacts on how the firing of groups of neurons is coordinated, but this remains largely unknown. Here we describe a method for applying ACh at different spatial locations within a single cortical circuit and measuring the resulting differences in population neural activity. We use two microdialysis probes implanted at opposite ends of a microelectrode array in barrel cortex of anesthetized rats. As a demonstration of the method, we applied ACh or neostigmine in different spatial locations via the microdialysis probes while we concomitantly recorded neural activity at 32 locations with the microelectrode array. First, we show that cholinergic changes in cortical state can vary dramatically depending on where the ACh was applied. Second, we show that cholinergic changes in cortical state can vary dramatically depending on where the state-change is measured. These results suggests that previous work with single-site recordings or single-site ACh application should be interpreted with some caution, since the results could change for different spatial locations.

Efficient calculation of heterogeneous non-equilibrium statistics in coupled firing-rate models

Understanding nervous system function requires careful study of transient (non-equilibrium) neural response to rapidly changing, noisy input from the outside world. Such neural response results from dynamic interactions among multiple, heterogeneous brain regions. Realistic modeling of these large networks requires enormous computational resources, especially when high-dimensional parameter spaces are considered. By assuming quasi-steady-state activity, one can neglect the complex temporal dynamics; however, in many cases the quasi-steady-state assumption fails. Here, we develop a new reduction method for a general heterogeneous firing-rate model receiving background correlated noisy inputs that accurately handles highly non-equilibrium statistics and interactions of heterogeneous cells. Our method involves solving an efficient set of nonlinear ODEs, rather than time-consuming Monte Carlo simulations or high-dimensional PDEs, and it captures the entire set of first and second order statistics while allowing significant heterogeneity in all model parameters.

Spontaneous behaviors drive multidimensional, brainwide activity

Neuronal populations in sensory cortex produce variable responses to sensory stimuli and exhibit intricate spontaneous activity even without external sensory input. Cortical variability and spontaneous activity have been variously proposed to represent random noise, recall of prior experience, or encoding of ongoing behavioral and cognitive variables. Recording more than 10,000 neurons in mouse visual cortex, we observed that spontaneous activity reliably encoded a high-dimensional latent state, which was partially related to the mouse's ongoing behavior and was represented not just in visual cortex but also across the forebrain. Sensory inputs did not interrupt this ongoing signal but added onto it a representation of external stimuli in orthogonal dimensions. Thus, visual cortical population activity, despite its apparently noisy structure, reliably encodes an orthogonal fusion of sensory and multidimensional behavioral information.

State-dependent geometry of population activity in rat auditory cortex

The accuracy of the neural code depends on the relative embedding of signal and noise in the activity of neural populations. Despite a wealth of theoretical work on population codes, there are few empirical characterizations of the high-dimensional signal and noise subspaces. We studied the geometry of population codes in the rat auditory cortex across brain states along the activation-inactivation continuum, using sounds varying in difference and mean level across the ears. As the cortex becomes more activated, single-hemisphere populations go from preferring contralateral loud sounds to a symmetric preference across lateralizations and intensities, gain-modulation effectively disappears, and the signal and noise subspaces become approximately orthogonal to each other and to the direction corresponding to global activity modulations. Level-invariant decoding of sound lateralization also becomes possible in the active state. Our results provide an empirical foundation for the geometry and state-dependence of cortical population codes.

Leveraging heterogeneity for neural computation with fading memory in layer 2/3 cortical microcircuits

Complexity and heterogeneity are intrinsic to neurobiological systems, manifest in every process, at every scale, and are inextricably linked to the systems' emergent collective behaviours and function. However, the majority of studies addressing the dynamics and computational properties of biologically inspired cortical microcircuits tend to assume (often for the sake of analytical tractability) a great degree of homogeneity in both neuronal and synaptic/connectivity parameters. While simplification and reductionism are necessary to understand the brain's functional principles, disregarding the existence of the multiple heterogeneities in the cortical composition, which may be at the core of its computational proficiency, will inevitably fail to account for important phenomena and limit the scope and generalizability of cortical models. We address these issues by studying the individual and composite functional roles of heterogeneities in neuronal, synaptic and structural properties in a biophysically plausible layer 2/3 microcircuit model, built and constrained by multiple sources of empirical data. This approach was made possible by the emergence of large-scale, well curated databases, as well as the substantial improvements in experimental methodologies achieved over the last few years. Our results show that variability in single neuron parameters is the dominant source of functional specialization, leading to highly proficient microcircuits with much higher computational power than their homogeneous counterparts. We further show that fully heterogeneous circuits, which are closest to the biophysical reality, owe their response properties to the differential contribution of different sources of heterogeneity.

Neuromodulation influences synchronization and intrinsic read-out

Background: The roles of neuromodulation in a neural network, such as in a cortical microcolumn, are still incompletely understood. Neuromodulation influences neural processing by presynaptic and postsynaptic regulation of synaptic efficacy. Neuromodulation also affects ion channels and intrinsic excitability.

Methods: Synaptic efficacy modulation is an effective way to rapidly alter network density and topology. We alter network topology and density to measure the effect on spike synchronization. We also operate with differently parameterized neuron models which alter the neuron's intrinsic excitability, i.e., activation function.

Results: We find that (a) fast synaptic efficacy modulation influences the amount of correlated spiking in a network. Also, (b) synchronization in a network influences the read-out of intrinsic properties. Highly synchronous input drives neurons, such that differences in intrinsic properties disappear, while asynchronous input lets intrinsic properties determine output behavior. Thus, altering network topology can alter the balance between intrinsically vs. synaptically driven network activity.

Conclusion: We conclude that neuromodulation may allow a network to shift between a more synchronized transmission mode and a more asynchronous intrinsic read-out mode. This has significant implications for our understanding of the flexibility of cortical computations.

Neuronal firing rates diverge during REM and homogenize during non-REM

Neurons fire at highly variable intrinsic rates and recent evidence suggests that low- and high-firing rate neurons display different plasticity and dynamics. Furthermore, recent publications imply possibly differing rate-dependent effects in hippocampus versus neocortex, but those analyses were carried out separately and with potentially important differences. To more effectively synthesize these questions, we analyzed the firing rate dynamics of populations of neurons in both hippocampal CA1 and frontal cortex under one framework that avoids the pitfalls of previous analyses and accounts for regression to the mean (RTM). We observed several consistent effects across these regions. While rapid eye movement (REM) sleep was marked by decreased hippocampal firing and increased neocortical firing, in both regions firing rate distributions widened during REM due to differential changes in high- versus low-firing rate cells in parallel with increased interneuron activity. In contrast, upon non-REM (NREM) sleep, firing rate distributions narrowed while interneuron firing decreased. Interestingly, hippocampal interneuron activity closely followed the patterns observed in neocortical principal cells rather than the hippocampal principal cells, suggestive of long-range interactions. Following these undulations in variance, the net effect of sleep was a decrease in firing rates. These decreases were greater in lower-firing hippocampal neurons but also higher-firing frontal cortical neurons, suggestive of greater plasticity in these cell groups. Our results across two different regions, and with statistical corrections, indicate that the hippocampus and neocortex show a mixture of differences and similarities as they cycle between sleep states with a unifying characteristic of homogenization of firing during NREM and diversification during REM.

Bridging large-scale neuronal recordings and large-scale network models using dimensionality reduction

A long-standing goal in neuroscience has been to bring together neuronal recordings and neural network modeling to understand brain function. Neuronal recordings can inform the development of network models, and network models can in turn provide predictions for subsequent experiments. Traditionally, neuronal recordings and network models have been related using single-neuron and pairwise spike train statistics. We review here recent studies that have begun to relate neuronal recordings and network models based on the multi-dimensional structure of neuronal population activity, as identified using dimensionality reduction. This approach has been used to study working memory, decision making, motor control, and more. Dimensionality reduction has provided common ground for incisive comparisons and tight interplay between neuronal recordings and network models.

Using computational models to predict in vivo synaptic inputs to interneuron specific 3 (IS3) cells of CA1 hippocampus that also allow their recruitment during rhythmic states

Brain coding strategies are enabled by the balance of synaptic inputs that individual neurons receive as determined by the networks in which they reside. Inhibitory cell types contribute to brain function in distinct ways but recording from specific, inhibitory cell types during behaviour to determine their contributions is highly challenging. In particular, the in vivo activities of vasoactive intestinal peptide-expressing interneuron specific 3 (IS3) cells in the hippocampus that only target other inhibitory cells are unknown at present. We perform a massive, computational exploration of possible synaptic inputs to IS3 cells using multi-compartment models and optimized synaptic parameters. We find that asynchronous, in vivo-like states that are sensitive to additional theta-timed inputs (8 Hz) exist when excitatory and inhibitory synaptic conductances are approximately equally balanced and with low numbers of activated synapses receiving correlated inputs. Specifically, under these balanced conditions, the input resistance is larger with higher mean spike firing rates relative to other activated synaptic conditions investigated. Incoming theta-timed inputs result in strongly increased spectral power relative to baseline. Thus, using a generally applicable computational approach we predict the existence and features of background, balanced states in hippocampal circuits.

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.

Circuit Models of Low-Dimensional Shared Variability in Cortical Networks

Trial-to-trial variability is a reflection of the circuitry and cellular physiology that make up a neuronal network. A pervasive yet puzzling feature of cortical circuits is that despite their complex wiring, population-wide shared spiking variability is low dimensional. Previous model cortical networks cannot explain this global variability, and rather assume it is from external sources. We show that if the spatial and temporal scales of inhibitory coupling match known physiology, networks of model spiking neurons internally generate low-dimensional shared variability that captures population activity recorded in vivo. Shifting spatial attention into the receptive field of visual neurons has been shown to differentially modulate shared variability within and between brain areas. A top-down modulation of inhibitory neurons in our network provides a parsimonious mechanism for this attentional modulation. Our work provides a critical link between observed cortical circuit structure and realistic shared neuronal variability and its modulation.