The Impact of Structural Heterogeneity on Excitation-Inhibition Balance in Cortical Networks

Top-down modulation in canonical cortical circuits with short-term plasticity

Cortical dynamics and computations are strongly influenced by diverse GABAergic interneurons, including those expressing parvalbumin (PV), somatostatin (SST), and vasoactive intestinal peptide (VIP). Together with excitatory (E) neurons, they form a canonical microcircuit and exhibit counterintuitive nonlinear phenomena. One instance of such phenomena is response reversal, whereby SST neurons show opposite responses to top-down modulation via VIP depending on the presence of bottom-up sensory input, indicating that the network may function in different regimes under different stimulation conditions. Combining analytical and computational approaches, we demonstrate that model networks with multiple interneuron subtypes and experimentally identified short-term plasticity mechanisms can implement response reversal. Surprisingly, despite not directly affecting SST and VIP activity, PV-to-E short-term depression has a decisive impact on SST response reversal. We show how response reversal relates to inhibition stabilization and the paradoxical effect in the presence of several short-term plasticity mechanisms demonstrating that response reversal coincides with a change in the indispensability of SST for network stabilization. In summary, our work suggests a role of short-term plasticity mechanisms in generating nonlinear phenomena in networks with multiple interneuron subtypes and makes several experimentally testable predictions.

A robust balancing mechanism for spiking neural networks

Dynamical balance of excitation and inhibition is usually invoked to explain the irregular low firing activity observed in the cortex. We propose a robust nonlinear balancing mechanism for a random network of spiking neurons, which works also in the absence of strong external currents. Biologically, the mechanism exploits the plasticity of excitatory-excitatory synapses induced by short-term depression. Mathematically, the nonlinear response of the synaptic activity is the key ingredient responsible for the emergence of a stable balanced regime. Our claim is supported by a simple self-consistent analysis accompanied by extensive simulations performed for increasing network sizes. The observed regime is essentially fluctuation driven and characterized by highly irregular spiking dynamics of all neurons.

Network-neuron interactions underlying sensory responses of layer 5 pyramidal tract neurons in barrel cortex

Neurons in the cerebral cortex receive thousands of synaptic inputs per second from thousands of presynaptic neurons. How the dendritic location of inputs, their timing, strength, and presynaptic origin, in conjunction with complex dendritic physiology, impact the transformation of synaptic input into action potential (AP) output remains generally unknown for in vivo conditions. Here, we introduce a computational approach to reveal which properties of the input causally underlie AP output, and how this neuronal input-output computation is influenced by the morphology and biophysical properties of the dendrites. We demonstrate that this approach allows dissecting of how different input populations drive in vivo observed APs. For this purpose, we focus on fast and broadly tuned responses that pyramidal tract neurons in layer 5 (L5PTs) of the rat barrel cortex elicit upon passive single whisker deflections. By reducing a multi-scale model that we reported previously, we show that three features are sufficient to predict with high accuracy the sensory responses and receptive fields of L5PTs under these specific in vivo conditions: the count of active excitatory versus inhibitory synapses preceding the response, their spatial distribution on the dendrites, and the AP history. Based on these three features, we derive an analytically tractable description of the input-output computation of L5PTs, which enabled us to dissect how synaptic input from thalamus and different cell types in barrel cortex contribute to these responses. We show that the input-output computation is preserved across L5PTs despite morphological and biophysical diversity of their dendrites. We found that trial-to-trial variability in L5PT responses, and cell-to-cell variability in their receptive fields, are sufficiently explained by variability in synaptic input from the network, whereas variability in biophysical and morphological properties have minor contributions. Our approach to derive analytically tractable models of input-output computations in L5PTs provides a roadmap to dissect network-neuron interactions underlying L5PT responses across different in vivo conditions and for other cell types.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population's mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn.

Reduced corticospinal drive to antagonist muscles of upper and lower limbs during hands-and-knees crawling in infants with cerebral palsy: Evidence from intermuscular EMG-EMG coherence

BACKGROUND

There is growing interest in understanding the central control of hands-and-knees crawling, especially as a significant motor developmental milestone for early assessment of motor dysfunction in infants with cerebral palsy (CP) who have not yet acquired walking ability. In particular, CP is known to be associated with walking dysfunctions caused by early damage and incomplete maturation of the corticospinal tract. However, the extent of damage to the corticospinal connections during crawling in infants with CP has not been fully clarified. Therefore, this study aimed to investigate the disparities in intermuscular EMG-EMG coherence, which serve as indicators of corticospinal drives to antagonist muscles in the upper and lower limbs during crawling, between infants with and without CP.

METHODS

This study involved 15 infants diagnosed with CP and 20 typically developing (TD) infants. Surface EMG recordings were obtained from two pairs of antagonist muscles in the upper limbs (triceps brachii (TB) and biceps brachii (BB)) and lower limbs (quadriceps femoris (QF) and hamstrings (HS)), while the infants performed hands-and-knees crawling at their self-selected velocity. Intermuscular EMG-EMG coherence was computed in two frequency bands, the beta band (15-30 Hz) and gamma band (30-60 Hz), which indicate corticospinal drive. Additionally, spatiotemporal parameters, including crawling velocity, cadence, duration, and the percentage of stance phase time, were calculated for comparison. Spearman rank correlations were conducted to assess the relationship between EMG-EMG coherence and crawling spatiotemporal parameters.

RESULTS

Infants with CP exhibited significantly reduced crawling velocity, decreased cadence, longer cycle duration, and a higher percentage of stance phase time compared to TD infants. Furthermore, CP infants demonstrated decreased coherence in the beta and gamma frequency bands (indicators of corticospinal drive) in both upper and lower limb muscles. Regarding limb-related differences in the beta and gamma coherence, significant disparities were found between upper and lower limb muscles in TD infants (p < 0.05), but not in infants with CP (p > 0.05). Additionally, significant correlations between coherence metrics and crawling spatiotemporal parameters were identified in the TD group (p < 0.05), while such correlations were not evident in the CP group.

CONCLUSIONS

Our findings suggest that the corticospinal drive may functionally influence the central control of antagonist muscles in the limbs during typical infant crawling. This functional role could be impaired by neurological conditions such as cerebral palsy. The neurophysiological markers of corticospinal drive, specifically intermuscular EMG-EMG coherence during crawling in infants with cerebral palsy, could potentially serve as a tool to assess developmental response to therapy.

On the physiological and structural contributors to the overall balance of excitation and inhibition in local cortical networks

Overall balance of excitation and inhibition in cortical networks is central to their functionality and normal operation. Such orchestrated co-evolution of excitation and inhibition is established through convoluted local interactions between neurons, which are organized by specific network connectivity structures and are dynamically controlled by modulating synaptic activities. Therefore, identifying how such structural and physiological factors contribute to establishment of overall balance of excitation and inhibition is crucial in understanding the homeostatic plasticity mechanisms that regulate the balance. We use biologically plausible mathematical models to extensively study the effects of multiple key factors on overall balance of a network. We characterize a network's baseline balanced state by certain functional properties, and demonstrate how variations in physiological and structural parameters of the network deviate this balance and, in particular, result in transitions in spontaneous activity of the network to high-amplitude slow oscillatory regimes. We show that deviations from the reference balanced state can be continuously quantified by measuring the ratio of mean excitatory to mean inhibitory synaptic conductances in the network. Our results suggest that the commonly observed ratio of the number of inhibitory to the number of excitatory neurons in local cortical networks is almost optimal for their stability and excitability. Moreover, the values of inhibitory synaptic decay time constants and density of inhibitory-to-inhibitory network connectivity are critical to overall balance and stability of cortical networks. However, network stability in our results is sufficiently robust against modulations of synaptic quantal conductances, as required by their role in learning and memory. Our study based on extensive bifurcation analyses thus reveal the functional optimality and criticality of structural and physiological parameters in establishing the baseline operating state of local cortical networks.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Logarithmically scaled, gamma distributed neuronal spiking

