Rhythm and Synchrony in a Cortical Network Model

A comprehensive data-driven model of cat primary visual cortex

Knowledge integration based on the relationship between structure and function of the neural substrate is one of the main targets of neuroinformatics and data-driven computational modeling. However, the multiplicity of data sources, the diversity of benchmarks, the mixing of observables of different natures, and the necessity of a long-term, systematic approach make such a task challenging. Here we present a first snapshot of a long-term integrative modeling program designed to address this issue in the domain of the visual system: a comprehensive spiking model of cat primary visual cortex. The presented model satisfies an extensive range of anatomical, statistical and functional constraints under a wide range of visual input statistics. In the presence of physiological levels of tonic stochastic bombardment by spontaneous thalamic activity, the modeled cortical reverberations self-generate a sparse asynchronous ongoing activity that quantitatively matches a range of experimentally measured statistics. When integrating feed-forward drive elicited by a high diversity of visual contexts, the simulated network produces a realistic, quantitatively accurate interplay between visually evoked excitatory and inhibitory conductances; contrast-invariant orientation-tuning width; center surround interactions; and stimulus-dependent changes in the precision of the neural code. This integrative model offers insights into how the studied properties interact, contributing to a better understanding of visual cortical dynamics. It provides a basis for future development towards a comprehensive model of low-level perception.

Firing rate models for gamma oscillations in I-I and E-I networks

Firing rate models for describing the mean-field activities of neuronal ensembles can be used effectively to study network function and dynamics, including synchronization and rhythmicity of excitatory-inhibitory populations. However, traditional Wilson-Cowan-like models, even when extended to include an explicit dynamic synaptic activation variable, are found unable to capture some dynamics such as Interneuronal Network Gamma oscillations (ING). Use of an explicit delay is helpful in simulations at the expense of complicating mathematical analysis. We resolve this issue by introducing a dynamic variable, u, that acts as an effective delay in the negative feedback loop between firing rate (r) and synaptic gating of inhibition (s). In effect, u endows synaptic activation with second order dynamics. With linear stability analysis, numerical branch-tracking and simulations, we show that our r-u-s rate model captures some key qualitative features of spiking network models for ING. We also propose an alternative formulation, a v-u-s model, in which mean membrane potential v satisfies an averaged current-balance equation. Furthermore, we extend the framework to E-I networks. With our six-variable v-u-s model, we demonstrate in firing rate models the transition from Pyramidal-Interneuronal Network Gamma (PING) to ING by increasing the external drive to the inhibitory population without adjusting synaptic weights. Having PING and ING available in a single network, without invoking synaptic blockers, is plausible and natural for explaining the emergence and transition of two different types of gamma oscillations.

A large-scale neuronal network modelling study: Stimulus size modulates gamma oscillations in the primary visual cortex by long-range connections

Stimulus size modulation of neuronal firing activity is a fundamental property of the primary visual cortex. Numerous biological experiments have shown that stimulus size modulation is affected by multiple factors at different spatiotemporal scales, but the exact pathways and mechanisms remain incompletely understood. In this paper, we establish a large-scale neuronal network model of primary visual cortex with layer 2/3 to study how gamma oscillation properties are modulated by stimulus size and especially how long-range connections affect the modulation as realistic neuronal properties and spatial distributions of synaptic connections are considered. It is shown that long-range horizontal synaptic connections are sufficient to produce dimensional modulation of firing rates and gamma oscillations. In particular, with increasing grating stimulus size, the firing rate increases and then decreases, the peak frequency of gamma oscillations decreases and the spectral power increases. These are consistent with biological experimental observations. Furthermore, we explain in detail how the number and spatial distribution of long-range connections affect the size modulation of gamma oscillations by using the analysis of neuronal firing activity and synaptic current fluctuations. Our results provide a mechanism explanation for size modulation of gamma oscillations in the primary visual cortex and reveal the important and unique role played by long-range connections, which contributes to a deeper understanding of the cognitive function of gamma oscillations in visual cortex.

Learning spiking neuronal networks with artificial neural networks: neural oscillations

First-principles-based modelings have been extremely successful in providing crucial insights and predictions for complex biological functions and phenomena. However, they can be hard to build and expensive to simulate for complex living systems. On the other hand, modern data-driven methods thrive at modeling many types of high-dimensional and noisy data. Still, the training and interpretation of these data-driven models remain challenging. Here, we combine the two types of methods to model stochastic neuronal network oscillations. Specifically, we develop a class of artificial neural networks to provide faithful surrogates to the high-dimensional, nonlinear oscillatory dynamics produced by a spiking neuronal network model. Furthermore, when the training data set is enlarged within a range of parameter choices, the artificial neural networks become generalizable to these parameters, covering cases in distinctly different dynamical regimes. In all, our work opens a new avenue for modeling complex neuronal network dynamics with artificial neural networks.

Rodents' visual gamma as a biomarker of pathological neural conditions

