Abstract
Experiments in visual cortex have shown that the firing rate of a neuron in response to the simultaneous presentation of a preferred and non-preferred stimulus within the receptive field is intermediate between that for the two stimuli alone (stimulus competition). Attention directed to one of the stimuli drives the response towards the response induced by the attended stimulus alone (selective attention). This study shows that a simple feedforward model with fixed synaptic conductance values can reproduce these two phenomena using synchronization in the gamma-frequency range to increase the effective synaptic gain for the responses to the attended stimulus. The performance of the model is robust to changes in the parameter values. The model predicts that the phase locking between presynaptic input and output spikes increases with attention.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Our retinas are constantly stimulated by an overwhelming amount of information and the brain faces the task of reducing a potentially overloading amount of information into a manageable flow that reflects both the current needs of the organism and the external demands placed on it. In order to solve this problem, the brain uses a strategy to select the relevant information and to suppress information which is not relevant. The focus on and selection of relevant information is referred to as “attention”. If just one single stimulus falls within the receptive field of a neuron, this stimulus can be attended or not, and in the latter case a stimulus outside the receptive field may be attended. Since higher cortical areas have large receptive fields (Smith et al. 2002), it is quite common that two (or even more) stimuli fall within the receptive field of a neuron. In that case one of them can be attended (selective attention) or none of them. In order to understand the neuronal substrate of attention, many single-unit studies in visual cortex have investigated how attended and unattended stimuli are encoded in the firing rate of neurons.
Neural correlates of attention have been studied using single-unit recordings in areas V1, V2, V4 and V5/MT in primate visual cortex. Several studies have shown that attention increases a neuron’s firing rate in response to a single stimulus in its receptive field (Treue and Maunsell 1999; Luck et al. 1997; Reynolds et al. 1999; McAdam and Maunsell 1999; Fries et al. 2001). When two stimuli are presented in the receptive field of the neuron, the firing rate lies between the firing rates elicited by each of the stimuli presented alone (Moran and Desimone 1985; Treue and Maunsell 1996, 1999; Luck et al. 1997; Chelazzi et al. 1998, 2001; Reynolds et al. 1999; Reynolds and Desimone 2003). This phenomenon is called stimulus competition, since populations of input neurons, encoding different stimuli, are thought to compete with one another to generate neuronal responses intermediate between the responses to the individual stimuli. When attention is directed to the neuron’s preferred stimulus, the neuron’s firing rate increases, whereas attention to the non-preferred stimulus decreases the firing rate (Chelazzi et al. 1998; Reynolds et al. 1999).
Several models have been proposed to reproduce these experimental observations regarding stimulus competition and selective attention. Reynolds et al. (1999) could explain their experimental results by assuming that the synaptic weights of an input representing one of the two stimuli increase five-fold when attention is directed towards that stimulus. However, it is not clear how synaptic efficacies could change five-fold at the time scale of attentional shifts.
Most approaches to come up with an explanation for stimulus competition and selective attention have focused on the effects of attention on the firing rate of neurons (see e.g. Tiesinga 2005; Deco and Rolls 2005; Buia and Tiesinga 2006; Mishra et al. 2006). In addition to firing rate, several studies have provided convincing evidence that selective attention also increases rhythmic synchronization among selected neuronal signals (Kreiter and Singer 1996; Fries et al. 2001; Schoffelen et al. 2005; Taylor et al. 2005; Womelsdorf et al. 2006). Several groups have published a model for neural implementation of attentional processes that attributes a possible role to the neuronal oscillatory activity in stimulus competition and/or selective attention (Tiesinga 2005; Buia and Tiesinga 2006; Mishra et al. 2006). Mishra et al. (2006) used gamma range correlations in the feedforward inhibitory inputs to the V4 neuron which are out of phase with the gamma band correlations within the excitatory input corresponding to the attended stimulus. Tiesinga (2005) used two asynchronous excitatory input populations and two stimulus-driven inhibitory input populations, which send 40 Hz spike volleys with some temporal dispersion to a V4 model neuron. In that study attention is modelled by changing the temporal dispersion or the relative phase between the volleys coming from the two inhibitory populations. Tiesinga (2005) used the crosscorrelation function as a measure for the synchronization between the responses of two V4 neurons. Since he did this only for the condition that two stimuli are presented in the same receptive field, it is difficult to compare the result with the experimental results of one stimulus within and one outside the receptive field of a neuron as measured by Fries et al. (2001). Another measure for the synchronization between two signals is the coherence function. We will use the coherence function as a measure in the frequency domain for the synchronization between the input and output of the excitatory neuron in our model for different conditions.
Since it is well known that the excitatory input in visual cortex from V1 to V2 and from V2 to V4 contains gamma frequency oscillations (Eckhorn et al. 1993; Frien et al. 1994; Maldonado et al. 2000), we have explored the possible role of gamma frequency oscillatory input in stimulus competition and selective attention. We tried to reproduce the experimental observations by a simple feedforward model. Therefore, the aim of the present study was to explore whether a simple feedforward model could explain the phenomenon of stimulus competition with a role for synchronous modulation of stimulus-related activity to implement the attentional bias. Our results show that a feedforward model, very similar to the gain modulation model of Reynolds et al. (1999) but with fixed synaptic weights, can explain stimulus competition. Assuming that attention is implemented by increased synchronization of multi-unit spike activity, the model can reproduce the results by Chelazzi et al. (1998) and Reynolds et al. (1999) on stimulus competition and selective attention. Although this model has a feedforward architecture, the underlying mechanism for changes in attention-related modulations of synchronous activity is not specified, this requires a role for some top–down feedback mechanism capable of enhancing synchrony.
2 Methods and theory
We will start this section with a description of our model and the input signals to the model. In the second part of this section we will describe the methods to calculate the coherence, the phase coherence and the phase locking value between synaptic input and spike output.
2.1 Model
Figure 1 shows the feedforward network, that we propose to explain stimulus-competition and selective attention. The output neuron Y receives excitatory spike-trains from two populations (X1 and X2) with 80 Poisson neurons each and also receives inhibitory input from a population of 40 inhibitory neurons, for brevity called interneurons, I. In this study X1 and X2 represent the population of neurons encoding the preferred and non-preferred stimulus, respectively. With two populations of 80 excitatory neurons and a population of 40 inhibitory neurons projecting to the output neuron, the ratio of excitatory versus inhibitory synapses is 80 vs. 20% in agreement with experimental observations (Beaulieu et al. 1992). The two excitatory populations of neurons also project to the interneuron population. There is a small time delay τ d of 2 ms between the spike times of the interneurons and the arrival times of these spikes at neuron Y. The interneurons and the output neuron Y have been implemented in NEURON, as Hodgkin–Huxley type neurons (see below).
2.2 Stimulus-related input signals
The outputs from X1 and X2 are Poisson trains of spikes with a time-dependent rate r i(t):
with i ε {1,2}, r the constant rate, η i bandpass filtered Gaussian white noise with 3dB points at 45 and 55 Hz, a quality factor Q of 5, zero mean and a variance of one, and with \(A^i_m\) the modulation amplitude of the Gaussian white noise (GWN) for population i. When the modulation amplitudes \(A_m^i\) are the same for the non-preferred and the preferred stimulus, the spike trains encoding the non-preferred and the preferred stimulus are statistically identical. The different responses of the output neuron to the two stimulus inputs are due to the differences in synaptic conductances of the projections of the two populations of Poisson neurons to the inhibitory neurons and output neuron (will be explained later). Since we are not aware of any hard physiological data about these synaptic conductances in the literature, the different projections of the preferred and non-preferred stimulus to the interneurons and to the output neuron are an assumption of the model.
Several studies have shown that attention to a visual stimulus results in increased coherence between the local field potential and the activity of neurons, especially in the γ-band range (Fries et al. 2001; Womelsdorf et al. 2006). In the visual system γ-band oscillations have been reported at frequencies in the range 40–80 Hz. Based on these findings we postulate that selective attention to a sensory stimulus is implemented as an increased amplitude A m for the neuronal activity encoding that stimulus. For the simulations of the responses of the output neuron Y to various input signals we used a time duration T of 8 092 s and time step dt of 0.1 ms. The spike trains of the two Poisson populations X 1 and X 2 were modulated by a constant mean rate r = 20 and with a modulation amplitude A m = 6 for a non-attended stimulus and A m = 8 for an attended stimulus (see Eq. (1)). If no input is presented to X 1 or X 2, r = 3 and A m = 0. In order to explore the role of the modulation amplitude on the results of this study, some simulations used a modulation amplitude of 12 and 16 for the unattended and attended stimulus, respectively.
2.3 Geometry and properties of the HH-type interneurons and output neuron
The interneurons and output neuron Y were implemented in the NEURON simulation environment (Hines and Carnevale 1997) as single-compartment Hodgkin–Huxley type neurons with an area of 34,636 \(\upmu\)m2, in agreement with Destexhe et al. (2001). The inhibitory interneurons contain two sets of 80 synapses, the output neuron Y has 40 inhibitory and two sets of 80 excitatory synapses. The synaptic conductivity g is modelled by the default alpha function in NEURON. In this study most results were obtained for modulation amplitude A m values of 6 or 8. In that case the excitatory synapses from the populations X 1 and X 2 onto the interneurons have a maximum conductance of \(g^{\rm int}_{\rm np}\) = 0.84 nS and \(g^{\rm int}_{\rm p}\) = 0.55 nS for the non-preferred and preferred stimulus input, respectively and a time constant τ e = 2 ms. For the excitatory synapses onto the output neuron Y, the following values are taken: \(g^{Y}_{\rm np}\) = 1.52 nS, \(g^{Y}_{\rm p}\) = 1.71 nS, τ e = 2 ms. For the synapses from the inhibiting interneurons to the output neuron Y we had g inh = 4.5 nS and τ i = 5 ms. For modulation amplitudes A m with values of 12 (‘no attention’) and 16 (‘with attention’) the synaptic conductance values were \(g^{\rm int}_{\rm np}\) = 0.84 nS, \(g^{\rm int}_{\rm p}\) = 0.55 nS, \(g^{Y}_{\rm np}\) = 1.52 nS, \(g^{Y}_{\rm p}\) = 1.71 nS and g inh = 3.8 nS. With these values for A m = 12 the output neuron in our model generates, in agreement with experimental data of Reynolds et al. (1999), a firing rate (f p) of about 20 sp/s in response to the ‘preferred’ stimulus condition and a firing rate (f np) of about 10 sp/s in response to the ‘non-preferred’ stimulus.
The somata of the Hodgkin–Huxley type neurons have passive and active cell properties. The passive properties are the leak reversal potential (− 80 mV), leak conductance (4.52 10 − 5 S/cm2) and membrane capacitance (1 \(\upmu\)F/cm2). The active properties refer to the voltage-dependent Na + current and the “delayed-rectifier” K + current. The parameter values for the voltage-dependent Na + and K + currents were as described by Traub and Miles (1991) (see Appendix 1).
The synaptic background activity of the Hodgkin–Huxley-like neurons (interneurons and output neuron) was approximated by conductance injections in the soma as described in Destexhe et al. (2001) (see Appendix 1). In agreement with Destexhe et al. (2001), we used the following parameter values for the output neuron: the reversal potentials of the excitatory and inhibitory inputs E e = 0 mV, E i = − 75 mV, the average conductances g eo = 12.1 nS, g i0 = 57.3 nS and the time constants τ e = 2.73 ms, τ i = 10.49 ms. The standard deviations of the conductances corresponding to the background activity of output neuron Y are given by σ e = 3.0 nS and σ i = 6.0 nS. For the interneurons the average conductances and the standard deviations of these conductances are 50% of the corresponding values of the output neuron.
In order to understand the responses of the interneurons, it is helpful to appreciate the relative size of the synaptic currents due to the background noise and due to stimulus related inputs. These synaptic currents due to the spike input are rough estimates, since the precise relation between spike input and synaptic current depends on the membrane potential of the neuron, and thereby also depends on other synaptic inputs that affect the membrane potential. Assuming that the mean membrane potential is near − 55 mV (i.e. halfway between the membrane potential at rest near − 75 mV and the threshold for firing) the mean current due to background activity for the interneurons is about 60% of the total excitatory input current. The remaining 40% comes from the mean excitatory input related to the preferred stimulus (16%) and to the non-preferred stimulus (24%). For the output neuron, the inhibitory stimulus related input is about 20% of the current due to the background activity, whereas the excitatory stimulus related input is about 85% of the background current. More details on these relative contributions and their effect on the relation between mean input current and firing rate is provided in Appendix 2.
2.4 Coherence estimate
One of the predictions that flows from our hypothesis (see Section 1) is that the output spike train is more coherent to the “attended” input spike train than to the “ignored” input spike train. To quantify this, we will use the coherence function, in addition to firing rate to investigate the effect of attention on the spike output of neuron Y. In order to distinguish between the effect of the non-preferred and the preferred stimulus on the spike output, the non-preferred and preferred stimulus are statistically uncorrelated (\(<\eta_i^{(t)}\eta_j^{(t)} > = \delta_{ij}\)). This is in agreement with Gray et al. (1989) and Kreiter and Singer (1996) who reported that correlations between neuron population activities encoding different stimuli are absent.
The coherence function γ(ω) reflects how much of the variations in the output y can be attributed to a linear filtering of the input signal x. The coherence function γ(ω) is defined by:
with C xy (ω) the Fourier transform of the cross covariance function (Marmarelis and Marmarelis 1978). The coherence takes values in the range between 0 (input and output are uncorrelated) and 1 (the output is equal to the input after convolution by a linear system). Since the neuron itself is not a linear system, the coherence between the bandpass filtered Gaussian white noise input of one of the two Poisson populations and the spike output of neuron Y will not reach the upper limit of one.
To estimate the coherence and its variance, we used the multi-taper method (Thomson 1982; Mitra and Pesaran 1999). The key idea behind the multi-taper method is that a physiological signal does not have discontinuities in the frequency spectrum and that the variance in the estimate of a signal can be reduced by smoothing in the frequency domain. The multi-taper method minimizes bias and variance of the estimate by using multiple orthonormal data tapers. We have used sine-tapers as described in Zeitler et al. (2006) with length N = 1.024 s and bandwidth W = 2.9 Hz. Since the number of tapers to be used is K = 2NW − 1 tapers, the values for N and W used in this study gave K = 5. The binwidth in the frequency domain is the Rayleigh frequency f r = 1/T = 1/(nfft/f s ) = 0.98 Hz, with sampling frequency f s (1,000 Hz) and where nfft (1,024) is the number of data points in the FFT . The input and output signals were both segmented in T/N non-overlapping time segments of 1,024 ms, with T the duration of the simulation.
2.5 Phase locking
A high value of the regular coherence Eq. (2) implies a strong relation between both amplitude and phase of input and output. Previous studies have shown that pairs of neuronal responses can undergo variations in relative amplitude even in the presence of tight phase coupling (Tass et al. 1998; Lachaux et al. 1999). For this reason, the phase coherence has been introduced, which only considers the variability in relative phase between two signals s1 and s2. In this study, the phase coherence is calculated by segmentation of the two signals s1(t) and s2(t), both segmented in T/N non-overlapping time segments of 1,024 ms. Each segment of the signal s1 and the corresponding segment of the second signal s2 form a pair. The phase difference Δϕ(f) at frequency f for each pair is given by:
where * refers to complex conjugate.
Figure 2(a) shows a typical polar plot of the phase differences between 150 stimulus–response pairs for the neuron model in Fig. 1. The full range of 360° was subdivided into 24 bins of 15° [15(j − 1), 15j] for j ε {1,..,24}. The number of phase differences falling into a bin, divided by the total number of phase differences in the unit circle, is the fraction of stimulus–response pairs with a phase difference in that bin. For each of the twenty bins, this fraction is represented by the length of the arrow, drawn in the middle of each bin [see Fig. 2(b)]. All fractions are connected by a line. When stimulus and response are not phase locked at all, the phase differences will be distributed uniformly over 360°. Complete phase locking with phase difference Φ corresponds to an arrow of unit length pointing in the direction Φ.
Lachaux et al. (1999) introduced a method to quantify the degree of phase-locking between two signals. The phase locking value (PLV) between the two periodically repeated signals measures the inter-trial variability of the phase difference between these two signals. In our simulations, we average the phase relation over all stimulus–response pairs of the M = T/N time segments:
This phase locking value measures the average variability of the phase difference and takes values between 0 (complete lack of phase-locking) and 1 (completely phased locked).
The phase locking value is a function of frequency. In order to reduce the variance of the phase locking value, we used the multi-taper method, with K = 5 sine-tapers to reduce the variance of the spectra S1(f) and S2(f) in Eq. (3). Since the phase locking value was very similar for all frequencies near 50 Hz, we determined the phase locking value for f = 50 Hz as this gave the best signal-to-noise ratio for the 50 Hz bandpass filtered Gaussian white noise input.
3 Results
In this section the simulation results will be described for the firing rate of the output neuron Y (Section 3.2) and the coherences between the spikes of the output neuron and each of the stimulus-related inputs to the populations of Poisson neurons (Section 3.3). We will conclude this section with the phase locking results (Section 3.4).
3.1 Input–output relation of an interneuron
The interneuron plays a crucial role to explain stimulus competition (Fig. 3). If an interneuron receives input of the preferred stimulus with firing rate f in (this implies that all 80 neurons encoding the preferred stimulus have a constant firing rate of f in and the 80 neurons encoding the (absent) non-preferred stimulus have a firing rate of 3 spikes/s) the interneuron starts to respond at relatively high input firing rates (dashed line). Since the synaptic projections to the interneuron of the neural activity encoding the non-preferred stimulus is stronger than for the neural activity related to the preferred stimulus, the relation shifts to the left for the non-preferred stimulus only (solid line, lower threshold for firing). The inset shows the population activity of the Gaussian white noise (GWN) modulated Poisson spike input. This explains why the mean firing rate of the interneuron increases with increasing amplitude of the GWN-modulated Poisson spike input. For the output neuron qualitatively simular relations hold for the excitatory stimulus-related input, except for the fact that the output neuron responds better to the preferred stimulus alone, than to the non-preferred stimulus alone (solid and dashed lines interchanged).
3.2 Simulation results for the firing rate
Figure 4 shows the firing rate of the output neuron for various stimulus conditions. The upper panel (a) shows the results for small modulation amplitudes (A m = 6 and 8), the lower panel (b) for larger modulation amplitudes (A m = 12 and 16). The results at the left part of the figure show the results for the ‘no attention’ condition, the right part of the figure the results for the stimulus conditions with one stimulus attended. We will first discuss Fig. 4(a).
As explained in Section 2, the statistical properties of the spike series, representing the non-preferred and preferred stimulus with no attention, are identical. The different effectiveness of the non-preferred and preferred stimulus is mainly due to the different conductance of the excitatory synapses from X 1 and X 2 to the output neuron. Since \(g_{\rm p}^{Y}\) (1.71 nS) > \(g_{\rm np}^{Y}\) (1.52 nS) the direct excitatory projections of the population representing the preferred stimulus to the output neuron induce more action potentials in the output neuron than that of the population of neurons representing the non-preferred stimulus.
The population activities representing the preferred (X 1) and non-preferred stimulus (X 2) also reach the output neuron via the interneurons. In case only one stimulus is offered, the interneurons have a low firing rate. This is shown in Fig. 5. For each stimulus separately, the induced firing rate of the interneurons is increasing as a function of the modulation amplitude A m . However, the firing rates in response to the preferred and non-preferred stimulus are rather small (range between 0 and 0.03 Sp/s and between 0 and 1.99 Sp/s for the preferred and non-preferred stimulus, respectively). Therefore, it is mainly the larger conductance of the excitatory synapses from population X 1 to Y which explains the higher firing rate of the output neuron to the activity of population X 1 (f p = 14.15 sp/s, SD = 0.05 sp/s) than to the population activity X 2 (f np = 8.09 sp/s, SD = 0.04 sp/s), see left side of Fig. 4(a).
If the two neuronal populations, representing the activity of the preferred and non-preferred stimulus, would project to the output neuron only via excitatory synapses, one would expect a summation of firing rates when the preferred and non-preferred stimulus are presented simultaneously. However, we find stimulus competition in the responses of the output neuron, which is in agreement with experimental single-unit recordings (Reynolds et al. 1999), which most likely reflect the activity of excitatory neurons by their greater number and larger extracellular spikes. The interneurons play a crucial role in stimulus competition. This is illustrated in Fig. 5, which shows the response of an interneuron for the three conditions: ‘non-preferred stimulus only’, ‘preferred stimulus only’ and ‘preferred and non-preferred stimulus simultaneously’ as a function of the modulation amplitudes A m of each stimulus. For all three conditions the response increases as a function of increasing A m which shows that the interneurons are sensitive to correlated input. As mentioned before, the firing rate of the interneuron is very small when only one stimulus is presented. Therefore, the inhibition is small. When two stimuli are presented simultaneously the firing rate of the interneurons increases more than linearly due to the sigmoidal relation between synaptic input and firing rate of neurons in general. For the preferred and non-preferred stimulus alone the interneurons operate at the bottom of the sigmoidal relation, whereas the combined input of the preferred and non-preferred input shifts the firing rate to the steep phase of the sigmoidal relation, see also Fig. 11 in Appendix 2. So stimulus competition is caused by the activity of inhibitory interneurons, which generate a much higher firing rate when two stimuli are presented simultaneously compared to the condition that only one stimulus is presented. This higher response causes more inhibition for the target neuron and thus explains why the firing rate of the output neuron to both stimuli (f both = 13.72 sp/s, SD = 0.05 sp/s) falls between the firing rates to the preferred and the non-preferred stimulus presented alone.
In summary, the responses to the preferred and non-preferred stimulus alone are mainly due to excitatory inputs and the difference in firing rates (f np < f p) is caused by the different synaptic conductances. Competition (f np < f both < f p) is the net effect of the two direct excitatory inputs plus the inhibition via the interneurons, which are mainly actively if both stimuli are offered simultaneously.
Based on experimental observations that have revealed larger amplitudes of γ-range activity during attention (Fries et al. 2001; Womelsdorf et al. 2006; Taylor et al. 2005), attention to the preferred or non-preferred stimulus is implemented by a larger amplitude A m of the band-pass filtered noise to the Poisson neurons. Increasing A m leads to more spikes in the bursts of the population activity. Since the interneurons and output neuron receive a background synaptic input, they are sensitive to synchronous input (Martinez 2006; Higley and Contreras 2005). This explains the higher firing rate of the output neuron to the preferred (\(f_{\rm p}^{\rm att}=\) 15.88 sp/s, SD = 0.05 sp/s) and non-preferred (\(f_{\rm np}^{\rm att}=\) 9.18 sp/s, SD = 0.04 sp/s) stimulus with attention, relative to the ‘no attention’ condition [see right sight of Fig. 4(a), which shows the responses to the attended stimuli].
A larger modulation depth causes larger excitatory spike volleys in the populations of Poisson neurons and results in higher firing rates of the 41 HH-like neurons. Since the larger modulation depth impacts also the interneurons, this increased modulation can increase or even decrease the firing rate of the output neuron depending on the net balance between excitatory and inhibitory input. When both stimuli are presented simultaneously and when the preferred stimulus is attended, the effect of the larger excitatory spike volleys encoding the attended preferred stimulus is larger than the effect of inhibition by the increased firing rate of the interneuron. Therefore, the resulting firing rate \(f_{\rm both}^{{\kern1pt}\rm att.pref}\) = 14.14 sp/s, SD = 0.05 sp/s [second line from top at the right side of Fig. 4(a)] is slightly larger than that in the condition of ‘no attention, both stimuli’ (f both = 13.72 sp/s, SD = 0.05 sp/s). If the non-preferred stimulus is attended instead of the preferred stimulus, the effect of larger excitatory spike volleys is smaller than the effect of inhibition by the increased firing rate of the interneurons. Therefore, the resulting firing rate \(f_{\rm both}^{{\kern1pt}\rm att.np} =\) 12.73 sp/s, SD = 0.05 sp/s [third line from top at the right side of Fig. 4(a)] is significantly lower than in the condition ‘no attention, both stimuli’.
In summary, the attended stimulus, presented alone, gives higher firing rates than the unattended stimulus alone due to the increased number of spikes in the population volleys. If both stimuli are presented and one is attended, the firing rate of the output neuron changes towards the firing rate elicited by that stimulus alone \(\left( f_{\rm both}^{\rm att.np}<f_{\rm both}< f_{\rm both}^{\rm att.pref}\right)\).
Figure 4(b) shows that stimulus competition and the stimulus selection effect can also occur for other values of the modulation amplitude A m .
Obviously, the performance of the model depends on the strength of the excitatory projections of the non-preferred and preferred stimulus \((g_{\rm np}^{Y}\) and \(g_{\rm p}^{Y})\), on their projections to the interneurons \((g_{\rm np}^{\rm int}\) and \(g_{\rm p}^{\rm int})\), and on the synaptic connection g inh of the interneurons to the output neuron. The results presented in Fig. 4(a) were obtained with a fixed set of parameter values. The results presented in Fig. 4(b) [with a modulation amplitude twice as large as in Fig. 4(a)] were obtained with the same parameter values except for g inh which was decreased to 3.8 nS. The increase in modulation amplitude gives rise to an increased excitatory drive to both the output neuron Y and the inhibitory neuron. Since the output neuron is inhibited by the interneurons, the change in firing rate of the output neuron Y related to the increased modulation amplitude depends on the relative amounts of background noise, excitatory input and the strength of inhibition by the interneurons (Table 1). The new value for g inh of 3.8 nS brings the firing rate of the output neuron to both stimuli halfway between that for the preferred and non-preferred stimulus only. Without reduction of g inh the firing rate to both stimuli would have been strongly biased towards the firing rate for the non-preferred stimulus only (Result not shown).
In order to investigate to what extent the results in Fig. 4(b) depend on the particular choice of synaptic conductances, we have analyzed the model for a range of values of the relevant five synaptic conductances, \(g_{\rm np}^{\rm int}\), \(g_{\rm p}^{\rm int}\), \(g_{\rm np}^{Y}\), \(g_{\rm p}^{Y}\) and g inh. As it is difficult to visualize a five dimensional parameter space, we have varied the synaptic conductances of the non-preferred and preferred stimulus to the interneurons (\(g_{\rm np}^{\rm int}\) and \(g_{\rm p}^{\rm int}\)), and tried to find the proper values for \(g_{\rm np}^{Y}\), \(g_{\rm p}^{Y}\) and g inh such, that the model reproduced the properties of stimulus competition and selective attention. In detail, we adjusted the values of \(g_{\rm np}^Y\), \(g_{\rm p}^Y\) and g inh for each pair of (\(g_{\rm np}^{\rm int}\), \(g_{\rm p}^{\rm int}\)) values such that the model had the following properties:
-
1)
the firing rate to the preferred and non-preferred stimulus alone should be in the range between 19.6 to 20.1 sp/s and 9.7 to 10.4 sp/s, respectively. As explained before this is implemented by the requirement that \(g_{\rm p}^Y>g_{\rm np}^Y\);
-
2)
the firing rate of the output neuron to the non-preferred and preferred stimulus presented together should be between the firing rates of the non-preferred and the preferred stimulus presented alone (stimulus competition);
-
3)
attention should give higher firing rates than without attention, when the non-preferred or preferred stimulus is presented alone;
-
4)
attention to either the non-preferred or preferred stimulus, presented simultaneously, changes the firing rate towards that for the attended non-preferred or preferred stimulus presented alone.
We found \(1.45~nS \leq g_{\rm np}^{Y} \leq 1.52~nS\), \(1.70~nS \leq g_{\rm p}^{Y} \leq 1.74~nS\) and 3.3 nS ≤ g inh ≤ 5.04 nS for the three values of the synaptic connections to the output neuron Y, which are not shown in Fig. 6(a).
The fitted ellipse in Fig. 6(a) shows the range of parameter values for \(g_{\rm np}^{\rm int}\) and \(g_{\rm p}^{\rm int}\) where the effects of competition and selective attention can be reproduced for A m = 12 (no attention) and A m = 16 (with attention). For the region with parameters left of the grey area, either the inhibition is too small to reproduce the effect of stimulus competition or the inhibition is too strong, such that attention to the preferred stimulus does not increase but decrease the firing rate of the output neuron Y. For the region with parameters at the lower right of the grey area the model fails on a third aspect: attention to the non-preferred stimulus only decreases rather than increases the firing rate of the output neuron. Outside the upper boundary either the condition \(f_{\rm both}^{\rm att.np} < f_{\rm both}\) or f np < f both is violated.
Figure 6(a) shows that stimulus competition and selective attention occur for different values of the pair (\(g_{\rm np}^{\rm int}, g_{\rm p}^{\rm int}\)). Variations in the parameter values \(g_{\rm np}^{\rm int}\) and \(g_{\rm p}^{\rm int}\) cause changes in firing rates of the interneurons. The fitted line in Fig. 6(b) shows the firing rate f Y of the output neuron in the ‘no attention, both stimuli’ condition as a function of the firing rate f I of the interneurons. The firing rate of the output neuron, when both the preferred and the non-preferred stimulus are presented falls between the firing rates of the preferred stimulus alone (20 sp/s) and of the non-preferred stimulus alone (10 sp/s). This firing rate is high (small) for low (high) firing rates of the inhibitory neurons.
In summary, our results show that i) the competition and attention effects as shown in Fig. 4 occur for a range of synaptic conductance values; ii) the firing rate in the ‘no attention, both stimuli’ condition takes values between f np and f p.
3.3 Simulation results for coherence estimate
Figure 7 shows the coherence between the response of the output neuron and the time-dependent rate to populations X 1 and X 2 when either the non-preferred (upper row) or the preferred (lower row) stimulus is presented. The left and right column show the results for the ‘no attention’ (A m =6) and ‘with attention’ (A m =8) condition, respectively. Each of the panels shows a peak at 50 Hz, corresponding to the frequency content of the band-pass filtered stimuli.
Both for the non-preferred and preferred stimulus, the peak value of the coherence is larger for the ‘with attention’ condition (0.50 and 0.60, respectively) than for the ‘no attention’ condition (0.40 and 0.50, respectively). The 95% confidence level corresponds roughly to the range of the mean value, plus or minus 0.04. The larger coherence for the ‘with attention’ condition relative to the ‘no attention’ condition is due to the fact that larger spike volleys in the input will cause more precise spike timing (less variability). The peak values of the coherence for the non-preferred stimulus (Fig. 7(a, b)) are smaller than those for the preferred stimulus (Fig. 7(c, d)). This is caused by two facts: the preferred stimulus has stronger excitatory synapses to the output neuron than the non-preferred stimulus (\(g_{\rm p}^Y>g_{\rm np}^Y\)) and will therefore cause spikes which are more precisely time-locked to the stimulus. The second reason is that the stronger synaptic projections of the non-preferred stimulus to the interneurons cause more frequent inhibitory post-synaptic potentials in the output neuron, which can delay or even prevent the non-preferred stimulus to elicit a spike in the output neuron, resulting in a smaller coherence peak value.
Figure 8 shows the coherence between the response of the output neuron and the input to population X 2 (non-preferred stimulus) (upper row) and to population X 1 (preferred stimulus) (lower row), respectively, when both stimuli are presented. The middle column [Fig. 8(b, e)] shows the results when both stimuli are presented simultaneously without attention (A m = 6). For the non-preferred and preferred stimulus the coherence estimate has a peak value of 0.35 and 0.40, respectively. These two peak values are smaller than for the condition when these stimuli were presented alone [Fig. 7(a, c)]. When the non-preferred and preferred stimulus are presented simultaneously, the spikes of the output neuron reflect the contribution of both stimuli. The effect of the non-preferred (preferred) stimulus on the spike responses acts as a noise term in the response to the preferred (non-preferred) stimulus, which explains the smaller coherence values in Fig. 8 compared to that in Fig. 7.
When the preferred or non-preferred stimulus is attended [Fig. 8(a, f)], this stimulus becomes more effective, causing a better locking of the spike response to that stimulus. This more precise locking of the neuron to the attended stimulus leads to a larger coherence value for the attended stimulus and a lower coherence for the non-attended stimulus [compare Fig. 8(a, f) with Fig. 8(c, d), respectively]. The larger coherence for the attended preferred stimulus (0.49 vs. 0.40) and for the attended non-preferred stimulus (0.43 vs. 0.35) is significant (p < 0.001). The tendency that the coherence for the non-attended stimulus decreases when the other stimulus is attended (0.37 vs. 0.40 for the preferred and 0.33 vs. 0.35 for the non-preferred stimulus) is significant (p < 0.001). The 95% confidence level of the values corresponds roughly to the range of the mean value, plus or minus 0.04.
The coherence results for modulation amplitudes, which are twice as large, are similar and therefore not shown.
In summary: by attending a stimulus, the peak value of the coherence between the attended input and the response is larger compared to the condition ‘no attention’. The coherence between the non-attended input and the response does significantly decrease compared to the ‘both stimuli’ condition.
3.4 Phase locking results
Figure 9 shows polar plots of the probability distributions of phase differences between stimulus and response. The solid line shows the results for the ‘no attention’, the dashed line for the ‘with attention’ condition. For the ‘non-preferred stimulus only’ condition Fig. 9(a) shows that there is clear phase locking between the stimulus and the response of the output neuron which increases with attention (dashed line). The narrower the ellipse, the better the signals are locked to a certain phase difference and the higher the phase locking value (PLV) will be. The increase of the PLV for the ‘attention’ condition is significant (0.69 ± 0.01 (‘non-preferred stimulus only’) versus 0.80 ± 0.01 (‘attended non-preferred stimulus only’), p < 0.001), where phase locking values are given as the mean plus or minus the standard deviation. Figure 9(b) shows similar results for the ‘preferred stimulus only’ condition (mean PLV 0.80 ± 0.01 and 0.88 ± 0.01, p < 0.001, for the ‘non-attended preferred stimulus only’ and ‘attended preferred stimulus only’ condition, respectively).
Figure 10(a, b) show the polar distribution of the phase relation between input and spike output, when the preferred and non-preferred stimulus are presented simultaneously. Figure 10(a, b) shows the phase relation between the output and the input to X 2 (non-preferred) and X 1 (preferred), respectively. The solid line shows the results for the condition ‘no attention, both stimuli’, the dashed line for the condition ‘with attention’. Both panels show that attention increases the phase locking between input and response. The phase locking values for the non-preferred (preferred) stimulus are significantly larger for the condition ‘with attention’ (0.73 ± 0.01, 0.79 ± 0.01, respectively, p < 0.001) than for the ‘no attention, both stimuli’ (0.63 ± 0.01, 0.70 ± 0.01, respectively, p < 0.001) condition. The mean PLV for the preferred (non-preferred) stimulus for the ‘no attention, both stimuli’ condition, PLV = 0.70 ± 0.01 (0.63 ± 0.01 ) is significantly different for the condition ‘with attention to the other stimulus’, PLV = 0.64 ± 0.01 (0.59 ± 0.01).
The phase locking value results for two times larger modulation amplitudes are similar and therefore not shown.
In our simulations, the results of the coherence function and phase locking values at 50 Hz are very similar: attention significantly increases the coherence and the phase locking value between the response and the attended stimulus.
4 Discussion
Many experimental and modelling studies have focussed on the neuronal implementation of attention (Bushnell et al. 1981; Spitzer et al. 1988; Motter 1993; McAdam and Maunsell 1999; Treue and Martinez-Trujillo 1999; Fries et al. 2001; Tiesinga 2005) and on stimulus competition (Moran and Desimone 1985; Desimone and Duncan 1995; Reynolds et al. 1999; Treue and Martinez-Trujillo 1999; Tiesinga 2005) at different levels of neuronal processing varying from brain areas (Corbetta and Shulman 2002) to single neurons (Deco and Rolls 2005; Tiesinga 2005). Most of these studies have focussed on firing rate to encode attended and unattended stimuli. However, it is well known that rhythmic neuronal activity, such as in β- and γ-oscillations, plays an important role in encoding sensory stimuli (see e.g. Kreiter and Singer 1996) and that attention affects the amplitude of the rhythmic neuronal oscillations. The latter is illustrated by the coherence between the local field potential and spike output, which provides a sensitive measure of local neuronal synchronization. Fries et al. (2001) found that for the ‘with attention’ condition, the coherence between the local field potential and the simultaneously recorded spike train was significantly larger with than without attention to the stimulus. Our results will be discussed in more detail below, starting with a comparison of the model responses with other models.
The architecture of our model is quite similar to the gain modulation model by Reynolds and coworkers (Reynolds et al. 1999). The main differences with respect to the gain modulation model are related to the nature of the neuronal input signals and to the neuronal implementation of attention. In the Reynolds model constant firing rates are used to encode the preferred and non-preferred stimuli and attention was implemented by a five-fold increase of the efficacy of the synapses that transmit the attended stimulus. This model left open the question of how synaptic efficacy can be modulated selectively for the attended stimulus input at such a short time scale. In agreement with experimental observations (Fries et al. 2001), our working hypothesis was that attention is implemented through enhanced gamma activity, which makes this input more effective in eliciting a spike in the output neuron, and thus increases the effective strength of the signal encoding the attended stimulus.
Our model is an alternative for the model proposed by Tiesinga (2005) which postulates a stimulus-related excitatory input without rhythmic oscillations and with top–down input from the frontal eye fields (FEF). The main difference between our model and that by Tiesinga (2005) is that we assume that attention is implemented in the γ-modulated stimulus-related neural input, whereas Tiesinga (2005) does not assume γ-modulated stimulus-related input. In the Tiesinga model the γ-oscillations are postulated to be induced by FEF input to the interneurons. Although it is well know that the FEF is involved in attention-related modulations of neuronal activity (Moore and Amstrong 2003), it is still a matter of debate how the FEF input affects the neuronal processing. Our model allows a role for top–down attentional modulation of the amplitude of the γ-oscillations in the input population activity representing the visual stimulus. It is a topic for future research to investigate the details of attention-related top–down mechanisms.
One of the values of models is that they can provide possible explanations for experimentally observed phenomena. When developing a model, one should always try to explain as many experimental findings with as few as possible model assumptions. In our model, we assumed that the stimulus-related neuronal activity has rhythmic oscillatory components. This assumption is supported by experimental observations which have revealed stimulus-related rhythmic activity in V1 (van der Togt et al. 2006; Roelfsema et al. 2004; Rols et al. 2001), V2 (Frien et al. 1994) and V4 (Fries et al. 2001; Taylor et al. 2005). Moreover, we assume that attention is implemented by increased amplitudes of the rhythmic excitatory activity. This is in agreement with experimental observations by Fries et al. (2001), Taylor et al. (2005), Womelsdorf et al. (2006), who reported that attention is related to an increased coherence between local field potentials and single-unit activity. The simple feed-forward model reproduces experimental data of stimulus competition and attention effects on firing rate (see e.g. Reynolds et al. (1999)). Moreover, our model predicts an increased peak value of the coherence due to attention, emphasizing the increased neuronal synchronization by attention. Our predictions concerning an increased coherence for attended stimuli and a decreased coherence in case the other stimulus within the receptive field is attended, are in agreement with what is found by Smiyukha et al. (2006). These authors placed two small stimuli close to each other, causing two spatially well separated foci of gamma-band activity in area V1 of a macaque. The corresponding foci in V4 were largely overlapping. Wavelet based analysis of correlations revealed strong synchronization of field potentials in the gamma-band between the site in V1, processing the attended shape, and the site in V4 responsive to both stimuli. Synchronization with activity in V4 is weak for other sites in V1, processing non-attended stimuli. This strong synchronization between the area in V1, which processes the attended stimulus, and the site in V4, is at least qualitatively similar to the increase in coherence between input and spike output in our model.
Recently, a model with an architecture very similar to our model was proposed by Mishra et al. (2006) to explain the phenomena of stimulus competition and selective attention. The neuronal mechanisms in their model to explain stimulus competition are feedforward inhibition, like in our model, and synaptic depression, which is effective for input frequencies of 40 Hz and above. Like in our model, each stimulus is represented by excitatory multi-unit activity. In their model the excitatory neuronal signals that encode the preferred and non-preferred stimuli are always in anti-phase. The phenomenon of selective attention in their model is achieved by imposing a phase shift of the response of the interneuron relative to the excitatory activity encoding the attended stimulus. This implies that the inhibition is more or less in anti-phase with the excitatory drive of the attended stimulus, but in phase with the excitatory drive of the unattended stimulus. Therefore, the excitatory input of the unattended stimulus is cancelled by inhibition from the interneuron. This works well when the excitatory drive and the inhibitory input from the interneurons is tuned at the same frequency (40 Hz in the paper by Mishra et al. (2006)) and more or less in anti-phase. However, experimental studies have found that rhythmic synchronization is broadly tuned and that the neuronal activity representing two different stimuli is uncorrelated (Gray et al. 1989; Kreiter and Singer 1996). Therefore, we decided to generate the neuronal signals for attended and unattended stimuli by band-pass filtering two independent noise signals.
Fries et al. (2001) showed that the amplitude of the input fluctuations of neurons in V4 is larger when the stimulus is attended than when the same stimulus is not attended. Therefore, we increased the amplitude of our input modulation by 33% to implement the effect of attention (from A m = 6 to 8). This increase in amplitude caused a 12% higher firing rate of the output neuron, a 20–25% increase in the coherence between input and output and a 10–16% larger phase locking value (PLV). We have also done the simulations for amplitudes of the input modulations, twice as large. Now, the 33% increase in amplitude of the input modulation, which implemented the effect of attention, caused a 15–25% higher firing rate of the output neuron, a 11% increase in the coherence between input and output and a 3% larger PLV. This indicates that all results for firing rate, coherence and phase coherence are qualitatively similar, independent of modulation amplitude, showing that our model is quite robust.
The results of this study were obtained for various modulation amplitudes of the stimulus-related input with the same set of parameters, except for the value of g inh which was reduced from 4.5 to 3.8 nS when the modulation amplitude was made twice as large. If we had kept the synaptic strength at 4.5 nS, the response of the output neuron to both stimuli would have been more biased towards the output for the non-preferred stimulus. As far as we know there have been no studies which have systematically investigated the effect of changes in modulation amplitude of excitatory drive to neurons in V2 and V4 on stimulus competition. Maybe a bias to the response to the non-preferred stimulus alone for larger modulation amplitudes is what will be observed. Another alternative might be that dynamic synapses (Tsodyks et al. 1998, 2000) reduce the effective synaptic strength of the projection of the inhibitory neurons to the output neuron when the increased modulation amplitude causes a larger increase of the firing rate of the inhibitory neurons. The latter seems a plausible mechanism which we saw as a justification to reduce the synaptic efficacy of the projections of the inhibitory neurons. The size of the reduction is certainly not critical to qualitatively reproduce the results in this study.
A similar robustness was found for variations in the synaptic strengths. As shown in Fig. 6(a) stimulus competition could be reproduced over a range from 0.7 to 1.0 nS for \(g_{\rm np}^{\rm int}\) and from 0.4 to 0.7 nS for \(g_{\rm p}^{\rm int}\). Changes in the parameter values lead to variations in the firing rate of the output neuron [Fig. 6(b)]. For the condition ‘no attention, both stimuli’ this range goes from about 15 sp/s to 18 sp/s depending on the firing rates of the inhibitory neurons. This range of attenuation of the firing rate of the output neuron to both stimuli compared to the firing rate for the preferred stimulus alone, is within the range reported by Reynolds et al. (1999) and Gawne and Martin (2002) for V4, by Miller et al. (1993) for the inferior temporal cortex and by Rolls and Tovee (1995) in the anterior part of the superior temporal sulcus.
We want to remark that the effect of competition in the experimental results in the literature is not always as large as in the paper shown by Reynolds et al. (1999). See for example the study of Gawne and Martin (2002). For a substantial fraction of the neurons these authors found that the firing rate to both stimuli was close to the highest firing rate to the stimuli presented separately.
Lachaux et al. (1999) showed that the coherence cannot distinguish phase and amplitude covariance. As an alternative they introduced the PLV to detect phase synchrony. For our simulation results we determined the coherence values (typically 0.33–0.60) as well as the phase locking values (typically 0.59–0.88). These two sets of values lead to the same qualitative conclusions: (1) the input and response of the output neuron Y are more synchronized if the stimulus, represented by the neuronal input, is attended compared to be not attended; (2) the input representing the not-attended stimulus is less synchronized with the response of the output neuron Y than the input which represents the simultaneously offered but attended stimulus. In addition, the polar plots of Figs. 9 and 10 show that the average phase difference between input and response of the output neuron Y are different for the stimulus conditions ‘non-preferred stimulus, only’, ‘preferred stimulus, only’ and ‘both stimuli’. The average phase difference for the ‘preferred stimulus, only’ condition [< ϕ p > ~ 125°, Fig. 9(b)] is smaller than for the ‘non-preferred stimulus, only’ condition [< ϕ np > ~ 133°, Fig. 9(a)]. This can be explained by the fact that the maximum conductance of the excitatory synapses to the output neuron is larger for the excitatory neuronal activity representing the preferred than for that representing the non-preferred stimulus. For the condition ‘both stimuli, no attention’ this average phase difference [< ϕ both > ~ 114°, solid line in Fig. 10(a, b)] is even smaller since the neuron receives more input, so that it can generate even faster a spike (in case it is not inhibited! The input-output relation for the output neuron without inhibition is comparable to that for the inhibitory neurons shown in Fig. 11).
In this study we have presented a feedforward model which can reproduce neuronal responses in visual cortex related to stimulus competition and selective attention effect, by: (1) using gamma-modulated population activities to represent the stimuli; (2) increasing the modulation depth of the population activity representing the attended stimulus; (3) determining the non-preferred and preferred stimulus response by using different values for each group of the various synaptical conductances of the interneuron and output neuron. Our numerically obtained firing-rate results are similar to experimental results reported by Reynolds et al. (1999), Miller et al. (1993) and Rolls and Tovee (1995).
References
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). Cerebral Cortex, 2(4), 295–309.
Buia, C., & Tiesinga, P. (2006). Attentional modulation of firing rate and synchrony in a model cortical network. Journal of Computational Neuroscience, 20(3), 247–264.
Bushnell, M. C., Goldberg, M. E., & Robinson, D. C. (1981). Behavioral enhancement of visual responses in monkey cerebral cortex. I modulation in posterior parietal cortex related to selective visual attention. Journal of Neurophysiology, 46(4), 755–772.
Chelazzi, L., Duncan, J., Miller, E. K., & Desimone, R. (1998). Responses of neurons in inferior temporal cortex during memory-guided visual search. Journal of Neurophysiology, 80, 2918–2940.
Chelazzi, L., Miller, E. K., Duncan, J., & Desimone, R. (2001). Responses of neurons in macaque area V4 during memory-guide visual search. Cerebral Cortex, 11, 761–772.
Corbetta, M., & Shulman, G. L. (2002). Control of goal-directed and stimulus-driven attention in the brain. Nature reviews, Neuroscience, 3, 201–215.
Deco, G., & Rolls, E. (2005). Neurodynamics of biased competition and cooperation for attention: A model with spiking neurons. Journal of Neurophysiology, 94, 295–313.
Desimone, R., & Duncan, J. (1995). Neural mechanisms of selective visual attention. Annual Review of Neuroscience, 18, 193–222.
Destexhe, A., & Paré, D. (1999) Impact of network activity on the integrative properties of neocortical pyramidal neurons in vivo. Journal of Neurophysiology, 81, 1531–1547.
Destexhe, A., Rudolph, M., Fellous, J.-M., & Sejnowski, T. J. (2001). Fluctuating synaptic conductances recreated in vivo-like activity in neocortical neurons. Neuroscience, 107(1), 13–24.
Eckhorn, R., Frien, A., Bauer, R., Woelbern, T., & Kehr, H. (1993). High frequency (60–90 Hz) oscillations in primary visual cortex of awake monkey. NeuroReport, 4, 243–246.
Frien, A., Eckhorn, R., Bauer, R., Woelbern, T., & Kehr, H. (1994). Stimulus specific fast oscillations at zero phase between areas V1 and V2 of awake monkey. NeuroReport, 5, 2273–2277.
Fries, P., Reynolds, J. H., Rorie, A. E., & Desimone, R. (2001). Modulation of oscillatory neuronal synchronization by selective visual attention. Science, 291, 1560–1563.
Gray, C. M., König, P., Engel, A. K., & Singer, W. (1989). Oscillatory responses in cat visual cortex exhibit intercolumnar synchronization which reflects global stimulus properties. Nature, 338, 334–337.
Gawne, T. J., & Martin, J. M. (2002). Responses of primate visual cortical V4 Neurons to simultaneously presented stimuli. Journal of Neurophysiology, 88, 1128–1135.
Higley, M. J., & Contreras, D. (2005). Balanced excitation and inhibition determine spike timing during frequency adaptation. Journal of Neuroscience, 26(2), 448–457.
Hines, M. L., & Carnevale, N. T. (1997). The NEURON simulation environment. Neural computation, 9, 1179–1209.
Kreiter, A. K., & Singer, W. (1996). Stimulus-dependent synchronization of neuronal responses in the visual cortex of the awake macaque monkey. Journal of Neuroscience, 16(7), 2381–2396.
Lachaux, J. P., Rodriguez, E., Martinerie, J., & Varela, F. (1999). Measuring phase synchrony in brain signals. Human Brain Mapping, 8, 194–208.
Luck, S. J., Chelazzi, L., Hillyard, S. A., & Desimone, R. (1997). Neural mechanisms of spatial selective attention in areas V1, V2 and V4 of macaque visual cortex. Journal of Neurophysiology, 77, 24–72.
Maldonado, P. E., Friedman-Hill, S., & Gray, C. (2000). Dynamics of striate cortical activity in the alert macaque: II. fast time scale synchronization. Cerebral Cortex, 10(11), 1117–1131.
Marmarelis, P. Z., & Marmarelis, V. Z. (1978). Analysis of physiological systems: The white-noise approach. New York: Plenum Press.
Martinez, D. (2006). Oscillatory synchronization requires precise and balanced feedback inhibition in a model of the insect antennal lobe. Neural Computation, 17, 2548–2570.
McAdam, C. J., & Maunsell, J. H. R. (1999). Effects of attention on orientation-tuning function of single neurons in macaque cortical area V4. Journal of Neuroscience, 19(1), 431– 441.
Miller, E. K., Gochin, P. M., & Gross, C. G. (1993). Suppression of visual responses of neurons in inferior temporal cortex of a awake macaque by addition of a second stimulus. Brain Research, 616, 25–29.
Mishra, J., Fellous, J.-M., & Sejnowski, T. J. (2006). Selective attention through phase relationship of excitatory and inhibitory input synchrony in a model cortical neuron. Neural Networks, 19, 1329–1346.
Mitra, P. P., & Pesaran, B. (1999). Analysis of dynamic brain imaging data. Biophysical Journal, 76(2), 691–708.
Moore, T., & Amstrong, K. M. (2003). Selective gating of visual signals by microstimulation of frontal cortex. Nature, 421, 370–373.
Moran, J., & Desimone, R. (1985). Selective attention gates visual processing in the extrastriate xortex. Science, 229, 782– 784.
Motter, B. C. (1993). Focal attention produces spatially selective processing in visual cortical areas V1, V2 and V4 in the presence of competing stimuli. Journal of Neurophysiology, 70(3), 909–919.
Reynolds, J., Chelazzi, L., & Desimone, R. (1999). Competitive mechanisms subserve attention in macaque areas V2 and V4. Journal of Neuroscience, 19(5), 1736–1753.
Reynolds, J. H., & Desimone, R. (2003). Interacting roles of attention and visual salience in V4. Neuron, 37, 853–863.
Roelfsema, P. R., Lamme, V. A. F., & Spekreijse, H. (2004). Synchrony and covariation of firing rates in the primary visual cortex during contour grouping. Nature Neuroscience, 7(9), 982–991.
Rolls, E. T., & Tovee, M. J. (1995). The responses of single neurons in the temporal visual cortical areas of the macaque when more than one stimulus is present in the receptive field. Experimental Brain Research, 103, 409–420.
Rols, G., Tallon-Baudry, C., Girard, P., Bertrand, O., & Bullier, J. (2001). Cortical mapping of gamma oscillations in areas V1 and V4 of the macaque monkey. Visual Neuroscience, 18(4), 527–540.
Schoffelen, J. M., Oostenveld, R., & Fries, F. (2005). Neuronal coherence as a mechanism of effective corticospinal interaction. Science, 308(5718), 111–113.
Smith, A. T., Singh, A. L., Williams, A. L., & Greenlee, M. W. (2002). Estimating receptive field size from fMRI data in human striate and extrastriate visual cortex. Cerebral Cortex, 11, 1182–1190.
Smiyukha, Y., Mandon, F., Galashan, F. O., Neitzel, S. D., & Kreiter, A. (2006). Attention-dependent switching of interareal synchronization between V4 neurons and different subpopulations of their V1 afferents. Soc. Neurosci. Abtract, 32, 11.2.
Spitzer, H., Desimone, R., & Moran, J. (1988). Increased attention enhances both behavioural and neuronal performance. Science, 240(4850), 338–340.
Tass, P., Rosenblum, M. G., Weule, J., Kurths, J., Pikovsky, A., Volkmann, J., et al. (1998). Detection of n:m phase locking from noisy data: Application to magnetoencephalography. Physical Review Letters, 81(15), 3291–3294.
Taylor, K., Mandon, S., Freiwald, W. A., & Kreiter, A. K. (2005). Coherent oscillatory activity in monkey area V4 predicts successfull allocation of attention. Cerebral Cortex, 15(9), 1424–1437.
Thomson, D. J. (1982). Spectrum estimation and harmonic analysis. Proceedings of the IEEE, 70, 1055–1096.
Tiesinga, P. H. E. (2005). Stimulus competition by inhibitory interference. Neural Computation, 17, 2421–2453.
Traub, R. D., & Miles, R. (1991). Neuronal networks of the hippocampus. Cambridge, UK: Cambridge Univ. Press.
Treue, S., & Martinez Trujillo, J. C. (1999). Feature-based attention influences motion processing Gain in macaque visual cortex. Nature, 399, 575–579.
Treue, S., & Maunsell, J. H. R. (1996). Attentional modulation of visual motion processing in cortical areas MT and MST. Nature, 382, 539–541.
Treue, S., & Maunsell, J. H. R. (1999). Effects of attention on the processing of motion in macaque visual cortical areas MT and MST. Journal of Neuroscience, 19(17), 7603–7616.
van der Togt, C., Kalitzin, S., Spekreijse, H., Lamme, V. A. F., & Super, H. (2006). Synchrony dynamics in monkey V1 predicts success in visual detection. Cerebral Cortex, 16(1), 136–148.
Tsodyks, M., Pawelzik, K., & Makram, H. (1998). Neural networks with dynamic synapses. Neural Computation, 10, 821–835.
Tsodyks, M., Uziel, A., & Makram, H. (2000). Synchrony generation in recurrent networks with frequency-dependent synapses. Journal of Neuroscience, 20(RC50), 1–5.
Womelsdorf, T., Fries, P., Mitra, P. P., & Desimone, R. (2006). Gamma-band synchronization in visual cortex predicts speed of change detection. Nature, 439, 733–736.
Zeitler, M., Fries, P., & Gielen, S. (2006). Assessing neuronal coherence with single-unit, multi-unit and local field potentials. Neural Computation, 18(9), 2256–2281.
Acknowledgements
We thank Michelle Lilith for fruitful discussions during and after the Advanced Course in Computational Neuroscience in Arcachon (France), where we started this project. This project was funded by the Netherlands Organisation for Scientific Research (NWO 051.02.050).
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Author information
Authors and Affiliations
Corresponding author
Additional information
Action Editor: Matthew Wiener
Appendices
Appendix 1
The parameter values for the voltage-dependent Na + and K + currents were described by (Traub and Miles 1991):
where V T = − 58 mV and V S = − 10 mV (Destexhe and Paré 1999). \(\bar{g}\) is the maximum conductance (\(\bar{g}_{\rm Na} = 361.2\, 10^{-4}\) S/cm2, \(\bar{g}_{Kd} = 70\, 10^{-4}\) S/cm2), m, h and n are the time-varying gate variables, E Na = 50 mV is the sodium reversal potential, E Kd = − 90 mV the potassium reversal potential, α is the forward and β the backward rate. The “delayed-rectifier” K + current was described by:
In our model both the interneurons and the output neuron receive background synaptic input as received by cortical neurons in vivo, represented by fluctuating background conductance injections in the soma. These conductances are produced by an Ornstein–Uhlenbeck process, as described by Destexhe et al. (2001):
where g 0 is the mean conductance, τ is the conductance time constant, σ 2 is the variance of the conductance and χ(t) is Gaussian white noise with zero mean and a standard deviation of 1. For the inhibitory background conductance of the output neuron we use g i0 = 57.3 nS, τ i = 10.49 ms and σ i = 6.0 nS with a reversal potential E i = − 75 mV. For the excitatory background conductance of the output neuron we use g e0 = 12.1 nS, τ e = 2.73 ms and σ e = 3.0 nS with a reversal potential E e = 0 mV.
For the interneurons the average conductances and the standard deviations of these conductances are 50% of the corresponding values of the output neuron.
For the implementation in NEURON of the passive and active properties and of the synaptic background we used parts of the code of example 5 of the NEURON tutorial from the Obidos 2004 course http ://www.neuron.yale.edu/ftp/neuron/contrib/obidos_tutorials/.
Appendix 2
In this appendix we provide a rough estimate of the currents injected as a noisy background and the currents due to the stimulus related inputs to the conductance based neurons in our model. Current and conductance are related by:
where G i,e is the total conductance, E i,e the reversal potential and V the membrane potential for the inhibitory (i) or excitatory (e) input. The total input current is the sum over all excitatory and inhibitory currents. The amount of current is time-dependent since it is a function of the fluctuating membrane potential and of the total amount of the conductance at time t. We will approximate the currents by taking time-averages for conductance and membrane potential. The mean value for the membrane potential V depends on the contribution of all excitatory and inhibitory inputs. From the simulations in NEURON we know that the average membrane potential < V > of the interneurons (output neuron Y) is about − 55 mV (− 60 mV, respectively), i.e., well between the rest membrane potential near − 75 mV and the threshold for action potential generation. The value of the reversal potentials are E i = − 75 mV and E e = 0 mV for both the interneuron and excitatory neuron, with these values the average currents caused by the noisy background in the interneurons is
The average input currents caused by the stimuli in the interneurons is
with N i,e the number of input spike trains, < f I > the average firing rate of the stimulus-related inhibitory input spike trains, < r > the average rate of the stimulus-related excitatory Poisson spike trains and τ i,e the rise time of the α-synapses. This means that the inhibitory input current of the interneurons is only due to the noisy background and the excitatory current comes for 58% from the noisy background and 42% comes from the stimulus related input.
The average input currents caused by the noisy background in the output neuron are
where g 0 is the mean conductance of the conductance injections in the soma as described in Appendix 1. The average input currents caused by the stimuli in output neuron Y are
with N i,e the number of input spike trains, < f I > the averaged inhibitory input spike train, < r > the averaged rate of the excitatory Poisson spike train, τ i,e the rise time of the α-synapses. The noisy background is responsible for 84% of the total inhibitory input current and the stimulus related input contributes 16%. The excitatory current contributes for 54% of the noisy background and 46% is due to the stimulus related input.
The input current for the ‘preferred stimulus only’ condition (p) consist of the background current plus the excitatory related to the preferred stimulus. For the ‘non-preferred stimulus only’ condition (np) the excitatory stimulus-related current to the inhibitory neuron is larger than the current for the preferred stimulus because of the larger value of the synaptic connections for the non-preferred stimulus input.
Rights and permissions
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zeitler, M., Fries, P. & Gielen, S. Biased competition through variations in amplitude of γ-oscillations. J Comput Neurosci 25, 89–107 (2008). https://doi.org/10.1007/s10827-007-0066-2
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10827-007-0066-2