Naturally log-scaled quantities abound in the nervous system. Distributions of these quantities have non-intuitive properties, which have implications for data analysis and the understanding of neural circuits. Here, we review the log-scaled statistics of neuronal spiking and the relevant analytical probability distributions. Recent work using log-scaling revealed that interspike intervals of forebrain neurons segregate into discrete modes reflecting spiking at different timescales and are each well-approximated by a gamma distribution. Each neuron spends most of the time in an irregular spiking 'ground state' with the longest intervals, which determines the mean firing rate of the neuron. Across the entire neuronal population, firing rates are log-scaled and well approximated by the gamma distribution, with a small number of highly active neurons and an overabundance of low rate neurons (the 'dark matter'). These results are intricately linked to a heterogeneous balanced operating regime, which confers upon neuronal circuits multiple computational advantages and has evolutionarily ancient origins.

On the physiological and structural contributors to the overall balance of excitation and inhibition in local cortical networks

Overall balance of excitation and inhibition in cortical networks is central to their functionality and normal operation. Such orchestrated co-evolution of excitation and inhibition is established through convoluted local interactions between neurons, which are organized by specific network connectivity structures and are dynamically controlled by modulating synaptic activities. Therefore, identifying how such structural and physiological factors contribute to establishment of overall balance of excitation and inhibition is crucial in understanding the homeostatic plasticity mechanisms that regulate the balance. We use biologically plausible mathematical models to extensively study the effects of multiple key factors on overall balance of a network. We characterize a network's baseline balanced state by certain functional properties, and demonstrate how variations in physiological and structural parameters of the network deviate this balance and, in particular, result in transitions in spontaneous activity of the network to high-amplitude slow oscillatory regimes. We show that deviations from the reference balanced state can be continuously quantified by measuring the ratio of mean excitatory to mean inhibitory synaptic conductances in the network. Our results suggest that the commonly observed ratio of the number of inhibitory to the number of excitatory neurons in local cortical networks is almost optimal for their stability and excitability. Moreover, the values of inhibitory synaptic decay time constants and density of inhibitory-to-inhibitory network connectivity are critical to overall balance and stability of cortical networks. However, network stability in our results is sufficiently robust against modulations of synaptic quantal conductances, as required by their role in learning and memory.

Predicting in vitro single-neuron firing rates upon pharmacological perturbation using Graph Neural Networks

Modern Graph Neural Networks (GNNs) provide opportunities to study the determinants underlying the complex activity patterns of biological neuronal networks. In this study, we applied GNNs to a large-scale electrophysiological dataset of rodent primary neuronal networks obtained by means of high-density microelectrode arrays (HD-MEAs). HD-MEAs allow for long-term recording of extracellular spiking activity of individual neurons and networks and enable the extraction of physiologically relevant features at the single-neuron and population level. We employed established GNNs to generate a combined representation of single-neuron and connectivity features obtained from HD-MEA data, with the ultimate goal of predicting changes in single-neuron firing rate induced by a pharmacological perturbation. The aim of the main prediction task was to assess whether single-neuron and functional connectivity features, inferred under baseline conditions, were informative for predicting changes in neuronal activity in response to a perturbation with Bicuculline, a GABA A receptor antagonist. Our results suggest that the joint representation of node features and functional connectivity, extracted from a baseline recording, was informative for predicting firing rate changes of individual neurons after the perturbation. Specifically, our implementation of a GNN model with inductive learning capability (GraphSAGE) outperformed other prediction models that relied only on single-neuron features. We tested the generalizability of the results on two additional datasets of HD-MEA recordings-a second dataset with cultures perturbed with Bicuculline and a dataset perturbed with the GABA A receptor antagonist Gabazine. GraphSAGE models showed improved prediction accuracy over other prediction models. Our results demonstrate the added value of taking into account the functional connectivity between neurons and the potential of GNNs to study complex interactions between neurons.

Circuit-selective cell-autonomous regulation of inhibition in pyramidal neurons by Ste20-like kinase

Maintaining an appropriate balance between excitation and inhibition is critical for neuronal information processing. Cortical neurons can cell-autonomously adjust the inhibition they receive to individual levels of excitatory input, but the underlying mechanisms are unclear. We describe that Ste20-like kinase (SLK) mediates cell-autonomous regulation of excitation-inhibition balance in the thalamocortical feedforward circuit, but not in the feedback circuit. This effect is due to regulation of inhibition originating from parvalbumin-expressing interneurons, while inhibition via somatostatin-expressing interneurons is unaffected. Computational modeling shows that this mechanism promotes stable excitatory-inhibitory ratios across pyramidal cells and ensures robust and sparse coding. Patch-clamp RNA sequencing yields genes differentially regulated by SLK knockdown, as well as genes associated with excitation-inhibition balance participating in transsynaptic communication and cytoskeletal dynamics. These data identify a mechanism for cell-autonomous regulation of a specific inhibitory circuit that is critical to ensure that a majority of cortical pyramidal cells participate in information coding.

A circuit mechanism for independent modulation of excitatory and inhibitory firing rates after sensory deprivation

The cortex is particularly vulnerable to perturbations during sensitive periods, such as the critical period when manipulating sensory experience can induce long-lasting changes in brain structure. Depriving rodents of vision in one eye (known as monocular deprivation [MD]) reduces network activity over two days, whereby inhibitory neurons decrease their firing rates one day after MD, while excitatory neurons are delayed by an additional day. We use spiking networks to mechanistically dissect the requirements for this independent firing-rate regulation after sensory deprivation. We find that in networks stabilized by recurrent inhibition, at least two interneuron subtypes (parvalbumin-expressing and somatostatin-expressing interneurons) are necessary to dynamically alter the circuit response after deprivation and generalize the result across sensory cortices.

Diverse interneuron subtypes shape sensory processing in mature cortical circuits. During development, sensory deprivation evokes powerful synaptic plasticity that alters circuitry, but how different inhibitory subtypes modulate circuit dynamics in response to this plasticity remains unclear. We investigate how deprivation-induced synaptic changes affect excitatory and inhibitory firing rates in a microcircuit model of the sensory cortex with multiple interneuron subtypes. We find that with a single interneuron subtype (parvalbumin-expressing [PV]), excitatory and inhibitory firing rates can only be comodulated—increased or decreased together. To explain the experimentally observed independent modulation, whereby one firing rate increases and the other decreases, requires strong feedback from a second interneuron subtype (somatostatin-expressing [SST]). Our model applies to the visual and somatosensory cortex, suggesting a general mechanism across sensory cortices. Therefore, we provide a mechanistic explanation for the differential role of interneuron subtypes in regulating firing rates, contributing to the already diverse roles they serve in the cortex.

Evaluating the extent to which homeostatic plasticity learns to compute prediction errors in unstructured neuronal networks