Neural gamma oscillations (indicatively 30-100 Hz) are ubiquitous: they are associated with a broad range of functions in multiple cortical areas and across many animal species. Experimental and computational works established gamma rhythms as a global emergent property of neuronal networks generated by the balanced and coordinated interaction of excitation and inhibition. Coherently, gamma activity is strongly influenced by the alterations of synaptic dynamics which are often associated with pathological neural dysfunctions. We argue therefore that these oscillations are an optimal biomarker for probing the mechanism of cortical dysfunctions. Gamma oscillations are also highly sensitive to external stimuli in sensory cortices, especially the primary visual cortex (V1), where the stimulus dependence of gamma oscillations has been thoroughly investigated. Gamma manipulation by visual stimuli tuning is particularly easy in rodents, which have become a standard animal model for investigating the effects of network alterations on gamma oscillations. Overall, gamma in the rodents' visual cortex offers an accessible probe on dysfunctional information processing in pathological conditions. Beyond vision-related dysfunctions, alterations of gamma oscillations in rodents were indeed also reported in neural deficits such as migraine, epilepsy and neurodegenerative or neuropsychiatric conditions such as Alzheimer's, schizophrenia and autism spectrum disorders. Altogether, the connections between visual cortical gamma activity and physio-pathological conditions in rodent models underscore the potential of gamma oscillations as markers of neuronal (dys)functioning.

Neuronal composition of processing modules in human V1: laminar density for neuronal and non-neuronal populations and a comparison with macaque

The neuronal composition of homologous brain regions in different primates is important for understanding their processing capacities. Primary visual cortex (V1) has been widely studied in different members of the catarrhines. Neuronal density is considered to be central in defining the structure-function relationship. In human, there are large variations in the reported neuronal density from prior studies. We found the neuronal density in human V1 was 79,000 neurons/mm3, which is 35% of the neuronal density previously determined in macaque V1. Laminar density was proportionally similar between human and macaque. In V1, the ocular dominance column (ODC) contains the circuits for the emergence of orientation preference and spatial processing of a point image in many mammalian species. Analysis of the total neurons in an ODC and of the full number of neurons in macular vision (the central 15°) indicates that humans have 1.3× more neurons than macaques even though the density of neurons in macaque is 3× the density in human V1. We propose that the number of neurons in a functional processing unit rather than the number of neurons under a mm2 of cortex is more appropriate for cortical comparisons across species.

Temporal characteristics of gamma rhythm constrain properties of noise in an inhibition-stabilized network model

Gamma rhythm refers to oscillatory neural activity between 30 and 80 Hz, induced in visual cortex by stimuli such as iso-luminant hues or gratings. The power and peak frequency of gamma depend on the properties of the stimulus such as size and contrast. Gamma waveform is typically arch-shaped, with narrow troughs and broad peaks, and can be replicated in a self-oscillating Wilson-Cowan (WC) model operating in an appropriate regime. However, oscillations in this model are infinitely long, unlike physiological gamma that occurs in short bursts. Further, unlike the model, gamma is faster after stimulus onset and slows down over time. Here, we first characterized gamma burst duration in local field potential data recorded from two monkeys as they viewed full screen iso-luminant hues. We then added different types of noise in the inputs to the WC model and tested how that affected duration and temporal dynamics of gamma. While the model failed with the often-used Poisson noise, Ornstein-Uhlenbeck noise applied to both the excitatory and the inhibitory populations replicated the duration and slowing of gamma and replicated the shape and stimulus dependencies. Thus, the temporal dynamics of gamma oscillations put constraints on the type and properties of underlying neural noise.

Non-Trivial Dynamics in the FizHugh-Rinzel Model and Non-Homogeneous Oscillatory-Excitable Reaction-Diffusions Systems

This article focuses on the qualitative analysis of complex dynamics arising in a few mathematical models in neuroscience context. We first discuss the dynamics arising in the three-dimensional FitzHugh-Rinzel (FHR) model and then illustrate those arising in a class of non-homogeneous FitzHugh-Nagumo (Nh-FHN) reaction-diffusion systems. FHR and Nh-FHN models can be used to generate relevant complex dynamics and wave-propagation phenomena in neuroscience context. Such complex dynamics include canards, mixed-mode oscillations (MMOs), Hopf-bifurcations and their spatially extended counterpart. Our article highlights original methods to characterize these complex dynamics and how they emerge in ordinary differential equations and spatially extended models.

Synchronization and oscillation behaviors of excitatory and inhibitory populations with spike-timing-dependent plasticity

