|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Laboratory of Auditory Neurophysiology, Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland
Submitted 5 June 2006; accepted in final form 17 October 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The way in which temporally modulated sound signals are encoded in the auditory thalamus serves as an interesting test case for investigating a coding transformation through the thalamus. Numerous studies of the inferior colliculus (IC), which is the auditory input to the thalamus, have shown that IC neurons represent time-varying sounds such as amplitude-modulated (AM) tones or noise by stimulus-synchronized discharges that are phase locked to the modulation envelope (Batra et al. 1989
; Krishna and Semple 2000
; Langner and Schreiner 1988
; Muller-Preuss et al. 1994
; for review see Joris et al. 2004
), including one study of multiunit responses to AM stimuli in awake squirrel monkeys (Muller-Preuss et al. 1994
). Although the firing rates of IC neurons also change in a modulation-frequencydependent manner, the range of modulation frequencies over which IC neurons respond vigorously is typically the same range over which they are well synchronized with the stimulus (Joris et al. 2004
; Krishna and Semple 2000
). This is not the case in the auditory cortex of awake marmoset monkeys, where there are two main response types. Synchronized responses are phase locked to the stimulus, similar to reported IC responses (Liang et al. 2002
; Lu et al. 2001a
). Cortical synchronized responses occur mainly for low repetition rates (typically <50 Hz or 20-ms interclick intervals). The other main response type represents decreasing interclick intervals (increasing click rates) with monotonic increases in firing rate (Lu et al. 2001a
). Thus temporal and rate representations are often dissociated in the auditory cortex but only rarely in the IC, although no studies have studied temporal processing in marmoset IC neurons. This transformation in response properties occurs either at the ICmedial geniculate body (MGB) synapse, the thalamocortical synapse, or intracortically. If nonsynchronized rate representations are prevalent in MGB neuronal responses, this would indicate that there is a conversion from a temporal to a rate code at the ICMGB synapse.
Data from brain slice studies in the rat MGB have demonstrated two response patterns to repetitive electrical stimuli that are reminiscent of the synchronized and nonsynchronized responses observed in the auditory cortex (Bartlett and Smith 2002
). One type produces large excitatory postsynaptic potentials (EPSPs) in MGB neurons and shows strong paired-pulse depression. At low stimulation frequencies, the synaptic responses produce highly entrained spiking in MGB neurons, similar to the synchronized responses observed in the auditory cortex in vivo. The other type produces small EPSPs and often shows paired-pulse facilitation. They generally respond poorly at low stimulation frequencies, instead requiring rapid stimulation frequencies for the EPSPs to summate and produce a suprathreshold response (Bartlett and Smith 2002
), similar to the nonsynchronized responses in cortex. Thus the synaptic physiology of MGB neurons demonstrates that many MGB neurons are capable of producing nonsynchronized responses in vivo, at least in rodents.
Using repetitive click stimuli, we investigated the temporal processing capabilities of MGB neurons in the awake marmoset monkey. We found a significant alteration in MGB response properties compared with reported IC responses and compared with auditory cortex properties obtained from neurons recorded from the same preparation.
| METHODS |
|---|
|
|
|---|
All experiments were performed at Johns Hopkins University in AAALAC approved facilities following protocols approved by the Institutional Animal Care and Use Committee (IACUC). The methods for preparing marmoset monkeys for electrophysiological recording and carrying out single-unit recordings in awake marmosets were previously reported (Bartlett and Wang 2005
; Lu et al. 2001a
,b
). Some modifications were necessary for recording from MGB neurons and are described here.
The recording electrode approached the MGB in a dorsolateral to ventromedial trajectory. The electrode entered approximately 3 mm lateral to the lateral sulcus near its caudal end at an angle of 57.5°, where 0° is horizontal and 90° is vertical. The penetrations were at an angle 0° relative to the coronal plane (i.e., no change in rostrocaudal plane as the electrode was advanced). In animals whose auditory cortices had been mapped for different experiments, the insertion site for MGB was known to be lateral to the high-frequency representation of A1. Auditory-responsive neurons were encountered approximately 69 mm after penetrating the cortical surface (usually 6.57.0 mm). The rostral half of the MGB is bordered laterally by the lateral geniculate nucleus. Therefore it was possible to determine that the MGB was nearby when one passed through a region of visually responsive neurons followed by 200700 µm without encountering any neurons.
The arousal state of the animal was monitored during experiments. Stimulus presentation and recording were halted or paused if the recording quality deteriorated or the animal was suspected to be drowsy when monitored through a camera placed inside the recording chamber (trials were resumed after the animal was confirmed to be awake with eyes open). Typically, one to three single units were fully studied in each recording session that lasted 35 h.
Acoustic stimuli
Stimuli were presented in free field from a speaker (B&W 601) located 70 cm in front of the animal and level with the animal's head (0° azimuth, 0° elevation). The interiors of the double-walled soundproof chambers (Industrial Acoustics) in which recordings took place were covered with 3-in. acoustic absorption foam (Sonex, Ilbruck, MN). Stimuli were created from custom programs using MATLAB software (The MathWorks, Natick, MA), generated through a D/A converter (Tucker-Davis Technologies, Alachua, FL) with a 100-kHz sampling rate, low-pass filtered at 50 kHz, attenuated with two serially linked attenuators (Tucker-Davis Technologies), and amplified by a power amplifier (Crown International, Elkhart, IN). The measured speaker output was within 4 dB of the intended output for frequencies from 100 Hz to 36 kHz, which encompasses the hearing range of marmosets, with a calibrated sound level of 90 dB SPL at 0 dB attenuation for 1-kHz tones.
The sound stimuli used were either 0.1-ms rectangular clicks or Gaussian clicks (Fig. 1, A and B). Gaussian clicks are tones, the amplitude envelope of which is a Gaussian function whose duration and bandwidth are determined by a variance parameter sigma (
). Gaussian clicks were used to produce a spectral peak at the neuron's characteristic frequency (CF) and to avoid the activation of inhibitory sidebands. Examples of the individual click waveforms and power spectra are shown for 1-s train Gaussian clicks presented at 40-ms interclick intervals (25-Hz click rate). Each Gaussian click in this example had a carrier frequency of 5.0 kHz. Gaussian clicks were previously shown to be more effective than rectangular clicks at evoking responses in most auditory cortical neurons (Lu et al. 2001a
). Gaussian clicks were set to the best frequency (BF) of the recorded neuron. Smaller sigma values produce shorter-duration clicks that have wider bandwidths, whereas the larger sigma values produce clicks with longer durations and narrower bandwidths (Fig. 1, A and B). To be consistent with the language of Lu et al. 2001a
, temporal modulation will be discussed mostly in terms of interclick intervals (ICIs) rather than click rates (1/ICI).
|
500 ms after the offset of the sound stimulus. There was at least 1 s of silence between each stimulus. Each click train stimulus was repeated five to ten times. Data analysis
MGB neurons (n = 126) from three animals were tested with click stimuli. Of those, 97 neurons produced significant responses, either as statistically significant increases in firing rate or as significant synchrony, and were analyzed further. Firing rate was computed for the entire stimulus duration. Spontaneous rate was computed as the mean rate of the 500-ms period that preceded each trial. The ability to synchronize to a click train was quantified by measuring the vector strength of the neuronal response at each ICI. Statistical significance was assessed using the Rayleigh statistic [2 x number of spikes x (vector strength)2], which also takes into account the number of spikes evoked by the stimulus (Lu et al. 2001a
; Mardia and Jupp 2000
). For each ICI tested, the Rayleigh statistic was computed for the first half and second half of the stimulus and the minimum value was used. This excluded some responses that were synchronized during only part of the stimulus. A threshold Rayleigh statistic value of 13.8 was considered significant (P < 0.001).
For each unit, the minimum interclick interval at which the neuron significantly synchronized with the stimulus, called the synchronization boundary, was computed. This is the lower bound at which there will be significant synchronization to the click train. A linear interpolation was made between the Rayleigh value at the minimum ICI for which the unit was synchronized and the Rayleigh value for the next shorter ICI at which the unit was nonsynchronized (Lu et al. 2001a
). The synchronization boundary was the ICI for which the Rayleigh value was 13.8. The maximum click rate at which the neuron was significantly synchronized was called Fmax, given by Fmax = 1/(synchronization boundary).
Many units produced nonsynchronized responses for short ICI stimuli. A rate boundary was computed for neurons that had statistically significantly driven firing rates (rank-sum test, P < 0.05 vs. spontaneous firing rate) in response to stimuli with ICI
5 ms. This is the upper boundary at which significant increases in firing rates are observed. A linear interpolation was made between the firing rate at the maximum ICI for which the unit was significantly driven and the firing rate for the next longer ICI at which the unit was not significantly driven (Lu et al. 2001a
). The rate boundary was the ICI for which the interpolated driven firing rate crossed the threshold of twice the SD of the spontaneous rate.
Anatomy
At the cessation of recording, electrolytic lesions were made in physiologically identified regions of MGB by passing 2- to 10-µA current through the recording electrode (610 s each polarity). Animals were euthanized by administering an initial intramuscular injection of ketamine followed by an intraperitoneal injection of Euthasol (100 mg/kg; Virbac Australia).
Animals were transcardially perfused with room-temperature PBS with heparinized phosphate buffer (pH
7.0) followed by 4% paraformaldehyde (EM Grade, Ted Pella, Redding, CA) in 0.1 M phosphate buffer (pH
7.0). After perfusion, the brain was removed, put in 30% sucrose solution, and frozen. Sections (30 µm) were cut on a freezing microtome through the entire extent of the MGB. Sections were processed for Nissl staining, parvalbumin immunoreactivity, and calbindin immunoreactivity using standard histological procedures (Jones and Hendry 1989
). Myelin staining used a modified Gallyas stain protocol (Pistorio et al. 2005).
MGB subdivisions were assigned based on parcellation schemes established in previous studies of the MGB in macaques (Hashikawa et al. 1995
; Jones 2003
; Kaas and Hackett 2000
; Molinari et al. 1995
), owl monkeys (Morel and Kaas 1992
), and marmosets (Aitkin et al. 1988
, 1993
). The locations of recorded neurons in the MGB were reconstructed based on the coordinates of the tracks relative to the tracks in which lesions were made and the depths at which the recordings occurred. All analyzed units were recorded from the dorsal (anterodorsal and posterodorsal) and ventral divisions of the MGB (Jones 2003
). A small number (n = 6) of units not included in the analysis were found in the medial division of MGB and the suprageniculate nucleus.
| RESULTS |
|---|
|
|
|---|
Response types
Response patterns to click stimuli typically fell into one of three main classes, shown in Fig. 2. Ninety-seven responses were analyzed from 87 units, with 10 units producing two different response types depending on the click sigma used (described later in detail). Synchronized (Synch) responses were observed in 16 units and consisted of discharges that were well entrained to the click train for long ICIs. At short ICIs, these neurons did not respond and were often inhibited after the onset response (14/16 units). The neuron in Fig. 2A produced highly synchronized discharges in response to rectangular click trains for ICIs from 10 to 150 ms, with the firing rate steadily decreasing from its peak at 10-ms ICI (100 Hz) as the ICI increased. This was a result of the decreased number of clicks in the stimulus with increasing ICI. Shorter ICIs from 4 to 7.5 ms resulted in a reduction of rate (Fig. 2A, solid line) and synchrony (Fig. 2A, dashed line) until the neuron produced an onset response only for 1- to 3-ms ICIs.
|
Neurons with Nonsynch responses were typically not driven at long ICIs and mainly generated high discharge rates in response to ICIs
10 ms. The neuron in Fig. 2B was significantly driven only for ICI
5 ms. In this example, firing rate was nonmonotonic with respect to ICI and peaked for 3-ms ICI stimuli. This differs from what was observed in the nonsynchronized responses of auditory cortex neurons, whose rates almost always increased monotonically with decreasing ICI (Lu et al. 2001a
). Of 41 neurons with Nonsynch responses, 14 were strongly nonmonotonic in their rate response, similar to the neuron in Fig. 2B (>50% change in firing rate between peak rate and rate in response to shortest ICI), 22/41 neurons were at least weakly nonmonotonic (>25% change in firing rate), and the firing rates of 17/41 neurons increased monotonically with decreasing ICI. Six nonsynchronized neurons were atypical in that they increased their firing rates with increasing ICI.
In contrast to purely synchronized or nonsynchronized responses, 40 neurons displayed responses that included two response regimes. These Mixed (Mixed) responses consisted of a stimulus-synchronized response at long ICIs and a nonsynchronized rate response at short ICIs. These response regimes were often separated by a range of ICIs over which the neuron was inhibited (20/40 units). In Fig. 2C, the neuron generated stimulus-synchronized responses and a vigorous firing rate for ICI from 7.5 to 100 ms. For ICIs from 3 to 5 ms, the neuron was clearly inhibited. Strong sustained firing occurred for ICIs from 1 to 2 ms. This response was nonsynchronized and the maximal firing rate exceeded that observed for longer ICIs when the response was synchronized. It should be noted that for the click sigma used for the stimulus (sigma = 0.1), the click waveforms were nonoverlapping. This meant that the depth of AM was 100%, even for a 1-ms ICI. The presence and prevalence of Mixed responses suggests that the intrinsic and synaptic mechanisms necessary to generate synchronized and nonsynchronized responses are present in single MGB neurons rather than forming separate neuronal populations.
Population averages of the firing rate, vector strength, and proportion of units synchronized are plotted in Fig. 3 for each response type as a function of ICI. Synch and Mixed responses were on average driven at long ICIs, whereas Nonsynch responses were not (Fig. 3A). The average firing rates of Synch responses dropped dramatically for ICI <7.5 ms and were usually inhibited for ICI from 1 to 5 ms. In contrast, the firing rates of Mixed and Nonsynch responses were maximal or near maximal for short ICIs. In fact, Nonsynch firing rates did not begin to increase until ICIs were <50 ms and the Nonsynch rate response range was complementary to the Synch rate response range. Both Synch and Mixed firing rates were significantly higher than Nonsynch rates for all ICIs from 12.5 to 150 ms (P < 0.05, KruskalWallis test), whereas Nonsynch and Mixed responses were significantly higher than Synch responses for 1- to 5-ms ICIs (P < 0.05, KruskalWallis test). For the Synch and Mixed response populations, the vector strength of synchronized responses was consistently high for ICIs as low as 5 ms (Fig. 2B). The proportion of units that were stimulus synchronized (Fig. 2C) decreased with decreasing ICI and became very low for ICI <3 ms for both Synch and Mixed populations. Overall, 56 neurons produced stimulus-synchronized discharges. Almost half (26/56 = 46%) of the Synch and Mixed neurons were still phase locked at 5-ms ICI, but only (9/56 = 17%) were still phase locked at 3-ms ICI, suggesting that this was approaching the upper bound for phase locking in the MGB. Only two units were phase locked at 2-ms ICI.
|
5 ms (Fig. 4A, open circles). By contrast, nonsynchronized firing rates were typically maximal at ICIs
5 ms and minimal at ICIs >10 ms (Fig. 4A, black squares). Depending on whether the maximum firing rate occurred in the synchronized or nonsynchronized portions of the Mixed response, the ICIs were widely distributed (Fig. 4A, gray triangles).
|
The nonsynchronized responses of auditory cortex neurons generally increased monotonically with decreasing ICI (Lu et al. 2001a
). By contrast, the responses of 34% (14/41) of Nonsynch responses and 50% (20/40) of Mixed responses exhibited a clear preference for ICI, such that responses to 1-ms ICI click trains were
50% lower than the maximum response (Fig. 4B).
There was no significant difference in the temporal precision of Synch or Mixed responses as assessed by the mean vector strength averaged over all ICIs (Fig. 4C) or the maximum vector strength (KruskalWallis, P > 0.05). There was also no significant difference between the ICIs that produced the highest value of the Rayleigh statistic value, which is a measure of the maximal sustained synchronized response (KruskalWallis, P > 0.05). Only 14/56 (25%) responses produced the maximum Rayleigh for ICIs
100 Hz, with most neurons in the 20- to 80-Hz range.
The different types of responses were not associated with systematic differences in the BFs of the neurons. Moreover, the nonsynchronized responses were not attributable to responses to click rates at or near BF in almost all cases. Figure 4D plots the normalized distribution of BFs derived from tone or noise responses for neurons with synchronized responses (combined Synch and Mixed populations, open white bars) and nonsynchronized responses (black bars). Synchronized and nonsynchronized responses were found in units over a wide range of BFs and there was no significant difference between the two BF distributions (KruskalWallis, P > 0.05). The BFs of the vast majority of units were well outside of the range at which any of the click rates used would impinge on the unit's frequency receptive field (synchronized units 2575% of distribution = 2.111.3 kHz, nonsynchronized units 2575% of distribution = 3.210.5 kHz). Only 11 units had BFs <1.5 kHz and only five units had BFs <1.0 kHz.
Although the type of response evoked was not dependent on sigma, there was a difference in the synchronization boundaries between the largest sigma used (0.4) and the other sigmas used (rectangular, 0.1, 0.2). Figure 4E shows the distribution of click sigmas that generated Synch (n = 16, open bars), Nonsynch (n = 41, black bars), and Mixed (n = 40, gray bars) responses. Each type of response was observed at all sigmas, although there was a slight tendency for synchronized responses to be produced by larger click sigmas (8/16 synchronized responses were produced by click sigma = 0.4). Figure 4F plots the distributions of the synchronization boundaries at each sigma. There were no significant differences between the synchronization boundaries for units tested with rectangular clicks, sigma = 0.1, or sigma = 0.2 (KruskalWallis, P > 0.05). However, all three groups had lower synchronization boundaries than those of units tested with sigma = 0.4 clicks (KruskalWallis, P < 0.05). This makes sense because for sigma = 0.4, individual clicks fused and the click train envelope was basically flat for ICI <5 ms. Note that if a unit was tested with two sigmas, the sigma associated with the shortest synchronization boundary was used.
Comparison of MGB and cortical responses
For each MGB neuron, the synchronization boundary and the rate boundary were computed (see METHODS). If a unit was tested with two or more different click sigmas or sound levels, the response that had the shortest synchronization boundary was chosen for analysis. Figure 5A shows the distribution of synchronization boundaries for the Synch and Mixed response types compared with the synchronization boundaries of auditory cortex neurons in awake marmosets from Lu et al. (2001a)
. MGB neurons had much shorter synchronization boundaries than those of auditory cortex neurons. Most commonly the MGB synchronization boundary was between 4 and 7 ms (median = 5.2 ms, 2575% = 3.311.4 ms, n = 56 neurons). This was approximately 2 octaves higher in terms of frequency than the median value for auditory cortex neurons in awake marmosets (Fig. 5A, solid line, median = 21.3 ms, 2575% = 12.454.3 ms; Lu et al. 2001a
). MGB Synch responses tended to have higher synchronization boundaries than Mixed responses, with 14/16 Synch responses having boundaries
5-ms ICI, compared with 18/40 Mixed responses (Fig. 5A, inset, P < 0.01,
2 test).
|
3 ms, the MGB firing rate of nonsynchronized neurons was near maximal. Then it decreased rapidly from 3 to 10 ms. Over this same ICI range, the MGB population vector strength of the combined Synch and Mixed populations increased rapidly. By 10-ms ICI, the MGB representations have clearly segregated (Fig. 5C), whereas the cortex representations are not clearly segregated until 20- to 30-ms ICIs. One question that we had was whether the nonsynchronized responses were initiated and generated by the auditory cortex. If this were the case, we would expect that the latency of the cortical nonsynchronized responses would be shorter than the MGB latencies. For each unit, the latency of the response was calculated. On average, the latencies of Synch and Mixed responses were significantly shorter than synchronized responses in the auditory cortex (KruskalWallis test, P = 0.02) (Fig. 6A). The median MGB synchronized response latency was 15.5 ms (n = 46), whereas the median cortex latency was 20.5 ms (n = 36). Nonsynchronized responses in MGB and auditory cortex occurred at much longer latencies than synchronized responses in both MGB and auditory cortex neurons (Fig. 6B). The median MGB nonsynchronized response latency was 49.8 ms (n = 39) and the distribution of MGB nonsynchronized response latencies overlapped extensively with auditory cortex latencies (median = 46.0 ms, n = 50). Therefore we cannot determine from these data whether nonsynchronized responses are initiated by MGB or auditory cortex neurons.
|
Shifts in response pattern arising from changes in click width
A critical issue in temporal processing by single units is whether the observed temporal processing properties are invariant with respect to the carrier. In other words, are neurons more sensitive to the periodicity of the temporal envelope, regardless of the details of the modulation waveform such as spectrum, duration, and envelope shape? To investigate this issue, we tested 38 neurons with two or more different click sigma over the full range of ICIs. Note that this differs from our initial testing for click sigma preference, described in METHODS, for which we sampled only three to four different ICIs. As shown in Fig. 1, altering the click sigma from rectangular to 0.1 up to 0.4 substantially decreases the bandwidth of the click stimulus while simultaneously increasing the duration of individual clicks and increasing the duration of each individual click (see Fig. 1).
A few neurons exhibited dramatic differences in their response patterns when the click sigma changed. An extremely striking example of these differences is shown in Fig. 7, A and B. This example was used to illustrate high sensitivity to click sigma, but such high sensitivity was observed in only a small minority of cases. Responses in Fig. 7, A and B were evoked by Gaussian click trains with the same CF and sound level and covering the same ICI range. The only difference was that
= 0.4 in Fig. 7A and
= 0.1 in Fig. 7B. When
= 0.4, the neuron was well synchronized for all stimuli with ICIs from 15 to 100 ms and it fired only weakly in response to shorter ICI stimuli (Fig. 7A). When
= 0.1, the neuron did not respond at all until the ICI was 10 ms. Even then, it fired vigorously for ICIs only from 3 to 7.5 ms (Fig. 7B), so this response pattern was complementary to the response pattern observed when
= 0.4. This can be seen clearly in Fig. 7C. The long ICIs that produced sustained firing and well-synchronized responses in Fig. 7A did not drive the unit at all in Fig. 7B. Note also that even when the firing rates were similar for 7.5-ms ICI stimuli, the latency of the response was much longer in Fig. 7B than in Fig. 7A.
|
Twenty-seven neurons produced synchronized responses for at least one sigma. For each of these neurons, the synchronization boundaries were computed for the smallest and largest click sigma values tested. Of 27 neurons only nine changed their Fmax (1/synchronization boundary) by less than one octave. Together, the temporal coding of 13/38 (9 synchronized + 4 nonsynchronized = 34%) neurons was therefore insensitive to the click stimuli used. However, 18/27 of units changed their Fmax by more than one octave. Thus 25/38 neurons (18 synchronized + 7 nonsynchronized = 66%) were sensitive to the spectrotemporal parameters of the click stimuli and 10/38 (26%) neurons changed their representation of modulation rate from synchronized to nonsynchronized when click sigma changed.
Despite this sensitivity to the click stimuli in some neurons, the sensitivity was predictable. We compared the synchronization boundaries for the 27 neurons that were synchronized for at least one click sigma. If a neuron generated synchronized responses for short-duration, spectrally broad clicks (rectangular clicks or sigma = 0.1), then it would usually also produce synchronized responses for long-duration, spectrally narrow clicks (sigma = 0.2 or 0.4). However, the converse was not true. The neurons that were most sensitive to click sigma generally synchronized only to sigma = 0.2 or 0.4 because these clicks were long enough in duration to evoke a response.
These observations are shown graphically in Fig. 7D. When the shorter synchronization boundary was observed for the smallest click sigma (Fig. 7D, open circles below y = x line, n = 14), the synchronization boundary was generally very short and did not change much when the click sigma changed. When the shorter synchronization boundary was observed for the largest click sigma (Fig. 7D, gray triangles above y = x line, n = 13), the synchronization boundary was generally longer and was quite sensitive to the click sigma. In fact, seven of ten neurons that were synchronized at one sigma and nonsynchronized at another click sigma were synchronized only for the largest click sigma. Neurons whose boundaries were shorter for the smallest sigma (Fig. 7D, open circles) had significantly shorter boundaries at the smallest and largest sigma values than neurons whose shorter boundaries were for the largest sigma values (Fig. 7D, gray triangles, P < 0.05, KruskalWallis). Furthermore, neurons whose boundaries were shorter for the smallest sigma (Fig. 7D, open circles) had significantly smaller differences between synchronization boundaries tested at different sigmas (i.e., closer to the y = x line, P < 0.05, KruskalWallis). These results suggest that some MGB neurons are highly sensitive to the spectrotemporal parameters of the overall stimulus, not just the ICI. Those neurons that synchronize poorly or not at all to short-duration, spectrally broad clicks appear to be particularly sensitive to parameters other than ICI. Conversely, some neurons appear able to phase lock at short ICI regardless of the click stimulus parameters, particularly those that synchronize well to short-duration clicks. These neurons may be specialized for rapid temporal processing.
| DISCUSSION |
|---|
|
|
|---|
We have demonstrated that MGB responses to temporally modulated stimuli differ from reported IC and auditory cortex responses to repetitive auditory stimuli. This implies that the MGB does not act as a simple relay and that there are sequential transformations of temporal coding from IC to MGB to cortex. Unlike reported IC responses, both stimulus-synchronized temporal and nonsynchronized rate representations were observed in MGB responses to repetitive click trains. Many individual MGB neurons were able to encode low click rates with a temporal representation and high click rates with a rate representation. With the appearance of nonsynchronized firing rate responses to encode high click rates, the MGB is the first auditory nucleus for which there is a complete dissociation of high synchrony and high firing rate. Some MGB neurons dramatically changed their discharge patterns when the click width changed, implying that the temporal processing capabilities of some MGB neurons depend on stimulus parameters other than the temporal modulation rate.
Comparison with previous studies of inferior colliculus
The representations of temporally modulated (sAM) stimuli were previously investigated extensively in the cat (Escabi and Schreiner 2002
; Langner and Schreiner 1988
), rat (Rees and Møller 1983
,1987
), bat (Burger and Pollak 1998
; Condon et al. 1994
), gerbil (Krishna and Semple 2000
), rabbit (Batra et al. 1989
), and squirrel monkey (Muller-Preuss et al. 1994
). Of these studies, only Batra et al. (1989
; single-unit) and Muller-Preuss et al. (1994
; multiunit) recorded in unanesthetized animals. IC neurons typically have band-pass or weakly modulated rate responses and band-pass or low-pass synchrony responses. Responsive IC neurons were always synchronized for some range of modulation frequencies, with synchronization boundaries in the range of 216 ms (Joris et al. 2004
). Thus the distribution of MGB synchronization boundaries in the present study is similar to what has been reported in IC (see Fig. 5A). The main difference between reported IC responses and the MGB responses in our study appears to be in the abundance of MGB neurons that did not respond significantly to low click rates and responded to high click rates only with nonsynchronized discharges. However, it is possible that some IC neurons are not phase locked to the stimulus and participate in the generation of nonsynchronized MGB responses (Escabi and Schreiner 2002
; Muller-Preuss et al. 1994
). Another possibility is that the abundance of nonsynchronized responses reflects a specialization found in primates but not in other species, although this remains to be tested.
Comparison with previous studies of MGB
Previous studies investigating the temporal processing abilities of MGB neurons observed primarily phase-locked responses (Creutzfeldt et al. 1980
; Miller et al. 2002
; Preuss and Muller-Preuss 1990
; Rouiller and deRibaupierre 1982
; Rouiller et al. 1981
; Vernier and Galambos 1957
), although neurons sometimes produced onset responses for high click rates. Our observation that the median Fmax of MGB neurons (192 Hz) is significantly higher than auditory cortex (47 Hz) confirms at the population level what was suggested by examples recorded from the awake guinea pig (Creutzfeldt et al. 1980
).
However, none of the previous studies reported the large number of nonsynchronized responses that we have observed. In some cases (Miller et al. 2002
; Rouiller et al. 1981
), this may have been the result of anesthetics that reduced sustained excitation produced by activation of N-methyl-D-aspartate (NMDA) receptors (Krystal et al. 2003
; Villars et al. 2004
). In studies using unanesthetized animals, either the range of repetition rates tested was below the range of most nonsynchronized responses (Preuss and Muller-Preuss 1990
) or the emphasis of the study was on synchronized responses (Creutzfeldt et al. 1980
; Vernier and Galambos 1957
).
It should be noted that one multiunit study in awake macaques found responses to clicks at very high rates (
800 Hz) in the putative thalamorecipient lamina in auditory cortex (Steinschneider et al. 1998). However, these responses occurred at sites whose BFs were matched to the preferred click rates, which differs from the nonsynchronized responses that we have reported. For the most part, the click rates that evoked nonsynchronized responses in our study were well outside of the frequency receptive field of the neurons (Fig. 4D). Only 11/70 neurons that produced nonsynchronized responses had BFs
1.5 kHz. Of those, only six of 11 had 50% response bandwidths that would produce a response to frequencies <1 kHz, which was the highest click rate tested.
The prevalent nonsynchronized responses that we observed did not arise from a fusion of click waveforms. Over 80% of nonsynchronized responses (34/41 nonsynchronized responses) produced significant rate responses at ICIs corresponding to distinct click waveforms and all of the other neurons produced significant responses at ICIs for which there was some fusion but still clear AM of the sound envelope (>30% modulation depth). It is possible that the use of Gaussian clicks biased the responses toward nonsynchronized responses because of their spectral peak and limited bandwidth relative to rectangular clicks. However, rectangular clicks also produced nonsynchronized responses in 13 neurons.
Comparison with previous studies of auditory cortex
Although MGB neurons exhibited synchronized and nonsynchronized responses, they differed from the synchronized and nonsynchronized responses observed in auditory cortex neurons using the same experimental preparation (awake marmoset) (Lu et al. 2001a
). The ICI ranges covered by the synchronized and nonsynchronized discharges were shifted to shorter ICIs (Fig. 5C). The rate boundaries of MGB neurons were also shifted to shorter ICIs (Fig. 5C). Nonsynchronized responses from the auditory cortex of anesthetized cats had rate boundaries similar to those reported here (Lu and Wang 2000
), suggesting that the nonsynchronized responses observed in the anesthetized auditory cortex reflected the responses of the thalamocortical inputs. More important than the simple shifts toward shorter ICI in the MGB was the prevalence of Mixed responses in the MGB compared with the IC or auditory cortex. This suggests that the MGB is beginning the transition between the synchronized and nonsynchronized responses that becomes more distinct in cortex.
Sensitivity to spectrotemporal parameters
Approximately two thirds (25/38) of the neurons tested with more than one click sigma exhibited a significant shift in their response, shifting either their Fmax or their preferred click rate by more than one octave (Fig. 7). Thus the temporal processing abilities of most neurons were stimulus dependent. Only one third of the MGB neurons were insensitive to the spectrotemporal stimulus parameters with respect to temporal processing. Neurons whose responses were stimulus dependent typically generated better synchronized responses for larger-click sigma values and poorly synchronized or nonsynchronized responses for smaller click sigma values (Fig. 7). For these neurons, this implies that despite the shorter duration and therefore increased temporal localization of the smaller sigma stimuli (Fig. 1A), synchronization was generally better in response to the longer-duration envelopes with sharper spectral peaks.
One possibility is that differences in the energy at or near the neuron's best frequency might be responsible for the stimulus dependency. However, similar to auditory cortex neurons (Lu et al. 2001a
), energy normalization did not change the temporal response properties of MGB neurons (n = 6 MGB neurons tested). Moreover, despite large changes in the bandwidths of stimuli with different sigma values, it seems unlikely that the sensitivity to click sigma was attributable to shifting the click spectrum into or out of the neuron's pure-tone frequency receptive field. Of the 38 neurons tested with more than one sigma, seven of 38 had BFs <1.5 kHz. None of these neurons switched from Nonsynch to Synch. All of the neurons whose behavior changed from Nonsynch to Synch had BFs 620 kHz, so the Nonsynch activity was almost definitely not the result of direct excitation (see Fig. 1B for example with 5.0-kHz carrier).
Another possibility is that the neurons require a minimum stimulus duration to produce sustained synchronization, making it more likely that better synchronization will occur for the larger-click sigma values. We conclude that many MGB neurons will respond to a repetitive stimulus if there is energy at their best frequency, but the exact nature of the response will depend not only on the stimulus repetition rate, but also on the carrier spectrum and envelope duration in a predictable manner. Similar sensitivity to the carrier's fine structure and envelope shape was also previously observed in IC (Rees and Møller 1987
) and auditory cortex (Lu et al. 2001b
).
Potential mechanisms for generating synchronized and nonsynchronized responses
MGB synchronized discharges originate from the synchronized inputs of IC neurons, but the Fmax of many MGB neurons is lower than that of IC neurons. Inhibition apparently regulates the maximum repetition rate at which many MGB neurons can synchronize (Fig. 2, A and C). Inhibition could arise from local MGB interneurons (Sherman 2004) or from the feedforward inhibitory projection from the IC (Peruzzi et al. 1997
), which can strongly affect the magnitude and timing MGB postsynaptic potentials (Bartlett and Smith 1999
, 2002
), similar to the way in which inhibition influences spike generation in auditory cortex (Wehr and Zador 2003).
Another inclusive possibility is that two classes of IC excitatory inputs to MGB neurons determine the response pattern. Data from rat brain slices demonstrated that one type of IC excitatory input generates large, short-latency inputs that exhibit strong paired-pulse depression (Bartlett and Smith 2002
). These inputs would reliably propagate IC synchronized discharge patterns at low repetition rates, but the sustained synchronized response would eventually disappear at high repetition rates, similar to what was observed in Synch and Mixed responses (Fig. 3). The other type of IC excitatory input generates smaller, longer-latency inputs that exhibit paired-pulse facilitation at short interstimulus intervals (Bartlett and Smith 2002
). Neurons receiving these inputs would be unlikely to respond to long ICI stimuli because multiple small inputs would be required to reach firing threshold. Instead, they would respond to short ICI stimuli when the small inputs could summate to generate a longer-latency suprathreshold response, which is supported by the differences in the latencies of synchronized and nonsynchronized responses (Fig. 6). Activation of NMDA receptors in these neurons would produce sustained excitation that would bolster weak responses and smear temporal precision, both of which were previously observed in MGB and auditory cortex neurons in vitro (Bartlett and Smith 2002
; Rose and Metherate 2005). This could result in nonsynchronized responses in MGB neurons.
MGB is not merely a relay station
The MGB is not a simple relay of IC inputs. The main difference between IC and cortex representations of repetition rate appears to be in the creation of separate synchronized and nonsynchronized responses. Then the MGB acts as a transition stage in this transformation, with the emergence of nonsynchronized responses and the prevalence of mixed responses. A transition stage from IC to MGB to cortex was recently observed for the integration of tone and noise responses (Las et al. 2005
). Figure 8 depicts the transformation of temporal representations suggested by our study by comparing the response types across regions using data from three studies that recorded from single units in awake animals. The IC study looked at single-unit responses in awake rabbits (Batra et al. 1989
) and the MGB and cortex studies were done using awake marmosets (Lu et al. 2001a
; present study). As a population, MGB neurons may create separate temporal and rate coding populations that could potentially represent temporal and spectral features, respectively. Complementary MGB rate and synchronization boundaries (Fig. 5C) are near perceptual boundaries of periodicity perception versus perception of fine structure, a transition from spectrotemporal to purely spectral features (Rosen 1992
). Rate responses at short ICIs in the absence of synchronized responses at longer ICIs is a response pattern that first appears in the MGB and is quite prevalent. Such a transformation of response pattern indicates that MGB neurons are not simple relays and contribute to the neural representations of sound features that serve as the inputs to auditory cortex.
|
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
Present address of E. L. Bartlett: Purdue University, Depts. of Biological Sciences and Biomedical Engineering, 206 S. Intramural Drive, BMED 2023, West Lafayette, IN 47907-1791.
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: E. L. Bartlett or X. Wang, Johns Hopkins University, Dept. of Biomedical Engineering, Traylor 412, 720 Rutland Avenue, Baltimore, MD 21205 (E-mail: ebartle{at}purdue.edu or xiaoqin.wang{at}jhu.edu)
| REFERENCES |
|---|
|
|
|---|
Aitkin LM, Kudo M, Irvine DR. Connections of the primary auditory cortex in the common marmoset, Callithrix jacchus jacchus. J Comp Neurol 269: 235248, 1988.[CrossRef][ISI][Medline]
Bartlett EL, Smith PH. Anatomic, intrinsic, and synaptic properties of dorsal and ventral division neurons in rat medial geniculate body. J Neurophysiol 81: 19992016, 1999.
Bartlett EL, Smith PH. Effects of paired-pulse and repetitive stimulation on neurons in the rat medial geniculate body. Neuroscience 113: 957974, 2002.[CrossRef][ISI][Medline]
Bartlett EL, Wang X. Long-lasting modulation by stimulus context in primate auditory cortex. J Neurophysiol 94: 83104, 2005.
Batra R, Kuwada S, Stanford TR. Temporal coding of envelopes and their interaural delays in the inferior colliculus of the unanesthetized rabbit. J Neurophysiol 61: 257268, 1989.
Burger RM, Pollak GD. Analysis of the role of inhibition in shaping responses to sinusoidally amplitude-modulated signals in the inferior colliculus. J Neurophysiol 80: 16861701, 1998.
Calford MB. The parcellation of the medial geniculate body of the cat defined by the auditory response properties of single units. J Neurosci 3: 23502364, 1983.[Abstract]
Condon CJ, White KR, Feng AS. Processing of amplitude-modulated signals that mimic echoes from fluttering targets in the inferior colliculus of the little brown bat, Myotis lucifugus. J Neurophysiol 71: 768784, 1994.
Creutzfeldt O, Hellweg FC, Schreiner C. Thalamocortical transformation of responses to complex auditory stimuli. Exp Brain Res 39: 87104, 1980.[ISI][Medline]
Escabi MA, Schreiner CE. Nonlinear spectrotemporal sound analysis by neurons in the auditory midbrain. J Neurosci 22: 41144131, 2002.
Hashikawa T, Molinari M, Rausell E, Jones EG. Patchy and laminar terminations of medial geniculate axons in monkey auditory cortex. J Comp Neurol 362: 195208, 1995.[CrossRef][ISI][Medline]
Hu B. Functional organization of lemniscal and nonlemniscal auditory thalamus. Exp Brain Res 153: 543549, 2003.[CrossRef][ISI][Medline]
Hu B, Senatorov V, Mooney D. Lemniscal and non-lemniscal synaptic transmission in rat auditory thalamus. J Physiol 479: 217231, 1994.[ISI][Medline]
Jones EG. Chemically defined parallel pathways in the monkey auditory system. Ann NY Acad Sci 999: 218233, 2003.
Jones EG, Hendry SHC. Differential calcium binding protein immunoreactivity distinguishes classes of relay neurons in monkey thalamic nuclei. Eur J Neurosci 1: 222246, 1989.[CrossRef][ISI][Medline]
Joris PX, Schreiner CE, Rees A. Neural processing of amplitude-modulated sounds. Physiol Rev 84: 541577, 2004.
Kaas JH, Hackett TA. Subdivisions of auditory cortex and processing streams in primates. Proc Natl Acad Sci USA 97: 1179311799, 2000.
Krishna BS, Semple MN. Auditory temporal processing: responses to sinusoidally amplitude-modulated tones in the inferior colliculus. J Neurophysiol 84: 255273, 2000.
Krystal JH, D'Souza DC, Mathalon D, Perry E, Belger A, Hoffman R. NMDA receptor antagonist effects, cortical glutamatergic function, and schizophrenia: toward a paradigm shift in medication development. Psychopharmacology (Berl) 169: 215233, 2003.