The brain is believed to operate in part by making predictions about sensory stimuli and encoding deviations from these predictions in the activity of "prediction error neurons." This principle defines the widely influential theory of predictive coding. The precise circuitry and plasticity mechanisms through which animals learn to compute and update their predictions are unknown. Homeostatic inhibitory synaptic plasticity is a promising mechanism for training neuronal networks to perform predictive coding. Homeostatic plasticity causes neurons to maintain a steady, baseline firing rate in response to inputs that closely match the inputs on which a network was trained, but firing rates can deviate away from this baseline in response to stimuli that are mismatched from training. We combine computer simulations and mathematical analysis systematically to test the extent to which randomly connected, unstructured networks compute prediction errors after training with homeostatic inhibitory synaptic plasticity. We find that homeostatic plasticity alone is sufficient for computing prediction errors for trivial time-constant stimuli, but not for more realistic time-varying stimuli. We use a mean-field theory of plastic networks to explain our findings and characterize the assumptions under which they apply.

The impact of neuron morphology on cortical network architecture

The neurons in the cerebral cortex are not randomly interconnected. This specificity in wiring can result from synapse formation mechanisms that connect neurons, depending on their electrical activity and genetically defined identity. Here, we report that the morphological properties of the neurons provide an additional prominent source by which wiring specificity emerges in cortical networks. This morphologically determined wiring specificity reflects similarities between the neurons’ axo-dendritic projections patterns, the packing density, and the cellular diversity of the neuropil. The higher these three factors are, the more recurrent is the topology of the network. Conversely, the lower these factors are, the more feedforward is the network’s topology. These principles predict the empirically observed occurrences of clusters of synapses, cell type-specific connectivity patterns, and nonrandom network motifs. Thus, we demonstrate that wiring specificity emerges in the cerebral cortex at subcellular, cellular, and network scales from the specific morphological properties of its neuronal constituents.

•

Neuronal network architectures reflect the morphologies of their constituents

•

Morphology predicts nonrandom connectivity from subcellular to network scales

•

Morphology predicts connectivity patterns consistent with those observed empirically

•

Neuron morphology is a major source for wiring specificity in the cerebral cortex

Computational synthesis of cortical dendritic morphologies

Neuronal morphologies provide the foundation for the electrical behavior of neurons, the connectomes they form, and the dynamical properties of the brain. Comprehensive neuron models are essential for defining cell types, discerning their functional roles, and investigating brain-disease-related dendritic alterations. However, a lack of understanding of the principles underlying neuron morphologies has hindered attempts to computationally synthesize morphologies for decades. We introduce a synthesis algorithm based on a topological descriptor of neurons, which enables the rapid digital reconstruction of entire brain regions from few reference cells. This topology-guided synthesis generates dendrites that are statistically similar to biological reconstructions in terms of morpho-electrical and connectivity properties and offers a significant opportunity to investigate the links between neuronal morphology and brain function across different spatiotemporal scales. Synthesized cortical networks based on structurally altered dendrites associated with diverse brain pathologies revealed principles linking branching properties to the structure of large-scale networks.

Chaotic and stochastic dynamics of epileptiform-like activities in sclerotic hippocampus resected from patients with pharmacoresistant epilepsy

The types of epileptiform activity occurring in the sclerotic hippocampus with highest incidence are interictal-like events (II) and periodic ictal spiking (PIS). These activities are classified according to their event rates, but it is still unclear if these rate differences are consequences of underlying physiological mechanisms. Identifying new and more specific information related to these two activities may bring insights to a better understanding about the epileptogenic process and new diagnosis. We applied Poincaré map analysis and Recurrence Quantification Analysis (RQA) onto 35 in vitro electrophysiological signals recorded from slices of 12 hippocampal tissues surgically resected from patients with pharmacoresistant temporal lobe epilepsy. These analyzes showed that the II activity is related to chaotic dynamics, whereas the PIS activity is related to deterministic periodic dynamics. Additionally, it indicates that their different rates are consequence of different endogenous dynamics. Finally, by using two computational models we were able to simulate the transition between II and PIS activities. The RQA was applied to different periods of these simulations to compare the recurrences between artificial and real signals, showing that different ranges of regularity-chaoticity can be directly associated with the generation of PIS and II activities.

Temporal lobe epilepsy (TLE) is the most prevalent type of epilepsy in adults and hippocampal sclerosis is the major pathophysiological substrate of pharmaco-refractory TLE. Different patterns of epileptiform-like activity have been described in human hippocampal sclerosis, but the standard analysis applied to characterize the activities usually do not consider the nonlinear features that epileptiform patterns exhibit. Here, using Poincaré map and Recurrence Quantitative Analysis we characterized the most prevalent type of epileptiform-like activities—interictal-like events (II) and periodic ictal spiking (PIS), recorded in vitro from resected hippocampi of pharmacoresistant patients with TLE—according to their levels of stochasticity, chaoticity and determinism. The II activities showed to be more chaotic with complex rhythmicity than PIS activities. The nonlinear dynamic differences between II and PIS leads us to conjecture that they are expressions of different seizure susceptibility. We also identified that each hippocampal subfield expresses II and PIS activities in a specific and different way. Finally, from the modulation of internal parameters of two computational models, we show the conversion of one type of activity into the other, showing how specific neuron networks synchronize over time, leading to II and PIS activities and then into a generalized seizure.

Synaptic Variability and Cortical Gamma Oscillation Power in Schizophrenia

OBJECTIVE

Cognitive impairments in schizophrenia are associated with lower gamma oscillation power in the prefrontal cortex (PFC). Gamma power depends in part on excitatory drive to fast-spiking parvalbumin interneurons (PVIs). Excitatory drive to cortical neurons varies in strength, which could affect how these neurons regulate network oscillations. The authors investigated whether variability in excitatory synaptic strength across PVIs could contribute to lower prefrontal gamma power in schizophrenia.

METHODS

In postmortem PFC from 20 matched pairs of comparison and schizophrenia subjects, levels of vesicular glutamate transporter 1 (VGlut1) and postsynaptic density 95 (PSD95) proteins were quantified to assess variability in excitatory synaptic strength across PVIs. A computational model network was then used to simulate how variability in excitatory synaptic strength across fast-spiking (a defining feature of PVIs) interneurons (FSIs) regulates gamma power.

RESULTS

The variability of VGlut1 and PSD95 levels at excitatory inputs across PVIs was larger in schizophrenia relative to comparison subjects. This alteration was not influenced by schizophrenia-associated comorbid factors, was not present in monkeys chronically exposed to antipsychotic medications, and was not present in calretinin interneurons. In the model network, variability in excitatory synaptic strength across FSIs regulated gamma power by affecting network synchrony. Finally, greater synaptic variability interacted synergistically with other synaptic alterations in schizophrenia (i.e., fewer excitatory inputs to FSIs and lower inhibitory strength from FSIs) to robustly reduce gamma power.

CONCLUSIONS

The study findings suggest that greater variability in excitatory synaptic strength across PVIs, in combination with other modest synaptic alterations in these neurons, can markedly lower PFC gamma power in schizophrenia.

Coherent oscillations in balanced neural networks driven by endogenous fluctuations

We present a detailed analysis of the dynamical regimes observed in a balanced network of identical quadratic integrate-and-fire neurons with sparse connectivity for homogeneous and heterogeneous in-degree distributions. Depending on the parameter values, either an asynchronous regime or periodic oscillations spontaneously emerge. Numerical simulations are compared with a mean-field model based on a self-consistent Fokker-Planck equation (FPE). The FPE reproduces quite well the asynchronous dynamics in the homogeneous case by either assuming a Poissonian or renewal distribution for the incoming spike trains. An exact self-consistent solution for the mean firing rate obtained in the limit of infinite in-degree allows identifying balanced regimes that can be either mean- or fluctuation-driven. A low-dimensional reduction of the FPE in terms of circular cumulants is also considered. Two cumulants suffice to reproduce the transition scenario observed in the network. The emergence of periodic collective oscillations is well captured both in the homogeneous and heterogeneous setups by the mean-field models upon tuning either the connectivity or the input DC current. In the heterogeneous situation, we analyze also the role of structural heterogeneity.