The effect of synaptic plasticity on the synchronization mechanism of the cerebral cortex has been a hot research topic over the past two decades. There are a great deal of literatures on excitatory pyramidal neurons, but the mechanism of interaction between the inhibitory interneurons is still under exploration. In this study, we consider a complex network consisting of excitatory (E) pyramidal neurons and inhibitory (I) interneurons interacting with chemical synapses through spike-timing-dependent plasticity (STDP). To study the effects of eSTDP and iSTDP on synchronization and oscillation behaviors emerged in an excitatory-inhibitory balanced network, we analyzed three different cases, a small-world network of purely excitatory neurons with eSTDP, a small-world network of purely inhibitory neurons with iSTDP and a small-world network with excitatory-inhibitory balanced neurons. By varying the number of inhibitory interneurons, and that of connected edges in a small-world network, and the coupling strength, these networks exhibit different synchronization and oscillation behaviors. We found that the eSTDP facilitates synchronization effectively, while iSTDP has no significant impact on it. In addition, eSTDP and iSTDP restrict the balance of the excitatory-inhibitory balanced neuronal network together and play a fundamental role in maintaining network stability and synchronization. They also can be used to guide the treatment and further research of neurodegenerative diseases.

Multi-band oscillations emerge from a simple spiking network

In the brain, coherent neuronal activities often appear simultaneously in multiple frequency bands, e.g., as combinations of alpha (8-12 Hz), beta (12.5-30 Hz), and gamma (30-120 Hz) oscillations, among others. These rhythms are believed to underlie information processing and cognitive functions and have been subjected to intense experimental and theoretical scrutiny. Computational modeling has provided a framework for the emergence of network-level oscillatory behavior from the interaction of spiking neurons. However, due to the strong nonlinear interactions between highly recurrent spiking populations, the interplay between cortical rhythms in multiple frequency bands has rarely been theoretically investigated. Many studies invoke multiple physiological timescales (e.g., various ion channels or multiple types of inhibitory neurons) or oscillatory inputs to produce rhythms in multi-bands. Here, we demonstrate the emergence of multi-band oscillations in a simple network consisting of one excitatory and one inhibitory neuronal population driven by constant input. First, we construct a data-driven, Poincaré section theory for robust numerical observations of single-frequency oscillations bifurcating into multiple bands. Then, we develop model reductions of the stochastic, nonlinear, high-dimensional neuronal network to capture the appearance of multi-band dynamics and the underlying bifurcations theoretically. Furthermore, when viewed within the reduced state space, our analysis reveals conserved geometrical features of the bifurcations on low-dimensional dynamical manifolds. These results suggest a simple geometric mechanism behind the emergence of multi-band oscillations without appealing to oscillatory inputs or multiple synaptic or neuronal timescales. Thus, our work points to unexplored regimes of stochastic competition between excitation and inhibition behind the generation of dynamic, patterned neuronal activities.

State-dependent modulation of thalamocortical oscillations by gamma light flicker with different frequencies, intensities, and duty cycles

Rhythmic light flickers have emerged as useful tools to modulate cognition and rescue pathological oscillations related to neurological disorders by entrainment. However, a mechanistic understanding of the entrainment for different brain oscillatory states and light flicker parameters is lacking. To address this issue, we proposed a biophysical neural network model for thalamocortical oscillations (TCOs) and explored the stimulation effects depending on the thalamocortical oscillatory states and stimulation parameters (frequency, intensity, and duty cycle) using the proposed model and electrophysiology experiments. The proposed model generated alpha, beta, and gamma oscillatory states (with main oscillation frequences at 9, 25, and 35 Hz, respectively), which were successfully transmitted from the thalamus to the cortex. By applying light flicker stimulation, we found that the entrainment was state-dependent and it was more prone to induce entrainment if the flicker perturbation frequency was closer to the endogenous oscillatory frequency. In addition, endogenous oscillation would be accelerated, whereas low-frequency oscillatory power would be suppressed by gamma (30–50 Hz) flickers. Notably, the effects of intensity and duty cycle on entrainment were complex; a high intensity of light flicker did not mean high entrainment possibility, and duty cycles below 50% could induce entrainment easier than those above 50%. Further, we observed entrainment discontinuity during gamma flicker stimulations with different frequencies, attributable to the non-linear characteristics of the network oscillations. These results provide support for the experimental design and clinical applications of the modulation of TCOs by gamma (30–50 Hz) light flicker.

The Use of Reduced Models to Generate Irregular, Broad-Band Signals That Resemble Brain Rhythms

The brain produces rhythms in a variety of frequency bands. Some are likely by-products of neuronal processes; others are thought to be top-down. Produced entirely naturally, these rhythms have clearly recognizable beats, but they are very far from periodic in the sense of mathematics. The signals are broad-band, episodic, wandering in amplitude and frequency; the rhythm comes and goes, degrading and regenerating. Gamma rhythms, in particular, have been studied by many authors in computational neuroscience, using reduced models as well as networks of hundreds to thousands of integrate-and-fire neurons. All of these models captured successfully the oscillatory nature of gamma rhythms, but the irregular character of gamma in reduced models has not been investigated thoroughly. In this article, we tackle the mathematical question of whether signals with the properties of brain rhythms can be generated from low dimensional dynamical systems. We found that while adding white noise to single periodic cycles can to some degree simulate gamma dynamics, such models tend to be limited in their ability to capture the range of behaviors observed. Using an ODE with two variables inspired by the FitzHugh-Nagumo and Leslie-Gower models, with stochastically varying coefficients designed to control independently amplitude, frequency, and degree of degeneracy, we were able to replicate the qualitative characteristics of natural brain rhythms. To demonstrate model versatility, we simulate the power spectral densities of gamma rhythms in various brain states recorded in experiments.

