Spontaneous Local Gamma Oscillation Selectively Enhances Neural Network Responsiveness

PLOS Computational Biology, Apr 2009

Synchronized oscillation is very commonly observed in many neuronal systems and might play an important role in the response properties of the system. We have studied how the spontaneous oscillatory activity affects the responsiveness of a neuronal network, using a neural network model of the visual cortex built from Hodgkin-Huxley type excitatory (E-) and inhibitory (I-) neurons. When the isotropic local E-I and I-E synaptic connections were sufficiently strong, the network commonly generated gamma frequency oscillatory firing patterns in response to random feed-forward (FF) input spikes. This spontaneous oscillatory network activity injects a periodic local current that could amplify a weak synaptic input and enhance the network's responsiveness. When E-E connections were added, we found that the strength of oscillation can be modulated by varying the FF input strength without any changes in single neuron properties or interneuron connectivity. The response modulation is proportional to the oscillation strength, which leads to self-regulation such that the cortical network selectively amplifies various FF inputs according to its strength, without requiring any adaptation mechanism. We show that this selective cortical amplification is controlled by E-E cell interactions. We also found that this response amplification is spatially localized, which suggests that the responsiveness modulation may also be spatially selective. This suggests a generalized mechanism by which neural oscillatory activity can enhance the selectivity of a neural network to FF inputs.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:


Spontaneous Local Gamma Oscillation Selectively Enhances Neural Network Responsiveness