Excitatory-Inhibitory Homeostasis and Diaschisis: Tying the Local and Global Scales in the Post-stroke Cortex

Maintaining a balance between excitatory and inhibitory activity is an essential feature of neural networks of the neocortex. In the face of perturbations in the levels of excitation to cortical neurons, synapses adjust to maintain excitatory-inhibitory (EI) balance. In this review, we summarize research on this EI homeostasis in the neocortex, using stroke as our case study, and in particular the loss of excitation to distant cortical regions after focal lesions. Widespread changes following a localized lesion, a phenomenon known as diaschisis, are not only related to excitability, but also observed with respect to functional connectivity. Here, we highlight the main findings regarding the evolution of excitability and functional cortical networks during the process of post-stroke recovery, and how both are related to functional recovery. We show that cortical reorganization at a global scale can be explained from the perspective of EI homeostasis. Indeed, recovery of functional networks is paralleled by increases in excitability across the cortex. These adaptive changes likely result from plasticity mechanisms such as synaptic scaling and are linked to EI homeostasis, providing a possible target for future therapeutic strategies in the process of rehabilitation. In addition, we address the difficulty of simultaneously studying these multiscale processes by presenting recent advances in large-scale modeling of the human cortex in the contexts of stroke and EI homeostasis, suggesting computational modeling as a powerful tool to tie the meso- and macro-scale processes of recovery in stroke patients.

Hyperexcitability and Homeostasis in Fragile X Syndrome

Fragile X Syndrome (FXS) is a leading inherited cause of autism and intellectual disability, resulting from a mutation in the FMR1 gene and subsequent loss of its protein product FMRP. Despite this simple genetic origin, FXS is a phenotypically complex disorder with a range of physical and neurocognitive disruptions. While numerous molecular and cellular pathways are affected by FMRP loss, there is growing evidence that circuit hyperexcitability may be a common convergence point that can account for many of the wide-ranging phenotypes seen in FXS. The mechanisms for hyperexcitability in FXS include alterations to excitatory synaptic function and connectivity, reduced inhibitory neuron activity, as well as changes to ion channel expression and conductance. However, understanding the impact of FMR1 mutation on circuit function is complicated by the inherent plasticity in neural circuits, which display an array of homeostatic mechanisms to maintain activity near set levels. FMRP is also an important regulator of activity-dependent plasticity in the brain, meaning that dysregulated plasticity can be both a cause and consequence of hyperexcitable networks in FXS. This makes it difficult to separate the direct effects of FMR1 mutation from the myriad and pleiotropic compensatory changes associated with it, both of which are likely to contribute to FXS pathophysiology. Here we will: (1) review evidence for hyperexcitability and homeostatic plasticity phenotypes in FXS models, focusing on similarities/differences across brain regions, cell-types, and developmental time points; (2) examine how excitability and plasticity disruptions interact with each other to ultimately contribute to circuit dysfunction in FXS; and (3) discuss how these synaptic and circuit deficits contribute to disease-relevant behavioral phenotypes like epilepsy and sensory hypersensitivity. Through this discussion of where the current field stands, we aim to introduce perspectives moving forward in FXS research.

Cortical Representation of Touch in Silico

With its six layers and ~ 12,000 neurons, a cortical column is a complex network whose function is plausibly greater than the sum of its constituents'. Functional characterization of its network components will require going beyond the brute-force modulation of the neural activity of a small group of neurons. Here we introduce an open-source, biologically inspired, computationally efficient network model of the somatosensory cortex's granular and supragranular layers after reconstructing the barrel cortex in soma resolution. Comparisons of the network activity to empirical observations showed that the in silico network replicates the known properties of touch representations and whisker deprivation-induced changes in synaptic strength induced in vivo. Simulations show that the history of the membrane potential acts as a spatial filter that determines the presynaptic population of neurons contributing to a post-synaptic action potential; this spatial filtering might be critical for synaptic integration of top-down and bottom-up information.

Emergence of Irregular Activity in Networks of Strongly Coupled Conductance-Based Neurons

Cortical neurons are characterized by irregular firing and a broad distribution of rates. The balanced state model explains these observations with a cancellation of mean excitatory and inhibitory currents, which makes fluctuations drive firing. In networks of neurons with current-based synapses, the balanced state emerges dynamically if coupling is strong, i.e., if the mean number of synapses per neuron K is large and synaptic efficacy is of the order of 1 / K . When synapses are conductance-based, current fluctuations are suppressed when coupling is strong, questioning the applicability of the balanced state idea to biological neural networks. We analyze networks of strongly coupled conductance-based neurons and show that asynchronous irregular activity and broad distributions of rates emerge if synaptic efficacy is of the order of 1/ log(K). In such networks, unlike in the standard balanced state model, current fluctuations are small and firing is maintained by a drift-diffusion balance. This balance emerges dynamically, without fine-tuning, if inputs are smaller than a critical value, which depends on synaptic time constants and coupling strength, and is significantly more robust to connection heterogeneities than the classical balanced state model. Our analysis makes experimentally testable predictions of how the network response properties should evolve as input increases.

Asynchronous and Coherent Dynamics in Balanced Excitatory-Inhibitory Spiking Networks

Dynamic excitatory-inhibitory (E-I) balance is a paradigmatic mechanism invoked to explain the irregular low firing activity observed in the cortex. However, we will show that the E-I balance can be at the origin of other regimes observable in the brain. The analysis is performed by combining extensive simulations of sparse E-I networks composed of N spiking neurons with analytical investigations of low dimensional neural mass models. The bifurcation diagrams, derived for the neural mass model, allow us to classify the possible asynchronous and coherent behaviors emerging in balanced E-I networks with structural heterogeneity for any finite in-degree K. Analytic mean-field (MF) results show that both supra and sub-threshold balanced asynchronous regimes are observable in our system in the limit N >> K >> 1. Due to the heterogeneity, the asynchronous states are characterized at the microscopic level by the splitting of the neurons in to three groups: silent, fluctuation, and mean driven. These features are consistent with experimental observations reported for heterogeneous neural circuits. The coherent rhythms observed in our system can range from periodic and quasi-periodic collective oscillations (COs) to coherent chaos. These rhythms are characterized by regular or irregular temporal fluctuations joined to spatial coherence somehow similar to coherent fluctuations observed in the cortex over multiple spatial scales. The COs can emerge due to two different mechanisms. A first mechanism analogous to the pyramidal-interneuron gamma (PING), usually invoked for the emergence of γ-oscillations. The second mechanism is intimately related to the presence of current fluctuations, which sustain COs characterized by an essentially simultaneous bursting of the two populations. We observe period-doubling cascades involving the PING-like COs finally leading to the appearance of coherent chaos. Fluctuation driven COs are usually observable in our system as quasi-periodic collective motions characterized by two incommensurate frequencies. However, for sufficiently strong current fluctuations these collective rhythms can lock. This represents a novel mechanism of frequency locking in neural populations promoted by intrinsic fluctuations. COs are observable for any finite in-degree K, however, their existence in the limit N >> K >> 1 appears as uncertain.