A Computational Model of Direction Selectivity in Macaque V1 Cortex Based on Dynamic Differences between On and Off Pathways

This paper is about neural mechanisms of direction selectivity (DS) in macaque primary visual cortex, V1. We present data (on male macaque) showing strong DS in a majority of simple cells in V1 layer 4Cα, the cortical layer that receives direct afferent input from the magnocellular division of the lateral geniculate nucleus (LGN). Magnocellular LGN cells are not direction-selective. To understand the mechanisms of DS, we built a large-scale, recurrent model of spiking neurons called DSV1. Like its predecessors, DSV1 reproduces many visual response properties of V1 cells including orientation selectivity. Two important new features of DSV1 are (1) DS is initiated by small, consistent dynamic differences in the visual responses of OFF and ON Magnocellular LGN cells, and (2) DS in the responses of most model simple cells is increased over those of their feedforward inputs; this increase is achieved through dynamic interaction of feedforward and intracortical synaptic currents without the use of intracortical direction-specific connections. The DSV1 model emulates experimental data in the following ways: (1) most 4Cα Simple cells were highly direction-selective but 4Cα Complex cells were not; (2) the preferred directions of the model's direction-selective Simple cells were invariant with spatial and temporal frequency (TF); (3) the distribution of the preferred/opposite ratio across the model's population of cells was very close to that found in experiments. The strong quantitative agreement between DS in data and in model simulations suggests that the neural mechanisms of DS in DSV1 may be similar to those in the real visual cortex.SIGNIFICANCE STATEMENT Motion perception is a vital part of our visual experience of the world. In monkeys, whose vision resembles that of humans, the neural computation of the direction of a moving target starts in the primary visual cortex, V1, in layer 4Cα that receives input from the eye through the lateral geniculate nucleus (LGN). How direction selectivity (DS) is generated in layer 4Cα is an outstanding unsolved problem in theoretical neuroscience. In this paper, we offer a solution based on plausible biological mechanisms. We present a new large-scale circuit model in which DS originates from slightly different LGN ON/OFF response time-courses and is enhanced in cortex without the need for direction-specific intracortical connections. The model's DS is in quantitative agreement with experiments.

Shape analysis of gamma rhythm supports a superlinear inhibitory regime in an inhibition-stabilized network

Visual inspection of stimulus-induced gamma oscillations (30–70 Hz) often reveals a non-sinusoidal shape. Such distortions are a hallmark of non-linear systems and are also observed in mean-field models of gamma oscillations. A thorough characterization of the shape of the gamma cycle can therefore provide additional constraints on the operating regime of such models. However, the gamma waveform has not been quantitatively characterized, partially because the first harmonic of gamma, which arises because of the non-sinusoidal nature of the signal, is typically weak and gets masked due to a broadband increase in power related to spiking. To address this, we recorded local field potential (LFP) from the primary visual cortex (V1) of two awake female macaques while presenting full-field gratings or iso-luminant chromatic hues that produced huge gamma oscillations with prominent peaks at harmonic frequencies in the power spectra. We found that gamma and its first harmonic always maintained a specific phase relationship, resulting in a distinctive shape with a sharp trough and a shallow peak. Interestingly, a Wilson-Cowan (WC) model operating in an inhibition stabilized mode could replicate this shape, but only when the inhibitory population operated in the super-linear regime, as predicted recently. However, another recently developed model of gamma that operates in a linear regime driven by stochastic noise failed to produce salient harmonics or the observed shape. Our results impose additional constraints on models that generate gamma oscillations and their operating regimes.

Gamma rhythm is not sinusoidal. Understanding these distortions could provide clues about the cortical network that generates the rhythm. Here, we use harmonic phase analysis to describe these waveforms quantitatively and show that gamma rhythm recorded from the primary visual cortex of macaques has a signature arch shaped waveform, with a sharp trough and a shallow peak, when visual stimuli such as full-screen plain hues and achromatic gratings are presented. This arch shaped waveform is observed over a wide range of stimuli, despite the variation in power and frequency of the rhythm. We then compare two population rate models that have been used to accurately describe the stimulus dependencies of gamma rhythm and show that this arch shaped waveform is obtained only in one of those models. Further, the waveform shape is dependent on the operating domain of the system. Therefore, shape analysis provides additional constraints on cortical models and their operating regimes.

Gamma rhythms in the visual cortex: functions and mechanisms

