Identifying and exploiting the anatomical origin of population rate oscillations in multi-layered spiking networks
© Bos et al. 2015
Published: 18 December 2015
Fast oscillations of the population firing rate in the high gamma range (50-200 Hz), where individual neurons fire slowly and irregularly, are observed in the living brain and in network models of leaky integrate-and-fire (LIF) neurons, that have also been studied analytically . However, a systematic approach identifying sub-circuits responsible for specific oscillations in a structured network of neural populations is currently unavailable.
We consider a large-scale, neural network consisting of 4 layers each composed of an excitatory and inhibitory population of LIF-neurons with connectivity determined by electrophysiological and anatomical studies . In simulations we observe a peak in the power spectrum around 83 Hz in all populations and low frequency oscillations with smaller power in a subset of the populations. Mapping the dynamics of the fluctuations to an effective linear rate model, using the recently derived transfer function for LIF-neurons with synaptic filtering , we derive the spectra of the population firing rates analytically.
Decomposing the noise-driven fluctuations into eigenmodes of the effective connectivity, we identify the modes responsible for peaks in the spectra. Applying perturbation theory, we quantify the influence of individual anatomical connections on the spectrum at given frequencies and identify a sub-circuitry, localized in the supra-granular and granular layer, generating the oscillation. These findings are in agreement with layer-specific local field potential measurements in the Macaque primary visual cortex, where gamma-frequency oscillations were mostly pronounced in layer 2,3 and 4B . We exploit this method i) to identify the connectivity loops responsible for the observed peaks and ii) to alter the circuitry in a targeted manner to control the position and amplitude of the peaks and the generation of slow frequency fluctuations. This requires removal and addition of only small numbers of synapses. The analytical framework moreover explains the suppression of higher frequencies by distributed delays and the amplification of population specific oscillatory input. Mapping the stimulus vector onto the eigenmodes of the system shows how the components of the input vector are processed in the network. Thus one can derive the sensitivity of the population rate dynamics to the direction and frequency of stimuli.
Our method finds application in the identification of the connectivity loops that determine emergent and externally driven global measures of activity observable in experiments as well as in engineering circuits that exhibit desired correlations on the population level.
We acknowledge funding by the Helmholtz Association: portfolio theme SMHB and Young Investigator's Group VH-NG-1028, and EU Grants 269921 (BrainScaleS) and 604102 (Human Brain Project). All network simulations were carried out with NEST (http://www.nest-initiative.org).
- Brunel N, Hakim V: Fast Global Oscillations in Networks of Integrate-and-Fire Neurons with Low Firing Rate. Neural Comput. 1998, 1621-1671.Google Scholar
- Potjans TC, Diesmann M: The cell-type specific cortical microcircuit: relating structure and activity in a full-scale spiking network model. Cereb Cortex. 2014, 24: 785-806.PubMedPubMed CentralView ArticleGoogle Scholar
- Schuecker J, Diesmann M, Helias M: Spectral properties of excitable systems subject to colored noise. arXiv:1411.0432 [cond-mat.stat-mech]Google Scholar
- Xing D, Yeh C-I, Burns S, Shapley MS: Laminar analysis of visually evoked activity in the primary visual cortex. PNAS. 2012, 109 (34): 13871-13876.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.