Neural mechanisms of interstimulus interval-dependent responses in the primary auditory cortex of awake cats
© Sakai et al; licensee BioMed Central Ltd. 2009
Received: 24 September 2008
Accepted: 10 February 2009
Published: 10 February 2009
Primary auditory cortex (AI) neurons show qualitatively distinct response features to successive acoustic signals depending on the inter-stimulus intervals (ISI). Such ISI-dependent AI responses are believed to underlie, at least partially, categorical perception of click trains (elemental vs. fused quality) and stop consonant-vowel syllables (eg.,/da/-/ta/continuum).
Single unit recordings were conducted on 116 AI neurons in awake cats. Rectangular clicks were presented either alone (single click paradigm) or in a train fashion with variable ISI (2–480 ms) (click-train paradigm). Response features of AI neurons were quantified as a function of ISI: one measure was related to the degree of stimulus locking (temporal modulation transfer function [tMTF]) and another measure was based on firing rate (rate modulation transfer function [rMTF]). An additional modeling study was performed to gain insight into neurophysiological bases of the observed responses.
In the click-train paradigm, the majority of the AI neurons ("synchronization type"; n = 72) showed stimulus-locking responses at long ISIs. The shorter cutoff ISI for stimulus-locking responses was on average ~30 ms and was level tolerant in accordance with the perceptual boundary of click trains and of consonant-vowel syllables. The shape of tMTF of those neurons was either band-pass or low-pass. The single click paradigm revealed, at maximum, four response periods in the following order: 1st excitation, 1st suppression, 2nd excitation then 2nd suppression. The 1st excitation and 1st suppression was found exclusively in the synchronization type, implying that the temporal interplay between excitation and suppression underlies stimulus-locking responses. Among these neurons, those showing the 2nd suppression had band-pass tMTF whereas those with low-pass tMTF never showed the 2nd suppression, implying that tMTF shape is mediated through the 2nd suppression. The recovery time course of excitability suggested the involvement of short-term plasticity. The observed phenomena were well captured by a single cell model which incorporated AMPA, GABAA, NMDA and GABAB receptors as well as short-term plasticity of thalamocortical synaptic connections.
Overall, it was suggested that ISI-dependent responses of the majority of AI neurons are configured through the temporal interplay of excitation and suppression (inhibition) along with short-term plasticity.
The perceptual quality of successive acoustic signals considerably varies depending on the inter-stimulus intervals (ISI). For example, when click signals are repetitively presented at ISI ≥ ~30 ms, individual signals are clearly heard as discrete events ; at ISI ≤ ~30 ms, those are perceptually fused together [2, 3]. This ISI boundary, denominated as "temporal-order threshold," has long been considered as an important indicator of temporal resolving capacity of the auditory system (reviewed by ). Another example of ISI-dependent perception is categorical perception of stop consonant-vowel syllables (CV syllables): if ISI between the consonant release and voicing onset (voice onset time [VOT]) is shorter than a critical value (VOT boundary), the consonant is perceived as "voiced"; if ISI exceeds this value, the consonant is perceived as "unvoiced" (eg.,/da/-/ta/continuum ). In many languages including English, the VOT boundary lies at 20–40 ms with some variance among place of articulation (reviewed by ). Monkeys , chinchillas  and birds  all place the VOT boundary at approximately the same value, indicating the categorical perception of CV syllables does not necessarily arise from a specific human speech mechanism but is based, at least partially, on general properties of the auditory system.
Case studies of patients with stroke lesions restricted to the bilateral primary auditory cortex (AI) reported that (1) their temporal-order threshold was elongated up to ~100 ms [10, 11] and (2) they were severally impaired in the categorical perception of CV syllables [11, 12]. These findings suggest that AI is critically involved in ISI-dependent differential perception regardless of whether the signals are phonetic or non-phonetic (reviewed by ).
The previous single unit study in un-anesthetized animals AI revealed that click trains produce qualitatively distinct response features depending on ISI: at ISI ≥ ~30 ms, stimulus-locking responses dominate; at ISI ≤ ~30 ms, responses occur only at the onset of the train . Similar finding was obtained for AI responses to CV syllables: at ISI (VOT) ≥ ~30 ms, stimulus-locking responses take place to both the consonant and vowel; at ISI (VOT) ≤ ~30 ms, responses occur only to the consonant . Since these neurophysiological ISI boundaries (~30 ms) match both the temporal-order threshold and VOT boundary (see above), it was suggested that the neural processes constraining AI stimulus-locking responses are also responsible for the perceptual boundaries of phonetic/non-phonetic acoustic signals .
The present study, by employing a single unit recording technique in un-anesthetized cats, thoroughly analyses how AI neurons respond to click trains of variable ISI. Then, by modeling the observed responses, we extract general principles governing various ISI-dependent behaviors of AI neurons especially stimulus-locking responses.
Response Features in the Click-train Paradigm
The results are based on 116 AI neurons that showed statistically significant excitatory responses to the click stimuli (see Methods). We classified those neurons into 2 types depending on whether they had the capacity for stimulus-locking responses to click trains (synchronization type: n = 72) or not (non-synchronization type: n = 44).
Response variables (mean ± S.D.) in the click-train paradigm
Best ISI (region β)
174 ± 112 ms
71.4 ± 19.6 ms
30.6 ± 25.9 ms
12.1 ± 13.5 ms
(5.7 ± 8.9 Hz)
(14 ± 51 Hz)
(33 ± 38 Hz)
(83 ± 74 Hz)
80.6 ± 23.8 ms
33.1 ± 28.7 ms
11.4 ± 16.9 ms
(13 ± 42 Hz)
(30 ± 34 Hz)
(87 ± 59 Hz)
14.9 ± 11.7 ms
(67 ± 85 Hz)
Effects of the Stimulus Level
To examine effects of the stimulus level, we adhered to the click-train paradigm at various stimulus levels (in pe-SPL; see Methods). We examined 34 synchronization neurons, firing activities of which could be isolated long enough for the detailed analysis. Among them, 29 neurons exhibited responses at 20 dB below the best SPL while 20 neurons did at 40 dB below the best SPL.
Accumulation Effects: Region β
It has been widely reported that repetitive stimulation exerts accumulation effects on auditory neurons especially those in the auditory cortex (reviewed by [16, 17]). We addressed whether and how such effects influenced the observed responses. In the present and the following section, we paid special attention to regions β and δ, where the causative relationship between a given stimulus and spikes can be clearly identified.
To address the durability of this response degeneration, we presented a single click (= probe stimulus) at 1.0, 1.8 or 3.6 s after the termination of the click trains that were delivered at best ISI (Fig. 4D). Figure 4E illustrates the normalized #spikes (see Methods) at each time point (Roman numerals correspond to the probe stimuli depicted in Fig. 4D). The value was still smaller than unity (broken line) at 1.0 s but became equivalent to unity at 1.8 and 3.6 s (P < 0.05; Fisher's PLSD test following ANOVA). This indicates that the response degeneration lasted 1.0–1.8 s after the termination of the trains.
The above features, both during and after the click trains, correspond to the phenomenon, so-called "frequency-dependent depression" .
Accumulation Effects: Region δ
At relatively long ISIs in region δ, spikes occasionally occurred after the onset responses (Fig. 1A–1). Such activities potentially hinder temporal precision in measuring the duration of the onset responses. Those activities were sufficiently suppressed at ISIs at or shorter than 0.9 multiples of the γ-δ border in all synchronization type neurons examined. This relative ISI value is uniformly adopted in the following analysis of region δ.
The onset response in region δ was typically followed by an unresponsive period (Fig. 5A). To estimate the recovery time course from this response degeneration, we presented a probe stimulus at 1.0, 1.8 or 3.6 s after the termination of the train (Fig. 5E). The normalized #spikes at each time point was illustrated as filled bars in Figure 5F (same format as Fig. 4E). For comparison, the data in Figure 4E (broken bars) was appended, which depicts the recovery time course after region β. The major findings are: (1) at 1.0 s, the value was smaller after region δ than after region β (P < 0.05, t-test); (2) at 1.8 s, the value after region δ was still below unity (P < 0.05) while the value after region β got equivalent to unity. Together, it appears that response degeneration is more profound and longer lasting after region δ than after region β. This fits well to the principle of "frequency-dependent depression": response degeneration grows larger and longer at higher stimulation rates . It is, thus, plausible that the unresponsive period in region δ was caused by such intense "frequency-dependent depression" as to abolish firing activities for a while.
Involvement of Post-activation Suppression in Stimulus-locking Responses
The neural processes of stimulus-locking responses (e.g., Fig. 4A) can be glimpsed if we pay special attention to the vector strength, the origin of the tMTF (see Methods). The vector strength measures the degree of temporal confinement of spikes against stimuli. It reaches a maximum (= 1.0) when spikes occur in exactly the same period with reference to the individual stimuli, and spikes (regardless of whether evoked or spontaneous) are completely absent in the remaining period. On the other hand, it reaches a minimum (= 0.0) when spikes occur entirely independently of the stimuli. It is, thus, quite conceivable that the capacity for stimulus-locking responses arise from neural processes that temporally confine spikes. To examine this, we conducted the single-click paradigm in which the dynamics of neural activity after single click presentation was qualified with referring to the spontaneous firing rate (excitation, suppression or spontaneous-level activities; see Methods). This analysis was performed on 35 neurons (non-synchronization type, n = 8; synchronization type, n = 27) that showed an appreciable amount of spontaneous firing rate for detecting suppression.
Sequence of response periods observed in the single-click paradigm
Modeling of AI Temporal Behavior
Consequently, stimulus-locking responses, if any, would be weak (Fig. 7C, broken curve, region δ). Second, when the 2nd stimulus is given during the fast-IPSP (Fig. 7B–2) or slow-IPSP (Fig. 7B-3), it may hardly elicit firing responses. Stimulus-locking responses would be negligibly small (Fig. 7C, broken curve, region γ and region α). Third, when the 2nd stimulus is given during the slow-EPSP (Fig. 7B–4), it may readily elicit firing responses. Those discharge clusters are expectedly separated by the fast-IPSP (ref., Fig. 7A). As a consequence, clear stimulus-locking responses would take place (Fig. 7C, broken curve, region β). Interestingly, the "predicted tMTF for paired stimuli"(broken curve) resembles the "observed tMTF on click trains" that was schematically drawn based on the data of the synchronization type with band-pass tMTF (solid curve; based on Table 1) except for the short ISI potion (shaded zone). This inconsistency may arise from the fact that if the stimuli are repetitively delivered at short ISI, impact of initial several stimuli leads to such intense "frequency-dependent depression" as to abolish firing responses to the following stimuli (Fig. 5A). Collectively it is suggested that (1) as a principle, the temporal interplay of the PSP components underlies AI stimulus-locking responses and (2) at short ISI, intense frequency-dependent depression abolishes stimulus-locking responses.
Nonetheless, the 4-receptor version does not necessarily parallel our physiological observation. First, this model predicts the occurrence of stimulus locking at much shorter ISI (e.g., 14.3 ms ISI; Fig. 8B–1) compared to physiological observation. In other words, the predicted β-γ border (4.5 ms; Fig. 8F, open arrow) was much shorter than the observed one (~30 ms, Table 1). Second, this model predicts "skipping" of firing responses at short ISI: for example, the spikes were expected to occur every other stimulus at 11.8 ms ISI (Fig. 8C–1) and every three stimuli at 5.6 ms ISI (Fig. 8D–1). These features, however, have been scarcely encountered in our physiological recording (but  in somatosensory cortex of anesthetized rats). In fact, the shortening of ISI led to the systematic reduction of responsivity to the latter clicks in a given train (e.g., Fig. 4A) while responsivity to the initial several clicks being relatively well preserved, as reported by plenty of single unit studies using periodic signals [14, 20, 23, 28, 29, 38]. These two discrepancies were diminished when we provided the external input with both frequency-dependent depression (FDD) and paired-pulse facilitation (PPF) ("full version"). Specifically, the capacity for stimulus locking at shorter ISI was considerably weakened (cf., Fig. 8B–2 to 8B–1, 8C–2 to 8C–1, and 8D–2 to 8D–1) with only the onset response being manifested at the γ-δ border (Fig. 8C–2) and in region δ (Fig. 8D–2). As a consequence, the β-γ border prolonged to 25 ms (Fig. 8F, filled arrow).
At longer ISI, in marked contrast to short ISI, the 4-receptor version and full version predicted similar response features. For instance, at 164 ms the initial two clicks elicited stimulus-locking responses while the following ones did not (Fig. 8E–1,-2), similarly to region α (Fig. 1A–1). This indicates that the frequency-dependent depression and paired-pulse facilitation are much less influential at longer ISI.
Next, we investigated the main constraints on responses at long ISI. First examined was the contribution of GABAB-receptor-mediated IPSP. When the conductance of the GABAB receptors (GABAB) was reduced by factor of 0.1 (Fig. 8E–3; conductance of the other receptors was kept constant, same as followings), responsivity to the initial several stimuli was enhanced (cf., Fig. 8E–2) thereby prolonging the α-β border (Fig. 8F, inset, dotted arrow). Contrariwise, when the conductance of the NMDA receptors (NMDA) was reduced, responsivity to the 2nd and latter stimuli decreased (Fig. 8E–4) enhancing low-cut effect. Manipulation of the conductance of the AMPA or GABAA receptors did not greatly influence responses at long ISI (data not shown). Taken together, it is suggested that the GABAB and NMDA act as main constraints on the temporal filtering at long ISI: if the former is reduced, the tMTF tends to be more low-pass whereas if the latter is reduced, the tMTF more band-pass.
By using click signals, we investigated neural mechanisms underlying ISI-dependent responses of the AI neurons which had the capacity for stimulus-locking responses (synchronization type; Fig. 1A and Fig. 1B). The β-γ border, i.e., the shorter cutoff ISI for stimulus-locking responses, lay at on average ~30 ms (Table 1) and was level tolerant over high SPLs (Fig. 3C). The time course of excitability during (Figs. 4C and 5D) and after (Fig. 5F) the click trains suggested the involvement of short-term plasticity of thalamocortical synaptic connections. Comparison between response features to the click trains and a single click (Table 2) led to the notion that the temporal interplay of excitation and suppression basically determines the capacity for stimulus-locking responses as well as tMTF shape. A single-cell dynamic model well replicated the physiological data (Figs. 8F and 8G) suggesting that ISI-dependent responses of the synchronization neurons are configured through the temporal interplay of the post-synaptic potentials (Fig. 7A) along with short-term plasticity of thalamocortical synaptic connections.
Perceptual Relevance of the Observed AI Responses
Case studies about AI-impaired patients indicated that AI is responsible for the temporal-order threshold (see Background). The mean value of the β-γ border (~30 ms; Table 1), i.e., shorter cutoff ISI for stimulus-locking responses (e.g., Figs. 1A–1 and 1A–2), agree well with the temporal-order threshold . The β-γ border was nearly invariant at the best SPL and 20 dB below it (Fig. 3C, dotted line and filled triangle in the right, respectively) in accordance with the level tolerance of the temporal-order threshold over high SPL . These findings strongly support the notion that the β-γ border serves as a neural correlate of the temporal-order threshold. Accordingly, it could be postulated that click trains are clearly heard as a series of "discrete" events as long as stimulus-locking responses dominate in AI. By taking count of the view that the temporal-order threshold and perceptual boundary of CV syllables at least partially share common neural processes (see Background), it is possible that the β-γ border also serves as the basis for the perceptual boundary of CV syllables.
The psychological studies have revealed that a total body of click trains of ISI ≥ ~30 ms produces two kinds of sensation: at ISI of ~30–200 ms, it leads to "rhythm" percept; at ISI ≥ ~200 ms, rhythm percept fades away while the sensation of "fluctuation" remains . Since the value of the α-β border (~175 ms; Table 1) as well as its level tolerance (Fig. 3C, open triangles) is consistent with the rhythm-fluctuation boundary (~200 ms), it is possible that region α dominantly represents "fluctuation" whereas region β does "rhythm" percept.
At ISI < ~30 ms where individual clicks are no longer clearly heard as discrete events, a total body of click trains leads to three kinds of sensation with partial overlap [31–33]. At ISI of ~3–30 ms, a buzz or rattle like sensation is produced, defined as "roughness"; at ISI of ~5–15 ms, tone quality of sensation dominates, whose perceived frequency is directly related to waveform periodicity (defined as "periodicity pitch"); at ISI < several ms, another mode of pitch sensation dominates, which depends on the fundamental frequency ("spectral pitch"). None of the ISI range of these sensations does not fit to that of region γ or region δ (Table 1) making it unlikely that our single unit data have direct relevance to these sensations. Steinschneider and colleagues  suggested that periodicity pitch may be represented in AI by oscillatory neuronal ensemble responses locking to temporal envelope, and spectral pitch by rate-place coding that is sensitive to both the fundamental frequency and other harmonics in the train. Further insight into neurophysiological bases of these sensations would be obtained by coordinated single/multiunit recordings and psychoacoustic experiments.
Comparison to Previous AI Studies in Un-anesthetized Animals regarding to Cell-type Classification
An accumulating body of single unit studies in un-anesthetized animals has investigated AI responses to periodic acoustic signals such as click trains [14, 23] and amplitude-modulated sounds [35–38]. Irrespective of methological differences (e.g., stimulus configuration, electrode properties and statistics), these studies appear to agree that AI neurons comprise largely two subsets: one subset responds predominantly at long ISI (≥ ~30 ms, occasionally extending to 10 ms or less) in a stimulus-locking manner (i.e., temporal code) while another subset does at short ISI (< ~10 ms) in a sustained manner (i.e., rate code).
The following three findings indicate that our synchronization neurons (eg., Fig. 1A) correspond to the subset that was reported to conduct temporal code. First, the shorter cutoff ISI for stimulus-locking responses in those studies (30–40 ms) was similar to the β-γ border of our synchronization neurons (~30 ms; Table 1). Second, at shorter ISI (10–30 ms), the studies reported that the neurons intermittently fired without stimulus locking. This feature resembles region γ responses of our synchronization neurons. Third, at much shorter ISI (≤ 10 ms), the studies reported that only the onset response was evident. This feature, as well as its cutoff ISI, is quite akin to region δ of our synchronization neurons.
Our non-synchronization neurons with high-pass rMTF showed non-stimulus-locking responses during the presence of the click trains (eg., Fig. 1D). Such responses occurred only at short ISI (< ~10 ms) with shorter ISI leading to larger driven rate (D-3). This feature closely resembles the responses of the subset that was reported to conduct rate code. On the other hand, our non-synchronization neurons with low-pass rMTF (eg., Fig. 1C) do not correspond to either subset mentioned above. They may belong to the "unclassified neurons" in Wang and colleagues' study , which were reported to respond in some manner to click signals without clearly defined stimulus-locking responses or non-stimulus-locking rate responses.
Comparison to Previous AI Studies of Neural Mechanisms underlying Stimulus-locking Responses
de Ribaupierre and colleagues  revealed that AI stimulus-locking responses are related to the temporal interplay of depolarization and hyperpolization. While this interplay potentially results from non-synaptically mediated after-hyperpolarization , growing evidence indicates that this interplay is based mainly on the sequence of EPSPs and IPSPs [25, 26]. In particular, Cox and colleagues  proved in rat AI slice preparations that EPSPs and IPSPs are mediated chiefly through AMPA/NMDA receptors and GABAA/GABAB receptors, respectively (Fig. 7A). To date, however, it has been unclear whether and how these PSP components are related to stimulus-locking responses of AI neurons.
On the other hand, there is physiological data suggesting that AI stimulus-locking responses are mediated through other neural processes than the interplay of PSP components. For example Wehr and Zador , by employing the whole-cell recording technique on ketamine-anesthetized rats AI, measured the excitatory and inhibitory synaptic conductances elicited by click pairs of variable ISI. They found that inhibitory conductances were too short-lived to account for suppression of spiking responses to the 2nd click which was delivered at ISIs ≥ several hundreds milliseconds. Eggermont , based on physiological data, proposed a model in which presynaptic facilitation and depression determine the low-pass characteristics of AI stimulus-locking responses.
Our single-cell dynamic model comprehensively integrates the above findings/suggestions in that the capacity for stimulus-locking responses (i.e., region β) is explained in terms of the temporal interplay of the PSP components along with short-term plasticity (Fig. 8A–5). This view is compatible with the pioneering work of Grothe  that proved the critical contribution of EPSPs and IPSPs to stimulus-locking responses of auditory brainstem neurons for encoding interaural time difference. Furthermore, our model can explain the neural processes that give rise to the other ISI-dependent AI responses such as in regions γ, δ and α (Figs. 8B–2, D–2, and 8E–2).
It has been proved, in a gap-in-noise detection paradigm where leading and trailing wideband noise had the same frequency content, that the minimum gap between the noises (i.e., ISI) for stimulus-locking responses was 30 ms for a 20-ms leading noise, 10 ms for a 50-ms leading noise, and reached an asymptote of 5 ms for a 200-ms leading noise [15, 42]. This implies that temporal resolving capacity of individual AI neurons is not fixed, but varies dynamically as depending on the duration of the leading signal. Since our model was based on the physiological data obtained using only 1-ms long clicks, this is best suited for neural processes that are triggered at the stimulus onset, not later in the stimulus.
The present study revealed that recovery period of AI spiking activities was in the order of seconds (Fig. 5F). This range is consistent with the values reported in awake (> 1s ) and ketamine-anesthetized animals (> 500 ms ) but is much longer than those measured in barbiturate-anesthetized animals (20–200 ms [44–46]). A number of observations indicate that barbiturate reduces spontaneous and evoked spiking activities [28, 47, 48]. This leads to a conjecture that frequency-dependent depression, which results chiefly from temporal exhaustion of readily releasable neurotransmitter pool , is much less potent and less durable under barbiturate anesthesia than the other conditions. Under barbiturate anesthesia, weight of the influence on the recovery period may shift from frequency-dependent depression, which lasts for seconds [17, 18], to IPSPs which extend for maximally several hundreds milliseconds.
Relative Contribution to AI Stimulus-locking Responses: Intra-AI Processing vs. Sub-AI Processing
There is a marked resemblance between the predicted tMTF (Fig. 7C, broken curve) and observed tMTF (solid curve) except for the short ISI portion (shaded area; for reason, see above). Remarkably, these tMTF curves derived from rather different experimental conditions: the former was based on the data obtained by electrical stimulation to rat AI slice preparations (Fig. 7A); the latter, during acoustic stimulation in un-anesthetized cats. By taking account of a large difference in the shorter cutoff ISI for stimulus-locking responses (i.e., β-γ border; Fig. 7C, filled triangle) between un-anesthetized (~30 ms; present study and ) and anesthetized animals (≥ 50 ms [28, 49, 50]), it seems likely that the neural processes of AI stimulus-locking responses are much less susceptible to the mechanical infringement of slicing the brain and species difference than the pharmacological effects of anesthetics. Some anesthetics, such as pentobarbital, are known to potentate GABAA-ergic inhibition . This may enhance the fast-IPSP and consequently diminish the slow-EPSP, especially in its early phase (ref., Fig. 7A). As a result, the transition point from the fast-IPSP to slow-EPSP will shift afterward (horizontal arrow), prolonging β-γ border (Fig. 7C, horizontal arrow). Slice preparations, on the other hand, do not suffer such artificial enhancement of GABAA-ergic inhibition but retain local circuitry. These considerations favor the idea that AI stimulus-locking responses are elaborated mainly through intra-AI processing rather than simple preservation of sub-AI processing. The above idea receives support from the previous studies that directly compared the best ISI within pairs of functionally connected medial geniculate body (MGB) neuron and AI neuron [37, 51]. In those studies, spiking activities of individual MGB and AI neurons were simultaneously recorded and the functional connection was confirmed if their activities showed a single cross-correlogram peak within 1–5 ms lag time, the MGB neuron leading the AI neuron, under both spontaneous and stimulus-driven conditions. Importantly, no rank correlation was revealed for the best ISI; MGB neurons with longer (shorter) best ISI did not preferentially connect with AI neurons with longer (shorter) best ISI. This suggests that the generally observed prolongation in the best ISI from MGB to AI [52, 53] cannot be simply attributed to a degradation of temporal resolution due to intrinsic membrane properties or synaptic delay but is rather due to more elaborated intra-AI processing.
Our physiological observation suggest that β-γ border, the shorter cutoff ISI for stimulus-locking responses of AI neurons, serves as a neural correlate of the temporal-order threshold and VOT boundary of CV syllables. The present modeling study supports the idea that the observed ISI-dependent responses are largely mediated through temporal interplay between EPSPs and IPSPs at the thalamocortical synapse along with its short-term plasticity.
The parameter values in our model were not directly measured in the current recording study but were determined referring to other published data (see Methods). Under this proviso, the weight of AMPA, GABAA, NMDA and GABAB receptors was arranged to better predict the observed phenomena. In fact, relative weight of the receptors may vary across neurons, and other cellular and/or network mechanisms may also contribute to the observed phenomena. Further insight into this issue would be glimpsed by measuring the membrane potential of the AI neurons during acoustic stimulation or analyzing the effect of selective receptor antagonists on their stimulus-response features.
Experiments were performed in a manner consistent with the Guidelines for Animal Experiments, University of Yamanashi, and the Guiding Principles for the Care and Use of Animals approved by the Council of the Physiological Society of Japan. Animal preparation, recording, and histology procedures were the same as in our previous report [54–57]. Briefly, the cats underwent surgery under pentobarbital sodium anesthesia. Aluminum cylinders and metal blocks were implanted for subsequent extracellular recording and restraining of the animal's head, respectively. At least two weeks were allowed for recovery before the recording. During the recording, the animal was placed in an electrically shielded and sound-attenuated room with its body wrapped in a cloth bag and its head restrained with a holding bar. The animals were kept awake throughout the recording period and were monitored with an online surveillance camera and electroencephalography. When drowsiness was suspected, the cat was awakened by gently tapping its body using a remote-controlled tapping tool or by briefly opening the door of the room. A glass microelectrode (tip diameter, 1.8–2.5 μ m; resistance, 2–3 MΩ; filled with 2 M NaCl) was inserted into AI. Tone bursts of variable frequencies and sound pressure levels (SPLs) were presented as search signals. Single unit activities were recorded and their occurrences were identified using a window discriminator. The spike-occurrence outputs were captured on a Pentium-based computer with a time resolution of 2 μ s as the digital input for data analysis. The animals sometimes voluntarily moved during recording sessions, creating artifacts in the recording. By carefully checking the online monitor screens of the animal and firing activities, these motion artifacts were marked in real time on the recording-computer while recordings were in progress. Data with artifacts were rejected. Daily recording sessions lasted < 6 hours, and the total duration of the experiment continued for 2–6 months per animal. At the termination of the experiment, some recording sites were marked with electrolytic lesions. The animals were then sacrificed with an overdose of pentobarbital sodium and perfused with 10% formalin. The brain was cut in transverse sections and stained with neutral red. The recording sites were reconstructed based on the electrolytic lesions and electrode tracks.
Sound generation and delivery
The sound signals were generated using user-written programs in MATLAB (Mathworks, Natick, MA) on a Pentium-based computer. The signals were fed into a 12-bit digital-to-analogue converter (BNC2090; National Instruments, Austin, TX) at a sampling interval of 100 kHz and to an eight-pole Chebyshev filter (P-86; NF Electric Instruments, Yokohama, Japan) with a high cutoff frequency of 20 kHz. The output was attenuated and sent to a low-output-impedance power amplifier (PMA2000III; Denon Electronic GmbH, Ratingen, Germany), and then the sound signals were presented from a speaker (K1000; AKG Acoustics, Wien, Austria) placed 2 cm away from the auricle contralateral to the recording site. We equalized and calibrated the sound delivery system between 128 and 16,000 Hz in 8 Hz steps, and the output varied by +/-1.5 dB. One set of stimuli was presented at variable intervals ranging 4.0–5.5 s. We employed (1) click-train paradigm and (2) single-click paradigm (see below).
Once single unit activities were isolated, we conducted a click-train paradigm: rectangular clicks (1-ms duration) were delivered in a train fashion (0.5-s duration; Fig. 1A–1, horizontal line). Since the neurons examined in the present study are included in our previous study for periodicity coding , click repetition rate was systematically varied at 0.5 (or 1.0, 2.0 or 4.0) Hz steps over 2.5–100 (or 2.5–480) Hz. Typical repetition rates were 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, ... 99.0, 99.5, 100 Hz (Fig. 1B–1, ordinate, right). Consequently, inter-stimulus intervals (ISI), which is the reciprocal of the repetition rate, varied 10.0, 10.1, 10.2, ... 200, 222, 250, 333, 400 ms (ordinate, left). The intensity of the clicks for 20 repetitions (peak equivalent sound pressure level [pe-SPL]) was set at 20–80 dB (10 dB steps; occasionally 5 dB steps). An identical set of signals was presented 5–20 times. The analytic time window was set at 50–500 ms after the onset of the click trains for eliminating nonspecific effect of the onset response that takes place regardless of ISI value (see below). This kind of compensation has been often adopted in the previous AI studies [14, 28, 58, 59].
where xi = cos θi and yi = sin θi, n: total number of spikes, and each spike is treated as a unit vector with a given phase 0–2 π assigned to the ISI of interest. The vector strength ranges between 0.0–1.0. A value of 0.0 indicates that spike occurrence is entirely independent of the signal periodicity, whereas a value of 1.0 indicates that all spikes occur at exactly the same phase as the signal. The significance of the vector strength was assessed using the Rayleigh Z-test  at the 5% significance level (Fig. 1A–2, dotted line). At a given SPL, Z value was filtered by a weighted average with its 5 neighbors in the ratio of 1:2:3:2:1 and was plotted against ISI (hereafter, "temporal modulation transfer function [tMTF]"). By comparing the tMTF obtained at different SPLs, we defined the value of the SPL and ISI to produce the maximum Z value as the "best SPL" and "best ISI" (arrow), respectively. Note, the following data were obtained at the best SPL, unless otherwise specified. If the maximum Z value of a given neuron exceed the significance level, we evaluated stimulus-locking responses taking place (solid bar) and classified the neuron as "synchronization type"; if not, classified as "non-synchronization type."
The second measure was based on the rate of firing activities. At the beginning, the driven rate at each ISI was calculated by subtracting the mean of spontaneous firing rate from firing rate for 50–500 ms after the onset of the click trains. Then, the driven rate was filtered by a weighted average with its 5 neighbors in the ratio of 1:2:3:2:1 and was plotted against ISI (hereafter, "rate modulation transfer function [rMTF]"; Fig. 1A–3). When this measure exceeded the threshold for excitation (= 2*SD of the mean of spontaneous firing rate; dotted line), we evaluated rate response taking place (solid bar).
To estimate the recovery time course of neural responsivity, we presented a single click (= probe stimulus) at 1.0, 1.8 or 3.6 s after the termination of the click trains and measured the number of evoked spikes (#spikes) (Fig. 4D). For simplicity, we set ISI of the click trains at either best ISI or 0.9 multiples of the γ-δ border (see Results).
To assess the involvement of suppressory processes in the stimulus-locking responses, we conducted a single-click paradigm: a single rectangular click (1-ms duration, at the best SPL) was presented 10–20 times at intervals of 4.0–5.5 s. The spike occurrence was examined with dot raster plots (Fig. 6A–1) and post-latency time histogram (Fig. 6B–2). For constructing the latter, we first obtained peri-stimulus time histogram (bin width = 2 ms). We defined the response latency as the beginning of three consecutive bins in which the firing rate exceeded the threshold for excitation (see above). Then, we obtained a post-latency time histogram for 500 ms after the latency (bin width = 10 ms). The threshold for "suppression" was set at the mean minus 2*SD of the spontaneous firing rate (dotted line, bottom). Note that this paradigm was executed only for part of the AI neurons (n = 35) that showed an appreciable amount of spontaneous firing rate to make the threshold for suppression > 0 (spikes/s).
Physiological data was presented as the mean ± SD. If necessary, #spikes was normalized to the value elicited by the 1st click in the train that was delivered at the best ISI (= control level; e.g., Fig. 4C). In general, we employed Student's t-test or a one-way repeated Analysis of variance (ANOVA) followed by (post-hoc) Fisher protected least-significant difference test (PLSD) for pairwise comparisons. The significance level was set at P < 0.05 against a null hypothesis of equal performance.
Minimal cortical models
where C m is the membrane capacitance (1 μ F/cm2). The right-hand side incorporates the intrinsic ionic currents (I ion ), synaptic currents (I syn ) and external input (I app ). In addition, the model includes noise current that comes from a presynaptic neuron (noise; λ = 500 Hz). The presynaptic neuron fires randomly with a uniform distribution in time. The usual method of integration was a fourth-order Runge-Kutta method with a time step of dt = 0.02.
where the values of the parameters were derived from the experimental data and other models [63, 64] (same as I k and I leak ) and are as follows: Na (maximal conductance of the sodium channel) = 120 mS/cm2, p = 3, q = 1, E Na = 50 mV.
where k = 36 mS/cm2, j = 4, E k = - 77 mV.
where g L = 0.3 mS/cm2 and E L = - 54.4 mV.
where k sn = 1 ms-1, τ NMDA = 120 ms, k xn = 1 ms-1 and ~ τ NMDA = 14 ms.
where represents the normalized G-protein concentration in the activated form, K d refers to the dissociation constant of G-protein binding to K+ channels, r GABAB is the fraction of the GABAB receptors in the activated form. GABA concentration in the synaptic cleft, [T GABAB ], was set to 1 mM for 1 ms when the membrane potential crossed zero (only rising phase). Fitting this model to whole-cell recorded GABAB currents  gave the following values: K d = 100 μ M4, K1 = 180 s-1, K2 = 34 s-1, K3 = 9 × 104 M-1s-1, K4 = 1.2 s-1 with n = 4 binding sites.
The external input (I app ) was modeled as 1-ms-long repetitive square-wave pulses (initial amplitude = 12.5 μ Acm-2) with or without frequency-dependent depression which results mainly from temporal exhaustion of readily releasable neurotransmitter pool at presynaptic terminals .
where λ d and λ r are the depletion and refilling time constant, respectively, while N c represents the maximum size of the readily releasable vesicles, and n d represents the number of stimuli required to deplete the vesicles. Parameter values are set such that λ r = 0.05, N c = 15 and n d = 14. These values were chosen from the range given in previous reports [75, 76] and were selected to fit to the present data that most probably reflects frequency-dependent depression (Fig. 4C, open circles). Since synaptic potency, which is defined as the average size of the synaptic response when transmitter release does occur, has been proven not to considerably vary despite the changes in release probability [75, 77], we considered that I app is linearly related to P(N).
here, int = 1000/f r . The best fit of this kinetic scheme to the present data (Fig. 5D, filled circle) gave C1 = 0.2,C2 = 1.3 and τ1 = 20 ms.
primary auditory cortex
analysis of variance
- CV syllable:
excitatory post-synaptic potential
- E-S sequence:
inhibitory post-synaptic potential
medial geniculate body
- PLSD test:
protected least-significant difference test
rate modulation transfer function
the number of spikes
sound pressure level
temporal modulation transfer function
voice onset time
We thank N. Yaguchi for technical assistance. The work was supported by a grant from the Ministry of Education, Science, Culture, Sports and Technology, Japan.
- Warren RM, Bashford JA: Perception of acoustic iterance: pitch and infrapitch. Percept Psychophys. 1981, 29: 395-402.View ArticlePubMedGoogle Scholar
- Hirsh IJ: Auditory perception of temporal order. J Acoust Soc Am. 1959, 31: 759-767. 10.1121/1.1907782.View ArticleGoogle Scholar
- Steinbüchel vN, Pöppel E: Temporal order threshold and language perception. Frontiers in knowledge based computing. Edited by: Bhatkar VP, Rege KM. 1991, New Delhi, India: Narosa Publishing House, 81-90.Google Scholar
- Fraisse P: Perception and estimation of time. Annu Rev Psychol. 1984, 35: 1-36. 10.1146/annurev.ps.35.020184.000245.View ArticlePubMedGoogle Scholar
- Lisker L, Abramson AS: A cross-language study of voicing in initial stops: acoustical measurements. Word. 1964, 20: 384-422.Google Scholar
- Stevens KN: Acoustic Phonetics. 1998, Cambridge, UK: MIT PressGoogle Scholar
- Kuhl PK, Padden DM: Enhanced discriminability at the phonetic boundaries for the voicing feature in macaques. Percept Psychophys. 1982, 32 (6): 542-550.View ArticlePubMedGoogle Scholar
- Kuhl PK, Miller JD: Speech perception by the chinchilla: identification function for synthetic VOT stimuli. J Acoust Soc Am. 1978, 63: 905-917. 10.1121/1.381770.View ArticlePubMedGoogle Scholar
- Dooling RJ, Okanoya K, Brown SD: Speech perception by budgerigars (Melopsittacus undulatus): the voiced-voiceless distinction. Percept Psychophys. 1989, 46: 65-71.View ArticlePubMedGoogle Scholar
- Albert ML, Bear D: Time to understand. A case study of word deafness with reference to the role of time in auditory comprehension. Brain. 1974, 97: 373-384. 10.1093/brain/97.1.373.View ArticlePubMedGoogle Scholar
- Auerbach SH, Allard T, Naeser M, Alexander MP, Albert ML: Pure word deafness. Analysis of a case with bilateral lesions and a defect at the prephonemic level. Brain. 1982, 105: 271-300. 10.1093/brain/105.2.271.View ArticlePubMedGoogle Scholar
- Miceli G: The processing of speech sounds in a patient with cortical auditory disorder. Neuropsychologia. 1982, 20: 5-20. 10.1016/0028-3932(82)90083-5.View ArticlePubMedGoogle Scholar
- Phillips DP, Farmer ME: Acquired word deafness, and the temporal grain of sound representation in the primary auditory cortex. Behav Brain Res. 1990, 40: 85-94. 10.1016/0166-4328(90)90001-U.View ArticlePubMedGoogle Scholar
- Lu T, Liang L, Wang X: Temporal and rate representations of time-varying signals in the auditory cortex of awake primates. Nat Neurosci. 2001, 4: 1131-1138. 10.1038/nn737.View ArticlePubMedGoogle Scholar
- Eggermont JJ: Neural correlates of gap detection and auditory fusion in cat auditory cortex. NeuroReport. 1995, 6: 1645-1648.View ArticlePubMedGoogle Scholar
- Amitai Y: Thalamocortical synaptic connections: efficacy, modulation, inhibition and plasticity. Rev Neurosci. 2001, 12: 159-173.View ArticlePubMedGoogle Scholar
- Zucker RS, Regehr WG: Short-term synaptic plasticity. Annu Rev Physiol. 2002, 64: 355-405. 10.1146/annurev.physiol.64.092501.114547.View ArticlePubMedGoogle Scholar
- Galarreta M, Hestrin S: Frequency-dependent synaptic depression and the balance of excitation and inhibition in the neocortex. Nat Neurosci. 1998, 1: 587-594. 10.1038/2882.View ArticlePubMedGoogle Scholar
- Kamiya H, Zucker RS: Residual Ca2+ and short-term synaptic plasticity. Nature. 1994, 371: 603-606. 10.1038/371603a0.View ArticlePubMedGoogle Scholar
- Eggermont JJ: Temporal modulation transfer functions for AM and FM stimuli in cat auditory cortex. Effects of carrier type, modulating waveform and intensity. Hear Res. 1994, 74: 51-66. 10.1016/0378-5955(94)90175-9.View ArticlePubMedGoogle Scholar
- Bao S, Chang EF, Woods J, Merzenich MM: Temporal plasticity in the primary auditory cortex induced by operant perceptual learning. Nat Neurosci. 2004, 7: 974-981. 10.1038/nn1293.View ArticlePubMedGoogle Scholar
- Garabedian CE, Jones SR, Merzenich MM, Dale A, Moore CI: Band-pass response properties of rat SI neurons. J Neurophysiol. 2003, 90: 1379-13791. 10.1152/jn.01158.2002.View ArticlePubMedGoogle Scholar
- de Ribaupierre F, Goldstein MH, Yeni-Komshian G: Intracellular study of the cat's primary auditory cortex. Brain Res. 1972, 48: 185-204. 10.1016/0006-8993(72)90178-3.View ArticlePubMedGoogle Scholar
- Metherate R, Ashe JH: Facilitation of an NMDA receptor-mediated EPSP by paired-pulse stimulation in rat neocortex via depression of GABAergic IPSPs. J Physiol. 1994, 481: 331-348.PubMed CentralView ArticlePubMedGoogle Scholar
- Ojima H, Murakami K: Intracellular characterization of suppressive responses in supragranular pyramidal neurons of cat primary auditory cortex in vivo. Cereb Cortex. 2002, 12: 1079-1091. 10.1093/cercor/12.10.1079.View ArticlePubMedGoogle Scholar
- Wehr M, Zador AM: Synaptic mechanisms of forward suppression in rat auditory cortex. Neuron. 2005, 47: 437-445. 10.1016/j.neuron.2005.06.009.View ArticlePubMedGoogle Scholar
- Cox CL, Metherate R, Weinberger NM, Ashe JH: Synaptic potentials and effects of amino acid antagonists in the auditory cortex. Brain Res Bull. 1992, 28: 401-410. 10.1016/0361-9230(92)90039-Z.View ArticlePubMedGoogle Scholar
- Gaese BH, Ostwald J: Temporal coding of amplitude and frequency modulation in the rat auditory cortex. Eur J Neurosci. 1995, 7: 438-450. 10.1111/j.1460-9568.1995.tb00340.x.View ArticlePubMedGoogle Scholar
- Eggermont JJ: Representation of spectral and temporal sound features in three cortical fields of the cat. Similarities outweigh differences. J Neurophysiol. 1998, 80: 2743-2764.PubMedGoogle Scholar
- Bruder GE: Dichotic listening: new developments and applications in clinical research. Ann NY Acad Sci. 1991, 620: 217-232. 10.1111/j.1749-6632.1991.tb51586.x.View ArticlePubMedGoogle Scholar
- Licklider JC: A duplex theory of pitch perception. Experientia. 1951, 7: 128-134. 10.1007/BF02156143.View ArticlePubMedGoogle Scholar
- Flanagan JL, Guttman N: On the pitch of periodic pulses. J Acoust Soc Am. 1960, 32: 1308-1319. 10.1121/1.1907900.View ArticleGoogle Scholar
- Guttman N, Flanagan JL: Pitch of high-pass-filtered pulse trains. J Acoust Soc Am. 1964, 36: 757-765. 10.1121/1.1919063.View ArticleGoogle Scholar
- Steinschneider M, Reser DH, Fishman YI, Schroeder CE, Arezzo JC: Click train encoding in primary auditory cortex of the awake monkey: evidence for two mechanisms subserving pitch perception. J Acoust Soc Am. 1998, 104: 2935-2955. 10.1121/1.423877.View ArticlePubMedGoogle Scholar
- Maruyama N, Saitoh K, Kudoh M: Unit responses of the cat's auditory cortex to synthesized formant. Proceedings of the Japan Academy. 1979, 55 (B): 413-417.View ArticleGoogle Scholar
- Bieser A, Müller-Preuss P: Auditory responsive cortex in the squirrel monkey: neural responses to amplitude-modulated sounds. Exp Brain Res. 1996, 108: 273-284. 10.1007/BF00228100.View ArticlePubMedGoogle Scholar
- Miller LM, Escabi MA, Read HL, Schreiner CE: Functional convergence of response properties in the auditory thalamocortical system. Neuron. 2001, 32: 151-160. 10.1016/S0896-6273(01)00445-7.View ArticlePubMedGoogle Scholar
- Liang L, Lu T, Wang X: Neural representations of sinusoidal amplitude and frequency modulations in the primary auditory cortex of awake primates. J Neurophysiol. 2002, 87: 2237-2261.PubMedGoogle Scholar
- Sherman SM, Koch C: The control of retinogeniculate transmission in the mammalian lateral geniculate nucleus. Exp Brain Res. 1986, 63: 1-20. 10.1007/BF00235642.View ArticlePubMedGoogle Scholar
- Eggermont JJ: Neural correlates of gap detection in three auditory cortical fields in the cat. J Neurophysiol. 1999, 81: 2570-2581.PubMedGoogle Scholar
- Grothe B: New roles for synaptic inhibition in sound localization. Nat Rev Neurosci. 2003, 4: 540-550. 10.1038/nrn1136.View ArticlePubMedGoogle Scholar
- Eggermont JJ: Neural Responses in Primary Auditory Cortex Mimic Psychophysical, Across-Frequency-Channel, Gap-Detection Thresholds. J Neurophysiol. 2000, 84: 1453-1463.PubMedGoogle Scholar
- Bartlett EL, Wang X: Long-lasting modulation by stimulus context in primate auditory cortex. J Neurophysiol. 2005, 94: 83-104. 10.1152/jn.01124.2004.View ArticlePubMedGoogle Scholar
- Brosch M, Schreiner CE: Time course of forward masking tuning curves in cat primary auditory cortex. J Neurophysiol. 1997, 77: 923-943.PubMedGoogle Scholar
- Reale RA, Brugge JF: Directional sensitivity of neurons in the primary auditory (AI) cortex of the cat to successive sounds ordered in time and space. J Neurophysiol. 2000, 84: 435-450.PubMedGoogle Scholar
- Nakamoto KT, Zhang J, Kitzes LM: Temporal nonlinearity during recovery from sequential inhibition by neurons in the cat primary auditory cortex. J Neurophysiol. 2006, 95: 1897-1907. 10.1152/jn.00625.2005.View ArticlePubMedGoogle Scholar
- Zurita P, Villa AE, de Ribaupierre Y, de Ribaupierre F, Rouiller EM: Changes of single unit activity in the cat's auditory thalamus and cortex associated to different anesthetic conditions. Neurosci Res. 1994, 19: 303-316. 10.1016/0168-0102(94)90043-4.View ArticlePubMedGoogle Scholar
- Franks NP, Lieb WR: Molecular and cellular mechanisms of general anaesthesia. Nature. 1994, 367: 607-614. 10.1038/367607a0.View ArticlePubMedGoogle Scholar
- Müller-Preuss P: On the mechanisms of call coding through auditory neurons in the squirrel monkey. Eur Arch Psychiatry Neurol Sci. 1986, 236: 50-55. 10.1007/BF00641059.View ArticlePubMedGoogle Scholar
- Schreiner CE, Urbas JV: Representation of amplitude modulation in the auditory cortex of the cat. II. Comparison between cortical fields. Hear Res. 1988, 32: 49-63. 10.1016/0378-5955(88)90146-3.View ArticlePubMedGoogle Scholar
- Creutzfeldt OD, Hellweg FC, Schreiner CE: Thalamocortical transformation of responses to complex auditory stimuli. Exp Brain Res. 1980, 39: 87-104. 10.1007/BF00237072.View ArticlePubMedGoogle Scholar
- Langner G: Periodicity coding in the auditory system. Hear Res. 1992, 60: 115-142. 10.1016/0378-5955(92)90015-F.View ArticlePubMedGoogle Scholar
- Joris PX, Schreiner CE, Rees A: Neural processing of amplitude-modulated sounds. Physiol Rev. 2004, 84: 541-577. 10.1152/physrev.00029.2003.View ArticlePubMedGoogle Scholar
- Chimoto S, Kitama T, Qin L, Sakayori S, Sato Y: Tonal response features of primary auditory cortex neurons in awake cats. Brain Res. 2002, 934: 34-42. 10.1016/S0006-8993(02)02316-8.View ArticlePubMedGoogle Scholar
- Qin L, Chimoto S, Sakai M, Wang J, Sato Y: Comparison between offset and onset responses of primary auditory cortex on-off neurons in awake cats. J Neurophysiol. 2007, 97: 3421-3431. 10.1152/jn.00184.2007.View ArticlePubMedGoogle Scholar
- Qin L, Sakai M, Chimoto S, Sato Y: Interaction of excitatory and inhibitory frequency-receptive fields in determining fundamental frequency sensitivity of primary auditory cortex neurons in awake cats. Cereb Cortex. 2005, 15: 1371-1383. 10.1093/cercor/bhi019.View ArticlePubMedGoogle Scholar
- Sakai M, Chimoto S, Qin L, Sato Y: Periodicity coding of primary auditory cortex neurons in awake cats: in relation to tonal response patterns. Soc Neurosci Abstr. 2005, 35: 621.1.Google Scholar
- Schulze H, Langner G: Periodicity coding in the primary auditory cortex of the Mongolian gerbil (Meriones unguiculatus): two different coding strategies for pitch and rhythm?. J Comp Physiol [A]. 1997, 181: 651-663. 10.1007/s003590050147.View ArticleGoogle Scholar
- Fishman YI, Volkov IO, Noh MD, Garell PC, Bakken H, Arezzo JC, Howard MA, Steinschneider M: Consonance and dissonance of musical chords: neural correlates in auditory cortex of monkeys and humans. J Neurophysiol. 2001, 86: 2761-2788.PubMedGoogle Scholar
- Goldberg JM, Brown PB: Response of binaural neurons of dog superior olivary complex to dichotic tonal stimuli: some physiological mechanisms of sound localization. J Neurophysiol. 1969, 32 (4): 613-636.PubMedGoogle Scholar
- Mardia KV, Jupp PE: Directional Statistics. 2000, New York, WileyGoogle Scholar
- Hodgkin A, Huxley A: A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol. 1952, 117: 500-544.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones SR, Pinto DJ, Kaper TJ, Kopell N: Alpha-frequency rhythms desynchronize over long cortical distances: a modeling study. J Comp Neurosci. 2000, 9: 271-291. 10.1023/A:1026539805445.View ArticleGoogle Scholar
- Kyriazi HT, Simons DJ: Thalamocortical response transformations in simulated whisker barrels. J Neurosci. 1993, 13: 1601-1615.PubMedGoogle Scholar
- Destexhe A, Mainen ZF, Sejnowski TJ: Kinetic models of synaptic transmission. Methods in neuronal modeling. Edited by: Koch C, Segev I. 1998, Cambridge, MA: MIT Press, 1-25.Google Scholar
- Stern P, Edwards FA, Sakmann B: Fast and slow components of unitary EPSCs on stellate cells elicited by focal stimulation in slices of rat visual cortex. J Physiol. 1992, 449: 247-278.PubMed CentralView ArticlePubMedGoogle Scholar
- Golomb D, Amitai Y: Propagating neuronal discharges in neocortical slices: computational and experimental study. J Neurophysiol. 1997, 78: 1199-1211.PubMedGoogle Scholar
- Wang XJ, Rinzel J: Spindle rhythmicity in the reticularis thalami nucleus: synchronization among mutually inhibitory neurons. Neuroscience. 1993, 53: 899-904. 10.1016/0306-4522(93)90474-T.View ArticlePubMedGoogle Scholar
- Jonas P, Spruston N: Mechanisms shaping glutamate-mediated excitatory postsynaptic currents in the CNS. Curr Opin Neurobiol. 1994, 4: 366-372. 10.1016/0959-4388(94)90098-1.View ArticlePubMedGoogle Scholar
- Golomb D, Shedmi A, Curtu R, Ermentrout GB: Persistent synchronized bursting activity in cortical tissues with low magnesium concentration: a modeling study. J Neurophysiol. 2006, 95: 1049-1067. 10.1152/jn.00932.2005.View ArticlePubMedGoogle Scholar
- Wang XJ, Buzsaki G: Gamma oscillation by synaptic inhibition in a hippocampal interneuronal network model. J Neurosci. 1996, 16: 6402-6413.PubMedGoogle Scholar
- Hansel D, Mato G: Asynchronous states and the emergence of synchrony in large networks of interacting excitatory and inhibitory neurons. Neural Comput. 2003, 15: 1-56. 10.1162/089976603321043685.View ArticlePubMedGoogle Scholar
- Destexhe A, Sejnowski TJ: G-protein activation kinetics and spill-over of GABA may account for differences between inhibitory responses in the hippocampus and thalamus. Proc Natl Acad Sci USA. 1995, 92: 9515-9519. 10.1073/pnas.92.21.9515.PubMed CentralView ArticlePubMedGoogle Scholar
- Otis TS, Mody I: Differential activation of GABAA and GABAB receptors by spontaneously released transmitter. J Neurophysiol. 1992, 67: 227-235.PubMedGoogle Scholar
- Dobrunz LE, Stevens CF: Heterogeneity of release probability, facilitation, and depletion at central synapses. Neuron. 1997, 18: 995-1008. 10.1016/S0896-6273(00)80338-4.View ArticlePubMedGoogle Scholar
- Murthy VN, Sejnowski TJ, Stevens CF: Heterogeneous release properties of visualized individual hippocampal synapses. Neuron. 1997, 18: 599-612. 10.1016/S0896-6273(00)80301-3.View ArticlePubMedGoogle Scholar
- Liang SL, Carlson GC, Coulter DA: Dynamic regulation of synaptic GABA release by the glutamate-glutamine cycle in hippocampal area CA1. J Neurosci. 2006, 26: 8537-8548. 10.1523/JNEUROSCI.0329-06.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Kretz R, Shapiro E, Kandel ER: Post-tetanic potentiation at an identified synapse in Aplysia is correlated with a Ca2+-activated K+ current in the presynaptic neuron: evidence for Ca2+ accumulation. Proc Natl Acad Sci USA. 1982, 79: 5430-5434. 10.1073/pnas.79.17.5430.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.