Gamma-band activity, peaking around 30–100 Hz in the local field potential's power spectrum, has been found and intensively studied in many brain regions. Although gamma is thought to play a critical role in processing neural information in the brain, its cognitive functions and neural mechanisms remain unclear or debatable. Experimental studies showed that gamma rhythms are stochastic in time and vary with visual stimuli. Recent studies further showed that multiple rhythms coexist in V1 with distinct origins in different species. While all these experimental facts are a challenge for understanding the functions of gamma in the visual cortex, there are many signs of progress in computational studies. This review summarizes and discusses studies on gamma in the visual cortex from multiple perspectives and concludes that gamma rhythms are still a mystery. Combining experimental and computational studies seems the best way forward in the future.

Multiple gamma rhythms carry distinct spatial frequency information in primary visual cortex

Gamma rhythms in many brain regions, including the primary visual cortex (V1), are thought to play a role in information processing. Here, we report a surprising finding of 3 narrowband gamma rhythms in V1 that processed distinct spatial frequency (SF) signals and had different neural origins. The low gamma (LG; 25 to 40 Hz) rhythm was generated at the V1 superficial layer and preferred a higher SF compared with spike activity, whereas both the medium gamma (MG; 40 to 65 Hz), generated at the cortical level, and the high gamma HG; (65 to 85 Hz), originated precortically, preferred lower SF information. Furthermore, compared with the rates of spike activity, the powers of the 3 gammas had better performance in discriminating the edge and surface of simple objects. These findings suggest that gamma rhythms reflect the neural dynamics of neural circuitries that process different SF information in the visual system, which may be crucial for multiplexing SF information and synchronizing different features of an object.

Extensive Structural Remodeling of the Axonal Arbors of Parvalbumin Basket Cells during Development in Mouse Neocortex

Parvalbumin-containing (PV+) basket cells are specialized cortical interneurons that regulate the activity of local neuronal circuits with high temporal precision and reliability. To understand how the PV+ interneuron connectivity underlying these functional properties is established during development, we used array tomography to map pairs of synaptically connected PV+ interneurons and postsynaptic neurons from the neocortex of mice of both sexes. We focused on the axon-myelin unit of the PV+ interneuron and quantified the number of synapses onto the postsynaptic neuron, length of connecting axonal paths, and their myelination at different time points between 2 weeks and 7 months of age. We find that myelination of the proximal axon occurs very rapidly during the third and, to a lesser extent, fourth postnatal weeks. The number of synaptic contacts made by the PV+ interneuron on its postsynaptic partner meanwhile is significantly reduced to about one-third by the end of the first postnatal month. The number of autapses, the synapses that PV+ interneurons form on themselves, however, remains constant throughout the examined period. Axon reorganizations continue beyond postnatal month 2, with the postsynaptic targets of PV+ interneurons gradually shifting to more proximal locations, and the length of axonal paths and their myelin becoming conspicuously uniform per connection. These continued microcircuit refinements likely provide the structural substrate for the robust inhibitory effects and fine temporal precision of adult PV+ basket cells.SIGNIFICANCE STATEMENT The axon of adult parvalbumin-containing (PV+) interneurons is highly specialized for fast and reliable neurotransmission. It is myelinated and forms synapses mostly onto the cell bodies and proximal dendrites of postsynaptic neurons for maximal impact. In this study, we follow the development of the PV+ interneuron axon, its myelination and synapse formation, revealing a rapid sequence of axonal reorganization, myelination of the PV+ interneuron proximal axon, and pruning of almost two-thirds of the synapses in an individual connection. This is followed by a prolonged period of axon refinement and additional myelination leading to a remarkable precision of connections in the adult mouse cortex, consistent with the temporal precision and fidelity of PV+ interneuron action.

Model Reduction Captures Stochastic Gamma Oscillations on Low-Dimensional Manifolds

Gamma frequency oscillations (25–140 Hz), observed in the neural activities within many brain regions, have long been regarded as a physiological basis underlying many brain functions, such as memory and attention. Among numerous theoretical and computational modeling studies, gamma oscillations have been found in biologically realistic spiking network models of the primary visual cortex. However, due to its high dimensionality and strong non-linearity, it is generally difficult to perform detailed theoretical analysis of the emergent gamma dynamics. Here we propose a suite of Markovian model reduction methods with varying levels of complexity and apply it to spiking network models exhibiting heterogeneous dynamical regimes, ranging from nearly homogeneous firing to strong synchrony in the gamma band. The reduced models not only successfully reproduce gamma oscillations in the full model, but also exhibit the same dynamical features as we vary parameters. Most remarkably, the invariant measure of the coarse-grained Markov process reveals a two-dimensional surface in state space upon which the gamma dynamics mainly resides. Our results suggest that the statistical features of gamma oscillations strongly depend on the subthreshold neuronal distributions. Because of the generality of the Markovian assumptions, our dimensional reduction methods offer a powerful toolbox for theoretical examinations of other complex cortical spatio-temporal behaviors observed in both neurophysiological experiments and numerical simulations.

A bio-syncretic phototransistor based on optogenetically engineered living cells