Self-organization of a doubly asynchronous irregular network state for spikes and bursts

Cortical pyramidal cells (PCs) have a specialized dendritic mechanism for the generation of bursts, suggesting that these events play a special role in cortical information processing. In vivo, bursts occur at a low, but consistent rate. Theory suggests that this network state increases the amount of information they convey. However, because burst activity relies on a threshold mechanism, it is rather sensitive to dendritic input levels. In spiking network models, network states in which bursts occur rarely are therefore typically not robust, but require fine-tuning. Here, we show that this issue can be solved by a homeostatic inhibitory plasticity rule in dendrite-targeting interneurons that is consistent with experimental data. The suggested learning rule can be combined with other forms of inhibitory plasticity to self-organize a network state in which both spikes and bursts occur asynchronously and irregularly at low rate. Finally, we show that this network state creates the network conditions for a recently suggested multiplexed code and thereby indeed increases the amount of information encoded in bursts.

Heterogeneous side effects of cortical inactivation in behaving animals

Cortical inactivation represents a key causal manipulation allowing the study of cortical circuits and their impact on behavior. A key assumption in inactivation studies is that the neurons in the target area become silent while the surrounding cortical tissue is only negligibly impacted. However, individual neurons are embedded in complex local circuits composed of excitatory and inhibitory cells with connections extending hundreds of microns. This raises the possibility that silencing one part of the network could induce complex, unpredictable activity changes in neurons outside the targeted inactivation zone. These off-target side effects can potentially complicate interpretations of inactivation manipulations, especially when they are related to changes in behavior. Here, we demonstrate that optogenetic inactivation of glutamatergic neurons in the superficial layers of monkey primary visual cortex (V1) induces robust suppression at the light-targeted site, but destabilizes stimulus responses in the neighboring, untargeted network. We identified four types of stimulus-evoked neuronal responses within a cortical column, ranging from full suppression to facilitation, and a mixture of both. Mixed responses were most prominent in middle and deep cortical layers. These results demonstrate that response modulation driven by lateral network connectivity is diversely implemented throughout a cortical column. Importantly, consistent behavioral changes induced by optogenetic inactivation were only achieved when cumulative network activity was homogeneously suppressed. Therefore, careful consideration of the full range of network changes outside the inactivated cortical region is required, as heterogeneous side effects can confound interpretation of inactivation experiments.

Optimal responsiveness and information flow in networks of heterogeneous neurons

Cerebral cortex is characterized by a strong neuron-to-neuron heterogeneity, but it is unclear what consequences this may have for cortical computations, while most computational models consider networks of identical units. Here, we study network models of spiking neurons endowed with heterogeneity, that we treat independently for excitatory and inhibitory neurons. We find that heterogeneous networks are generally more responsive, with an optimal responsiveness occurring for levels of heterogeneity found experimentally in different published datasets, for both excitatory and inhibitory neurons. To investigate the underlying mechanisms, we introduce a mean-field model of heterogeneous networks. This mean-field model captures optimal responsiveness and suggests that it is related to the stability of the spontaneous asynchronous state. The mean-field model also predicts that new dynamical states can emerge from heterogeneity, a prediction which is confirmed by network simulations. Finally we show that heterogeneous networks maximise the information flow in large-scale networks, through recurrent connections. We conclude that neuronal heterogeneity confers different responsiveness to neural networks, which should be taken into account to investigate their information processing capabilities.

Balanced networks under spike-time dependent plasticity

The dynamics of local cortical networks are irregular, but correlated. Dynamic excitatory–inhibitory balance is a plausible mechanism that generates such irregular activity, but it remains unclear how balance is achieved and maintained in plastic neural networks. In particular, it is not fully understood how plasticity induced changes in the network affect balance, and in turn, how correlated, balanced activity impacts learning. How do the dynamics of balanced networks change under different plasticity rules? How does correlated spiking activity in recurrent networks change the evolution of weights, their eventual magnitude, and structure across the network? To address these questions, we develop a theory of spike–timing dependent plasticity in balanced networks. We show that balance can be attained and maintained under plasticity–induced weight changes. We find that correlations in the input mildly affect the evolution of synaptic weights. Under certain plasticity rules, we find an emergence of correlations between firing rates and synaptic weights. Under these rules, synaptic weights converge to a stable manifold in weight space with their final configuration dependent on the initial state of the network. Lastly, we show that our framework can also describe the dynamics of plastic balanced networks when subsets of neurons receive targeted optogenetic input.

Animals are able to learn complex tasks through changes in individual synapses between cells. Such changes lead to the coevolution of neural activity patterns and the structure of neural connectivity, but the consequences of these interactions are not fully understood. We consider plasticity in model neural networks which achieve an average balance between the excitatory and inhibitory synaptic inputs to different cells, and display cortical–like, irregular activity. We extend the theory of balanced networks to account for synaptic plasticity and show which rules can maintain balance, and which will drive the network into a different state. This theory of plasticity can provide insights into the relationship between stimuli, network dynamics, and synaptic circuitry.

Learning excitatory-inhibitory neuronal assemblies in recurrent networks

Understanding the connectivity observed in the brain and how it emerges from local plasticity rules is a grand challenge in modern neuroscience. In the primary visual cortex (V1) of mice, synapses between excitatory pyramidal neurons and inhibitory parvalbumin-expressing (PV) interneurons tend to be stronger for neurons that respond to similar stimulus features, although these neurons are not topographically arranged according to their stimulus preference. The presence of such excitatory-inhibitory (E/I) neuronal assemblies indicates a stimulus-specific form of feedback inhibition. Here, we show that activity-dependent synaptic plasticity on input and output synapses of PV interneurons generates a circuit structure that is consistent with mouse V1. Computational modeling reveals that both forms of plasticity must act in synergy to form the observed E/I assemblies. Once established, these assemblies produce a stimulus-specific competition between pyramidal neurons. Our model suggests that activity-dependent plasticity can refine inhibitory circuits to actively shape cortical computations.

Neuronal heterogeneity modulates phase synchronization between unidirectionally coupled populations with excitation-inhibition balance

Several experiments and models have highlighted the importance of neuronal heterogeneity in brain dynamics and function. However, how such a cell-to-cell diversity can affect cortical computation, synchronization, and neuronal communication is still under debate. Previous studies have focused on the effect of neuronal heterogeneity in one neuronal population. Here we are specifically interested in the effect of neuronal variability on the phase relations between two populations, which can be related to different cortical communication hypotheses. It has been recently shown that two spiking neuron populations unidirectionally connected in a sender-receiver configuration can exhibit anticipated synchronization (AS), which is characterized by a negative phase lag. This phenomenon has been reported in electrophysiological data of nonhuman primates and human EEG during a visual discrimination cognitive task. In experiments, the unidirectional coupling could be accessed by Granger causality and can be accompanied by either positive or negative phase difference between cortical areas. Here we propose a model of two coupled populations in which the neuronal heterogeneity can determine the dynamical relation between the sender and the receiver and can reproduce phase relations reported in experiments. Depending on the distribution of parameters characterizing the neuronal firing patterns, the system can exhibit both AS and the usual delayed synchronization regime (DS, with positive phase) as well as a zero-lag synchronization regime and phase bistability between AS and DS. Furthermore, we show that our network can present diversity in their phase relations maintaining the excitation-inhibition balance.

Hopf Bifurcation in Mean Field Explains Critical Avalanches in Excitation-Inhibition Balanced Neuronal Networks: A Mechanism for Multiscale Variability