Citation: Paik S-B, Kumar T, Glaser DA ( Spontaneous Local Gamma Oscillation Selectively Enhances Neural Network Responsiveness Se-Bum Paik 0 Tribhawan Kumar 0 Donald A. Glaser 0 Karl J. Friston, University College London, United Kingdom 0 1 Department of Physics, University of California Berkeley, Berkeley, California, United States of America, 2 Department of Molecular & Cell Biology, University of California Berkeley , Berkeley, California , United States of America Synchronized oscillation is very commonly observed in many neuronal systems and might play an important role in the response properties of the system. We have studied how the spontaneous oscillatory activity affects the responsiveness of a neuronal network, using a neural network model of the visual cortex built from Hodgkin-Huxley type excitatory (E-) and inhibitory (I-) neurons. When the isotropic local E-I and I-E synaptic connections were sufficiently strong, the network commonly generated gamma frequency oscillatory firing patterns in response to random feed-forward (FF) input spikes. This spontaneous oscillatory network activity injects a periodic local current that could amplify a weak synaptic input and enhance the network's responsiveness. When E-E connections were added, we found that the strength of oscillation can be modulated by varying the FF input strength without any changes in single neuron properties or interneuron connectivity. The response modulation is proportional to the oscillation strength, which leads to self-regulation such that the cortical network selectively amplifies various FF inputs according to its strength, without requiring any adaptation mechanism. We show that this selective cortical amplification is controlled by E-E cell interactions. We also found that this response amplification is spatially localized, which suggests that the responsiveness modulation may also be spatially selective. This suggests a generalized mechanism by which neural oscillatory activity can enhance the selectivity of a neural network to FF inputs. - Understanding the responsiveness of a cortical neural network is a fundamental requirement for any study of sensory information processing in the brain. Several experiments show that various factors can affect the neuronal response property and information flow in nervous systems: In the primary visual cortex, spiking responses of neurons can be enhanced by slow cortical oscillation [1]. The spike transfer function of thalamo-cortical neurons is modulated by noisy synaptic background activity [2]. Gain of neuronal responses is modulated by background synaptic input [3]. Even at the single cell level, cellular responsiveness is significantly influenced by the presence of voltage fluctuations [4]. It was shown recently that neuronal oscillations can increase response gain and decrease reaction time as a mechanism of attention selection [5]. Cortical neurons commonly show synchronous or oscillatory patterns of activity, which is thought to be important for cortical functions of information flow [6]. In particular, synchronous gamma frequency oscillations (30,70 Hz) have been observed in various neural circuits [7,8,9], and they are thought to provide a temporal structure for information processing in the brain [10]. This gamma-band synchronization can be generated within local networks by coupling between GABAergic I- (inhibitory) interneurons and E- (excitatory) neurons [11,12,13,14], and is related to cognitive functions [15,16], and information delivery [15]. This population activity also has been studied in numerical simulations and mathematical models [12,17,18,19]. Previous analyses have shown that cortical oscillations are generated in networks with appropriate connectivity and can be correlated with the firing phases of E- and I- neurons, but the effect of these oscillations on the neural network responsiveness to external inputs remains elusive. In this research, using a large network model of Hodgkin-Huxley type E- and I- neurons, we study how spontaneous cortical oscillation - particularly in the gamma frequency band - modulates the response property of a neural network. We examine the cortical responsiveness to external FF inputs at the single-spike level because the inputoutput response function for a single input spike is a fundamental feature of neural networks for information processing. A recent study emphasizing the importance of single spike level analysis showed that a significant amount of visual information can be delivered by the very first spike emitted by a neural population [20]. We found that spontaneous cortical oscillation activity noticeably changed the cortical input-output response function. For example, weak inputs that are normally missed in the responses of single neurons, were significantly enhanced by cortical oscillations in the network. This response modulation was similar to the observed effect of the membrane potential oscillation reported in a previous experimental study [21]. More importantly, we found that this cortical response modulation by the oscillation activity In the nervous system, information is delivered and processed digitally via voltage spikes transmitted between cells. A neural system is characterized by its input/output spike signal patterns. Generally, a network of neurons shows a very different response pattern than that of a single neuron. In some cases, a neural network generates interesting population activities, such as synchronized oscillations, which are thought to modulate the response properties of the network. However, the exact role of these neural oscillations is unknown. We investigated the relationship between the oscillatory activity and the response modulation in neural networks using computational simulation modeling. We found that the response of the system is significantly modified by the oscillations in the network. In particular, the responsiveness to weak inputs is remarkably enhanced. This suggests that the oscillation can differentially amplify sensory information depending on the input signal conditions. We conclude that a neural network can dynamically modify its response properties by the selective amplification of sensory signals due to oscillation activity, which may explain some experimental observations and help us to better understand neural systems. was controlled by external feed-forward (FF) input strength variation. This means that the cortical network can self-regulate by differentially amplifying its FF inputs according to their strength, without physiological changes in single neuron properties or structural modulation of interneuron connectivity. We show that this differential amplification results from the modulation of gamma oscillation by cortical E-E neuron interaction. We suggest that this is an important example in which the modulation of gamma oscillation by cortical E-E interactions [12] can significantly change the population responsiveness. We also found that this cortical amplification effect was restricted spatially to an oscillation active region, which enables the spatially-selective tuning of responsiveness to given FF input. We find that this spatial localization is determined by the range of anatomical interneuron connectivity. This is consistent with recent experimental findings concerning the effects of cortical oscillations [1,5,22], and points to aspects of this oscillation which effectively enhance the response selectivity of a neural network to FF inputs. Cortical Network Model We examined a variety of neural network activity patterns produced under different conditions. We used a cortical network model in which E- and I- neurons interact with each other via lateral synaptic connections. We constructed isotropic local cortical connections, using physiological and earlier modeling data [23,24]. We varied the strength of each type of cortical interaction (E-E, E-I, I-E and I-I) in order to study different cortical connectivity conditions. Artificially generated random spikes were injected into the cells in the center area (diameter,450 mm, ,500 neurons: E- 377, I- 113) of this network model (1 mm by 1 mm, consisting of ,3300 neurons: E- 2500, I- 841) to simulate localized FF spike input to the network. The actual spike pattern for each neuron was generated by a Poisson process with constant mean firing rate, and FF input strength (The amount of intracellular conductance fluctuation caused by a single FF input spike) was varied within the range 5,100 mS/cm2, as a control parameter. By performing many simulations of different cortical parameters with FF input rates in the range 5,40 spike/s, we observed several different types of cortical activity patterns. Cortical E- and I-Cells Activity: Generation of Gamma Oscillation Gamma oscillation can be generated by interactions among Eand I- cells: The E- cells synchronize the I- cells, and the gammamodulated I- cells drive E- cells to generate gamma frequency rhythms in the network [11,14,25]. Our simulations agreed with earlier studies that gamma oscillations are generated by applying E-I and I-E cortical connections; I- cells were synchronized by E-I connections first (Figure 1A), then I-E connections generated gamma rhythm in E-cells via periodic inhibitory activity (Figure 1B) as in the PING model [11,14]. I- cells fire at higher rate than E- cells, just as fast-spiking cells fire at higher rates than regular-spiking cells [10]. The relative firing phase of E- and I- cell (Figure 1A and 1B) also showed a similar phase relation as reported in the previous experiments [25,26,27]. E- cells fire 3,5 ms before I- cells fire, with an approximately 70-degree phase difference in a gamma cycle. We extend the previous studies by also explicitly considering EE and I-I cortical connections. In ref. 25, Morita et al. showed that gamma- modulated cortical excitation increases the firing rate of E- cells and decreases the E-I firing phase difference. Based on these observations, they predicted that the gamma- modulated E-E cell interaction will suppress the cortical oscillation. In this study, we found that the gamma- modulated E-E coupling does not extinguish the cortical oscillation. When E-E connections were turned on, the firing rate of E- cells was increased and the E-I phase difference was diminished (Figure 1D and 1E), as shown in ref. 25. But these changes did not actually suppress the cortical oscillation. Instead, they caused a transition of operating mode such that the oscillation frequency dropped to a low gamma or near beta rhythm [28]. This transition could not be observed using the methods reported in ref 25, because they have only static data of pre- and post- synaptic activities to examine static equilibrium conditions. The approximately zero E-I firing phase difference is an important feature of E-E coupling, and is distinguishable from that of the case involving no or little E-E coupling (,70-degrees is usually observed in experiments [26,27]), so we call the former situation E-E interaction modulated E-I phase, in contrast to the normal gamma oscillation E-I phase. We will analyze below how this phase modulation results from E-E coupling. The difference between the firing rates of E- and I- cells is also diminished to some degree by E-E interactions, and the population oscillation frequency is also lowered [28,29]. We maintain this non-zero E-E interaction throughout the following simulations. We also allow I-I connections, which reduce the firing rate of Icells to some extent (Figure 1C and 1D). In this study, we do not examine the role of I-I connection in detail, and the I-I connection strength was not varied. Control of Spontaneous Cortical Network Activity Strength When the cortical interactions (E-E, E-I, I-E and I-I) are completely turned off and each neuron is driven by FF inputs only, there is no correlated behavior in the network and the average network firing rate simply follows the instantaneous input firing rate (Figure 2A, FF). On the other hand, when cortical interactions (E-E, E-I, I-E and I-I) are introduced, the network exhibits some synchronized patterns depending on FF input and cortical connectivity parameters. For example, with a moderate FF input strength (,35 mS/cm2) and rate (,40 spikes/s), the neural population shows an oscillation pattern in its firing rate for a broad range of cortical connectivity parameters. The crosscorrelogram among cortical spikes shows a clear oscillatory pattern for both E- and I- cells in this case (Figure 2B, OA). Generally, the amplitude of the instantaneous output firing rate of the neural population depends on the FF input firing rate. The frequency of the oscillation is somewhat affected by the FF input strength and cortical connectivity parameters, but the oscillation frequencies are mostly within the gamma band range (25,50 Hz), as in earlier experimental observations and theoretical models [8,9,30,31]. For some parameter range of cortical connections, this gamma oscillation can be very small. When the FF input spike rate was low (,10 spikes/s) with the other parameters unchanged, the oscillating firing pattern in E- cells became barely detectable even though the oscillation in I- cells was maintained to some extent (Figure 2C, OI). We examined the phase of the firing pattern of E-I cells in three cases (FF, OA, and OI). The E-I spike correlogram in the OA network showed that the effect of interactions among E-E cells is significant because the E- and Icells fire with little phase difference (Figure 2B (iii)), similar to Figure 1D and 1E (E-E interaction modulated E-I phase), while the other two cases (FF and OI) showed no clear phase relation (Figure 2A (iii) and 2C (iii)). Since OA and OI networks have the same parameters except for the FF input firing rate, we conclude that FF input firing rate is a crucial factor in determining E-E interaction strength. When E-E cortical connections are very strong, extremely sharp cortical spike synchronization is generated, resulting in spatially propagating waveform patterns [32,33,34]. The amplitude of the instantaneous output firing rate was almost constant, independent of FF input rate. In this case, the neural response was not controlled very much by FF inputs at each moment but mostly by the spatial cortical bursting pattern. We observed that the frequency of this periodically propagating pattern is in the beta oscillation range induced by E-E cortical interactions [28]. We do not examine the generation or the effects of periodically propagating patterns any further here. In the following simulations, we chose cortical conditions such that the system didnt enter this phase for the FF input strengths tested. To examine the responsiveness of the neural network, we chose a range of parameters that provided moderate and stable oscillatory behavior on application of FF input spikes. We tested this condition using sinusoidal time-varying input rates with a peak amplitude range of 0,60 spikes/s and peak frequency of 5,10 Hz. The network rapidly and reliably restored its oscillating state whenever the FF input firing rate became greater than some level (,20 spikes/s), and the oscillations diminished significantly and very rapidly when the input rate fell below that level. Throughout this research, we did not change any individual neuronal property (e.g. ionic channel parameters). Network Activity Increases the Cortical Response to Weak FF Inputs and Modifies the Response Function We compared the neural responsiveness for the following three states of network activity: (i) network with no cortical connectivity (and no spontaneous network activity) (FF), (ii) laterally connected cortical network with spontaneous oscillation activity (OA) and (iii) the same network as in (ii) but with little or no oscillation activity (OI). We generated random FF spike inputs by a Poisson process and provided this input to the center area (diameter,450 mm) of the 1 mm by 1 mm network. All the response properties were measured within this small center area in order to avoid finite size effects from the network boundary. Neurons outside this area received no FF input. FF input strength was varied from 5 to 100 mS/cm2, and FF spike rate was kept constant at 40 spikes/s for (i) and (ii). For (iii) the oscillation inactive case, the FF input rate was lowered to 10 spikes/s in order to maintain minimal oscillations while still providing enough input spikes to allow measurement of the network responsiveness. Figure 3A shows FF input and cortical output spike trains, the membrane voltage, and the internal current fluctuation of a sample neuron with and without network oscillation activity. Each FF input spike induces a synaptic conductance fluctuation in a cortical neuron and the corresponding intracellular current fluctuation. When cortical interactions are turned off, the response of each neuron depends only on the direct FF input (FF network, Figure 3A (i)). When the FF input strength is weak (25 mS/cm2), a single spike input could not produce a large enough current fluctuation to generate an output spike. Only when two or more inputs are temporally paired within a short time interval (,,10 ms), can the conductance fluctuations from each spike overlap to generate an output spike (Figure 3A (i) ***), as found in the measurements of correlated inputs in a previously reported experiment [35]. In this case, there is a threshold FF input strength (Sthresh) that determines whether each single FF input spike can generate an output spike or not. We calculated the responsiveness of network neurons, using a cross-correlation method (Figure 3B) [35,36]. The responsiveness of the cortical network was defined as (net peak integral)/(number of FF input spikes), where net peak integral is the total area of the maximum response peak above the background activity in each crosscorrelogram. Measuring this quantity at each FF input strength produces the generalized response function of the network (Figure 4A). Since any subthreshold single input produces no response when cortical connections are turned off (Figure 3B unpaired), the response to a single FF input spike is like a step function, similar to the measured thalamo-cortical transfer function in the absence of noisy background activity [2], providing the cell with a simple spike relay capability with limited encoding abilities. When the cortical connections are turned on, each neuron receives cortical synaptic inputs from other neurons as well as FF inputs. In the presence of the spontaneous network activity of synchronized oscillating patterns, the cortical spike inputs that each neuron receives are also oscillating (Figure 3A (ii)). The current fluctuation due to cortical interactions is added to that by FF inputs, and as a result, a single-spike sub-threshold FF input can produce an output spike response with the help of this additional cortical activity (Figure 3A (ii) *). This input amplification depends on the phase of the cortical activity. When FF input timing is out of phase with the cortical oscillation, it fails to produce an output spike for lack of additional cortical current (Figure 3A (ii) **) just as in the FF network. The general response function of the network to a single unpaired input spike is plotted in (Figure 4A), as a function of FF input strength. Different from the step-like FF case, the slope of the response function changes more gradually, with a plateau near the FF threshold input strength (Sthresh,30 mS/cm2). This broader and more gradual change of response function can provide richer encoding capabilities [2]. Next, we examine how responsiveness changes when the oscillation is inactive while the connectivity of neural population is kept the same. We lower the FF input rate to 10 spike/s, so that the spontaneous oscillation is almost absent. All the other parameters including cortical connectivity are kept the same so that each synaptic interaction via spike delivery in the network can give the same contribution as before. This time, neurons do not exhibit enhanced responsiveness for weak inputs; each neuron still experiences some conductance change by cortical interaction, but its contribution is negligible. The network response character is similar to that of the FF network (Figure 3A). Any unpaired weak inputs cannot generate a cortical response, losing its information. We found the response function of neurons is almost the same as that of the FF network (Figure 4A) when the cortical connections are turned off. Oscillation Strength is Self-Modifiable: Cortical Modulation Is Controlled by FF Input Strengths In Figure 4A, the absolute difference in responsiveness between the oscillatory network and the FF network is large when FF input strength is weak (near the FF response threshold, Sinput,30 mS/ cm2). This difference becomes smaller as the FF input strength increases. Finally the two response functions become equal at very strong inputs (Sinput.80 mS/cm2). In other words, the cortical amplification due to oscillation activity is large for weak FF inputs, and becomes insignificant as inputs become strong. We examined how this cortical response modulation is related to the change of oscillation strength. Figure 5A shows values of the network oscillation as a function of FF input strength. For weak FF inputs (Sinput = 25, 35 mS/cm2), the cortical oscillation is strong but it becomes weaker for a medium input strength (Sinput = 50 mS/cm2), and almost disappears for a strong input (Sinput = 80 mS/cm2). This change in the oscillation strength is almost proportional to the extent of the cortical response modulation which can be roughly defined as the difference between the response function and that of the FF network (Figure 5B). As shown by examining the E-I firing phase pattern (Figure 5A (iii)), the oscillation modulation is significantly affected by the E-E cell interaction. When the FF input is strong, the E-I spike correlogram shows normal E-I phase difference of the gamma oscillation (Figure 5A ***). On the other hand, when the FF input is weak, the relative E-I firing phase is near zero (E-E interaction modulated), showing that the network is affected by strong E-E interaction (Figure 5A (iii), Sinput = 25, 35 mS/cm2). For some range of FF input strength, two different peaks coexist in the E-I correlerogram (Figure 5A * and **). Whenever the E-I firing phase is significantly affected by the E-E interaction (,70-degreesR,0degrees), the cortical oscillation becomes strong and the responsiveness of the network to weak FF input is enhanced. When the EE connections were turned off (Figure 6), there was no E-I firing phase modulation (Figure 6 (iii)), and the cortical oscillation was far less affected by the FF input strength (Figure 6 (i)). For weak FF input, the cortical oscillation almost disappeared (Figure 6 (i) input strength 25 mS/cm2) and consequently the responsiveness did not exhibit any enhancement. We examined how the E-I firing phase is modulated from ,70degrees to ,0-degrees by the excitatory interactions in the E-E couplings. We turned on only excitatory cortical connections (E-E and E-I) and measured the firing patterns of E- and I- cells (Figure 7). In this condition, the network generated the periodically propagating waves near beta rhythm that we observed in the earlier part of this study (when E-E connections were relatively stronger than other types of connections). In each wave cycle, a small number of E- cells fired due to the FF inputs (Figure 7A i). Then they excited nearby E- and I- cells through the E-E and the E-I connections. The stimulated E- and I- cells fired almost simultaneously, or I- cells fired slightly before E- cells fire (Figure 7A ii and iii) because E-I connections are stronger than EE connections. The firing of the E- cells continued for a while due to the chain reaction of E-E couplings and the I- cells occasionally produced spike doublets by this sustained excitation (Figure 7A iv and 7B). These inhibitory spike doublets are frequently observed during long range synchronization processes in neural networks [9], usually along with the delayed excitation from distant cells. In our study the ERERI route could provide the delayed or sustained excitation. As a result of these E- and Ifiring patterns, the E-I firing phase was ,0-degrees on average. When the average excitation to each I- cell is strong enough to produce I- spike doublets, the E-I phase correlogram has two peaks around 0-degrees (Figure 5A * and **). If the excitation to Icells is not enough, or the inhibition of E- cells is strong so that the sustained E- cell activity cannot drive I- cells to produce the second spike of spike doublets, the E-I phase correlogram does not show clear peak separation (Figure 5A (iii) input strength 25 mS/cm2). In any case, the average E-I phase difference is close to 0-degrees, clearly different from ,70-degrees normal gamma phase distribution. The cortical response modification caused by this oscillation can be explained as an effect of membrane potential oscillations [21]. In addition, the modulation of the E-E interaction by the FF input strength can be understood as follows: as shown in Figure 3A, when the FF input strength is weak, each E- neuron can respond to individual FF input only with the help of the collective cortical activity of E- cells whose contribution is intrinsically periodic (oscillatory). The probability of a cortical response to an FF input spike strongly depends on the phase of oscillation of the E- cells (Figure 3A * and **), and this dependence is significantly strengthened by the E-E coupling loop. In the beginning of a gamma oscillation cycle, only a small number of E- cells fire together, but they trigger an E-E coupling loop which drives more E- cells to fire simultaneously. As a result, the peak firing rate of E- cells in a gamma cycle is quite high, causing much higher spiking probability near the oscillation peak. Therefore, (i) the cortical responsiveness is dependent on the oscillation phase, and (ii) the gain or the cortical amplification is proportional to the oscillation strength. When an FF input is strong enough to produce an output spike in each E- cell, there is no significant dependence of the E- cell response on the periodic cortical activity and E-E couplings, and each individual E- cell responds directly to its FF input pattern, independently of the network activity. Since the FF input spike train was generated by a random Poisson process, the cortical response pattern is also random and not oscillatory. Although the average output firing rate generally increases with increasing FF input strength, the cortical response modulation actually decreases. The oscillations are depressed by increasing the FF input strength while keeping the FF input rate constant. The average number of FF input spikes does not decrease but the oscillation is depressed, and the response function converges to the FF result. Although cortical oscillation is initially established by the E-I and the I-E interactions, the E-E cell interaction is crucial for responsiveness modulation because its strength strongly depends on the FF input strength. Since the cortical response modulation is controlled by the FF input strength in the system via the self modification of network oscillation, this seems to be a very effective automatic gain control system that does not require any synaptic adaptation or learning mechanism [37]. Next, we examine how the oscillation activity affects the response delay of the network. Generally, the response time (time interval between an FF input spike and a corresponding output spike) is relatively long (,10 ms) for weak inputs and becomes shorter as input strength increased for all cortical states (Figure 4B). But there is a significant difference in the average response time depending on the network activity state. Figure 4B shows that the average response delay is shorter during spontaneous oscillations than for the other two cortical activity states considered, especially when the FF input is weak. As input strength increases, the response of the oscillating network becomes similar to that of the oscillation inactive network, but still faster than that of the FF network. For a sufficiently strong FF input, the response time delay Figure 6. Oscillation modulation without E-E interaction. (i) (ii) The cortical oscillation strength is weaker, and far less modified by the FF input, compared with Figure 5A. The oscillation pattern in E- cells disappears for both weak and strong FF inputs. (iii) The relative E-I firing phase does not change for weak inputs in this case, unlike that shown in Figure 5A (iii). doi:10.1371/journal.pcbi.1000342.g006 curve of the oscillating network and the oscillation depressed network were almost the same, with a delay of about 4 ms, agreeing with the experimentally known value for monosynaptic connections [35]. This is still faster than that measured for the FF network. Cortical interaction adds some positive current to each neuron, and this additional current causes the membrane voltage to reach the action potential threshold faster. In the presence of cortical oscillation, the amount of cortical current added is larger in the positive phase of the oscillation and on average, the larger the oscillation, the greater the average positive current that is added to the cell. That leads to faster responses than observed for the simple FF network. As the FF input strength increases above 35 mS/cm2, the oscillation amplitude decreases (Figure 5A), and the net average positive current added is less. The response time difference between oscillatory and FF networks decreases accordingly (Figure 5B). EE cell interactions control the cortical oscillation strength, and the extent of response delay modulation and the response modulation of the network are in turn proportional to the cortical oscillation strength. Cortical Response Modification Is Localized: The FF Input Can Be Selectively Amplified in Particular Regions To examine the spatial localization of the cortical gain resulting from the spontaneous oscillation of network activity we used a network model four times the size (2 mm by 2 mm, 13000 neurons) of that used in the studies described above. In the center region (diameter,450 mm, including ,500 neurons, Figure 8A G1), the cells were activated as before with a FF input rate of 40 spikes/s, a rate that was also used to set up the spontaneous network activity oscillations for the simulations described above. The surrounding neurons (Figure 8A G2,G11) received a signal with an FF input rate of 10 spikes/s, a rate at which spontaneous oscillations are hardly evident previously. The cortical neural connectivity is the same everywhere so that the center and the surrounding neurons could interact with each other. There are differences in the responsiveness of the network between the central region that shows strong spontaneous oscillations and the surrounding regions with little to no oscillations. In other words, the responsiveness modification by the spontaneous oscillation can be localized. A control simulation where all cortical neurons received an FF input of 10 spike/s in both the center and the surrounding regions is an approximation of the neuronal property at infinite distance from the center region (region G in Figure 8A). For measurement purposes the network is divided into circular annuli (Figure 8A G2,G11). Each annulus contained 400 to 500 neurons. Figure 8A shows the network response function and the response delay time for single unpaired FF spike inputs in each region. Figure 8B shows the measured values in all the regions for FF input strengths of 25 and 35 mS/cm2 which are slightly smaller and slightly larger, respectively, than the FF response threshold value (Sthresh = 30 mS/cm2, Figure 4A). When the FF input strength was weaker than Sthresh, the responsiveness and the response delay pattern gradually moved from the oscillation-active center region (G1) to approach the control oscillation suppressed case (G). Subsets of the Figure 8A graphs show that the quantities measured in the surrounding regions G2,G6 are continuously distributed in the interval between the values obtained for G1 and G. Values measured in surrounding regions G7,G11 showed negligible differences from those measured for G, and are not plotted in the figure. The effect of oscillatory network activity is restricted to a central local region about 500 mm in radius. When the FF input strength is larger than Sthresh but not very large (30 mS/cm2,Sinput,50 mS/cm2), the spatial attenuation of the cortical activity effect is more apparent (Figure 5B, Sinput = 35 mS/ cm2). In this case, all surrounding regions, even including the nearest region G2, are clearly separated from G1, and are close to G. In this case, the cortical oscillation activity effect is restricted to the central activated region (G1). For very strong FF inputs (50 mS/cm2,Sinput), all response properties converge to control group (G) behavior as expected since under that condition oscillatory behavior is barely evident even in the central region. It seems reasonable to expect that the spatial localization of gamma oscillation is dependent upon the range of single neuron synaptic connectivity so that the shorter synaptic connection range, the smaller the surrounding area that is affected by the oscillation in the center. To verify this expectation, we reduced the excitatory and the inhibitory synaptic interaction range from the initial value (radius of 200 mm for E- cells, 100 mm for I- cells), keeping the ratio of E- and I- range the same. In Figure 8C, the range of oscillation effect is proportional to the single cell synaptic connection range, as expected. This suggests that the effective range of gamma oscillations is strongly dependent upon the details of the anatomical connectivity of neurons in experimental observations. Thus the effect of cortical oscillation is fairly well localized for weak and moderate FF input strengths. Neural response properties are modified only within or near the area in which spontaneous oscillation is activated. This suggests that spatially selective cortical response modification is possible. Spontaneous cortical oscillations can be locally induced by spatially localized thalamic inputs, and the cortical response character can also be selectively tuned within a limited region. As shown above, cortical responsiveness to a single unpaired FF input spike is enhanced by synchronized gamma frequency oscillations with the help of E-E neuronal interaction, particularly for weak FF input strengths. The cortical response modulation is proportional to the oscillation strength which is controlled by the network itself depending on FF input strengths. This cortical effect is spatially localized fairly tightly depending upon the range of cell connectivity, suggesting that each cortical area can be tuned selectively by well localized FF inputs. These findings are relevant to previous experimental and simulational results, and improve our understanding of the response character of the visual pathway. Gamma Oscillation with E-E Cells Interaction It was previously thought that cortical E-E activity interaction is not essential for gamma rhythm generation but can modify the oscillation frequency and the phase of cell firing pattern [11,28,29]. A recent experimental study showed that gamma rhythms in E-E cells activity can change the E- and I- cells firing phase distribution, and suggested that it may suppress the gamma oscillation [25]. We have shown that recurrent E-E interaction significantly modulates the oscillation frequency, the firing phase distribution of E- and I- cells, and the oscillation strength, without extinguishing the cortical oscillation. Moderate levels of E-E interaction generally strengthen the oscillation, causing the ,0degrees E-I firing phase and the lower oscillation frequency (near beta range). As a result, it modulates the cortical response function that is clearly relevant to encoding/decoding of information. The fact that the effect of E-E interaction is self regulatory for FF input strength variation suggests a useful mechanism for the cortical gain control, without neuronal feedback loops from the visual cortex to earlier visual stages. This suggests a general mechanism by which the same types of neurons in different cortical layers can have different properties due to the different coupling within each layer. In previous studies, it was reported that the upper and the lower layers of the cortex can have different oscillation characters and phase response properties [38,39,40]. Our observations about the phase and the frequency modulation of cortical oscillation by E-E coupling, suggest that different neuronal properties in different cortical layers may originate from the different lateral connectivity (especially E-E coupling) strength in each layer. For example, neurons in the different hippocampus regions (CA1 and CA3) show noticeably different firing phase distribution and correlation activity patterns in gamma oscillations [27]. Considering that the E-E couplings are significant in CA3 [41] while they are sparse in CA1 [42], this may be a situation in which the E-E coupling property plays an important role in the modulation of neural activity, as suggested above. Responsiveness Tuning with Fast and Slow Cortical Oscillations Previous experimental work has shown that (i) the visual responsiveness of the cortical network is significantly enhanced by slow cortical oscillation [1] and (ii) the phase of slow theta rhythm (4,8 Hz) oscillation modulates the high frequency gamma (80,150 Hz) band oscillation power [22]. Here we provide a clue to a mechanism for modification of neural response properties by various types of synchronized cortical network activities. As shown above in the results section, when the gamma frequency oscillation is generated, it enhances the neuronal responsiveness to weak FF input. If this gamma power is modulated by slower (theta or lower frequency) rhythm, then the network responsiveness could depend on the phase of this slow oscillation. This suggests a simple and consistent basis for the modulation of high frequency oscillation power by lower frequency activity. Some previous experimental work has shown that the power and the frequency of gamma oscillations can be modulated by various factors such as the operation of fast spiking interneurons and some neuromodulators [10]. In our simulation, the strength of the gamma oscillations can be significantly modified by changing the strength or rate of the FF input, with the help of E-E interactions but without changing any individual neuronal properties or network connectivity features. In addition, our results on the gamma oscillation modulation mechanism do not require modifying the FF input correlation pattern, learning/adaptation behaviors [37] in cortical synapses, or cortico-thalamic feedbacks [43]. If the slow frequency oscillations affect the FF input strength or the input rate locally within the network, the gamma oscillations will be readily modulated. The consequent modulation of the responsiveness will depend on the phase of the low frequency oscillation [1]. Controlled and Selective Responsiveness Modification This proposed responsiveness tuning mechanism does not require any accompanied background activity control. Therefore it is simpler than those gain control models suggested in the previous reports [2,3,4] that are mostly dependent on the modulation of the background activity. An important advantage of the present model is that the response modification can be dynamically selective according to the FF input strength and the FF input projection range variation. Since the cortical gain is modulated by FF input strength, the cortical amplification is selective to FF input strength. The system effectively determines the optimized gain via modulations of oscillation power spontaneously, and can avoid unnecessary adjustments even without any delayed feedbacks to thalamus or thalamo-cortical neurons. Moreover, the tuning is spatially localized to distances of less than about 500 mm for weak FF inputs, and less than 50 mm for strong inputs. This is comparable or even smaller than the size of the receptive field of single neuron in the mammalian primary visual cortex [44]. Therefore responsiveness modulation can be spatially selective, and this is more effective than mechanisms proposed in previous studies [2,3,4] where the cortical modulations were generally global and the gain optimization could not be achieved easily. We also suggest that this mechanism might be applicable to the functional modulation of the relative weight between thalamic inputs versus cortical inputs to the visual cortex neurons [45]. Attended or Awake State Animal Behavior with Oscillations In some experiments with monkeys, when attention is directed, visual sensitivity increased due to the increased synchronization between the visually evoked potentials and the stimulus [46]. In another report, neurons activated by the attended stimulus showed increased gamma frequency synchronization [47]. Considering the response enhancement by gamma oscillation in our model together with these experimental measurements, spontaneous gamma band activity seems to play an important role for regulating the information flow in the visual nervous system, as suggested previously [6,16]. Our findings support these experimental observations, and may suggest a new mechanism for attention modulation that is different from that of other models [48,49,50]. Cortical Neural Network Modeling This neural network model consists of a two-dimensional layer of coupled neurons, each modeled as a Hodgkin-Huxley neuron with Na+, K+ and Cl2 ion channels and E- and I- synaptic conductance channels. 75% of the neurons are E- and 25% are Ias suggested by experimental data [51], and adopted in other publications [24]. We used networks of two sizes in this research: (i) 1 mm by 1 mm, including ,3300 neurons for responsiveness simulation and (ii) 2 mm by 2 mm, including ,13400 neurons for locality simulations. The membrane potential of the jth neuron, vj , is determined by where s is E or I depending upon whether the neuron is E- or I-, respectively, C is the membrane capacitance, and gL is the leakage conductance. gjsE is the synaptic conductance of the jth neuron, Eor I- as specified by s, providing the cortical inputs from the neurons in its spatial neighborhood, and gjsI is the synaptic conductance of the jth neuron providing the I- input from the neurons in its spatial neighborhood. We used the commonly accepted biophysical parameters [52,53]: the capacitance C = 1026 Fcm22, the leakage reversal potential VL = 270 mV, the Na+ equilibrium potential VNa = 55 mV, the K+ equilibrium potential VK = 280 mV, the E- reversal potential Vexc = 0 mV, the I- reversal potential Vinh = 280 mV, the leakage conductance gL = 50*1026 Scm22, and Hodgkin-Huxley Na+ and K+ conductances gNa = 120*1023 Scm22, gK = 36 *1023 Scm22. The Hodgkin-Huxley ion channel conductance takes the form [54]: GNa~gNam3h, GK ~gK n4 where m, h and n denote the channel activation variables. anv~0:0110{v bnv~0:125 exp{v=80 amv~0:125{v {1 , bmv~4 exp{v=18 ahv~0:07 exp{v=20, bhv~1 A jth neurons synaptic conductance is given by: Djs,sk, GE t{tlk and GI t{tlk are the spatial, the E- temporal, and the I- temporal kernels of cortical interaction, respectively, which describe the contribution of lth spike from kth neuron to the jth neuron. For cortico-cortical connection, the spatial kernel in the synaptic conductance equation takes the form: Djs,sk~Css exp { ~rj{~rk where ~rj and ~rk are the jth and kth neurons spatial positions respectively. The decay constant, lE is 200 mm (for Econnections) and lI is 100 mm (for I-). The temporal kernel in the equation is set to be and the time constants t1, t2 in milliseconds were chosen as (3, 1) for E- and (7, 1) for I- synapses where Css9 and Cts9 are normalization constants chosen so that that the sum of the contributions of the two kernels would sum to unity. We assume spatially isotropic local connections with a range of 200 mm in radius for E- and 100 mm in radius for I- synapses. Wss9 are strengths of synaptic connections for the neuron pair of type (s, s9). If all Wss9 = 0, the network is then equivalent to the simple FF model. When the cortical synaptic connections were turned on, these values ratios were set to (WEE, WIE, WEI, WII) = (0.03: 0.06: 0.12: 0.12) throughout the simulation. This condition was achieved from the parameter search simulations shown in the first part of results section. The contribution to the E- conductance by the FF input spikes was given by: Sinput is the weighting factor for FF input synaptic strength and gjinput was varied within 5,100 mS/cm2, throughout the simulations reported here. The temporal kernel GE t{tjl has the same form as the E- temporal cortical kernel given above. The spike timings, tlj , of input were generated by Poisson processes. Computer Simulation and Data Analysis Techniques All of our simulations were coded using the GENESIS 2.3 environment [55], and performed with a Pentium IV PC system. Simulation outputs were analyzed using Matlab R2006b scripts. Acknowledgments We are grateful to Yang Dan and Michael DeWeese for helpful comments on an earlier version of the manuscript. 1. Haider B , et al. ( 2007 ) Enhancement of visual responsiveness by spontaneous local network activity in vivo . J Neurophysiol 97 : 4186 - 4202 . 2. Wolfart J , et al. ( 2005 ) Synaptic background activity controls spike transfer from thalamus to cortex . Nat Neurosci 8 : 1760 - 1767 . 3. Chance FS , Abbott LF , Reyes AD ( 2002 ) Gain modulation from background synaptic input . Neuron 35 : 773 - 782 . 4. Ho N , Destexhe A ( 2000 ) Synaptic background activity enhances the responsiveness of neocortical pyramidal neurons . J Neurophysiol 84 : 1488 - 1496 . 5. Lakatos P , et al. ( 2008 ) Entrainment of neuronal oscillations as a mechanism of attentional selection . Science 320 : 110 - 113 . 6. Salinas E , Sejnowski TJ ( 2001 ) Correlated neuronal activity and the flow of neural information . Nat Rev Neurosci 2 : 539 - 550 . 7. Gray CM ( 1994 ) Synchronous oscillations in neuronal systems: mechanisms and functions . J Comput Neurosci 1 : 11 - 38 . 8. Whittington MA , Traub RD , Jefferys JGR ( 1995 ) Synchronized oscillations in interneuron networks driven by metabotropic glutamate receptor activation . Nature 373 : 612 - 615 . 9. Traub RD , Whittington MA , Stanford IM , Jefferys JGR ( 1996 ) A mechanism for generation of long-range synchronous fast oscillations in the cortex . Nature 383 : 621 - 624 . 10. Bartos M , Vida I , Jonas P ( 2007 ) Synaptic mechanisms of synchronized gamma oscillations in inhibitory interneuron networks . Nat Rev Neurosci 8 : 45 - 56 . 11. Whittington MA , Traub RD , Kopell N , Ermentrout B , Buhl EH ( 2000 ) Inhibition-based rhythms: experimental and mathematical observations on network dynamics . Int J Psychophysiol 38 : 315 - 336 . 12. Traub RD , Jefferys JGR , Whittington MA ( 1997 ) Simulation of gamma rhythms in networks of interneurons and pyramidal cells . J Comput Neurosci 4 : 141 - 150 . 13. Battaglia D , Brunel N , Hansel D ( 2007 ) Temporal decorrelation of collective oscillations in neural networks with local inhibition and long-range excitation . Phys Rev Lett 99 : 238106 . 14. Bogers C , Epstein S , Kopell NJ ( 2005 ) Background gamma rhythmicity and attention in cortical local circuits: a computational study . Proc Natl Acad Sci U S A 102 : 7002 - 7007 . 15. Womelsdorf T , et al. ( 2007 ) Modulation of neuronal interactions through neuronal synchronization . Science 316 : 1609 - 1612 . 16. Schoffelen J , et al. ( 2005 ) Neuronal coherence as a mechanism of effective corticospinal interaction . Science 308 : 111 - 113 . 17. Tiesinga PHE , Fellous JM , Jose JV, Sejnowski TJ ( 2001 ) Computational model of carbachol-induced delta, theta, and gamma oscillations in the hippocampus . Hippocampus 11 : 251 - 274 . 18. Bojak I , Liley DTJ ( 2007 ) Self-organized 40 Hz synchronization in a physiological theory of EEG . Neurocomputing 70 : 2085 - 2090 . 19. Coombes S , et al. ( 2007 ) Modeling electrocortical activity through improved local approximations of integral neural field equations . Phys Rev E 76 : 051901 . 20. Gollisch T , Meister M ( 2008 ) Rapid neural coding in the retina with relative spike latencies . Science 319 : 1108 - 1111 . 21. Volgushev M , Chistiakova M , Singer W ( 1998 ) Modification of discharge patterns of neocortical neurons by induced oscillations of the membrane potential . Neuroscience 83 : 15 - 25 . 22. Canolty RT , et al. ( 2006 ) High gamma power is phase-locked to theta oscillations in human neocortex . Science 313 : 1626 - 1628 . 23. Yoshioka T , Blasdel G , Lebitt J , Lund J ( 1996 ) Relation between patterns of intrinsic lateral connectivity, ocular dominance, and cytochrome oxidasereactive regions in macaque monkey striate cortex . Cereb Cortex 6 : 297 - 310 . 24. McLaughlin D , Shapley R , Shelley M , Wielaard DJ ( 2000 ) A Neuronal network model of macaque primary visual cortex (V1): orientation selectivity and dynamics in the input layer 4Calpha . Proc Natl Acad Sci U S A 97 : 8087 - 8092 . 25. Morita K , Kalra R , Aihara K , Robinson HPC ( 2008 ) Recurrent synaptic input and the timing of gamma-frequency-modulated firing of pyramidal cells during neocortical ''UP'' states . J Neurosci 28 : 1871 - 1881 . 26. Hasenstaub A , et al. ( 2005 ) Inhibitory postsynaptic potentials carry synchronized frequency information in active cortical networks . Neuron 47 : 423 - 435 . 27. Csicsvari J , Jamieson B , Wise KD , Buzsaki G ( 2003 ) Mechanisms of gamma oscillations in the hippocampus of the behaving rat . Neuron 37 : 311 - 322 . 28. Whittington MA , Traub RD , Faulkner HJ , Stanford IM , Jefferys JGR ( 1997 ) Recurrent excitatory postsynaptic potentials induced by synchronized fast cortical oscillations . Proc Natl Acad Sci U S A 94 : 12198 - 12203 . 29. Traub RD , Jefferys JGR , Whittington MA ( 1997 ) Simulation of gamma rhythms in networks of interneurons and pyramidal cells . J Comput Neurosci 4 : 141 - 150 . 30. Bauer R , Brosch M , Eckhorn R ( 1995 ) Different rules of spatial summation from beyond the receptive field for spike rates and oscillation amplitudes in cat visual cortex . Brain Res 669 : 291 - 297 . 31. Ermentrout GB , Kopell N ( 1998 ) Fine structure of neural spiking and synchronization in the presence of conduction delays . Proc Natl Acad Sci U S A 95 : 1259 - 1264 . 32. Winfree AT ( 2001 ) The Geometry of Biological Time . New York : Springer. pp 411 - 442 . 33. Richardson KA , Schiff SJ , Gluckman BJ ( 2005 ) Control of traveling waves in the mammalian cortex . Phys Rev Lett 94 : 028103 . 34. Benucci A , Frazor RA , Carandini M ( 2007 ) Standing waves and traveling waves distinguish two circuits in visual cortex . Neuron 55 : 103 - 117 . 35. Usrey WM , Alonso J , Reid RC ( 2000 ) Synaptic interactions between thalamic inputs to simple cells in cat visual cortex . J Neurosci 20 : 5461 - 5467 . 36. Perkel DH , Gerstein GL , Moore GP ( 1967 ) Neuronal spike trains and stochastic point processes . II. Simultaneous spike trains . Biophys J 7 : 419 - 440 . 37. Dan Y , Poo M ( 2006 ) Spike timing-dependent plasticity: from synapse to perception . Physiol Rev 86 : 1033 - 1048 . 38. Buffalo EA , Fries P , Desimone R ( 2004 ) Layer-specific attentional modulation in early visual areas . Soc Neurosci Abstr 30 : 717.6. 39. Sejnowski TJ , Paulsen O ( 2006 ) Network oscillations: Emerging computational principles . J Neurosci 26 : 1673 - 1676 . 40. Tsubo Y , Takada M , Reyes AD , Fukai T ( 2007 ) Layer and frequency dependencies of phase response properties of pyramidal neurons in rat motor cortex . Eur J Neurosci 25 : 3429 - 3441 . 41. Amaral DG , Ishizuka N , Claiborne B ( 1990 ) Neurons, numbers and the hippocampal network . Prog Brain Res 83 : 1 - 11 . 42. Knowles W , Schwartzkroin PA ( 1981 ) Local circuit synaptic interactions in hippocampal brain slices . J Neurosci 1 : 318 - 322 . 43. Andolina IM , Jones HE , Wang W , Sillito AM ( 2007 ) Corticothalamic feedback enhances stimulus response precision in the visual system . Proc Natl Acad Sci U S A 104 : 1685 - 1690 . 44. Jones JP , Palmer LA ( 1987 ) The two-dimensional spatial structure of simple receptive fields in cat striate cortex . J Neurophysiol 58 : 1187 - 1211 . 45. Nauhaus I , Busse L , Carandini M , Ringach DL ( 2009 ) Stimulus contrast modulates functional connectivity in visual cortex . Nat Neurosci 12 : 70 - 76 . 46. Kim YJ , et al. ( 2007 ) Attention induces synchronization-based response gain in steady-state visual evoked potentials . Nat Neurosci 10 : 117 - 125 . 47. Fries P , et al. ( 2001 ) Modulation of oscillatory neuronal synchronization by selective visual attention . Science 291 : 1560 - 1563 . 48. Buia C , Tiesinga P ( 2006 ) Attentional modulation of firing rate and synchrony in a model cortical network . J Comput Neurosci 20 : 247 - 264 . 49. Tiesinga PHE , Fellous JM , Salinas E , Jose JV, Sejnowski TJ ( 2004 ) Synchronization as a mechanism for attentional gain modulation . Neurocomputing 58 - 60 : 641 - 646 . 50. Buehlmann A , Deco G ( 2008 ) The neuronal basis of attention: rate versus synchronization modulation . J Neurosci 28 : 7679 - 7686 . 51. Beaulieu C , Kisvarday Z , Somogyi P , Cynader M , Cowey A ( 1992 ) Quantitative distribution of GABA-immunopositive and -immunonegative neurons and synapses in the monkey striate cortex (area 17) . Cereb Cortex 2 : 295 - 309 . 52. Koch C ( 1999 ) Biophysics of Computation . Oxford, UK: Oxford University Press. pp 146 - 151 . 53. Kandel ER , Schwartz JH , Jessell TM ( 2000 ) Principles of Neural Science. New York : McGraw-Hill . pp 125 - 139 . 54. Hodgkin A , Huxley A ( 1952 ) A quantitative description of membrane current and its application to conduction and excitation in nerve . J Physiol (London) 117 : 500 - 544 . 55. Bower J , Beeman D ( 2003 ) The Book of Genesis , Internet Edition . pp 44 - 47 .

This is a preview of a remote PDF: http://www.ploscompbiol.org/article/fetchObject.action?uri=info%3Adoi%2F10.1371/journal.pcbi.1000342&representation=PDF

Se-Bum Paik, Tribhawan Kumar, Donald A. Glaser. Spontaneous Local Gamma Oscillation Selectively Enhances Neural Network Responsiveness, PLOS Computational Biology, 2009, DOI: 10.1371/journal.pcbi.1000342