Human eyes rely on photosensitive receptors to convert light intensity into action potentials for visual perception, and thus bio-inspired photodetectors with bioengineered photoresponsive elements for visual prostheses have received considerable attention by virtue of superior biological functionality and better biocompatibility. However, the current bioengieered photodetectors based on biological elements face a lot of challenges such as slow response time and lack of effective detection of weak bioelectrical signals, resulting in difficulty to perform imaging. Here, we report a human eye-inspired phototransistor by integrating optogenetically engineered living cells and a graphene-based transistor. The living cells, engineered with photosensitive ion channels, channelrhodopsin-2 (ChR2), and thus endowed with the capability of transducing light intensity into bioelectrical signals, are coupled with the graphene layer of the transistor and can regulate the transistor's output. The results show that the photosensitive ion channels enable the phototransistor to output stronger photoelectrical currents with relatively fast response (~25 ms) and wider dynamic range, and demonstrate the transistor owns optical and biological gating with a significant large on/off ratio of 197.5 and high responsivity of 1.37 mA W^-1. An artificial imaging system, which mimics the pathway of human visual information transmission from the retina through the lateral geniculate nucleus to the visual cortex, is constructed with the transistor and demonstrate the feasibility of imaging using the bioengineered cells. This work shows a potential that optogenetically engineered cells can be used to develop novel visual prostheses and paves a new avenue for engineering bio-syncretic sensing devices.

Malleability of gamma rhythms enhances population-level correlations

An important problem in systems neuroscience is to understand how information is communicated among brain regions, and it has been proposed that communication is mediated by neuronal oscillations, such as rhythms in the gamma band. We sought to investigate this idea by using a network model with two components, a source (sending) and a target (receiving) component, both built to resemble local populations in the cerebral cortex. To measure the effectiveness of communication, we used population-level correlations in spike times between the source and target. We found that after correcting for a response time that is independent of initial conditions, spike-time correlations between the source and target are significant, due in large measure to the alignment of firing events in their gamma rhythms. But, we also found that regular oscillations cannot produce the results observed in our model simulations of cortical neurons. Surprisingly, it is the irregularity of gamma rhythms, the absence of internal clocks, together with the malleability of these rhythms and their tendency to align with external pulses - features that are known to be present in gamma rhythms in the real cortex - that produced the results observed. These findings and the mechanistic explanations we offered are our primary results. Our secondary result is a mathematical relationship between correlations and the sizes of the samples used for their calculation. As improving technology enables recording simultaneously from increasing numbers of neurons, this relationship could be useful for interpreting results from experimental recordings.

Conduction Velocity Along the Local Axons of Parvalbumin Interneurons Correlates With the Degree of Axonal Myelination

Parvalbumin-containing (PV+) basket cells in mammalian neocortex are fast-spiking interneurons that regulate the activity of local neuronal circuits in multiple ways. Even though PV+ basket cells are locally projecting interneurons, their axons are myelinated. Can this myelination contribute in any significant way to the speed of action potential propagation along such short axons? We used dual whole cell recordings of synaptically connected PV+ interneurons and their postsynaptic target in acutely prepared neocortical slices from adult mice to measure the amplitude and latency of single presynaptic action potential-evoked inhibitory postsynaptic currents. These same neurons were then imaged with immunofluorescent array tomography, the synapses between them identified and a precise map of the connections was generated, with the exact axonal length and extent of myelin coverage. Our results support that myelination of PV+ basket cells significantly increases conduction velocity, and does so to a degree that can be physiologically relevant.

The Generation and Modulation of Distinct Gamma Oscillations with Local, Horizontal, and Feedback Connections in the Primary Visual Cortex: A Model Study on Large-Scale Networks

Gamma oscillation (GAMMA) in the local field potential (LFP) is a synchronized activity commonly found in many brain regions, and it has been thought as a functional signature of network connectivity in the brain, which plays important roles in information processing. Studies have shown that the response property of GAMMA is related to neural interaction through local recurrent connections (RC), feed-forward (FF), and feedback (FB) connections. However, the relationship between GAMMA and long-range horizontal connections (HC) in the brain remains unclear. Here, we aimed to understand this question in a large-scale network model for the primary visual cortex (V1). We created a computational model composed of multiple excitatory and inhibitory units with biologically plausible connectivity patterns for RC, FF, FB, and HC in V1; then, we quantitated GAMMA in network models at different strength levels of HC and other connection types. Surprisingly, we found that HC and FB, the two types of large-scale connections, play very different roles in generating and modulating GAMMA. While both FB and HC modulate a fast gamma oscillation (around 50-60 Hz) generated by FF and RC, HC generates a new GAMMA oscillating around 30 Hz, whose power and peak frequency can also be modulated by FB. Furthermore, response properties of the two GAMMAs in a network with both HC and FB are different in a way that is highly consistent with a recent experimental finding for distinct GAMMAs in macaque V1. The results suggest that distinct GAMMAs are signatures for neural connections in different spatial scales and they might be related to different functions for information integration. Our study, for the first time, pinpoints the underlying circuits for distinct GAMMAs in a mechanistic model for macaque V1, which might provide a new framework to study multiple gamma oscillations in other cortical regions.