Cortical neural circuits display highly irregular spiking in individual neurons but variably sized collective firing, oscillations and critical avalanches at the population level, all of which have functional importance for information processing. Theoretically, the balance of excitation and inhibition inputs is thought to account for spiking irregularity and critical avalanches may originate from an underlying phase transition. However, the theoretical reconciliation of these multilevel dynamic aspects in neural circuits remains an open question. Herein, we study excitation-inhibition (E-I) balanced neuronal network with biologically realistic synaptic kinetics. It can maintain irregular spiking dynamics with different levels of synchrony and critical avalanches emerge near the synchronous transition point. We propose a novel semi-analytical mean-field theory to derive the field equations governing the network macroscopic dynamics. It reveals that the E-I balanced state of the network manifesting irregular individual spiking is characterized by a macroscopic stable state, which can be either a fixed point or a periodic motion and the transition is predicted by a Hopf bifurcation in the macroscopic field. Furthermore, by analyzing public data, we find the coexistence of irregular spiking and critical avalanches in the spontaneous spiking activities of mouse cortical slice in vitro, indicating the universality of the observed phenomena. Our theory unveils the mechanism that permits complex neural activities in different spatiotemporal scales to coexist and elucidates a possible origin of the criticality of neural systems. It also provides a novel tool for analyzing the macroscopic dynamics of E-I balanced networks and its relationship to the microscopic counterparts, which can be useful for large-scale modeling and computation of cortical dynamics.

Nonlinear stimulus representations in neural circuits with approximate excitatory-inhibitory balance

Balanced excitation and inhibition is widely observed in cortex. How does this balance shape neural computations and stimulus representations? This question is often studied using computational models of neuronal networks in a dynamically balanced state. But balanced network models predict a linear relationship between stimuli and population responses. So how do cortical circuits implement nonlinear representations and computations? We show that every balanced network architecture admits stimuli that break the balanced state and these breaks in balance push the network into a "semi-balanced state" characterized by excess inhibition to some neurons, but an absence of excess excitation. The semi-balanced state produces nonlinear stimulus representations and nonlinear computations, is unavoidable in networks driven by multiple stimuli, is consistent with cortical recordings, and has a direct mathematical relationship to artificial neural networks.

Anatomy and Physiology of Macaque Visual Cortical Areas V1, V2, and V5/MT: Bases for Biologically Realistic Models

The cerebral cortex of primates encompasses multiple anatomically and physiologically distinct areas processing visual information. Areas V1, V2, and V5/MT are conserved across mammals and are central for visual behavior. To facilitate the generation of biologically accurate computational models of primate early visual processing, here we provide an overview of over 350 published studies of these three areas in the genus Macaca, whose visual system provides the closest model for human vision. The literature reports 14 anatomical connection types from the lateral geniculate nucleus of the thalamus to V1 having distinct layers of origin or termination, and 194 connection types between V1, V2, and V5, forming multiple parallel and interacting visual processing streams. Moreover, within V1, there are reports of 286 and 120 types of intrinsic excitatory and inhibitory connections, respectively. Physiologically, tuning of neuronal responses to 11 types of visual stimulus parameters has been consistently reported. Overall, the optimal spatial frequency (SF) of constituent neurons decreases with cortical hierarchy. Moreover, V5 neurons are distinct from neurons in other areas for their higher direction selectivity, higher contrast sensitivity, higher temporal frequency tuning, and wider SF bandwidth. We also discuss currently unavailable data that could be useful for biologically accurate models.

Spatially extended balanced networks without translationally invariant connectivity

Networks of neurons in the cerebral cortex exhibit a balance between excitation (positive input current) and inhibition (negative input current). Balanced network theory provides a parsimonious mathematical model of this excitatory-inhibitory balance using randomly connected networks of model neurons in which balance is realized as a stable fixed point of network dynamics in the limit of large network size. Balanced network theory reproduces many salient features of cortical network dynamics such as asynchronous-irregular spiking activity. Early studies of balanced networks did not account for the spatial topology of cortical networks. Later works introduced spatial connectivity structure, but were restricted to networks with translationally invariant connectivity structure in which connection probability depends on distance alone and boundaries are assumed to be periodic. Spatial connectivity structure in cortical network does not always satisfy these assumptions. We use the mathematical theory of integral equations to extend the mean-field theory of balanced networks to account for more general dependence of connection probability on the spatial location of pre- and postsynaptic neurons. We compare our mathematical derivations to simulations of large networks of recurrently connected spiking neuron models.

Modeling Hippocampal CA1 Gabaergic Synapses of Audiogenic Rats

Wistar Audiogenic Rats (WARs) are genetically susceptible to sound-induced seizures that start in the brainstem and, in response to repetitive stimulation, spread to limbic areas, such as hippocampus. Analysis of the distribution of interevent intervals of GABAergic inhibitory postsynaptic currents (IPSCs) in CA1 pyramidal cells showed a monoexponential trend in Wistar rats, suggestive of a homogeneous population of synapses, but a biexponential trend in WARs. Based on this, we hypothesize that there are two populations of GABAergic synaptic release sites in CA1 pyramidal neurons from WARs. To address this hypothesis, we used a well-established neuronal computational model of a CA1 pyramidal neuron previously developed to replicate physiological properties of these cells. Our simulations replicated the biexponential trend only when we decreased the release frequency of synaptic currents by a factor of six in at least 40% of distal synapses. Our results suggest that almost half of the GABAergic synapses of WARs have a drastically reduced spontaneous release frequency. The computational model was able to reproduce the temporal dynamics of GABAergic inhibition that could underlie susceptibility to the spread of seizures.

Inhibitory Plasticity: From Molecules to Computation and Beyond

Synaptic plasticity is the cellular and molecular counterpart of learning and memory and, since its first discovery, the analysis of the mechanisms underlying long-term changes of synaptic strength has been almost exclusively focused on excitatory connections. Conversely, inhibition was considered as a fixed controller of circuit excitability. Only recently, inhibitory networks were shown to be finely regulated by a wide number of mechanisms residing in their synaptic connections. Here, we review recent findings on the forms of inhibitory plasticity (IP) that have been discovered and characterized in different brain areas. In particular, we focus our attention on the molecular pathways involved in the induction and expression mechanisms leading to changes in synaptic efficacy, and we discuss, from the computational perspective, how IP can contribute to the emergence of functional properties of brain circuits.

Impact of higher order network structure on emergent cortical activity

Synaptic connectivity between neocortical neurons is highly structured. The network structure of synaptic connectivity includes first-order properties that can be described by pairwise statistics, such as strengths of connections between different neuron types and distance-dependent connectivity, and higher order properties, such as an abundance of cliques of all-to-all connected neurons. The relative impact of first- and higher order structure on emergent cortical network activity is unknown. Here, we compare network structure and emergent activity in two neocortical microcircuit models with different synaptic connectivity. Both models have a similar first-order structure, but only one model includes higher order structure arising from morphological diversity within neuronal types. We find that such morphological diversity leads to more heterogeneous degree distributions, increases the number of cliques, and contributes to a small-world topology. The increase in higher order network structure is accompanied by more nuanced changes in neuronal firing patterns, such as an increased dependence of pairwise correlations on the positions of neurons in cliques. Our study shows that circuit models with very similar first-order structure of synaptic connectivity can have a drastically different higher order network structure, and suggests that the higher order structure imposed by morphological diversity within neuronal types has an impact on emergent cortical activity.

