Dimensionally-reduced visual cortical network model predicts network response and connects system- and cellular-level descriptions

Dimensional reduction of a V1 ring model with simple and complex cells

In this paper, we extend a framework for constructing low-dimensional dynamical systems models of mammalian primary visual cortex to a cortical network model that incorporates the full nonlinear effects of complex cells. The procedure consists of capturing the essential dynamics in a low-dimensional subspace using empirical methods, then recasting the equations in the reduced vector space. Previously, we considered visual cortical network models consisting of only simple cells with nearly linear responses to external stimuli. Here we show that fully nonlinear effects can be incorporated by examining the dimensional reduction of an idealized ring model of V1 with both simple and complex cells. We found it expedient to divide the subspace into four separate neuronal populations: excitatory simple, excitatory complex, inhibitory simple and inhibitory complex. In order to reproduce the fluctuation-driven dynamics in this reduced space, we incorporated (1) white noises with different intensities into individual neuronal populations, and (2) firing rate estimates to capture the probability of firing due to subthreshold fluctuations. With a more accurate, fitted connectivity, our modified dimensional reduced models can reproduce the firing rates, circular variances and modulation ratios observed in the original ring model.

A principled dimension-reduction method for the population density approach to modeling networks of neurons with synaptic dynamics

The population density approach to neural network modeling has been utilized in a variety of contexts. The idea is to group many similar noisy neurons into populations and track the probability density function for each population that encompasses the proportion of neurons with a particular state rather than simulating individual neurons (i.e., Monte Carlo). It is commonly used for both analytic insight and as a time-saving computational tool. The main shortcoming of this method is that when realistic attributes are incorporated in the underlying neuron model, the dimension of the probability density function increases, leading to intractable equations or, at best, computationally intensive simulations. Thus, developing principled dimension-reduction methods is essential for the robustness of these powerful methods. As a more pragmatic tool, it would be of great value for the larger theoretical neuroscience community. For exposition of this method, we consider a single uncoupled population of leaky integrate-and-fire neurons receiving external excitatory synaptic input only. We present a dimension-reduction method that reduces a two-dimensional partial differential-integral equation to a computationally efficient one-dimensional system and gives qualitatively accurate results in both the steady-state and nonequilibrium regimes. The method, termed modified mean-field method, is based entirely on the governing equations and not on any auxiliary variables or parameters, and it does not require fine-tuning. The principles of the modified mean-field method have potential applicability to more realistic (i.e., higher-dimensional) neural networks.

Mapping Functional Connectivity between Neuronal Ensembles with Larval Zebrafish Transgenic for a Ratiometric Calcium Indicator

The ability to map functional connectivity is necessary for the study of the flow of activity in neuronal circuits. Optical imaging of calcium indicators, including FRET-based genetically encoded indicators and extrinsic dyes, is an important adjunct to electrophysiology and is widely used to visualize neuronal activity. However, techniques for mapping functional connectivities with calcium imaging data have been lacking. We present a procedure to compute reduced functional couplings between neuronal ensembles undergoing seizure activity from ratiometric calcium imaging data in three steps: (1) calculation of calcium concentrations and neuronal firing rates from ratiometric data; (2) identification of putative neuronal populations from spatio-temporal time-series of neural bursting activity; and then, (3) derivation of reduced connectivity matrices that represent neuronal population interactions. We apply our method to the larval zebrafish central nervous system undergoing chemoconvulsant-induced seizures. These seizures generate propagating, central nervous system-wide neural activity from which population connectivities may be calculated. This automatic functional connectivity mapping procedure provides a practical and user-independent means for summarizing the flow of activity between neuronal ensembles.