Contrast response in a comprehensive network model of macaque V1

The response to contrast is one of the most important functions of the macaque primary visual cortex, V1, but up to now there has not been an adequate theory for it. To fill this gap in our understanding of cortical function, we built and analyzed a new large-scale, biologically constrained model of the input layer, 4Cα, of macaque V1. We called the new model CSY2. We challenged CSY2 with a three-parameter family of visual stimuli that varied in contrast, orientation, and spatial frequency. CSY2 accurately simulated experimental data and made many new predictions. It accounted for 1) the shapes of firing-rate-versus-contrast functions, 2) orientation and spatial frequency tuning versus contrast, and 3) the approximate contrast-invariance of cortical activity maps. Post-analysis revealed that the mechanisms that were needed to produce the successful simulations of contrast response included strong recurrent excitation and inhibition that find dynamic equilibria across the cortical surface, dynamic feedback between L6 and L4, and synaptic dynamics like inhibitory synaptic depression.

DNN-assisted statistical analysis of a model of local cortical circuits

In neuroscience, computational modeling is an effective way to gain insight into cortical mechanisms, yet the construction and analysis of large-scale network models—not to mention the extraction of underlying principles—are themselves challenging tasks, due to the absence of suitable analytical tools and the prohibitive costs of systematic numerical exploration of high-dimensional parameter spaces. In this paper, we propose a data-driven approach assisted by deep neural networks (DNN). The idea is to first discover certain input-output relations, and then to leverage this information and the superior computation speeds of the well-trained DNN to guide parameter searches and to deduce theoretical understanding. To illustrate this novel approach, we used as a test case a medium-size network of integrate-and-fire neurons intended to model local cortical circuits. With the help of an accurate yet extremely efficient DNN surrogate, we revealed the statistics of model responses, providing a detailed picture of model behavior. The information obtained is both general and of a fundamental nature, with direct application to neuroscience. Our results suggest that the methodology proposed can be scaled up to larger and more complex biological networks when used in conjunction with other techniques of biological modeling.

Functional Clusters of Neurons in Layer 6 of Macaque V1

Layer 6 appears to perform a very important role in the function of macaque primary visual cortex, V1, but not enough is understood about the functional characteristics of neurons in the layer 6 population. It is unclear to what extent the population is homogeneous with respect to their visual properties or if one can identify distinct subpopulations. Here we performed a cluster analysis based on measurements of the responses of single neurons in layer 6 of primary visual cortex in male macaque monkeys (Macaca fascicularis) to achromatic grating stimuli that varied in orientation, direction of motion, spatial and temporal frequency, and contrast. The visual stimuli were presented in a stimulus window that was also varied in size. Using the responses to parametric variation in these stimulus variables, we extracted a number of tuning response measures and used them in the cluster analysis. Six main clusters emerged along with some smaller clusters. Additionally, we asked whether parameter distributions from each of the clusters were statistically different. There were clear separations of parameters between some of the clusters, particularly for f1/f0 ratio, direction selectivity, and temporal frequency bandwidth, but other dimensions also showed differences between clusters. Our data suggest that in layer 6 there are multiple parallel circuits that provide information about different aspects of the visual stimulus.SIGNIFICANCE STATEMENT The cortex is multilayered and is involved in many high-level computations. In the current study, we have asked whether there are subpopulations of neurons, clusters, in layer 6 of cortex with different functional tuning properties that provide information about different aspects of the visual image. We identified six major functional clusters within layer 6. These findings show that there is much more complexity to the circuits in cortex than previously demonstrated and open up a new avenue for experimental investigation within layers of other cortical areas and for the elaboration of models of circuit function that incorporate many parallel pathways with different functional roles.

Excitatory/Inhibitory Responses Shape Coherent Neuronal Dynamics Driven by Optogenetic Stimulation in the Primate Brain