Unstructured network topology begets order-based representation by privileged neurons

How spiking activity reverberates through neuronal networks, how evoked and spontaneous activity interacts and blends, and how the combined activities represent external stimulation are pivotal questions in neuroscience. We simulated minimal models of unstructured spiking networks in silico, asking whether and how gentle external stimulation might be subsequently reflected in spontaneous activity fluctuations. Consistent with earlier findings in silico and in vitro, we observe a privileged subpopulation of 'pioneer neurons' that, by their firing order, reliably encode previous external stimulation. We also confirm that pioneer neurons are 'sensitive' in that they are recruited by small fluctuations of population activity. We show that order-based representations rely on a 'chain' of pioneer neurons with different degrees of sensitivity and thus constitute an emergent property of collective dynamics. The forming of such representations is greatly favoured by a broadly heterogeneous connection topology-a broad 'middle class' in degree of connectedness. In conclusion, we offer a minimal model for the representational role of pioneer neurons, as observed experimentally in vitro. In addition, we show that broadly heterogeneous connectivity enhances the representational capacity of unstructured networks.

Network and cellular mechanisms underlying heterogeneous excitatory/inhibitory balanced states

Recent work has explored spatiotemporal relationships between excitatory (E) and inhibitory (I) signaling within neural networks, and the effect of these relationships on network activity patterns. Data from these studies have indicated that excitation and inhibition are maintained at a similar level across long time periods and that excitatory and inhibitory currents may be tightly synchronized. Disruption of this balance-leading to an aberrant E/I ratio-is implicated in various brain pathologies. However, a thorough characterization of the relationship between E and I currents in experimental settings is largely impossible, due to their tight regulation at multiple cellular and network levels. Here, we use biophysical neural network models to investigate the emergence and properties of balanced states by heterogeneous mechanisms. Our results show that a network can homeostatically regulate the E/I ratio through interactions among multiple cellular and network factors, including average firing rates, synaptic weights and average neural depolarization levels in excitatory/inhibitory populations. Complex and competing interactions between firing rates and depolarization levels allow these factors to alternately dominate network dynamics in different synaptic weight regimes. This leads to the emergence of distinct mechanisms responsible for determining a balanced state and its dynamical correlate. Our analysis provides a comprehensive picture of how E/I ratio changes when manipulating specific network properties, and identifies the mechanisms regulating E/I balance. These results provide a framework to explain the diverse, and in some cases, contradictory experimental observations on the E/I state in different brain states and conditions.

Hippocampal µ-opioid receptors on GABAergic neurons mediate stress-induced impairment of memory retrieval

Stressful life events induce abnormalities in emotional and cognitive behaviour. The endogenous opioid system plays an essential role in stress adaptation and coping strategies. In particular, the µ-opioid receptor (μR), one of the major opioid receptors, strongly influences memory processing in that alterations in μR signalling are associated with various neuropsychiatric disorders. However, it remains unclear whether μR signalling contributes to memory impairments induced by acute stress. Here, we utilized pharmacological methods and cell-type-selective/non-cell-type-selective μR depletion approaches combined with behavioural tests, biochemical analyses, and in vitro electrophysiological recordings to investigate the role of hippocampal μR signalling in memory-retrieval impairment induced by acute elevated platform (EP) stress in mice. Biochemical and molecular analyses revealed that hippocampal μRs were significantly activated during acute stress. Blockage of hippocampal μRs, non-selective deletion of μRs or selective deletion of μRs on GABAergic neurons (μR_GABA) reversed EP-stress-induced impairment of memory retrieval, with no effect on the elevation of serum corticosterone after stress. Electrophysiological results demonstrated that stress depressed hippocampal GABAergic synaptic transmission to CA1 pyramidal neurons, thereby leading to excitation/inhibition (E/I) imbalance in a μR_GABA-dependent manner. Pharmaceutically enhancing hippocampal GABA_A receptor-mediated inhibitory currents in stressed mice restored their memory retrieval, whereas inhibiting those currents in the unstressed mice mimicked the stress-induced impairment of memory retrieval. Our findings reveal a novel pathway in which endogenous opioids recruited by acute stress predominantly activate μR_GABA to depress GABAergic inhibitory effects on CA1 pyramidal neurons, which subsequently alters the E/I balance in the hippocampus and results in impairment of memory retrieval.

Jensen's force and the statistical mechanics of cortical asynchronous states

Cortical networks are shaped by the combined action of excitatory and inhibitory interactions. Among other important functions, inhibition solves the problem of the all-or-none type of response that comes about in purely excitatory networks, allowing the network to operate in regimes of moderate or low activity, between quiescent and saturated regimes. Here, we elucidate a noise-induced effect that we call “Jensen’s force” –stemming from the combined effect of excitation/inhibition balance and network sparsity– which is responsible for generating a phase of self-sustained low activity in excitation-inhibition networks. The uncovered phase reproduces the main empirically-observed features of cortical networks in the so-called asynchronous state, characterized by low, un-correlated and highly-irregular activity. The parsimonious model analyzed here allows us to resolve a number of long-standing issues, such as proving that activity can be self-sustained even in the complete absence of external stimuli or driving. The simplicity of our approach allows for a deep understanding of asynchronous states and of the phase transitions to other standard phases it exhibits, opening the door to reconcile, asynchronous-state and critical-state hypotheses, putting them within a unified framework. We argue that Jensen’s forces are measurable experimentally and might be relevant in contexts beyond neuroscience.

Firing rate distributions in spiking networks with heterogeneous connectivity

Mean-field theory for networks of spiking neurons based on the so-called diffusion approximation has been used to calculate certain measures of neuronal activity which can be compared with experimental data. This includes the distribution of firing rates across the network. However, the theory in its current form applies only to networks in which there is relatively little heterogeneity in the number of incoming and outgoing connections per neuron. Here we extend this theory to include networks with arbitrary degree distributions. Furthermore, the theory takes into account correlations in the in-degree and out-degree of neurons, which would arise, e.g., in the case of networks with hublike neurons. Finally, we show that networks with broad and positively correlated degrees can generate a large-amplitude sustained response to transient stimuli which does not occur in more homogeneous networks.

Dynamical modeling of multi-scale variability in neuronal competition

Variability is observed at multiple-scales in the brain and ubiquitous in perception. However, the nature of perceptual variability is an open question. We focus on variability during perceptual rivalry, a form of neuronal competition. Rivalry provides a window into neural processing since activity in many brain areas is correlated to the alternating perception rather than a constant ambiguous stimulus. It exhibits robust properties at multiple scales including conscious awareness and neuron dynamics. The prevalent theory for spiking variability is called the balanced state; whereas, the source of perceptual variability is unknown. Here we show that a single biophysical circuit model, satisfying certain mutual inhibition architectures, can explain spiking and perceptual variability during rivalry. These models adhere to a broad set of strict experimental constraints at multiple scales. As we show, the models predict how spiking and perceptual variability changes with stimulus conditions.

Correlated states in balanced neuronal networks

Understanding the magnitude and structure of interneuronal correlations and their relationship to synaptic connectivity structure is an important and difficult problem in computational neuroscience. Early studies show that neuronal network models with excitatory-inhibitory balance naturally create very weak spike train correlations, defining the "asynchronous state." Later work showed that, under some connectivity structures, balanced networks can produce larger correlations between some neuron pairs, even when the average correlation is very small. All of these previous studies assume that the local network receives feedforward synaptic input from a population of uncorrelated spike trains. We show that when spike trains providing feedforward input are correlated, the downstream recurrent network produces much larger correlations. We provide an in-depth analysis of the resulting "correlated state" in balanced networks and show that, unlike the asynchronous state, it produces a tight excitatory-inhibitory balance consistent with in vivo cortical recordings.

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.