Coherent neuronal dynamics play an important role in complex cognitive functions. Optogenetic stimulation promises to provide new ways to test the functional significance of coherent neural activity. However, the mechanisms by which optogenetic stimulation drives coherent dynamics remain unclear, especially in the nonhuman primate brain. Here, we perform computational modeling and experiments to study the mechanisms of optogenetic-stimulation-driven coherent neuronal dynamics in three male nonhuman primates. Neural responses arise from stimulation-evoked, temporally dynamic excitatory (E) and inhibitory (I) activity. Spiking activity is more likely to occur during E/I imbalances. Thus the relative difference in the driven E and I responses precisely controls spike timing by forming a brief time interval of increased spiking likelihood. Experimental results agree with parameter-dependent predictions from the computational models. These results demonstrate that optogenetic stimulation driven coherent neuronal dynamics are governed by the temporal properties of E/I activity. Transient imbalances in excitatory and inhibitory activity may provide a general mechanism for generating coherent neuronal dynamics without the need for an oscillatory generator.SIGNIFICANCE STATEMENT We examine how coherent neuronal dynamics arise from optogenetic stimulation in the primate brain. Using computational models and experiments, we demonstrate that coherent spiking and local field potential activity is generated by stimulation-evoked responses of excitatory and inhibitory activity in networks, extending the growing literature on neuronal dynamics. These responses create brief time intervals of increased spiking tendency and are consistent with previous observations in the literature that balanced excitation and inhibition controls spike timing, suggesting that optogenetic-stimulation-driven coherence may arise from intrinsic E/I balance. Most importantly, our results are obtained in nonhuman primates and thus will play a leading role in driving the use of causal manipulations with optogenetic tools to study higher cognitive functions in the primate brain.

A case study in the functional consequences of scaling the sizes of realistic cortical models

Neuroscience models come in a wide range of scales and specificity, from mean-field rate models to large-scale networks of spiking neurons. There are potential trade-offs between simplicity and realism, versatility and computational speed. This paper is about large-scale cortical network models, and the question we address is one of scalability: would scaling down cell density impact a network’s ability to reproduce cortical dynamics and function? We investigated this problem using a previously constructed realistic model of the monkey visual cortex that is true to size. Reducing cell density gradually up to 50-fold, we studied changes in model behavior. Size reduction without parameter adjustment was catastrophic. Surprisingly, relatively minor compensation in synaptic weights guided by a theoretical algorithm restored mean firing rates and basic function such as orientation selectivity to models 10-20 times smaller than the real cortex. Not all was normal in the reduced model cortices: intracellular dynamics acquired a character different from that of real neurons, and while the ability to relay feedforward inputs remained intact, reduced models showed signs of deficiency in functions that required dynamical interaction among cortical neurons. These findings are not confined to models of the visual cortex, and modelers should be aware of potential issues that accompany size reduction. Broader implications of this study include the importance of homeostatic maintenance of firing rates, and the functional consequences of feedforward versus recurrent dynamics, ideas that may shed light on other species and on systems suffering cell loss.

With the vast numbers of neurons in the cerebral cortex, models in neuroscience are, for practical reasons, often orders of magnitude smaller than the actual network. We examine in this article the scalability of cortical networks. We find that function and dynamics in a network depend on network size. For illustration, we use a previously constructed realistic model of monkey visual cortex. Reducing the number of cells in the model, we find that small changes in synaptic weights can help maintain firing rates. However, model characteristics change fundamentally in the reduced models. Neurons have abnormal intracellular dynamics. The model becomes dominated by feedforward inputs and is less capable of functions requiring neuronal interaction. Modelers need to be aware of the potential issues with reduced cortical network models.

Gamma Oscillations in the Basolateral Amygdala: Biophysical Mechanisms and Computational Consequences

The basolateral nucleus of the amygdala (BL) is thought to support numerous emotional behaviors through specific microcircuits. These are often thought to be comprised of feedforward networks of principal cells (PNs) and interneurons. Neither well-understood nor often considered are recurrent and feedback connections, which likely engender oscillatory dynamics within BL. Indeed, oscillations in the gamma frequency range (40 − 100 Hz) are known to occur in the BL, and yet their origin and effect on local circuits remains unknown. To address this, we constructed a biophysically and anatomically detailed model of the rat BL and its local field potential (LFP) based on the physiological and anatomical literature, along with in vivo and in vitro data we collected on the activities of neurons within the rat BL. Remarkably, the model produced intermittent gamma oscillations (∼50 − 70 Hz) whose properties matched those recorded in vivo, including their entrainment of spiking. BL gamma-band oscillations were generated by the intrinsic circuitry, depending upon reciprocal interactions between PNs and fast-spiking interneurons (FSIs), while connections within these cell types affected the rhythm’s frequency. The model allowed us to conduct experimentally impossible tests to characterize the synaptic and spatial properties of gamma. The entrainment of individual neurons to gamma depended on the number of afferent connections they received, and gamma bursts were spatially restricted in the BL. Importantly, the gamma rhythm synchronized PNs and mediated competition between ensembles. Together, these results indicate that the recurrent connectivity of BL expands its computational and communication repertoire.

How well do reduced models capture the dynamics in models of interacting neurons?

This paper introduces a class of stochastic models of interacting neurons with emergent dynamics similar to those seen in local cortical populations. Rigorous results on existence and uniqueness of nonequilibrium steady states are proved. These network models are then compared to very simple reduced models driven by the same mean excitatory and inhibitory currents. Discrepancies in firing rates between network and reduced models are investigated and explained by correlations in spiking, or partial synchronization, working in concert with "nonlinearities" in the time evolution of membrane potentials. The use of simple random walks and their first passage times to simulate fluctuations in neuronal membrane potentials and interspike times is also considered.