Cortical microcircuits are highly inhomogeneous dynamical systems whose information processing capacity is determined by the characteristics of its heterogeneous components and their complex interactions. The high degree of variability that characterizes macroscopic population dynamics, both during ongoing, spontaneous activity and active processing states reflects the underlying complexity and heterogeneity which has the potential to dramatically constrain the space of functions that any given circuit can compute, leading to richer and more expressive information processing systems. In this study, we identify different tentative sources of heterogeneity and assess their differential and cumulative contribution to the microcircuit’s dynamics and information processing capacity. We study these properties in a generic Layer 2/3 cortical microcircuit model, built and constrained by multiple sources of experimental data, and demonstrate that heterogeneity in neuronal properties and microconnectivity structure are important sources of functional specialization, greatly improving the circuit’s processing capacity, while capturing various important features of cortical physiology.

Rich-club connectivity, diverse population coupling, and dynamical activity patterns emerging from local cortical circuits

Experimental studies have begun revealing essential properties of the structural connectivity and the spatiotemporal activity dynamics of cortical circuits. To integrate these properties from anatomy and physiology, and to elucidate the links between them, we develop a novel cortical circuit model that captures a range of realistic features of synaptic connectivity. We show that the model accounts for the emergence of higher-order connectivity structures, including highly connected hub neurons that form an interconnected rich-club. The circuit model exhibits a rich repertoire of dynamical activity states, ranging from asynchronous to localized and global propagating wave states. We find that around the transition between asynchronous and localized propagating wave states, our model quantitatively reproduces a variety of major empirical findings regarding neural spatiotemporal dynamics, which otherwise remain disjointed in existing studies. These dynamics include diverse coupling (correlation) between spiking activity of individual neurons and the population, dynamical wave patterns with variable speeds and precise temporal structures of neural spikes. We further illustrate how these neural dynamics are related to the connectivity properties by analysing structural contributions to variable spiking dynamics and by showing that the rich-club structure is related to the diverse population coupling. These findings establish an integrated account of structural connectivity and activity dynamics of local cortical circuits, and provide new insights into understanding their working mechanisms.

To integrate essential anatomical and physiological properties of local cortical circuits and to elucidate mechanistic links between them, we develop a novel circuit model capturing key synaptic connectivity features. We show that the model explains the emergence of a range of connectivity patterns such as rich-club connectivity, and gives rise to a rich repertoire of cortical states. We identify both the anatomical and physiological mechanisms underlying the transition of these cortical states, and show that our model reconciles an otherwise disparate set of key physiological findings on neural activity dynamics. We further illustrate how these neural dynamics are related to the connectivity properties by analysing structural contributions to variable spiking dynamics and by showing that the rich-club structure is related to diverse neural population correlations as observed recently. Our model thus provides a framework for integrating and explaining a variety of neural connectivity properties and spatiotemporal activity dynamics observed in experimental studies, and provides novel experimentally testable predictions.

Homeostatic scaling of active zone scaffolds maintains global synaptic strength

Synaptic terminals grow and retract throughout life, yet synaptic strength is maintained within stable physiological ranges. To study this process, we investigated Drosophila endophilin (endo) mutants. Although active zone (AZ) number is doubled in endo mutants, a compensatory reduction in their size homeostatically adjusts global neurotransmitter output to maintain synaptic strength. We find an inverse adaptation in rab3 mutants. Additional analyses using confocal, STED, and electron microscopy reveal a stoichiometric tuning of AZ scaffolds and nanoarchitecture. Axonal transport of synaptic cargo via the lysosomal kinesin adapter Arl8 regulates AZ abundance to modulate global synaptic output and sustain the homeostatic potentiation of neurotransmission. Finally, we find that this AZ scaling can interface with two independent homeostats, depression and potentiation, to remodel AZ structure and function, demonstrating a robust balancing of separate homeostatic adaptations. Thus, AZs are pliable substrates with elastic and modular nanostructures that can be dynamically sculpted to stabilize and tune both local and global synaptic strength.

Suppressive Traveling Waves Shape Representations of Illusory Motion in Primary Visual Cortex of Awake Primate

How does the brain link visual stimuli across space and time? Visual illusions provide an experimental paradigm to study these processes. When two stationary dots are flashed in close spatial and temporal succession, human observers experience a percept of apparent motion. Large spatiotemporal separation challenges the visual system to keep track of object identity along the apparent motion path, the so-called "correspondence problem." Here, we use voltage-sensitive dye imaging in primary visual cortex (V1) of awake monkeys to show that intracortical connections within V1 can solve this issue by shaping cortical dynamics to represent the illusory motion. We find that the appearance of the second stimulus in V1 creates a systematic suppressive wave traveling toward the retinotopic representation of the first. Using a computational model, we show that the suppressive wave is the emergent property of a recurrent gain control fed by the intracortical network. This suppressive wave acts to explain away ambiguous correspondence problems and contributes to precisely encode the expected motion velocity at the surface of V1. Together, these results demonstrate that the nonlinear dynamics within retinotopic maps can shape cortical representations of illusory motion. Understanding these dynamics will shed light on how the brain links sensory stimuli across space and time, by preformatting population responses for a straightforward read-out by downstream areas.SIGNIFICANCE STATEMENT Traveling waves have recently been observed in different animal species, brain areas, and behavioral states. However, it is still unclear what are their functional roles. In the case of cortical visual processing, waves propagate across retinotopic maps and can hereby generate interactions between spatially and temporally separated instances of feedforward driven activity. Such interactions could participate in processing long-range apparent motion stimuli, an illusion for which no clear neuronal mechanisms have yet been proposed. Using this paradigm in awake monkeys, we show that suppressive traveling waves produce a spatiotemporal normalization of apparent motion stimuli. Our study suggests that cortical waves shape the representation of illusory moving stimulus within retinotopic maps for a straightforward read-out by downstream areas.

Balanced Active Core in Heterogeneous Neuronal Networks

It is hypothesized that cortical neuronal circuits operate in a global balanced state, i.e., the majority of neurons fire irregularly by receiving balanced inputs of excitation and inhibition. Meanwhile, it has been observed in experiments that sensory information is often sparsely encoded by only a small set of firing neurons, while neurons in the rest of the network are silent. The phenomenon of sparse coding challenges the hypothesis of a global balanced state in the brain. To reconcile this, here we address the issue of whether a balanced state can exist in a small number of firing neurons by taking account of the heterogeneity of network structure such as scale-free and small-world networks. We propose necessary conditions and show that, under these conditions, for sparsely but strongly connected heterogeneous networks with various types of single-neuron dynamics, despite the fact that the whole network receives external inputs, there is a small active subnetwork (active core) inherently embedded within it. The neurons in this active core have relatively high firing rates while the neurons in the rest of the network are quiescent. Surprisingly, although the whole network is heterogeneous and unbalanced, the active core possesses a balanced state and its connectivity structure is close to a homogeneous Erdös-Rényi network. The dynamics of the active core can be well-predicted using the Fokker-Planck equation. Our results suggest that the balanced state may be maintained by a small group of spiking neurons embedded in a large heterogeneous network in the brain. The existence of the small active core reconciles the balanced state and the sparse coding, and also provides a potential dynamical scenario underlying sparse coding in neuronal networks.