JN Ad Instruments
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 93: 3224-3247, 2005. First published December 15, 2004; doi:10.1152/jn.00445.2004
0022-3077/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
93/6/3224    most recent
00445.2004v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Grubb, M. S.
Right arrow Articles by Thompson, I. D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Grubb, M. S.
Right arrow Articles by Thompson, I. D.

Visual Response Properties of Burst and Tonic Firing in the Mouse Dorsal Lateral Geniculate Nucleus

Matthew S. Grubb and Ian D. Thompson

University Laboratory of Physiology, Oxford, United Kingdom

Submitted 29 April 2004; accepted in final form 11 December 2004


 ABSTRACT
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Thalamic relay cells fire action potentials in two modes: burst and tonic. Previous studies in cats have shown that these two modes are associated with significant differences in the visual information carried by spikes in the dorsal lateral geniculate nucleus (dLGN). Here we describe the visual response properties of burst and tonic firing in the mouse dLGN. Extracellular recordings of activity in single geniculate cells were performed under halothane and nitrous oxide anesthesia in vivo. After confirming that the criteria used to isolate burst spikes from these recordings identify firing events with properties described for burst firing in other species and preparations, we show that burst firing in the mouse dLGN occurs during visual stimulation. We then compare burst and tonic firing across a wide range of visual response characteristics. While the two firing modes do not differ with respect to spatial summation or spatial frequency tuning, they show significant differences in the temporal domain. Burst spikes are phase advanced relative to their tonic counterparts. Burst firing is also more rectified, possesses sharper temporal frequency tuning, and prefers lower temporal frequencies than tonic firing. In addition, contrast-response curves are more step-like for burst responses. Finally, we present analyses that describe the stimulus detection abilities and spike timing reliability of burst and tonic firing.


 INTRODUCTION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
At its most abstract and fundamental level, the role of the dorsal lateral geniculate nucleus (dLGN) is to feed the cerebral cortex with visual information. Most of this information comes from the retina the afferents of which provide the "driving" influence that determines the large majority of geniculate single cell response properties (e.g., Hubel and Wiesel 1961Go; Sherman and Guillery 1996Go, 2002Go; Usrey et al. 1999Go). However, retinal afferents do not represent the large majority of inputs to the dLGN. Dominating the dLGN's afferent population are the 85–90% of synapses that provide "modulating" influences on geniculate processing (e.g., Sherman and Guillery 1996Go, 2002Go). Although these modulatory inputs cannot change fundamental receptive field (RF) properties established by retinal drivers, they can produce subtle alterations in the response properties of dLGN neurons, such increasing or decreasing the magnitude of visually evoked responses (e.g., Uhlrich et al. 1995Go, 2002Go). Thus due to modulatory inputs onto dLGN cells, the thalamic relay of visual information from retina to cortex becomes a dynamic operation.

One important way in which modulatory inputs to the dLGN can alter the retinocortical transfer of visual information is by changing the firing mode of geniculate relay cells. These neurons, like all thalamic relay cells, can fire action potentials in two different modes: "burst" and "tonic." As demonstrated by in vitro recordings (e.g., Jahnsen and Llinas 1984Go; Kim et al. 2001Go; Ramcharan et al. 2000aGo; Sherman 2001Go; Zhan et al. 1999Go), firing mode is controlled by a T-type calcium channel that produces IT, a voltage- and time-dependent current. IT is inactivated at relatively depolarized membrane potentials (>60 mV). In this state, depolarizing inputs stimulating a geniculate relay cell will produce conventional, tonic action potentials. If the cell's membrane potential is relatively hyperpolarized for >100 ms, however, IT is de-inactivated. Cell stimulation then produces a depolarizing waveform carried predominantly by calcium ions, known as a low-threshold spike (LTS). This LTS is usually accompanied by a burst of sodium-based action potential firing. Modulatory geniculate inputs that produce (long-lasting) depolarizations or hyperpolarizations of cell membrane potentials can therefore switch relay cell responses from burst to tonic firing, and vice versa.

One of the most important situations in which such burst-tonic switches occur is the transition between sleep and wakefulness. Thalamocortical relay neurons in naturally sleeping animals display rhythmic bouts of bursting ( McCarley et al. 1983Go; Weyand et al. 2001Go), which in vitro studies have established are due to cyclic interactions between IT and a depolarizing, hyperpolarization-activated current Ih ( McCormick and Bal 1997Go). With the transition from slow-wave sleep to wakefulness, modulatory inputs from the brain stem produce a progressive depolarization of thalamocortical relay cells. This inactivates IT, decreasing burst firing and promoting tonic responses ( McCarley et al. 1983Go; McCormick and Bal 1997Go).

Because rhythmic burst firing is most prominent during slow-wave sleep ( McCarley et al. 1983Go; Weyand et al. 2001Go) and because during slow-wave sleep the responsiveness of dLGN neurons to RF stimulation is decreased ( Livingstone and Hubel 1981Go), it might seem that burst firing in thalamocortical relay cells is unimportant as a carrier of sensory information. However, studies in the past decade have begun to uncover the possibility that both burst and tonic firing modes represent specific signals in sensory processing. Lu et al. (1992)Go used intracellular recordings in vivo to demonstrate that burst responses occur in the cat dLGN during periods of visual stimulation. These recordings also showed that burst spikes could be accurately identified purely on the basis of temporal patterns in spike firing ( Lu et al. 1992Go), allowing less techni-cally demanding extracellular recordings to confirm that bursting, like tonic firing, occurs during sensory stimulation in anesthetized (e.g., Guido et al. 1992Go) and awake behaving (e.g., Fanselow et al. 2001Go; Guido and Weyand 1995Go; Ramcharan et al. 2000bGo; Swadlow and Gusev 2001Go; Weyand et al. 2001Go) animals. Furthermore, while burst and tonic dLGN spikes carry similar amounts of visual information ( Reinagel et al. 1999Go), recordings in the cat dLGN have shown that the type of visual information transmitted by the two firing modes can be significantly different. Among other differences, burst spikes in the cat follow temporal changes in a sinusoidal stimulus less faithfully ( Guido et al. 1992Go; Lu et al. 1992Go), offer better stimulus detection ( Guido et al. 1995Go), are more tightly tuned for TF ( Mukherjee and Kaplan 1995Go), and have more reliable timing ( Guido and Sherman 1998Go) than their tonic counterparts. Changes in firing mode in thalamic relay cells might therefore affect the information relayed by those cells to cortex.

Although studies in the cat have identified several differences in the stimulus coding properties of burst and tonic spikes in the dLGN ( Guido and Sherman 1998Go; Guido et al. 1992Go, 1995Go; Lu et al. 1992Go; Mukherjee and Kaplan 1995Go) (see DISCUSSION), the comprehensive set of data gathered in our laboratory concerning the response properties of mouse dLGN neurons ( Grubb and Thompson 2003Go) offers an opportunity to compare burst and tonic firing across a large range of visual response characteristics. As well as providing a full and detailed description of the functional differences between these two firing modes, such a comparison will also offer an indication as to whether such differences are conserved across mammalian species. We know that the properties of the thalamic IT current and its accompanying burst activity are similar across different mammals (guinea pig: Jahnsen and Llinas 1984Go; McCormick and Feeser 1990Go; cat: Zhan et al. 1999Go; monkey: Ramcharan et al. 2000aGo; mouse: Kim et al. 2001Go), but it is unknown whether the distinct visual response characteristics of burst and tonic firing ( Guido and Sherman 1998Go; Guido et al. 1992Go, 1995Go; Lu et al. 1992Go; Mukherjee and Kaplan 1995Go) are a common feature of mammalian visual processing at thalamic level. The mouse dLGN data presented here represent a first step toward addressing this issue.

Perhaps more importantly, though, studying the visual coding properties of burst and tonic firing in the mouse thalamus opens up the study of firing mode function to genetic manipulation techniques. In mammals, these techniques are most advanced in mice, where they can be invaluable tools for investigating neuronal function at the molecular level. In particular, because the switch between burst and tonic thalamic firing modes can be determined by a single ion channel, investigations of firing mode function might especially benefit from approaches involving genetic manipulation. Mice that lack the T-type calcium channel underlying IT have been generated ( Kim et al. 2001Go), but before the consequences of this mutation on visual processing at the thalamic or cortical level can be investigated, we need to know the probable contribution of this channel in the normal visual system. We can start to address this issue by investigating the basic visual response properties of burst and tonic firing in the wild-type mouse dLGN. In addition, genetic manipulations might be extremely useful in dissecting the molecular components underlying the control of response mode in thalamic relay cells. Various modulatory systems have been proposed to switch thalamic relay cells from burst to tonic mode or vice versa (e.g., Godwin et al. 1996Go; Lu et al. 1993Go; Uhlrich et al. 2002Go). Manipulating the receptors or other molecular machinery involved in this modulation could be extremely informative concerning the circuitry that controls or does not control thalamocortical information transmission.

As an essential precursor to such studies, we demonstrate here that burst-like events can be identified in extracellularly recorded responses of mouse dLGN cells and that these events have some properties expected of thalamic burst firing. We go on to show that burst firing in the mouse dLGN occurs during visual stimulation. Finally, we investigate the visual response characteristics of burst and tonic spikes in the mouse dLGN across a broad range of stimulus parameters.


 METHODS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Animals

Animals were adult (>3 mo of age) pigmented C57Bl/6J mice of either sex weighing 20–36 g (Harlan Olac, Oxon, UK). Mice were housed under a 12:12-h cycle of light and dark. All procedures were conducted under the auspices of United Kingdom Home Office project and personal licenses held by the authors.

Electrophysiology

We carried out electrophysiological recordings as described previously ( Grubb and Thompson 2003Go; Grubb et al. 2003Go). Mice were first anesthetized with an intraperitoneal injection of 25% fentanyl and fluanisone (Hypnorm, Janssen Animal Health, Bucks, UK) and 25% midazolam, Down, UK), administered at 2.7 µl/g. Surgical anesthesia was maintained throughout the setup procedures by subsequent applications of this induction dose. To prevent the accumulation of bronchial secretions, 5 µl atropine sulfate (600 µg/ml, Animalcare, Yorks, UK) was given subcutaneously. A tracheotomy was then performed ( Schwarte et al. 2000Go), and a small plastic tube (1 mm OD, 0.5 mm ID, SIMS Portex, Kent, UK) inserted and secured by ligation around the trachea. Application of a spot of cyanoacrylate glue added extra stability by attaching the tube to the skin of the upper thorax.

The animal was situated in a custom-built head holder and secured using a bite bar and ear bars before the tracheal tube was connected to a respiratory pump (MiniVent 845, Hugo Sachs Elektronik, March-Hugstetten, Germany). Pump rate (usually 120–140 strokes/min) and stroke volume (usually 150–200 µl) were adjusted to maintain the expired carbon dioxide level, monitored with a Capstar-100 CO2 analyzer (CWE, Ardmore, PA), at 2.5–4%. Once on the pump, the mice breathed a 1:3 mixture of oxygen and nitrous oxide, along with 1–1.5% halothane for maintained anesthesia. The mice were not paralyzed at any stage, allowing continued assessment of anesthetic depth via the paw pinch response. Throughout the experiment, heart rate was monitored via an ECG (~5 Hz was normal). Core body temperature, measured with a rectal probe, was maintained at ~37°C by combining a high ambient temperature with heat from a thermostatically coupled blanket (NP 50-7061-R, Harvard Apparatus, Kent, UK).

A circular craniotomy ~0.5 mm in diameter was made unilaterally above the presumptive location of the dLGN, ~2.5 mm posterior and 2 mm lateral of the Bregma suture ( Paxinos and Franklin 2001Go). A small durotomy was made in the center of this exposed zone and, after any whiskers occluding the eyes had been trimmed, a tungsten-in-glass recording electrode with a 5- to 10-µm exposed tip (Alan Ainsworth, Northants, UK) was lowered vertically into the brain using a microdrive. Changes in background activity while advancing through the brain were used to locate the electrode relative to the dLGN ( Grubb and Thompson 2003Go). Once in the dLGN, the electrode was advanced or retracted in 5-µm steps. Signals were amplified locally by a FET headstage, further amplified by a Neurolog preamp (NL104; typically 10–20 k), then filtered (Neurolog NL125, low-frequency ~150 Hz, high-frequency ~3,000 Hz) for display on an oscilloscope. A Schmidt trigger was used to isolate responses of single neurons by adjusting the signal threshold; the trigger was also connected to a storage oscilloscope that enabled us to monitor spike shapes and to ensure that no spikes were seen in the refractory period. To isolate the responses of single units, electrode position was adjusted until spikes were deemed to have the same shape, amplitude, and rise time. Responses of optic tract fibers—typically large and fast, with high levels of spontaneous firing and a lack of burst activity—were easily distinguished from cell soma responses. Spike events were then time-stamped with a resolution of 100 µm and were linked to the timing of visual stimulus changes using in-house software written by D. Smyth.

Although we did not apply contact lenses, optical quality in our setup was good and did not degrade over the course of experiments. As described previously ( Grubb and Thompson 2003Go), no significant correlations were found between time from experiment start (≤13 h) and any measures of spatial tuning in dLGN neurons. Furthermore, within animals, there were no significant differences in these measures between cells recorded in the first versus the second half of experiments (paired t-test or Wilcoxon paired test, depending on normality, P > 0.05). Eye movements under halothane anesthesia were rare and small ( Grubb and Thompson 2003Go; cf. Dräger 1975Go; Gordon and Stryker 1996Go; Wagor et al. 1980Go). As described previously ( Grubb and Thompson 2003Go), remapping of RFs 0.5–1 h after their initial localization revealed a mean shift (±SE) in position of only 3.6 ± 0.5° (n = 11). In addition, over a time scale of ~10 min, response features such as phase-locked raster plots (Fig. 5), sinusoidal dependence of cell responses on stimulus phase (Fig. 6), and regular shifts in response phase with increasing temporal frequency (Fig. 5), all suggest that RF position changed very little. Moreover, in a study where properties of burst and tonic firing were always compared within individual tuning curve experiments, small eye movements were unlikely to affect the present findings.



View larger version (30K):
[in this window]
[in a new window]
 
FIG. 5. Response latency and response phase of burst and tonic spikes in the mouse dLGN. A: burst spikes tend to precede tonic spikes in responses to sinusoidal grating stimuli. Plots show responses of a mouse dLGN cell to a drifting sinusoidal grating at optimal spatial frequency, optimal temporal frequency and 100% contrast. In raster plots (left) {bullet}, 1 spike. In poststimulus time histograms (PSTHs; right) spikes are binned across all stimulus cycles according to their poststimulus time (pst). The "all spikes" PSTH possesses 2 peaks. Separating responses into burst and tonic components shows that these 2 peaks correspond to these 2 different firing modes. Burst spikes (middle) are concentrated at the start of the cell's response to the stimulus, whereas tonic spikes (bottom) tend to occur later in the stimulus cycle. B: dissecting possible factors underlying the precedence of burst spikes. In these plots of temporal frequency (TF) vs. F1 response phase, calculated separately for a given cell's burst and tonic spikes, the slope of the line relating the 2 variables can be taken as a measure of response latency and the y intercept of this line can give an indication of relative response phase (e.g., Guido et al. 1992Go; Mukherjee and Kaplan 1995Go). Plots 1 and 2 were typical of mouse dLGN cells. At low TFs, the response phase of burst spikes is higher than that of tonic spikes, representing a phase advance (i.e., burst spikes occur first). At higher TFs, the difference between burst and tonic response phase becomes small, such that the slope of the burst spikes line is greater than that of the tonic spikes line. In such cells, therefore burst spikes possess greater latency but are still phase advanced relative to their tonic counterparts. Plot 3 shows the 1 cell (of 30) in which tonic spikes were phase advanced relative to burst firing. Note the strong nonlinearity of the TF-response phase relationship in this case. C: burst spikes have longer latencies than tonic spikes. In the scatter plot each {bullet}, 1 cell (n = 30). Most dots lie beneath the (- - -) line of unity, suggesting that burst latency is on average longer than tonic latency. Inset: the medians of the 2 distributions are significantly different (burst latency median: 106 ms; tonic: 90 ms; Wilcoxon matched pairs test, P = 0.001). b, burst; t, tonic. D: despite their longer latency, burst spikes are almost always phase advanced relative to tonic spikes. The histogram plots the distribution of response phase advances across the mouse dLGN sample (n = 30). Only one cell, the cell in B3, possessed a negative phase advance indicative of earlier tonic firing. In all other cells, burst spikes were evoked by earlier phases of the sinusoidal stimulus than tonic spikes. The mean of the distribution, at 1 radian, is marked ({blacktriangledown}).

 


View larger version (35K):
[in this window]
[in a new window]
 
FIG. 6. Linearity of burst and tonic firing in mouse dLGN neurons. A: assessing linearity of spatial summation in burst and tonic firing. Plots show F1 and F2 response amplitudes during presentation of stationary counterphased sinusoidal gratings at various spatial phases. For each example cell, separate plots are shown for spikes that formed part of bursts (burst firing, top) and those that did not (tonic firing, bottom). By convention, F1 responses are plotted as negative when their phase differs by 90–270° from that of the maximum F1 response; this explains why the burst and tonic F1 plots in 1 and 2 look like mirror images of each other. Plots 1 and 2: typical of mouse dLGN cells with burst and tonic firing both displaying linear spatial summation. In both response modes, F1 responses vary sinusoidally with stimulus phase, cross the 0 line at 2 "null" positions and are larger than corresponding F2 responses (cf. Hochstein and Shapley 1976Go). Relatively large and modulating F2 responses here reflect substantial response rectification, especially for burst spikes; they are not indicative of frequency-doubled responses. Plot 3: the only cell in our sample to display substantial nonlinearity in spatial summation. Here this nonlinearity, evidenced by a dominant and modulating F2 response, is observed in both response modes but is most prominent in tonic firing. Linearity values for all 3 cells: cell 1: burst 0.12; tonic 0.16; cell 2: burst 0.40; tonic 0.15; cell 3: burst 0.56; tonic 1.64. B: no difference between burst and tonic firing in linearity of spatial summation. {bullet}, a single null test experiment, plotted according to Hochstein and Shapley's (1976)Go index of linearity for both burst and tonic firing. Linearity values >1 reflect nonlinear summation, here shown only by the tonic responses of the cell in A3 (linearity = 1.64). Most {bullet} are distributed evenly across the - - - line of unity, with the inset showing that the mean linearity values for burst and tonic firing are not significantly different (burst means ± SE: 0.42 ± 0.02; tonic: 0.37 ± 0.04; paired t-test, P = 0.17). Error bars show SE. C: burst firing is more rectified than tonic firing. The scatter plot shows all tuning curve experiments (not just null tests), each plotted as a single dot according to the rectification index (see METHODS; Guido et al. 1992) of its burst and tonic firing components. The majority of the dots lie below the - - - line of unity, and the inset shows that the medians of the 2 response modes are significantly different (burst median: 0.46; tonic: 0.36; Wilcoxon matched pairs test, P < 0.0001).

 
Stimulus presentation

Visually responsive neurons were first identified using a hand-held ophthalmoscope. Monochrome visual stimuli were then presented on a cathode ray tube monitor under the control of a Visual Stimulus Generator 2/4 graphics card (Cambridge Research Systems, Cambridge, UK); this had a pseudo-15-bit analog output allowing precise control of the luminance of each pixel on the display and correction for expansive luminance nonlinearities. Stimuli could be presented with 256 linearly spaced gray levels. The display comprised 800 x 600 pixels and at viewing distances between 135 and 400 mm subtended 53–111 x 41–95°. With maximum and minimum luminances of 99 and 2cd/m2, respectively, the display was the main source of illumination in an otherwise dim room. The frame rate of the monitor was 160 Hz.

The display was initially placed with a given cell's presumptive RF at its center; mapping with reverse correlation ( Grubb and Thompson 2003Go) then confirmed this location, or directed movement of the monitor in order that the display covered the entire RF area. Only cells the RFs of which could be localized in this way were subjected to subsequent analysis of response characteristics. Stimuli were presented with both eyes open, but responses in all cells were most likely driven entirely by the contralateral eye: not only are all wild-type mouse dLGN cells monocularly driven ( Grubb et al. 2003Go), but in these investigations, the display was often placed in positions that could be viewed only by the contralateral eye. Furthermore, although it is still possible that interactions between the two eyes could modulate the responses of single mouse dLGN neurons, such interactions are unlikely to affect present analyses in which burst and tonic response properties were compared within cells.

For quantitative investigation of the visual response properties of burst and tonic firing, sinusoidal monochromatic vertical gratings were presented covering the full extent of the display. These gratings were drifting in all experiments except the null test, in which they were stationary and sinusoidally contrast-modulated. In all experiments, all grating parameters were kept constant except the one under study; gratings varying in this parameter were presented in a pseudorandom sequence. Grating parameters in spatial frequency (SF) and temporal frequency (TF) tuning assessments were varied systematically and logarithmically to produce 8 stimuli ranging between ~0.01 and ~1 c/° inclusive for SF assessments, and between 0.5 and 16 Hz inclusive for TF assessments. In contrast and null phase test assessments, grating parameters were varied systematically and linearly to produce 12 stimuli ranging between 1 and 100% inclusive for contrast assessments and between –180 and +150° inclusive for null phase tests. Every stimulus repeat also contained a blank stimulus at mean luminance from which spontaneous activity levels were calculated. Each presentation of a grating stimulus in a given repeat occurred for seven cycles; of these only responses from the last five cycles were used to discard any flash responses to the initial presentation of the stimulus. The exception was assessing TF tuning. In these experiments, stimuli ≤1 Hz were presented for seven cycles (with the 1st 2 discarded) as usual. However, to ensure capture of enough responses to high TF stimuli, stimuli >1 Hz were presented for seven seconds, with the first 2 s being discarded to allow for flash responses. At least two repeats of each 7-cycle (or 7 s) stimulus were presented to every cell, although in the majority of cases four repeats were utilized.

Histology

Each experiment was ended with an intraperitoneal injection of 0.3 ml pentobarbitone sodium (Sagatal, Rhône Mérieux, Essex, UK). The mouse was then transcardially perfused with phosphate-buffered saline (Sigma, Dorset, UK) followed by 4% paraformaldehyde (TAAB Laboratories, Berks, UK) in 0.1 M phosphate buffer. The brain was removed, sunk in 30% sucrose for cryoprotection, and sectioned coronally at 30–50 µm on a freezing microtome. Sections were then mounted from distilled water onto gelatinized slides. Subsequent staining for cell bodies with cresyl violet allowed the identification of small (6 µA for 6 s) electrolytic lesions made at precise locations on successful penetrations. Locating these lesions made it possible to reconstruct electrode tracks and localize recorded units. Only cells unequivocally located within the dLGN were included in subsequent physiological analyses.

Data analysis

IDENTIFYING BURST AND TONIC RESPONSE COMPONENTS.  Spikes forming part of bursts triggered by low-threshold calcium spikes (LTS bursts) were identified from temporal patterns in the firing of dLGN cells. From intracellular recordings in the cat dLGN, Lu et al. (1992)Go developed criteria that could reliably identify bursts driven by LTS activity; we applied the same criteria here. Specifically, the first action potential in a burst had a preceding silent period of ≥100 ms and a following interspike interval (ISI) of ≤4 ms. Any subsequent spikes with preceding ISIs of ≤4 ms were designated as additional action potentials in a burst. All spikes that did not form part of LTS bursts were designated as part of tonic firing.

RHYTHMICITY.  Rhythmicity in spontaneous bursting was assessed after Weyand et al. (2001)Go. For each 5-s presentation of a 51 cd/m2 blank screen during which at least four bursts were fired, we calculated the SD of interburst intervals, the latter being defined as the time between the first spikes of consecutive bursts. This SD is a direct measure of burst rhythmicity—perfect rhythmic bursting should produce a SD of 0. To assess the statistical significance of each SD measure, we employed a Monte Carlo approach ( Weyand et al. 2001Go). For each blank stimulus presentation, 999 simulations were run in which an identical number of bursts were randomly distributed within the 5-s interval. Standard deviations of interburst intervals were calculated for each simulation and were ranked alongside the original, measured value. The rank of the value of the observed SD then gave a statistical measure of burst rhythmicity: significant rhythmicity (P < 0.05) was detected if the observed value fell within the 50 lowest simulated values.

FOURIER ANALYSIS.  Fourier analysis was applied separately to a given cell's tonic and burst spikes. Burst or tonic responses to a given grating stimulus were first binned according to their poststimulus time. One hundred twenty eight bins were used per stimulus cycle, producing a poststimulus time histogram (PSTH) for each stimulus. A fast Fourier transform of these data then converted these PSTH data into amplitudes (in spikes/s) and phases (in degrees or radians) of harmonics of the stimulus' modulation frequency ( Hochstein and Shapley 1976Go).

RESPONSE LINEARITY.  Linearity of spatial summation was assessed using a linearity measure first applied by Hochstein and Shapley (1976)Go. Over the many stimulus phases presented in a single modified null test, the mean second harmonic (F2) response amplitude was divided by the maximum fundamental harmonic (F1) response amplitude. Linearity values <1 are indicative of linear spatial summation, while linearity values greater than one are indicative of nonlinear spatial summation.

Another measure of response linearity, "rectification," was calculated over a given tuning curve as follows

(1)
where F2 and F1 are the amplitudes of the second and first harmonics, respectively, in responses fired to a stimulus s. Although rectification might be best measured by comparing the amplitude of F1 response components with the combined amplitudes of all higher response harmonics (F2, F3, F4, and so on), this particular measure was employed here because of its previous use as an indicator of response nonlinearity in cat dLGN cells ( Guido et al. 1992Go; Lu et al. 1992Go).

CURVE FITTING.  All curve fits were carried out using a least-squares minimalization algorithm ( Grubb and Thompson 2003Go).

Plots of stimulus SF versus cell F1 response amplitude were fitted with a difference of Gaussians curve ( Grubb and Thompson 2003Go; Rodieck 1965Go; So and Shapley 1981Go)

(2)
where R is F1 response amplitude, {upsilon} is SF, b is baseline response, kc is the area under the RF center's Gaussian function, ks is the relative area under the RF surround's Gaussian function, and rc and rs are the radii of the center and surround Gaussian functions, respectively, at the point each mechanism has reached 1/e of its peak.

Plots of stimulus TF versus cell F1 response amplitude were fitted with an atheoretical function comprising two half-Gaussians ( Grubb and Thompson 2003Go)

(3a)

(3b)
where R is F1 response amplitude, {omega} is TF, p is peak TF, a is the response amplitude at optimum TF, s is the Gaussian spread, b1 is the baseline on the low-frequency side of the curve, and b2 is the baseline on the high-frequency side of the curve. Plots of stimulus TF versus cell F1 response phase were simply fitted with a straight line.

Plots of stimulus Michelson contrast versus cell F1 response amplitude were fitted using a hyperbolic function ( Albrecht and Hamilton 1982Go; Grubb and Thompson 2003Go) of the form

(4)
where R is F1 response amplitude, c is contrast, Rmax is maximum response, h is the contrast at which R reaches half of Rmax, and n is rate of change.

SIGNAL DETECTION ANALYSIS.  Receiver operating characteristic (ROC) curves were derived from probability distributions of spike counts obtained during a sampling episode {tau} of visually driven ({tau}s) and spontaneous ({tau}n) activity ( Guido et al. 1995Go). {tau}s comprised half the temporal period of the stimulus and was centered on the phase of the F1 response component. {tau}n comprised an equal period during the presentation of a blank screen (51 cd/m2). For all spikes fired to a given stimulus, a criterion domain was set corresponding to the full range of spike counts (including 0) occurring for {tau}s and {tau}n. ROC curves were then obtained by plotting, for each criterion level, the probability P(false alarm) of obtaining a criterion response during {tau}n versus the probability P(hit) of obtaining a response during {tau}s. To produce a nonparametric measure of discrimination performance ( Green and Swets 1966Go; Macmillan and Creelman 1991Go), the area underneath each ROC curve was then calculated using simple geometry.

MEASURES OF SPIKE TIMING RELIABILITY.  Spike timing reliability was assessed using two measures. The first, "reliability," was applied to all burst or tonic spikes fired to a given stimulus, and represented the mean deviation of spike times from their F1 phase, normalized by the stimulus period

(5)
where s is the time of a given spike i, tF1 is the time of the F1 phase, n is the number of spikes, and p is the stimulus period.

First-spike reliability was applied only to the first burst or tonic spikes fired during a particular stimulus cycle and was calculated as the SD of first-spike latencies. Latencies were taken from the start of a half-cycle period centered on the F1 component of either burst or tonic firing, as appropriate

(6)
where f is the time of a given first spike i, t0 is the start time of the half-cycle period, and n is the number of first spikes. This measure was only calculated when n ≥ 10.

Statistical analysis

Sample sizes in this study were large enough (n > 30) to reliably assess their Normality using a Kolmogorov-Smirnov test. Those deemed Normal were described using the means ± SE and were compared using parametric tests. Those deemed non-normal were described using the median and were compared using nonparametric tests. Paired tests were used to compare burst and tonic spikes fired by the same cell to the same stimuli. All comparison tests were two-tailed, with the level of significance set at 0.05.


 RESULTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Identifying burst firing in the responses of mouse dLGN cells

Visual response properties of burst and tonic firing were assessed in 102 mouse dLGN cells, all with classic geniculate ON- or OFF-center receptive fields (see Grubb and Thompson 2003Go). Responses were recorded from cells with RFs at visual eccentricities ranging from 0 to 92° (means ± SE: 44 ± 3°). Because the response properties of mouse dLGN cells change very little with increasing visual eccentricity ( Grubb and Thompson 2003Go), we did not limit our analysis to cells with RFs in particular regions of visual space.

Burst firing in the responses of cells to sinusoidal grating stimuli was identified on the basis of temporal patterns of spikes. The first spikes in bursts were preceded by ≥100 ms of silence and were followed by another spike within 4 ms. Subsequent spikes in bursts followed at ISIs of ≤4 ms ( Lu et al. 1992Go) (see METHODS). Figure 1 illustrates the way in which the responses of mouse dLGN cells to drifting sinusoidal grating stimuli were analyzed using a double interspike interval (ISI) plot, revealing distinct temporal patterns of firing associated with burst responses. For each spike, subsequent ISI is plotted against preceding ISI. The distribution of points depends on the temporal properties of the stimulus and, more importantly, the balance between burst and tonic firing. The significance of the latter can be seen by comparing Fig. 1, A and C, in which example cells fire almost all of their spikes as part of burst or tonic firing, respectively. Figure 1B illustrates the more common occurrence in which both burst and tonic spikes contribute to a cell's responses.



View larger version (37K):
[in this window]
[in a new window]
 
FIG. 1. Separating burst and tonic components in the responses of mouse dorsal lateral geniculate nucleus (dLGN) cells. Each interspike interval (ISI) plot shows spikes fired during a particular tuning curve experiment. Each spike, represented by a single dot, is plotted according to the period between it and the spike preceding it (ISI pre, x axis) and the period between it and the spike following it (ISI post, y axis). Note the logarithmic axes. In the "all spikes" plot (left), spikes are clustered into prominent groups. In box 1, spikes are preceded by ≥100 ms of silence and are followed by another spike within 4 ms. These spikes are the 1st spikes of bursts. In box 2, spikes possess preceding and proceding ISIs of no more than 4 ms. Most of these spikes are the middle spikes of bursts. In box 3, spikes have preceding ISIs of no more than 4 ms and are followed by variable periods of inactivity. Most of these spikes are the last spikes of bursts. Spikes falling within other areas of the plot are all included in tonic firing. The line marked TF displays the temporal period of the stimulus. The "burst spikes" plot (center) shows just those spikes that are part of bursts. The "tonic spikes" plot (right) shows all spikes that are not part of bursts. A: a cell dominated by burst firing. Most spikes fall within boxes 1–3. Note that the x-axis position of spikes in box 1 and the y-axis position of spikes in box 3 are largely determined by the period of the stimulus. The y-axis position of cluster 1 and the x-axis position of cluster 3, however, are determined by features of thalamic burst firing. B: a cell with prominent burst and tonic firing. This cell's spikes show temporal patterns in addition to the clusters in boxes 1–3. Region 4 outlines last spikes in bursts that are followed by a relatively long spike-free period. These burst spikes are probably the last stimulus-driven spikes fired in a given stimulus cycle. In contrast, region 5 outlines last spikes in bursts that are quickly followed by tonic firing. Region 6, where tonic spikes are followed by a long dead time almost equal to the stimulus period, outlines last spikes in stimulus cycles. Region 7, where tonic spikes are preceded by a long dead time slightly shorter than the stimulus period, outlines the first spikes fired in stimulus cycles that did not evoke any burst responses. C: a cell dominated by tonic firing. Box 2 is empty and only very few 2-spike bursts occurred.

 
In the double ISI plots of Fig. 1, distinct clusters of spike ISI distributions can be seen. Some of the outlined zones are constrained by the features of burst firing; others are determined by the temporal structure of the visual stimulus. Zones 1–3 are all conditioned by the properties of bursts. Zone 1 contains the first spikes in bursts: these spikes are preceded by ≥100 ms of silence and are followed by another spike within 4 ms. Spikes within bursts, with pre- and post-ISIs of ≤4 ms, form the majority of spikes in zone 2, while the x-axis position of the cluster in zone 3 is due to these spikes being mostly the last spikes in bursts—they follow within 4 ms of a previous spike but are themselves followed by >4 ms of nonspiking "dead" time. In the example cell in Fig. 1A, dominated by burst responses, these zones contain almost all of the spikes fired. However, the temporal structure of the stimulus also affects the positioning of clusters in double ISI plots: assuming a neuron fires spikes to every stimulus cycle, maximum pre- and post-ISIs are constrained by the grating's temporal frequency (TF). In the example cell in Fig. 1C, dominated by tonic responses, the majority of spikes fall outside zones 1–3 but are constrained by the stimulus' TF. The example cell in Fig. 1B, in contrast, displays a complex ISI pattern indicative of important contributions from both burst and tonic firing. This produces a broad distribution of post-ISIs in zone 3: while many of these last spikes in bursts are the last fired in a cycle and are thus constrained largely by stimulus TF (zone 4), many others are followed by tonic firing and thus possess much shorter post-ISIs (zone 5). Similarly, within spikes fired as part of tonic firing, spikes clustered in zone 6 are the last spikes fired in a cycle, and those clustered in zone 7 are the first fired in a cycle.

Firing in mouse dLGN cells thus conforms to particular temporal patterns, some of which correspond well to those predicted by the properties of intracellularly recorded LTS bursts in the cat dLGN ( Lu et al. 1992Go). Separating mouse dLGN cell responses into burst and tonic components on the basis of the preceding criteria is therefore nonarbitrary: it is effective in classifying preexisting temporal patterns of neuronal firing. However, this classification is not quite perfect. Figure 1 shows that "tonic" spikes can sometimes occupy ISI clusters (zones 2 and 3) dominated by burst firing. Indeed, Lu et al. (1992)Go concede that their criteria are conservative, allocating to tonic firing ~5–10% of spikes that actually form part of bursts.

Characteristics of burst firing in the mouse dLGN

Bursting appears to be a well conserved thalamic phenomenon. As assessed in vitro, the properties of the IT current and its accompanying burst activity are similar across thalamic nuclei ( Jahnsen and Llinas 1984Go; Ramcharan et al. 2000aGo) and across different mammalian species (guinea pig: Jahnsen and Llinas 1984Go; McCormick and Feeser 1990Go; cat: Zhan et al. 1999Go; monkey: Ramcharan et al. 2000aGo; mouse: Kim et al. 2001Go). Furthermore, it was shown in the previous section that distinct patterns of firing in the mouse dLGN can be identified by criteria used to identify burst spikes in other species (cat: Lu et al. 1992Go; monkey: Ramcharan et al. 2000bGo). However, we wanted to be more confident that the bursting criteria used here were identifying the same LTS-based phenomenon studied in other species. We therefore investigated whether certain features of burst firing reported to occur in those species also existed in the mouse dLGN.

Bursts are longer if preceded by greater periods of silence

It is known from in vitro recordings in the cat ( Zhan et al. 1999Go) and monkey ( Ramcharan et al. 2000aGo) that the number of spikes in a burst is not related to the strength of the depolarizing input stimulus but rather to the hyperpolarization of the membrane potential when that depolarizing stimulus arrives. In vivo, membrane hyperpolarization should depend on the degree of inhibition on a cell. Such inhibition will vary over time, but might have most chance to hyperpolarize a cell when that cell has been silent for an extended period. One might therefore predict that, although nonspiking "dead time" and membrane hyperpolarization need not go hand in hand, the number of spikes in bursts in vivo might be greater when bursts are preceded by longer periods of silence. In fact, across all bursts recorded in all experiments in all cells here, a small but significant positive relationship was noted between burst length and the duration of the preburst silent period (Spearman r = 0.16, P < 0.0001, n = 42,400). These data were summarized by taking the median preburst silent period for different burst lengths (Fig. 2A). For bursts consisting of ≤5 spikes, which represented >99% of all bursts recorded, there is a clear increase in the median period of preburst silence with increasing burst length. For larger bursts, this relationship breaks down. This could be due to the relatively small sample of long bursts. Alternatively, Rowe and Fischer (2001)Go used paired retinal and geniculate recordings to show that the first spikes in thalamic bursts are usually stereotypical events triggered by a single retinal spike but that later spikes in long bursts can be triggered directly by subsequent retinal action potentials. In other words, long bursts in vivo may not be long because of a large initial membrane hyperpolarization ( Zhan et al. 1999Go), but because they are triggered by two or more retinal action potentials closely spaced in time. If the latter is true, the relationship between preburst dead time and burst length might be expected to break down for the very longest bursts recorded.



View larger version (19K):
[in this window]
[in a new window]
 
FIG. 2. Characteristics of burst firing in the mouse dLGN. A: bursts tend to be longer when preceded by longer periods of silence. The plot shows the median duration of the nonspiking silent period prior to bursts of different lengths. Numbers in brackets show the number of observations, n, for each burst length. For bursts ≤5 spikes in length, which represent >99% of all bursts observed, there is a clear gradual increase in the duration of the initial ISI with increasing burst length. B: bursts tend to be longer when the 1st intraburst ISI is short. The median initial intraburst ISI decreases steadily as burst length increases. All conventions as in A. C: interspike intervals increase in length during bursts. In C1, the analysis is carried out across bursts. Median durations are shown for increasing orders of intraburst ISIs: for example, point 2 shows the median duration of all ISIs occurring between the 2nd and 3rd spikes in bursts. Intraburst ISIs tended to be shorter if they occurred closer to the start of the burst. Only intraburst ISIs occurring toward the end of long bursts did not follow this trend, but these observations account for only <1% of all intraburst ISIs analyzed. In C2, the analysis is carried out within bursts. Bars show median differences between successive intraburst ISIs: bar 1, for example, shows the median difference between 2nd and 1st intraburst ISIs across all bursts of length ≥3. The important point here is that all the median values are >0. Any given intraburst ISI therefore tends to be shorter than the intraburst ISI that follows it. D: bursts are stereotyped events. Across all tuning curve experiments, the histogram shows the distribution of burst length CVs. The CV, calculated here as the SD of burst length divided by mean burst length in any given experiment, should be ~1 if burst lengths vary by chance ( Reinagel et al. 1999Go). That all CV values are <1, with a median of 0.25 (marked with the clear arrow), reflects the fact that most bursts in the mouse dLGN comprise few spikes.

 
Longer bursts have shorter initial intraburst ISIs

It has been previously reported that burst length in the cat ( McCarley et al. 1983Go) and human ( Radhakrishnan et al. 1999Go) thalamus is correlated with the duration of the first within-burst ISI. This was confirmed for the mouse dLGN here. Across all bursts there was a strong negative relationship between the first ISI within a burst and burst length (Spearman r = –0.41, P < 0.0001, n = 42,400). Taking median initial intraburst ISI values for each burst length also reveals that bursts are longer when the first intraburst ISI is shorter (Fig. 2B).

ISI duration increases within bursts

Another widely reported feature of thalamic bursting is that ISIs become progressively longer within bursts (e.g., McCarley et al. 1983Go; Radhakrishnan et al. 1999Go; Reinagel et al. 1999Go). This feature was also observed here. Grouping data across bursts in the mouse dLGN, there was a significant positive correlation between within-burst ISI number and ISI duration (Spearman r = 0.22, P < 0.0001, n = 77,384). Taking median ISI duration values for each ISI position also shows a gradual increase in duration from the average first through to the average fourth ISI in a burst (Fig. 2C1). Later within-burst ISIs do not continue the trend but represent <1% of the ISIs analyzed and, as noted in the preceding text, might be expected to be less stereotyped than earlier within-burst ISIs ( Rowe and Fischer 2001Go). A more powerful demonstration of increasing ISI duration within bursts comes, however, from analyzing differences in within-burst ISIs within individual bursts. Figure 2C2 shows that the median difference in duration between a within-burst ISI and its preceding within-burst ISI is always positive. All median difference values here are significantly greater than zero (Wilcoxon rank sum test, P < 0.05). In other words, any given intraburst ISI tends to be shorter than the intraburst ISI that follows it.

Bursts are stereotyped events

Following Reinagel et al. (1999)Go, the coefficient of variation (CV) of burst length was calculated for each experiment. Defined as the SD divided by the mean of any given sample, the CV should be ~1 if sample values vary by chance. Here, most bursts (47%) consist of only two spikes. This stereotyped length is reflected in low CV values in every tuning curve experiment (Fig. 2D). The median CV value across all tuning curve experiments was 0.25, strikingly similar to mean values reported for cat X (0.25) or Y (0.23) cells ( Reinagel et al. 1999Go).

Rhythmicity of spontaneous burst firing

Spontaneous bursting in the dLGN neurons of naturally sleeping animals often occurs rhythmically ( McCarley et al. 1983Go; Weyand et al. 2001Go), with rhythmicity established by cyclic interactions between two depolarizing geniculate currents IT and Ih ( McCormick and Bal 1997Go). We assessed whether spontaneous bursting is also rhythmic in the mouse dLGN under conditions of halothane anesthesia. Our analysis was based on 5-s periods of spontaneous activity in which the stimulus display consisted of a blank screen at mean luminance (51 cd/m2); such blank presentations occurred once for every repeat of sinusoidal grating stimuli in a given tuning experiment (see METHODS). Rhythmicity was assessed by analyzing interburst intervals during such 5-s spontaneous epochs. For each epoch that contained at least four bursts, we calculated the SD of interburst intervals, with this SD measure being directly related to bursting rhythmicity: perfect rhythmic bursting would produce a SD of zero. We then assessed the statistical significance of each SD rhythmicity measure using Monte Carlo simulations that ranked the observed rhythmicity value among 999 simulated values (see METHODS).

Examples of burst timing in the responses of a single mouse dLGN cell during several 5-s spontaneous epochs are presented in Fig. 3A. These plots show that evenly-spaced bursts produce significant rhythmicity values and also demonstrate that spontaneous bursting in the mouse dLGN could alternate between rhythmic and nonrhythmic behavior over time in the same neuron. The histogram in Fig. 3B presents P values for all 234 spontaneous epochs (involving 45 cells) that contained at least four bursts. Forty four of these epochs (19%, involving 22 cells) were accompanied by significantly rhythmic bursting (P < 0.05). Of the 22 cells showing rhythmic spontaneous bursting, most (n = 12) did so in only a single epoch, although 3 cells had five rhythmic bursting epochs each. Rhythmicity in spontaneous bursting therefore occurs ~20% of the time in the dLGN cells of a halothane-anesthetized mouse and is restricted to around half of geniculate neurons under such conditions, although a minority of neurons display rhythmic bursting rather often. This level of rhythmicity is considerably lower than the ~47% of rhythmic episodes observed in the naturally sleeping cat ( Weyand et al. 2001Go). The discrepancy could arise because our spontaneous epochs were 5 s long, whereas Weyand et al. based their analysis on "burst bouts" of ~2 s. However, when we also analyzed the rhythmicity of such 2-s burst bouts in our mouse preparation, we observed a reduction in rhythmic episodes: 48 of 438 bouts (11%) displayed significantly rhythmic firing. dLGN cells in anesthetized mice therefore do burst rhythmically but compared with geniculate neurons in naturally sleeping animals do so only rather rarely.



View larger version (14K):
[in this window]
[in a new window]
 
FIG. 3. Rhythmicity of spontaneous bursting in the mouse dLGN. A: spontaneous bursting in a single neuron can alternate between rhythmic and nonrhythmic behavior. Plots show examples of burst timing in spontaneous spikes fired by a single mouse dLGN cell during 7 successive 5-s presentations of a blank gray screen. {bullet}, the time of the 1st spike in a burst. P values shown above each plot are a measure of rhythmicity (see METHODS): values <0.05 represent significantly rhythmic firing. B: spontaneous bursting was rhythmic ~19% of the time. The histogram shows P values calculated from all 234 blank screen presentations (involving 45 cells) that were accompanied by ≥4 bursts. Significant rhythmicity in spontaneous bursting, represented by P values <0.05, was observed in 44 presentations (involving 22 cells).

 
Burst firing is produced by visual stimulation

Neurons in the mouse dLGN fire spikes that can be identified as part of bursts (Fig. 1) and that can be shown to possess features characteristic of burst firing in other species and preparations (Fig. 2). But do these burst spikes encode visual information? A first, crude step toward answering this question comes from a comparison of bursting prevalence during periods of spontaneous and visually driven activity. Every investigation of a dLGN cell's responses to a certain stimulus parameter (known as a "tuning experiment") involved presenting a pseudorandom series of sinusoidal grating stimuli which varied in the parameter under study (spatial frequency, contrast, etc.; see METHODS) and also included the presentation of a blank stimulus at mean luminance. For each such tuning experiment, we calculated the percentage of spikes forming part of bursts (burst %) over all presentations of all grating stimuli (visually driven activity), and over all periods of blank screen presentation (spontaneous activity). We recognize that this approach does not compare matched periods of visually driven and spontaneous spiking, but our use of a percentage measure of bursting prevalence should compensate for this.

Although some cells' visual firing was dominated by burst firing (maximum burst % = 78%; see Fig. 1A), in most experiments bursts occupied a small fraction of all spikes fired (Fig. 4A1). However, burst percentage values for periods of spontaneous activity were usually even lower with many cells not firing any burst spikes in the absence of visual stimulation (Fig. 4A2). Indeed, across all tuning experiments (n = 323 in 102 cells), burst percentage values were significantly higher for visually driven versus spontaneous activity (visual median: 19.8%; spontaneous median: 11.1%; Wilcoxon matched pairs test, P = 0.002; Fig. 4B). Because this analysis was applied across all visual stimuli in a given experiment, not just those producing the highest firing rates, this is a rather clear demonstration that bursting is more prevalent when cells are responding to visual stimuli. This begs the question: exactly what visual information does burst firing carry?



View larger version (23K):
[in this window]
[in a new window]
 
FIG. 4. Mouse dLGN cells fire more burst spikes during visual stimulation than during periods of spontaneous activity. A: distribution of burst percentage values across all tuning experiments (see RESULTS) during visual stimulation (A1) and spontaneous firing (A2). The measure simply represents the percentage of spikes that formed part of bursts. During visual stimulation (A1), tonic firing was predominant in most cells. The sample distribution is skewed toward low burst percentage values, with the median, marked with the clear arrow, at 19.8%. During periods of spontaneous firing to a blank gray screen (A2), however, the sample is even more heavily left-skewed. In most experiments, <5% of all spikes fired spontaneously were part of bursts. The clear arrow marks the median at 11.1%. B: within tuning experiments, burst spikes contribute more to visual responses than to spontaneous firing. {bullet}, a single tuning experiment (n = 323 in 102 cells). While in many cases the burst percentage for spontaneous firing was greater than its visually driven counterpart ({bullet} lying above the - - - of unity), in more cases, including many in which spontaneous activity did not contain any burst spikes, burst percentage was higher during visual stimulation ({bullet} lying below the unity line). As shown in the inset bar plot, the medians of the 2 distributions differed significantly (visual stimulation: 19.8%; spontaneous firing: 11.1%; Wilcoxon matched pairs test, P = 0.002). vs, visual stimulation; sa, spontaneous activity.

 
Burst spikes are fired in response to earlier stimulus phases

Investigations into the relative timing of burst and tonic spikes during presentation of sinusoidal gratings suggest that burst spikes might carry important information about stimulus onset. Figure 5A shows the burst and tonic responses of a mouse dLGN cell to a sinusoidal grating of optimum spatial and temporal frequencies. As often reported by others (e.g., Guido et al. 1992Go, 1995Go; Lu et al. 1992Go) burst spikes, when they occur, tend to occur at the start of a response to a given stimulus cycle. Tonic spikes tend to occur a little later. This precedence of burst firing can have two possible mechanisms: burst spikes can have shorter latencies than tonic spikes or they can be triggered by earlier phases of the stimulus. The relative contribution of each of these candidate mechanisms can be assessed using plots of stimulus temporal frequency (TF) versus cell F1 response phase. The slope of a line relating these two variables can be used as a measure of response latency ( Grubb and Thompson 2003Go; Hawken et al. 1996Go). In addition, the y intercept of the line can give an indication of relative response phase (e.g., Guido et al. 1992Go; Mukherjee and Kaplan 1995Go). In mouse dLGN cells, burst responses to low TFs are usually phase advanced relative to tonic responses to the same stimuli. At high TFs, though, response phases for the two firing modes tend to be similar (Fig. 5B). This makes the TF-F1 phase line steeper for burst firing, meaning that burst spikes possess longer latencies than tonic spikes. Across all TF experiments, burst latencies were significantly longer than the latencies of tonic spikes fired to the same stimuli (burst median: 106 ms; tonic: 90 ms; Wilcoxon matched pairs test, P = 0.001; Fig. 5C). However, burst spikes still occur earlier than tonic spikes because they are triggered by an earlier phase of the stimulus. Phase advances for burst spikes, taken by subtracting tonic from burst y-intercept values in each experiment, were almost all positive, with the sample median being significantly different from zero (more positive values reflect earlier phases; 1-sample t-test, P < 0.0001; Fig. 5D). Despite their longer latency, burst spikes usually carry the first stimulus information relayed to cortex.

Response linearity

When all spikes are taken into consideration, almost all mouse dLGN cells sum inputs in a linear manner across their RFs ( Grubb and Thompson 2003Go). It is unlikely that such linear behavior could have a substantially nonlinear component, so it was reassuring to note that in almost all mouse dLGN cells both burst and tonic spikes show linear spatial summation. The plots in Fig. 6A display typical results of modified null tests ( Hochstein and Shapley 1976Go) applied separately to the burst and tonic components of cells' responses. These tests involve displaying stationary sinusoidal gratings the contrast of which is sinusoidally modulated in time, at various spatial phases across a cell's RF. Null tests were always initially applied with stimuli at a cell's optimal spatial frequency (SF) but were often repeated at higher SFs because nonlinearities of spatial summation can become more prevalent as SF increases (e.g., Hochstein and Shapley 1976Go). Twenty-eight cells were tested at optimal SF, 12 at twice optimal SF, and 3 at three times optimal SF. Temporal frequency of these stimuli was 1 Hz (11 cells) or optimal (17 cells), maximum stimulus contrast was 70%. By convention, fundamental (F1) response amplitudes are represented as negative in null test plots when their phase differs by 90–270° from that of the maximum response; this explains why burst and tonic plots from the same cell sometimes look like mirror images of each other.

In cells that sum inputs linearly across their RFs, responses to counterphased stimuli should be largest at the fundamental (F1) stimulus harmonic and the amplitude of this harmonic should vary sinusoidally with stimulus phase ( Hochstein and Shapley 1976Go). In addition, a linearly summating cell should possess two null phases where excitatory and inhibitory inputs are perfectly balanced across the RF ( Enroth-Cugell and Robson 1966Go). These null phases are represented by zero crossings in a null test plot ( Hochstein and Shapley 1976Go). In almost all cases in our mouse dLGN sample (27/28 cells, 42/43 null tests, Fig. 6A, 1 and 2), the F1 components of both burst and tonic firing were largest, varied sinusoidally in amplitude with stimulus phase, and crossed the zero line at 2 null locations. Such behavior was accompanied by linear "linearity" values (see METHODS) ( Hochstein and Shapley 1976Go) of <1.

One cell showed evidence of SF-dependent nonlinear spatial summation but did so in both response modes (Fig. 6A3). In both burst and tonic firing, although appearing to possess two null phases, the F2 response component of this neuron to stimuli at twice optimal SF was generally larger than the F1 component and varied sinusoidally with stimulus phase. Such a response pattern is characteristic of a RF with overlapping, nonantagonistic On and Off subregions ( Lennie and Perry 1981Go). While this cell displayed evidence of such ON-OFF nonlinearities in both its burst and tonic firing (Fig. 6A3), a high maximum F1 response in its burst firing meant that only its tonic component possessed a linearity value greater than unity, at 1.64.

Across all null test experiments, linearity values of burst spikes were not significantly different from those of tonic spikes fired to the same stimuli (burst means ± SE: 0.42 ± 0.02; tonic: 0.37 ± 0.04; paired t-test, P = 0.17; Fig. 6B). Both burst firing and tonic firing in the responses of mouse dLGN cells therefore sum spatial inputs in an equally linear manner.

Although burst and tonic firing do not differ in terms of their linearity of spatial summation, they do differ on a very different measure of response linearity, rectification. The response of a linear neuron to a sinusoidal grating stimulus should be sinusoidal, with a dominant F1 component (e.g., Shapley and Lennie 1985Go). However, even when a neuron displays wholly linear spatial summation its response to a sinusoid may not be sinusoidal because of response rectification: when spontaneous activity is low and visual responses strong and transient, cells cannot decrease their responses to the nonpreferred half cycle of a given grating, and may not precisely follow the temporal form of the grating's preferred half cycle. Such response rectification acts to increase all nonfundamental harmonics in a cell's responses. Because burst spikes tend to be high-frequency and transient, one might expect their rectification to be greater. Indeed, over all tuning experiments (n = 332 in 102 cells), the "rectification" of burst spikes (Eq. 1) was significantly greater than that of tonic spikes fired to the same stimuli (burst median: 0.46; tonic: 0.36; Wilcoxon matched pairs test, P < 0.0001; Fig. 6C). Thus although both burst and tonic spikes display linear spatial summation, burst spikes are significantly more nonlinear in terms of response rectification.

Spatial frequency tuning

As well as being identical with regard to spatial summation, burst and tonic spikes in the mouse dLGN are no different with regard to their SF tuning properties. SF tuning was assessed by fitting separate Difference of Gaussians (DoG) curves (Eq. 2) to the F1 amplitudes of burst and tonic responses to sinusoidal gratings of various SFs. Temporal frequency of these stimuli was 1 Hz, contrast was 70%. The DoG equation ( Rodieck 1965Go, So and Shapley 1981Go) models the center and surround mechanisms of a cell as symmetrical, antagonistic Gaussian functions, and can provide measures of the strengths and sizes of a given cell's center and surround RF regions. Because a DoG curve also provides a smooth fit to raw data, it can also be used to precisely calculate a cell's SF optimum (or peak) and cutoff (taken as the high SF at which response amplitude decayed to 1% of its maximum). However, due to the large size of RFs in the mouse dLGN, coupled with the limited size of our stimulus display at appropriate viewing distances, we were unable to accurately describe a cell's full low-SF roll-off (see Grubb and Thompson 2003Go).

As shown in Fig. 7A, DoG curves fitted to burst and tonic responses sometimes differed in their relative amplitudes (though not consistently in favor of either burst or tonic firing), but were usually similar in shape. We included in our analysis all cells (n = 61) in which the DoG function accounted for ≥80% of the variance of both burst and tonic raw data. In this sample, there were no significant differences in spatial tuning between the two response modes. Burst and tonic spikes fired to the same stimuli did not differ in terms of kc, the strength of the RF center (burst median: 7.59; tonic: 13.42; Wilcoxon matched pairs test, P = 0.43), ks, the relative strength of the RF surround (burst median: 0.98; tonic: 0.981; Wilcoxon matched pairs test, P = 0.06), rc, the radius of the RF center (burst median: 5.0°; tonic: 4.5°, Wilcoxon matched pairs test, P = 0.38), rs, the radius of the RF surround (burst median: 13.6°; tonic: 17.1°, Wilcoxon matched pairs test, P = 0.96), peak SF (burst means ± SE: 0.033 ± 0.003c/°; tonic: 0.032 ± 0.002c/°; paired t-test, P = 0.64; Fig. 7B), or SF cutoff (burst means ± SE: 0.17 ± 0.01; tonic: 0.16 ± 0.01; paired t-test, P = 0.56; Fig. 7C).



View larger version (24K):
[in this window]
[in a new window]
 
FIG. 7. Spatial frequency tuning of burst and tonic spikes in the mouse dLGN. A: examples of burst and tonic spatial frequency (SF) tuning curves in 3 neurons. {bullet} and {circ}, show F1 response amplitudes to drifting sinusoidal gratings of various SFs; —, the best fit to these data of a Difference of Gaussians model (see METHODS). Note the different y-axis scales for the 3 examples, and the logarithmic x axis. Across these examples, and the mouse dLGN sample as a whole, the shape of SF tuning curves tended to be very similar for burst and tonic components of a given cell's responses. In some cells, burst responses were stronger than tonic responses (A1). More commonly, tonic firing was stronger than burst firing (A2). In other cells, burst and tonic tuning curves were similar in amplitude as well as shape (A3). Spatial tuning parameters for all 3 cells: cell 1: burst firing: kc 18.4, ks 0.67, rc 2.8°, rs 4.0°, peak 0.063 c/°, cutoff 0.27c/°; tonic firing: kc 6.7, ks 0.70, rc 2.9°, rs 6.0°, peak 0.064 c/°, cutoff 0.25 c/°; cell 2: burst firing: kc 7.4, ks 0.69, rc 3.8°, rs 17.0°, peak 0.031 c/°, cutoff 0.18 c/°; tonic firing: kc 12.4, ks 1.0, rc 2.6°, rs 22.9°, peak 0.029 c/°, cutoff 0.26 c/°; cell 3: burst firing: kc 5.2, ks 0.96, rc 3.8°, rs 12.4°, peak 0.041 c/°, cutoff 0.19 c/°; tonic firing: kc 4.3, ks 1.20, rc 3.5°, rs 14.0°, peak 0.04 c/°, cutoff 0.20 c/°. B: burst and tonic firing modes do not have different peak SFs. {bullet}, single cell, and are just as prevalent above the - - - of unity as below it. Inset: the means of the 2 distributions do not differ significantly (burst means ± SE: 0.033 ± 0.003 c/°; tonic: 0.032 ± 0.002 c/°; paired t-test, P = 0.64). Error bars show SE. C: spatial acuity is no different for burst and tonic spikes. {bullet}, 1 cell each and are evenly spaced about the - - - of unity. Inset: no significant difference between the means of the burst and tonic distributions (burst means ± SE: 0.17 ± 0.01; tonic: 0.16 ± 0.01; paired t-test, P = 0.56). All conventions as in B.

 
Temporal frequency tuning

Burst and tonic spikes differ with respect to their temporal tuning. Temporal frequency (TF) tuning was assessed by presenting drifting sinusoidal gratings of various TFs at optimal SF and 70% contrast. F1 response amplitudes of burst and tonic firing to these stimuli were then separately fitted with two-half-Gaussian equations (Eq. 3). These functions tended to be narrower, peak at lower TFs, and fall off at lower high TFs for burst versus tonic firing (Fig. 8A). We included in our analysis all cells (n = 42) in which the two-half-Gaussians equation accounted for ≥80% of the variance of both burst and tonic raw data. In this sample, burst peak TFs were significantly lower than those of tonic spikes fired to the same stimuli (burst means ± SE: 3.9 ± 0.2 Hz; tonic: 4.5 ± 0.3 Hz; paired t-test, P = 0.007; Fig. 8B). Burst high50 values, taken as the high TF at which response amplitude dropped to half of its maximum, were also significantly lower than their tonic counterparts (burst means ± SE: 6.0 ± 0.3 Hz; tonic: 8.0 ± 0.5 Hz; paired t-test, P < 0.0001; Fig. 8C). To assess the tightness of TF tuning, a bandwidth measure was introduced. This measured, in octaves, the width of a TF tuning curve at half-maximum height. The bandwidth of burst spikes was significantly lower than that for tonic firing (burst means ± SE: 1.6 ± 0.1 octaves; tonic: 2.3 ± 0.2 octaves; paired t-test, P = 0.002; Fig. 8D). However, the low-frequency roll-off of burst and tonic TF tuning curves was not different. The measure low50, the low TF at which a curve reached half its maximum height, did not differ significantly between burst and tonic activity (burst means ± SE: 2.2 ± 0.2 Hz; tonic: 1.8 ± 0.3 Hz; paired t-test, P = 0.25).



View larger version (26K):
[in this window]
[in a new window]
 
FIG. 8. Temporal frequency tuning of burst and tonic spikes in the mouse dLGN. A: examrples of burst and tonic TF tuning for 3 different cells (A, 1–3). {bullet} and {circ}, F1 response amplitudes to drifting sinusoidal gratings of various TFs; —, the best fit to these data of 2-Gaussian-halves equations (see METHODS). Note the different y-axis scales for each example and the logarithmic x axis. Peak responses can occur at lower TFs for burst firing (A2). In all 3 examples, cells burst spikes are more tightly tuned for TF and possess lower temporal acuity than their tonic counterparts. Temporal tuning parameters for these example cells: cell 1: burst firing: peak 5.4 Hz, high50 7.2 Hz, low50 3.4 Hz, bandwidth 1.1 octaves; tonic firing: peak 5.3 Hz, high50 9.6 Hz, low50 0.9 Hz, bandwidth 3.4 octaves; cell 2: burst firing: peak 3.5 Hz, high50 5.2 Hz, low50 1.8 Hz, bandwidth 1.6 octaves; tonic firing: peak 5.3 Hz, high50 12.7 Hz, low50 2.0 Hz, bandwidth 2.7 octaves; cell 3: burst firing: peak 2.6 Hz, high50 4.3 Hz, low50 1.2 Hz, bandwidth 1.9 octaves; tonic firing: peak 2.8 Hz, high50 8.3 Hz, low50 0.6 Hz, bandwidth 3.7 octaves. B: burst spikes prefer lower TFs. {bullet}, 1 cell, more concentrated above the (- - -) unity line. Inset: that the means of the burst and tonic distributions are significantly different (burst means ± SE: 3.9 ± 0.2 Hz; tonic: 4.5 ± 0.3 Hz; paired t-test, P = 0.007). Error bars show SE. C: tonic firing can signal higher TFs than burst firing. {bullet}, cells, which are most concentrated above the unity line. Inset: the mean tonic high50 value is significantly greater than the corresponding burst firing mean (burst means ± SE: 6.0 ± 0.3 Hz; tonic: 8.0 ± 0.5 Hz; paired t-test, P < 0.0001). All conventions as in B. D: tonic firing signals a wider range of TFs than burst firing. {bullet}, 1 cell and predominantly occur above the line of unity. That the tonic firing mean bandwidth is significantly greater than the burst firing mean bandwidth is confirmed by the inset (burst means ± SE: 1.6 ± 0.1 octaves; tonic: 2.3 ± 0.2 octaves; paired t-test, P = 0.002). All conventions as in B.

 
It is possible that some of the preceding differences in the temporal tuning of burst and tonic spikes are due not to fundamental differences between the two firing modes but instead to the stringent criteria used to identify burst spikes. Bursts as we define them must be preceded by a silent period of ≥100 ms (see METHODS). For this reason, spikes triggered by LTS activity might go undetected as bursts when, as is the case at high TFs, the period of sinusoidal grating stimuli approaches or is <100 ms. This is unlikely to affect any other data presented in this paper because stimuli in all non-TF tuning experiments were presented at 1 Hz or peak TF—always with a stimulus period well above 100 ms. However, it could affect our TF tuning results. We therefore investigated the effects of relaxing the requisite preburst silent period to 50 ms for stimuli with TF ≥8 Hz, following a precedent set by Sherman and colleagues ( Guido et al. 1992Go; Lu et al. 1992Go). Surprisingly, this had remarkably little effect on burst-tonic differences in temporal tuning. With bursts preceded by ≥50-ms silence when stimulus TFs were ≥8 Hz, we still saw significant differences between burst and tonic spikes in high50 (burst means ± SE: 5.9 ± 0.3 Hz; tonic: 7.6 ± 0.4 Hz; paired t-test, P = 0.0003) and bandwidth (burst means ± SE: 1.6 ± 0.1 octaves; tonic: 2.3 ± 0.2 octaves; paired t-test, P = 0.002), and there remained no burst-tonic difference in low50 (burst means ± SE: 2.0 ± 0.2 Hz; tonic: 1.6 ± 0.2 Hz; paired t-test, P = 0.16). In addition, response timing results, obtained using data collected during TF tuning experiments (see Fig. 5) were also unchanged: bursts were still phase-advanced (mean advance ± SE: 1.3 ± 0.2 rad; 1-sample t-test vs. 0, P < 0.0001) with a longer latency (burst median: 107 ms; tonic: 90 ms; Wilcoxon matched pairs test, P = 0.0008). The only different result noted when the preburst silent period was relaxed at high TFs was that peak TF was no longer significantly different between burst and tonic firing (burst means ± SE: 4.0 ± 0.3 Hz; tonic: 4.4 ± 0.3 Hz; paired t-test, P = 0.11). Note, however, that the trend for a higher TF peak in tonic firing remained. It therefore appears that the large majority, perhaps all, of burst-tonic differences in temporal tuning stem from fundamental differences in stimulus response, rather than from the ≥100-ms silence criterion used to identify burst firing in the present study.

Contrast response characteristics

The contrast response functions of burst and tonic spikes were assessed by presenting drifting sinusoidal gratings of various contrasts at optimal SF. TF of these stimuli was either 1 Hz (18 cells) or optimal (14 cells). F1 response amplitudes of burst and tonic firing were then separately fitted with a hyperbolic function (Eq. 4) ( Albrecht and Hamilton 1982Go) and measures were extracted from the portion of this curve between 0 and 100% contrast ( Grubb and Thompson 2003Go). In most cases, both burst and tonic response amplitude grew rather gradually with increasing stimulus contrast (Fig. 9A1) or at least changed in a similar manner as contrast grew stronger (Fig. 9A3). In a few cases, though, a step-like function was found only in a cell's burst responses (Fig. 9A2). Contrast gain was calculated as the slope of a tangent to the curve where response amplitude was 20% of its value at 100% contrast. We included in our analysis all cells (n = 32) in which the hyperbolic function accounted for ≥80% of the variance of both burst and tonic raw data. In this sample, a trend was observed toward higher gain values for tonic firing. This trend, however, was not significant (burst median: 0.26 spikes · s–1 · %–1; tonic: 0.28 spikes · s–1 · %–1; Wilcoxon matched pairs test, P = 0.06; Fig. 9B). Neither was a significant difference found between burst and tonic spikes in terms of c50, the contrast at which response amplitude reached 50% of its value at 100% contrast (burst means ± SE: 33.5 ± 3.2%; tonic: 36.8 ± 3.2%; paired t-test, P = 0.33; Fig. 9C). However, because in vitro experiments have suggested that the contrast response function for burst spikes might be more step-like than that of tonic spikes ( McCormick and Feeser 1990Go; Ramcharan et al. 2000aGo; Zhan et al. 1999Go), the two firing modes were also compared with respect to the fitted curve parameter n (see Eq. 4). The larger this rate-of-change variable, the more step-like a curve is. n was significantly higher for burst spikes than for tonic spikes fired to the same stimuli (burst median: 3.34; tonic: 1.86; Wilcoxon matched pairs test, P = 0.0005; Fig. 9D).



View larger version (24K):
[in this window]
[in a new window]
 
FIG. 9. Contrast response characteristics of burst and tonic firing in the mouse dLGN. A: examples of contrast response functions for burst and tonic spikes in 3 different cells (A, 1–3). {bullet} and {circ}, F1 response amplitudes of burst and tonic spikes to each contrast; —, the best hyperbolic fit to these data. Note the different y-axis scales. Contrast response curves for both burst and tonic firing tended to show gradual increases in F1 amplitude with increasing contrast (A, 1 and 3). General curve shape and initial curve gradient were often very similar for the 2 modes of firing (A, 1 and 3). Only rarely were sharp step functions seen solely in the bursting data (A2). Contrast response parameters for the 3 units: cell 1: burst firing: gain 0.33 spikes · s–1 · %–1, c50 38%, n 3.7; tonic firing: gain 0.43 spikes · s–1 · %–1, c50 47%, n 2.5; cell 2: burst firing: gain 0.49 spikes · s–1 · %–1 c50 13%, n 3.1; tonic firing: gain 0.06 spikes · s–1 · %–1, c50 56%, n 1.2; cell 3: burst firing: gain 0.28 spikes · s–1 · %–1, c50 15%, n 2.6; tonic firing: gain 0.19 spikes · s–1 · %–1, c50 13.7%, n 2.6. B: no significant difference between contrast gain for burst and tonic firing. Although {bullet}, each representing 1 cell, appear most concentrated above the - - -, the inset bar plot shows that the burst and tonic median gains did not differ significantly (burst median: 0.26 spikes · s–1 · %–1; tonic: 0.28 spikes · s–1 · %–1; Wilcoxon matched pairs test, P = 0.06). C: c50 values are no different for burst and tonic spikes. {bullet} are scattered rather widely but equally about the unity line. Inset: the means of the 2 distributions do not differ significantly (burst means ± SE: 33 ± 3%; tonic: 37 ± 3%; paired t-test, P = 0.99). Error bars show SE. D: contrast response curves are more step-like for burst firing. Note the logarithmic axes in the scatter plot, in which each {bullet} represents 1 cell. {bullet} are clearly concentrated below the unity line, and the inset bar plot shows that the median n values (see Eq. 4) of the 2 firing modes are significantly different (burst median: 3.3; tonic: 1.9; Wilcoxon matched pairs test, P = 0.0005). Conventions as in B.

 
Signal detection

The more step-like contrast response functions of burst spikes (Fig. 9D), as well as their strong rectification (Fig. 6C), are consistent with the hypothesis that burst firing in an individual dLGN cell provides more information about whether a given stimulus is present than about that stimulus' fine visual details (e.g., Sherman and Guillery 1996Go). In this hypothesis, burst spikes are primarily important for detecting visual stimuli. The relationship between burst firing and stimulus detection can be directly tested using receiver operating characteristic (ROC) analysis ( Green and Swets 1966Go; Guido et al. 1995Go; MacMillan and Creelman 1991Go). This approach assesses how well a cell's responses discriminate visual stimulation from periods of spontaneous activity, without ever employing a fixed detection criterion. This is achieved through the construction of ROC curves (Fig. 10A), which for a given stimulus plot the probability of observing different response sizes during spontaneous activity, P (false alarm), against the probability of observing those response sizes during visual stimulation, P (hit) (see METHODS). The curve always runs from (0,0) to (1,1), and if real signal detection is occurring all points lie above the unity line. Area underneath the ROC curve (ROC area) can then be used as a nonparametric measure of stimulus detection ( Green and Swets 1966Go; Guido et al. 1995Go; MacMillan and Creelman 1991Go).



View larger version (12K):
[in this window]
[in a new window]
 
FIG. 10. Levels of bursting correlate with stimulus detection in the responses of mouse dLGN cells. A: examples of receiver operating characteristic (ROC) analysis applied to all spikes (burst and tonic) in the responses of 3 individual mouse dLGN cells. Each plot is generated from the responses of a cell to a single sinusoidal grating stimulus. In each case here, the stimulus was nonoptimal, in that, of all stimuli presented during a tuning experiment, it did not produce the largest-amplitude F1 response. For a given spike count, 1 point in an ROC curve plots the probability P (hit) of attaining this count during a presentation of a visual stimulus vs. the probability P (false alarm) of attaining this count during an equal period of spontaneous activity. The area underneath the curve, expressed as a value between 0 and 1, then gives a nonparametric measure of a cell's ability to distinguish visual from nonvisual stimulation ( Guido et al. 1995Go) (see METHODS). ROC area values for each curve are shown in the bottom right of each plot, along with a burst index representing the percentage of stimulus cycles that produced ≥1 burst response. In these 3 example stimulus presentations, it is clear that higher levels of bursting are coupled with better stimulus detectability. B: burst index values are positively and significantly correlated with ROC area values (Spearman's r = 0.34, P = 0.0005, n = 102). Although both measures were calculated for every stimulus presented to a given mouse dLGN cell (see A), the plot here shows a single burst index and ROC area for each cell, produced by averaging over all sinusoidal grating stimuli presented to that neuron. The positive and significant relationship is also observed when correlations are calculated over every sinusoidal grating stimulus, over all optimal stimuli, over all nonoptimal stimuli, and over all stimuli that produce ROC areas >0.5, representing significant signal detection (see text).

 
Given our experimental approach, we were unable to compare ROC curves directly for burst and tonic spikes fired to the same stimuli. This can be readily done when responses are recorded intracellularly ( Guido et al. 1995Go) but cannot be carried out with confidence when data are obtained from extracellular recordings. This is because the construction of ROC curves requires a direct comparison of spontaneously and visually driven activity, a comparison that is considerably complicated by recording periods in which no spikes are fired. During intracellular recordings, it is a simple matter to identify the underlying response mode of a neuron based on its membrane potential. An intracellularly recorded lack of spontaneous firing, for example, can then be easily categorized as a lack of either burst or tonic firing and can be compared with a visually driven response (or lack of response) that occurred while the cell was in the same firing mode. When responses are recorded extracellularly, though, a response failure during, for example, visual stimulation, cannot similarly be compared with the appropriate period of spontaneously driven activity. With extracellular recording techniques therefore, we were not able to effectively compare spontaneous and visually driven activity in burst and tonic firing. Furthermore, we felt we could not simply ignore all stimulus cycles that produced a lack of burst or tonic activity because such cycles constituted a rather large proportion of all stimulus presentations: 75% of all stimulus cycles did not evoke a burst response, whereas 22% of all stimulus cycles did not evoke a single tonic spike.

However, we were able to assess the contribution of burst firing to stimulus detection more indirectly by asking whether the latter improved as the former became more prominent (cf. Guido et al. 1995Go). For each sinusoidal grating stimulus presented to each cell, we therefore calculated ROC area based on the total complement (burst + tonic) of the cell's responses. Assigning response failures to a particular firing mode was therefore not an issue. We also calculated a "burst index" for each stimulus defined as the proportion of all stimulus cycles generating at least one burst response ( Guido et al. 1995Go). The example ROC curves in Fig. 10A show that ROC area, and thus signal detection ability, tended to be larger when the burst index was higher. Indeed, over all stimuli presented to all cells, we observed a positive and significant correlation between ROC area and burst index (Spearman r = 0.59; P < 0.0001; n = 3,040). This relationship was similarly positive and significant when assessed only for the optimal stimulus in each tuning experiment (Spearman r = 0.54; P < 0.0001; n = 339), only for the nonoptimal stimuli in each tuning experiment (Spearman r = 0.58; P < 0.0001; n = 2701), or only for those stimuli associated with ROC areas >0.5, representing significant stimulus detection (Spearman r = 0.54; P < 0.0001; n = 2347). The correlation between ROC area and burst index was also positive and significant if, instead of using one ROC area and one burst index per stimulus, we averaged ROC area and burst index values across all stimuli presented to each cell (Spearman r = 0.34; P =0.0005; n =102; Fig. 10B). Indirect evidence therefore suggests that burst firing is related to greater signal detection in the responses of mouse dLGN cells.

Reliability in spike timing

If burst firing tends to occur first in a dLGN cell's responses to visual stimuli (Fig 5), does it provide reliable information concerning the time of stimulus onset? Although reliability in the absolute timing of spikes from stimulus onset cannot be assessed using the cyclic sinusoidal grating stimuli presented here, this question can be answered by taking a fixed point in the stimulus cycle and assessing reliability in the relative timing of spikes from that point. Reliability in burst and tonic spike timing was first assessed by calculating the mean deviation of spike times from the appropriate F1 response phase, normalized by the stimulus period (Eq. 6, Fig. 11A ). Note that higher timing reliability here is reflected in lower values of the reliability measure. Over all grating stimuli, burst spikes occupy less of the stimulus period than tonic spikes (burst reliability median: 0.07; tonic: 0.13; Wilcoxon matched pairs test, P < 0.0001; Fig. 11B). An individual burst spike thus informs the brain more about stimulus timing than an individual tonic spike fired to the same grating. In addition, this advantage in timing reliability for burst spikes increases with increasing response size. While both burst and tonic spike times become less variable as overall F1 response size moves closer to optimal (burst Spearman r = –0.33, P < 0.0001, n = 2,178; tonic Spearman r = –0.53, P < 0.0001, n = 2,178), the ratio of burst to tonic reliability also becomes lower (i.e., burst spikes are relatively less spread than tonic spikes) with increasing response optimality (Spearman r = –0.15, P < 0.0001, n = 2,178).



View larger version (27K):
[in this window]
[in a new window]
 
FIG. 11. Reliability of burst and tonic spike timing in mouse dLGN cells. A: assessing reliability in the timing of burst and tonic firing. Raster plots, in which each {circ} represents 1 spike, show response timing for all burst and tonic spikes fired to a drifting sinusoidal grating (A1) or for just the 1st burst and tonic spikes fired to each cycle of the stimulus (A2). Note the close superposition of burst spikes in A1 that cause this plot to look very similar to the 1st-spike plot in A2. In A1, the vertical line represents the phase of the burst or tonic F1 component, as appropriate. Reliability values, displayed at bottom left, represent the mean deviation of spike times from this line, normalized by the stimulus period. Note that higher timing reliability is reflected in lower values of this measure. In A2, —, the start of the half stimulus period centered on the phase of the burst or tonic F1 component, as appropriate. First-spike reliability values, shown at bottom right, represent the SD of 1st-spike time distances from this line (see METHODS). Again, higher timing reliability is reflected in lower values of this measure. As reflected in these values, burst spikes show more reliability in timing when all spikes of a particular firing mode are taken into consideration. However, when just the first spikes fired to each stimulus cycle are analyzed, there is very little difference in timing reliability between burst and tonic firing. B: over all spikes, burst timing is more reliable. {bullet}, which each represent a single sinusoidal grating stimulus, are concentrated above the - - - unity line. Tonic spikes thus occur over a larger portion of a cycle than burst spikes fired to the same stimulus. This is confirmed in the inset, which shows that the burst and tonic firing median reliability values are significantly different (burst median: 0.07; tonic: 0.13; Wilcoxon matched pairs test, P < 0.0001). C: the 1st tonic spikes fired to each stimulus cycle are more reliable in their timing then the 1st burst spikes fired to the same stimuli. Burst and tonic 1st-spike reliability values appear well correlated with each other, and {bullet} appear rather evenly distributed about the unity line. Nevertheless, the inset shows that median 1st-spike reliability value for tonic firing is significantly smaller than that for burst spikes (burst median: 70 ms; tonic: 64 ms; Wilcoxon matched pairs test, P < 0.0001). Conventions as in B. D: burst and tonic spike timing becomes more reliable with increasing stimulus contrast. The plot shows median reliability values for each stimulus contrast, averaged over all contrast response experiments. E: burst and tonic 1st-spike timing also becomes more reliable with increasing stimulus contrast. Conventions as in D.

 
In terms of the timing of the first spikes fired to grating cycles, however, tonic firing is more reliable than burst activity. First-spike reliability was measured by taking the SD of first-spike time latencies from the start of a half-cycle period centered on the burst or tonic F1 response phase, as appropriate (Eq. 7, Fig. 11A). Again, higher timing reliability is reflected in lower values of this reliability measure. Although this measure produced very similar reliability estimates for tonic and burst firing (e.g., Fig. 11A), over all grating stimuli tonic first spikes were less jittery in time than burst first spikes fired to the same stimuli (burst median: 70 ms; tonic: 64 ms; Wilcoxon matched pairs test, P < 0.0001; Fig. 11C). This reliability advantage of tonic firing, however, decreased as stimuli became more optimal. First-spike reliability values for both burst and tonic firing decreased as overall F1 response amplitudes approached those produced by optimal stimuli (burst Spearman r = –0.45, P < 0.0001, n = 1,283; tonic Spearman r = –0.48, P < 0.0001, n = 1,283), but the ratio of burst to tonic first-spike reliability values also decreased (i.e., burst 1st-spike reliability became relatively better) with increasing stimulus optimality (Spearman r = –0.19, P < 0.0001, n = 1,283). So, although tonic first spikes provide more accurate information concerning stimulus "onset," this advantage over burst first spikes becomes less marked for more optimal stimuli.

This relationship between response strength and reliability in spike timing was investigated in more detail for contrast response characteristics because of a specific prediction arising from an in vitro electrophysiological study. Zhan et al. (1999)Go found that the latency of the first burst action potential caused by a low-threshold calcium spike was longer after weaker depolarizing current injections. Presumably assuming that longer latencies are associated with greater variability in response timing, these authors predicted that increasing stimulus contrast in vivo would be met with increased reliability in burst spike timing. Indeed, the present data show that this is the case. Across all stimuli in all contrast experiments, stimulus contrast is negatively and significantly correlated with burst reliability (Spearman r = –0.22, P < 0.0001, n = 652) and burst first-spike reliability (Spearman r = –0.29, P < 0.0001, n = 405; see Fig. 11, D and E). And within contrast experiments, burst reliability and first-spike reliability values were both higher for the highest contrast for which they could be measured versus the lowest contrast for which they could be measured (reliability: low contrast means ± SE: 0.11 ± 0.01; high contrast: 0.07 ± 0.007; paired t-test, P < 0.0001; 1st-spike reliability: low contrast means ± SE: 126 ± 15 ms; high contrast: 100 ± 15 ms; paired t-test, P = 0.01). However, increasing reliability with increasing stimulus contrast is not a feature restricted only to burst firing. Across all contrast stimuli, the variability of tonic spike timing also decreased with increasing contrast (reliability: Pearson r = –0.43, P < 0.0001, n = 691; 1st-spike reliability: Pearson r = –0.34, P < 0.0001, n = 792; see Fig. 11, D and E). Within contrast experiments, too, the highest contrasts for which reliability values could be measured produced significantly more reliable tonic spike timing than the lowest contrasts for which these values could be obtained (reliability: low contrast median: 0.20; high contrast: 0.11; Wilcoxon matched pairs test, P < 0.0001; first-spike reliability: low contrast means ± SE: 154 ± 11 ms; high contrast: 70 ± 7 ms; paired t-test, P < 0.0001). Thus although the reliability of burst spike timing increases as stimulus contrast increases, as predicted by Zhan et al. (1999)Go, tonic firing reliability becomes better with increasing contrast as well.

Analyses treating bursts as unitary events

In the preceding analyses, visual response properties are compared between spikes that form part of bursts and spikes that do not form part of bursts. We extended these analyses by asking what happens if instead of treating bursts as events made up of component spikes, bursts are analyzed as unitary events. In other words, are the same burst-tonic functional differences (and similarities) observed if we analyze the response properties, not of burst spikes, but of bursts per se? To address this issue, bursts were treated as unitary events, occurring at the poststimulus time of their first spike. Using this approach, we found an almost identical pattern of burst-tonic functional relationships to that outlined in the spike-based analyses described in the preceding text. As was the case when burst and tonic spikes were compared, burst events, in terms of bursts per second, were not significantly different (P > 0.05) from tonic spikes with respect to linearity of spatial summation, spatial frequency tuning, and contrast c50. Also in line with the preceding analyses, bursts as unitary events had higher visual than spontaneous firing rates, had advanced response phases but longer latencies than tonic firing, were more rectified than tonic spikes, had smaller ROC areas than tonic spikes, had more step-like contrast response curves than tonic firing, and, in their TF responses, had lower peak, high50, and bandwidth values (all P < 0.05). The only difference between analyses based on burst spikes versus those based on bursts per se occurred for the TF tuning measure low50—the low TF at which responses fell to half their maximal height. While burst spikes show a nonsignificant trend toward higher low50 values than tonic spikes (see preceding text), bursts as unitary events do possess significantly higher low50 values compared with tonic firing (burst median: 1.8 Hz; tonic: 1.4 Hz; Wilcoxon test, P = 0.03). This very small difference between two separate ways of assessing burst response characteristics suggests that the burst-tonic functional differences outlined in the preceding text are rather robust phenomena.

Burst-tonic differences in ON-center and OFF-center cells

In our previous description of visual response properties in the mouse dLGN ( Grubb and Thompson 2003Go), we looked for evidence of parallel processing of visual information by different cell types and found it only in functional differences between cells that respond primarily to increases in stimulus luminance (ON center) and cells that respond primarily to decreases in stimulus luminance (OFF center). It was therefore worth asking whether the burst-tonic differences in visual processing described in the preceding text occurred predominantly and/or consistently in either of the two main RF types in the dLGN.

For those response property parameters, which were normally distributed or which could be rendered normal by simple logarithmic transforms, we employed repeated-measures ANOVA tests, in which firing mode (burst/tonic) was the within-subject variable and center type (ON-OFF) was the between-subject variable. Any significance in the firing mode x center type interaction would then imply that a given burst-tonic difference was significantly stronger in cells of one or other center type. This ANOVA analysis was possible for 16 response parameters: latency, phase, linearity, rectification, kc, rc, rs, SF peak, SF cutoff, TF peak, high50, bandwidth, low50, contrast gain, c50, and contrast n. In most cases, the firing mode x center type interaction was not significant (P > 0.05). In four cases, however, the interaction was significant (P < 0.05), suggesting that burst-tonic differences in these parameters were stronger in either ON- or OFF-center neurons. Burst-tonic differences in response phase and rectification were significantly larger in ON-center cells, while differences between the two firing modes in latency and high50 were significantly greater in OFF-center cells. Still, the fact that burst-tonic comparisons were not significantly different between ON- and OFF-center cells for the majority of visual response measures, coupled with the observation that the four significant effects described in the preceding text were split evenly between ON- and OFF-center neurons, show that the burst-tonic differences described in the present paper are not due primarily to firing mode distinctions within cells of one single center type.

This conclusion is backed by consideration of those parameters that could not be rendered normally distributed by simple transforms. These response measures were compared by simply noting the results of paired burst-tonic comparisons applied to ON- and OFF-center data separately. We found positive and significant correlations between ROC area and burst index in both ON- and OFF-center neurons, observed that spike timing reliability is significantly better in the burst spikes fired by both ON- and OFF-center cells and noted that both center types possess no significant difference between burst and tonic spikes with respect to the RF surround parameter ks. Only in the remaining 2 response measures did burst-tonic differences appear to be dominated by cells of one particular center type. Like our dLGN sample as a whole, ON-center cells showed a significantly greater percentage of burst spikes in visually driven versus spontaneous firing. OFF-center cells, however, did not. Conversely, while OFF-center cells mirrored our overall results in that their burst first-spike reliability was significantly worse than that of tonic spikes, this effect was not significant in ON-center cells. Still, one burst-tonic difference due mainly to the responses of ON-center cells, and one due mainly to the responses of OFF-center cells, along with the ON-OFF similarities described in the preceding text, suggests that, overall, burst-tonic differences in the visual response properties of mouse dLGN cells are not due predominantly to the responses of cells of a single center type.

Why are burst spikes functionally different from tonic spikes?

Bursts possess two properties that distinguish them from the large majority of tonic spikes: high-frequency firing and a long preceding silent period. The possibility exists that the differences between burst and tonic spikes outlined in the preceding text might be due primarily to one or other of these distinguishing features—for example, differences in temporal tuning could be caused purely by the long silent period required for burst spikes to occur rather than being caused by LTS activity per se. To investigate this possibility, we compared the visual response properties of burst spikes with three types of tonic responses, all identified, like burst spikes themselves (see METHODS), on the basis of temporal patterns in spike firing. High-frequency tonic (HFT) spikes were separated by short interspike intervals of ≤4 ms, but, unlike bursts, groups of HFT spikes were not preceded by ≥100 ms of nonspiking dead time. In contrast, gap spikes were preceded by ≥100 ms of silence but unlike bursts had following ISIs of >4 ms. The final group, TX spikes, comprised the remainder of nonburst firing, and represented normal tonic activity (Fig. 12A). By comparing the visual response properties of these different types of activity, we hoped to glean some understanding of the underlying reasons for burst-tonic functional differences. For instance, the 4-way analysis might find that in the majority of cases where burst spikes differ from tonic spikes, HFT spikes differ from TX spikes in the same direction as the burst-tonic difference, but gap spikes do not. This might then indicate that the functional distinctions between burst and tonic responses outlined above are due mainly to the high-frequency firing component of bursts, rather than the requisite preceding silent period.



View larger version (38K):
[in this window]
[in a new window]
 
FIG. 12. Assessing possible causes of burst-tonic functional differences in the mouse dLGN. A: in 4-way analyses, burst firing was compared with 3 distinct subgroups of tonic spikes. High-frequency tonic (HFT) spikes were separated by short interspike intervals of ≤4 ms, but, unlike bursts, each cluster of HFT spikes was not preceded by ≥100 ms of nonspiking dead time. "Gap" spikes were preceded by ≥100 ms of silence, but unlike bursts had following ISIs of >4 ms. TX spikes comprised the remainder of nonburst firing, and represented "normal" tonic activity. ISI plots (see Fig. 1) were generated for each group of spikes. In the example shown, all spikes were produced by a single mouse dLGN cell during a contrast tuning experiment (see METHODS). SF and TF were optimal. Please note that, although the large number of high-frequency spikes fired by the example cell means that certain spikes appear to be shared between burst and HFT firing, these spike groups actually represent mutually exclusive response subsets. The line marked TF displays the temporal period of the stimulus (here 6 Hz). All other conventions as in Fig. 1. B: average values for burst (B), TX, HFT, and gap (G) spikes are shown for the 9 response measures in which significant burst-tonic differences were found in our original analyses. In cases where all 4 samples were normally distributed, bar plots show the means ± SE (2, 4, 5, 6); in all other cases, plots show medians (1, 3, 7, 8, 9). The P values at top left show the result of paired comparisons across the 4 spike types (repeated-measures ANOVA or Friedman tests, depending on normality). When P < 0.05, there was significant variation among groups. The results of individual post hoc comparisons (Tukey or Dunn, depending on normality) are shown above the appropriate bars in cases where the overall repeated-measures test produced significant results. Only the results of B-TX, TX-HFT, and TX-G comparisons are shown because these are the most informative and are discussed in the text. ns, nonsignificant; *, P < 0.05; ** P < 0.01; *** P < 0.001. It is clear that relationships between burst, TX, HFT, and gap spikes are complex and do not follow similar patterns across all response parameters. The relationships between burst and TX spikes are not always mirrored in either the relationships between TX and HFT spikes, or the relationships between TX and gap spikes. Neither the high-frequency firing component of burst spikes nor the requisite preceding silence of bursts can therefore be wholly responsible for the burst-tonic functional differences described in the preceding text.

 
The results of the 4-way comparisons are shown in Fig. 12B. For each visual response measure, we employed a repeated-measures ANOVA or Friedman test (depending on normality) followed by appropriate posttests (Tukey or Dunn, depending on normality) to assess the significance of differences between the four spike types. Our main observation is that the burst-tonic differences reported in the main body of paper are not as a whole attributable to either the high-frequency firing or the preceding silence features of burst firing because the relationships between HFT, gap, and TX spikes are not consistent. There were nine visual response measures reported in the preceding text in which burst and tonic firing differed significantly. In three of these (rectification, bandwidth, 1st-spike timing reliability), both HFT and gap responses were significantly different from TX spikes in the same direction as the burst-TX difference (Fig. 12B, 3, 6, and 9). This might suggest that the original burst-tonic differences in these parameters were due to both high-frequency firing and preburst silence. However, there were two cases (latency, phase) in which gap spikes and TX spikes differed significantly and in the same direction as the burst-TX difference, but no significant difference was seen between HFT and TX firing (Fig. 12B, 1 and 2). Trends in the contrast n parameter were also in this direction, although sample sizes were small and differences were insignificant (Fig. 12B7). In these parameters, the original burst-tonic differences appear to be mainly due to burst firing's feature of preceding silence. This is immediately contradicted, though, by the case of spike timing reliability, in which HFT and TX spikes differ in same way as burst and TX spikes, but where the gap-TX difference is significant in the opposite direction (Fig. 12B8). For this parameter, high-frequency firing seems to be the main contributor to the burst-tonic difference reported above. Finally, there were two parameters in which the four-way comparisons revealed a rather puzzling result. TF peak and high50 measures did not differ between burst and TX spikes; instead, the original burst-tonic difference appeared to be made entirely by gap spikes' contribution to tonic sample averages, since gap responses differed significantly from all other firing types (Fig. 12B, 4 and 5).

These results show that it was not possible to attribute burst-tonic functional differences, as a whole, to a single feature of burst firing. Perhaps assessing the firing properties of different subgroups of tonic spikes is not the best approach to this problem, but it was a sensible first step. Maybe intracellular recordings, by allowing a more accurate classification of burst, gap, and HFT spikes based on membrane hyperpolarizations and IT activity, might shed more light on the issue. Until such methods are employed, however, it is perhaps most straightforward to conclude that the collection of response properties that differentiate burst spikes from their tonic counterparts is produced by a complex combination of both high-frequency firing and a requisite preceding silent period.


 DISCUSSION
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
We present a comprehensive comparison of the visual response properties of burst and tonic firing in mouse dLGN cells. Once burst-like events in the firing of mouse dLGN neurons had been identified on the basis of temporal patterns of spike firing, various features of these events were shown to match those previously described for thalamic bursting in other species and preparations. Burst spikes were shown to occur more often during visual stimulation than during periods of spontaneous activity. Burst and tonic spikes within individual mouse dLGN neurons were then compared with respect to their response timing, linearity, SF and TF tuning, contrast response characteristics, signal detection properties, and spike time reliability. This analysis uncovered a number of functional distinctions between the two firing modes.

Comparison with previous studies

INCIDENCE OF BURSTING.  Criteria previously used to identify burst firing from extracellular recordings in the cat dLGN (e.g., Guido et al. 1992Go; Lu et al. 1992Go) were used successfully here in mouse dLGN neurons to classify firing events that already showed distinct temporal patterns (Fig. 1) and that possessed properties similar to those reported for burst firing in other species and preparations (Fig. 2). Coupled with in vitro thalamic recordings showing that IT and burst spikes are very similar across different mammalian species (guinea pig: Jahnsen and Llinas 1984Go; McCormick and Feeser 1990Go; cat: Zhan et al. 1999Go; monkey: Ramcharan et al. 2000aGo; mouse: Kim et al. 2001Go), these data suggest that burst firing is well conserved across the mammalian order. It will require intracellular recordings from mouse dLGN cells in vivo, however, for us to know for sure that the action potentials classified as burst firing in this study really were riding on a low-threshold calcium spike.

In terms of burst prevalence, ~20% of spikes form part of bursts in the anesthetized mouse dLGN during presentation of sinusoidal grating stimuli (Fig. 4). This is comparable to the value of 20–25% reported for X cells during visual stimulation in the anesthetized cat ( Lesica and Stanley 2004Go; Reinagel et al. 1999Go), but much lower than the Y cell value under the same conditions (60%) ( Reinagel et al. 1999Go), and much higher than the 1–2% reported during visual stimulation in the awake monkey dLGN ( Ramcharan et al. 2000bGo). Bursting was likely facilitated by the anesthetized preparation employed here (cf. Guido and Weyand 1995Go; Reinagel et al. 1999Go). If this is the case, experiments in anesthetized animals may not produce results that can translate directly to behavioral situations but do present the best opportunity to compare the response properties of many burst and many tonic spikes fired by the same cells to the same sensory stimuli.

RELATIVE SPIKE TIMING.  Burst spikes in the mouse dLGN have longer latencies, but because they are triggered by earlier stimulus phases tend to occur before tonic spikes in responses to sinusoidal gratings (Fig. 5). The situation in the cat dLGN is almost identical. Cells that fire more burst spikes have longer latencies and larger phase advances ( Mukherjee and Kaplan 1995Go), while stimulating the histaminergic tuberomammillary nucleus increases tonic firing and produces a phase lag in responses to grating stimuli ( Uhlrich et al. 2002Go). Furthermore, a Fourier analysis of in vitro responses to sinusoidal inputs shows that the line relating TF to F1 response phase is steeper (indicative of longer latency), with a higher y intercept (indicative of a phase advance), for burst versus tonic spikes ( Smith et al. 2000Go). And although Guido et al. (1992)Go report that burst and tonic spikes in the cat dLGN have similar latencies, with burst spikes phase advanced, close inspection of their raw data reveals that their plots of TF versus F1 response phase are strikingly similar to those shown in Fig. 5 here: burst spikes are phase advanced at low but not high TFs. Had line fitting procedures not differed between the two studies, with Guido et al. (1992)Go employing response-weighted fits that decrease the influence of weak responses to high TFs, both might have reported longer latencies for burst spikes.

The phase advance and the longer latency of burst spikes could be consequences of the LTS that produces burst firing ( Jahnsen and Llinas 1984Go). By definition, the LTS is low threshold, and there is in vitro evidence that burst spikes can be elicited by smaller stimulation amplitudes than tonic spikes (e.g., Zhan et al. 1999Go). In vivo, smaller stimulation strengths, in the form of low rates of firing in retinal afferents, should occur at lower stimulus contrasts and therefore at earlier phases of the stimulating half cycle (dark or bright) of a sinusoidal grating. Burst spikes would thus be the earliest spikes triggered by a given cycle of such a stimulus. Long latencies, on the other hand, could be a consequence of the LTS's slow dynamics. In vitro, burst spikes require longer durations of stimulation than tonic spikes ( McCormick and Feeser 1990Go), and in vivo, retinal spikes preceding burst firing occur much longer before geniculate activity than retinal spikes preceding tonic firing ( Rowe and Fischer 2001Go).

Linearity

Burst and tonic spikes in the mouse dLGN do not differ with respect to linearity of spatial summation, but burst firing is more rectified than tonic activity (Fig. 6). This is the first study to use a null test ( Enroth-Cugell and Robson 1966Go; Hochstein and Shapley 1976Go) to assess linearity of spatial summation in the two thalamic firing modes. Because both burst and tonic spikes display linear summation in linear mouse dLGN cells, it would be extremely interesting to know whether the ON-OFF doubling nonlinearity of cat Y cells (e.g., So and Shapley 1979Go) is carried by burst spikes, tonic spikes, or both.

In contrast, rectification nonlinearities have been well studied previously. As in the mouse dLGN, burst spikes in the cat dLGN in vivo ( Guido et al. 1992Go; Lu et al. 1992Go; Mukherjee and Kaplan 1995Go) and in vitro ( Smith et al. 2000Go) follow the temporal form of a sinusoidal stimulus less closely than tonic spikes.

Spatial and temporal tuning

Burst and tonic firing in the mouse dLGN do not differ in terms of SF tuning (Fig. 7). No other studies have directly compared SF tuning in the two firing modes, although Guido et al. (1992)Go did compare the percentage of burst spikes fired to different stimuli across cells and found no variation with SF. In addition, stimulation of the tuberomammillary nucleus, which increases tonic firing in the cat dLGN, produced very little change in the SF tuning properties of geniculate neurons ( Uhlrich et al. 2002Go). Those small but significant changes that were produced—increased RF center strength (kc) and decreased RF center radius (rc)—could well have been produced by nonburst related effects of the stimulation but were nevertheless mirrored by nonsignificant trends in the mouse dLGN data. In reverse correlation studies of the spatial RF structure of burst and tonic spikes, initial results are inconsistent: whereas Rivadulla et al. (2003)Go find RF center size to be smaller for burst responses, Allitto et al. (2003)Go do not but do observe increased surround strength during burst firing. There are therefore no large or consistent differences in the spatial response properties of burst and tonic spikes.

Burst spikes in the mouse dLGN prefer lower TFs and show sharper TF tuning than tonic spikes (Fig. 8). Because burst firing can only occur after ≥100 ms of hyperpolarization, one might expect it to have difficulty in following higher TFs. Indeed, in vitro studies in the guinea pig ( McCormick and Feeser 1990Go) and cat ( Smith et al. 2000Go) also report that burst spikes stop responding at lower cutoff TFs than tonic spikes. The present finding of sharper TF tuning (i.e., lower bandwidth) in burst firing is also well supported by previous data. Tonic TF tuning in slices taken from the cat dLGN is more broadband than its burst counterpart ( Smith et al. 2000Go), while a higher incidence of bursting was correlated with more band-pass TF tuning in the intact cat dLGN ( Mukherjee and Kaplan 1995Go). There is some slight disagreement between the present data and the latter in vivo cat study, however. While burst spikes here displayed sharper TF tuning mainly by dint of reduced responses at high TFs, cat dLGN cells showed sharpened tuning associated with bursting because of reduced responses at low TFs ( Mukherjee and Kaplan 1995Go). This may reflect species or methodological differences between the two studies but could also be caused by random fluctuations between different samples. There was certainly a trend toward higher low50 values for burst firing in the mouse dLGN.

Contrast response characteristics

Comparison of contrast response characteristics showed that burst and tonic firing modes do not differ in terms of contrast gain or c50. However, the contrast-response curves of burst firing are more step-like (reflected in higher values of the parameter n) than those of tonic spikes (Fig. 9). This latter difference might have been predicted on the basis of in vitro investigations. Due to the all-or-nothing nature of the LTS, the amplitude of which depends not on the strength of depolarizing inputs but rather on the relative hyperpolarization of the membrane potential when those inputs arrive, the burst response of a dLGN cell to increasing current injection at a constant membrane potential resembles a step function. Tonic firing, in contrast, increases almost linearly with increasing stimulation strength ( McCormick and Feeser 1990Go; Ramcharan et al. 2000aGo; Zhan et al. 1999Go). Burst contrast-response functions observed in vivo may not be entirely step-like, and in vivo tonic functions are certainly not always completely linear (see Fig. 9A). Nevertheless, the relatively more step-like relationship between burst responses and stimulus contrast might be expected based on the properties of the LTS and represents a significant functional difference between burst and tonic firing.

Stimulus detection

The preceding discussion shows that differences in burst and tonic visual coding in the mouse dLGN usually reflect similar differences reported in the cat dLGN. The functions of the two firing modes therefore seem rather well conserved across these two mammalian species. This conclusion is bolstered by our analysis of stimulus detection: in the mouse, as in the cat ( Guido et al. 1995Go), better signal detection is associated with more prominent burst firing in a dLGN cell's responses (Fig. 10). It appears, however, that bursting may not be associated with better stimulus detection in all thalamic nuclei (or maybe in all species), given that a greater incidence of burst firing is not related to better stimulus detection in the guinea pig auditory thalamus ( Massaux et al. 2004Go).

Reliability of spike timing

In the mouse dLGN, the relative timing of all burst spikes is more reliable than that of all tonic spikes, but tonic firing is more reliable in terms of the first spikes fired to a stimulus (Fig. 11). A single burst spike therefore carries more information concerning the timing of a stimulus event than a single tonic spike. Combined with the fact that burst spikes occur earlier than tonic spikes in responses to a given stimulus, this may make bursting ideal for signaling to cortex when a stimulus appeared, even if the first spikes in bursts are not quite so reliable in their timing as the first spikes in tonic responses. This latter result contradicts a finding in the cat dLGN, where the reliability of first-spike timing was better for burst than for tonic responses ( Guido and Sherman 1998Go). This contradiction may reflect true species differences, possible differences in halothane anesthetic depth (see preceding text), or differences in the stimuli employed in the two studies. Unlike the flashed spots employed by Guido and Sherman (1998)Go, the sinusoidal stimuli used here do not possess a single onset time and therefore might not be entirely appropriate for assessing spike timing properties. Still, the increased variability in the timing of burst first spikes in the mouse dLGN does appear to fit well with a different cat study, in which paired retinal and dLGN recordings showed burst soikes to be far more loosely and variably coupled to retinal inputs than tonic spikes (Rowe and Fischer 2001).

Although the main focus of this paper was simply to describe the different types of visual information that burst and tonic spikes in the mouse dLGN might send to visual cortex, it is worth briefly considering some implications of the mechanism by which burst and tonic stimulus selectivity arises in the first place. Tonic firing in thalamocortical relay cells is largely determined by excitatory retinal inputs with the correspondence rather tight between presynaptic spikes in retinal ganglion cells and postsynaptic spikes in geniculate neurons ( Rowe and Fischer 2001Go). The properties of the IT current underlying burst responses, however, mean that although burst spikes are normally triggered by excitatory retinal inputs, their stimulus selectivity should be largely determined by local inhibitory circuitry. Not only is membrane hyperpolarization needed to de-inactivate IT and allow burst firing to occur (e.g., Jahnsen and Llinas 1984Go), the degree of membrane hyperpolarization also determines the amplitude of burst responses: the more hyperpolarized the membrane, the more burst spikes result from a constant excitatory trigger ( Zhan et al. 1999Go). Measures of F1 response amplitude, such as those widely employed to assess stimulus selectivity in this paper and elsewhere, depend on both the probability of a response occurring to a particular stimulus cycle and the size of that response when it does occur. In burst firing, both of these factors will be largely determined by the degree of hyperpolarization present in dLGN cells before they receive an excitatory triggering input from the retina. It follows that burst spike response selectivity during the presentation of drifting grating stimuli is determined in large part by the stimulus half-cycle that inhibits the RF center rather than the stimulus half-cycle that excites the RF center and actually causes a spike response to occur. Because the excitatory "push" and inhibitory "pull" mechanisms operating on dLGN neurons are extremely well matched, at least spatially ( Hirsch 2003Go), this might not render burst and tonic stimulus selectivity too different. Indeed, we show that differences in burst and tonic visual response properties in the mouse dLGN, although significant, are rather small. However, this similarity might arise because of the cyclic nature of the stimuli used. Perhaps nonperiodic stimuli, such as natural scenes (cf. Lesica and Stanley 2004Go), that do not provide balanced and rhythmic durations of inhibitory pull and excitatory push will produce large differences in the visual response properties of burst and tonic firing and will allow us to better understand the functions of the two thalamic response modes.

What do bursts do?

The analyses presented here extend our knowledge of the visual encoding properties of burst and tonic firing in the dLGN and will be important when future studies investigate the molecular basis of thalamic processing using mutant mice (see INTRODUCTION), but do they shed any light on the fundamental functions of the two firing modes? They certainly confirm that tonic spikes are used to extract detailed information about visual stimulus attributes. Where tuning properties differ between burst and tonic firing, tonic spikes, with less rectification, broader TF tuning, better responsivity to high TFs, and more linear responses to increasing contrast, are better placed to signal information about a wider range of visual stimulus attributes. But it is not at all clear from the present analysis what the role of burst firing might be.

The prevailing hypothesis of burst function revolves around previous reports of greater rectification and signal detection in burst firing (e.g., Guido et al. 1992Go; Sherman 2001Go; Sherman and Guillery 1996Go). Burst firing is seen as a "wake-up call" for cortex, detecting stimuli but then leaving any detailed analysis to tonic firing. Indeed, the increase in stimulus detection that occurs with increased burst prominence in mouse dLGN cells (Fig. 10), along with the phase advance of burst firing (Fig. 5) ( Guido et al. 1992Go; Mukherjee and Kaplan 1995Go) and the reliable timing of burst spiking as a whole (Fig. 11) suggest that burst spikes may be ideal for signaling to cortex ( Swadlow and Gusev 2001Go; Swadlow et al. 2002Go) that a stimulus has occurred and that it started a certain time ago. Furthermore, given that the visual encoding differences between tonic and burst spikes are not huge (Figs. 58), this bursting detection signal could also contain a great deal of information concerning stimulus attributes, something for cortex to start to process before more detailed tonic spike-carried information arrives. In awake cats, burst spikes are fired at the very start of stimulus presentations and visual fixation periods ( Guido and Weyand 1995Go; Weyand et al. 2001Go), suggesting that bursting may indeed signal the presence of new stimuli to cortex. In addition, burst spikes fired by X cells in the anesthetized cat dLGN tend to occur when excitatory features follow the prolonged presence of inhibitory features in natural movie stimuli ( Lesica and Stanley 2004Go). In natural scenes, therefore bursts might encode the appearance of new objects. However, burst firing is also more prevalent during "passive" viewing in awake cats ( Weyand et al. 2001Go) and, when occurring before stimulus onset, increases the probability of a subsequent stimulus-induced response in rat somatosensory thalamus ( Fanselow et al. 2001Go). This suggests that burst spikes may not be detecting stimuli per se but are instead priming the brain to do so. It would be extremely useful to know whether bursting during stimulus presentation correlates with the animal's detection of that stimulus.

Another view of firing mode function, based on the different temporal tuning properties of burst and tonic spikes ( Mukherjee and Kaplan 1995Go) (Fig. 8), argues that the thalamus is a "tuneable filter" for visual information. Various modulatory inputs to the dLGN, by altering the probability of burst firing in relay cells (e.g., Godwin et al. 1996Go; Uhlrich et al. 2002Go), could alter the tuning of information sent to visual cortex. Given the present data this could certainly work in the contrast (Fig. 9) and TF (Fig. 8) domains, although it is unlikely to apply to spatial information (Fig. 7). However, it is unclear to what extent bursting tuning properties are functional adaptations and to what extent they are simply artifacts of burst properties that might have evolved for different reasons. It certainly seems strange that so much circuitry would be invested into changing temporal tuning in the dLGN when such a large amount of temporal information is discarded at the geniculocortical synapse (e.g., Hawken et al. 1996Go). And while increasing bursting might tighten temporal tuning in the dLGN, this would necessarily occur at the expense of detailed information about stimulus contrast (Fig. 9) and stimulus temporal form (increased rectification; Fig. 6). Still, until a definitive test of this hypothesis is formulated, it remains a valid suggestion.

Finally, there is a view that burst firing might be vital not for signaling specific stimulus attributes but for synchronizing activity across thalamic neurons. Sillito et al. (1994)Go showed that synchronous activity in dLGN cells is dependent on cortical feedback. In a model of these data, Kirkland and Gerstein (1998)Go show that this cortically induced synchronization is critically dependent on burst firing. The reliability of burst timing (Fig. 11) might be useful in this regard, although direct tests of the hypothesis, and a better understanding of whether bursts are linked more strongly to changes in sensory stimuli or changes in internal state, await multi-unit thalamic recordings.


 ACKNOWLEDGMENTS
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
The authors thank D. Smyth for computational assistance and C. Akerman for helpful discussions. M. S. Grubb was a student on the Wellcome Trust 4-year PhD program in Neuroscience.


 FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: I. D. Thompson University Laboratory of Physiology, Parks Road, Oxford, OX1 3PT UK (E-mail: ian.thompson{at}physiol.ox.ac.uk)


 REFERENCES
 
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGMENTS
 REFERENCES
 
Albrecht DG and Hamilton DB. Striate cortex of monkey and cat: contrast response function. J Neurophysiol 48: 217–237, 1982.[Free Full Text]

Allito HJ, Weyand TG, and Usrey WM. Corticothalamic interactions and the role of thalamic bursts in visual processing. Soc Neurosci Abstr 29: 485.16, 2003.

Dräger UC. Receptive fields of single cells and topography in mouse visual cortex. J Comp Neurol 160: 269–290, 1975.[CrossRef][Web of Science][Medline]

Enroth-Cugell C and Robson JG. The contrast sensitivity of retinal ganglion cells of the cat. J Physiol 187: 517–552, 1966.

Fanselow EE, Sameshima K, Baccala LA, and Nicolelis MA. Thalamic bursting in rats during different awake behavioral states. Proc Natl Acad Sci USA 98: 15330–15335, 2001.[Abstract/Free Full Text]

Godwin DW, Vaughan JW, and Sherman SM. Metabotropic glutamate receptors switch visual response mode of lateral geniculate nucleus cells from burst to tonic. J Neurophysiol 76: 1800–1816, 1996.[Abstract/Free Full Text]

Gordon JA and Stryker MP. Experience-dependent plasticity of binocular responses in the primary visual cortex of the mouse. J Neurosci 16: 3274–3286, 1996.[Abstract/Free Full Text]

Green DM and Swets JA. Signal Detection Theory and Psychophysics. New York: Wiley, 1966.

Grubb MS, Rossi FM, Changeux JP, and Thompson ID. Abnormal functional organization in the dorsal lateral geniculate nucleus of mice lacking the {beta}2 subunit of the nicotinic acetylcholine receptor. Neuron 40: 1161–1172, 2003.[CrossRef][Web of Science][Medline]

Grubb MS and Thompson ID. Quantitative characterization of visual response properties in the mouse dorsal lateral geniculate nucleus. J Neurophysiol 90: 3594–3607, 2003.[Abstract/Free Full Text]

Guido W, Lu SM, and Sherman SM. Relative contributions of burst and tonic responses to the receptive field properties of lateral geniculate neurons in the cat. J Neurophysiol 68: 2199–2211, 1992.[Abstract/Free Full Text]

Guido W, Lu SM, Vaughan JW, Godwin DW, and Sherman SM. Receiver operating characteristic (ROC) analysis of neurons in the cat's lateral geniculate nucleus during tonic and burst response mode. Vis Neurosci 12: 723–741, 1995.[Web of Science][Medline]

Guido W and Sherman SM. Response latencies of cells in the cat's lateral geniculate nucleus are less variable during burst than tonic firing. Vis Neurosci 15: 231–237, 1998.[CrossRef][Web of Science][Medline]

Guido W and Weyand T. Burst responses in thalamic relay cells of the awake behaving cat. J Neurophysiol 74: 1782–1786, 1995.[Abstract/Free Full Text]

Hawken MJ, Shapley RM, and Grosof DH. Temporal-frequency selectivity in monkey visual cortex. Vis Neurosci 13: 477–492, 1996.[Web of Science][Medline]

Hirsch JA. Synaptic physiology and receptive field structure in the early visual pathway of the cat. Cereb Cortex 13: 63–69, 2003.[Abstract/Free Full Text]

Hochstein S and Shapley RM. Quantitative analysis of retinal ganglion cell classifications. J Physiol 262: 237–264, 1976.[Abstract/Free Full Text]

Hubel DH and Wiesel TN. Integrative action in the cat's lateral geniculate body. J Physiol 155: 385–398, 1961.[Free Full Text]

Jahnsen H and Llinas R. Electrophysiological properties of guinea-pig thalamic neurons: an in vitro study. J Physiol 349: 205–226, 1984.[Abstract/Free Full Text]

Kim D, Song I, Keum S, Lee T, Jeong MJ, Kim SS, McEnery MW, and Shin HS. Lack of the burst firing of thalamocortical relay neurons and resistance to absence seizures in mice lacking alpha(1G) T-type Ca(2+) channels. Neuron 31: 35–45, 2001.[CrossRef][Web of Science][Medline]

Kirkland KL and Gerstein GL. A model of cortically induced synchronization in the lateral geniculate nucleus of the cat: a role for low-threshold calcium channels. Vision Res 38: 2007–2022, 1998.[CrossRef][Web of Science][Medline]

Lennie P and Perry VH. Spatial contrast sensitivity of cells in the lateral geniculate nucleus of the rat. J Physiol: 31: 569–579, 1981.

Lesica NA and Stanley GB. Encoding of natural scene movies by tonic and burst spikes in the lateral geniculate nucleus. J Neurosci 24: 10731–10740, 2004.[Abstract/Free Full Text]

Livingstone MS and Hubel DH. Effects of sleep and arousal on the processing of visual information in the cat. Nature 291: 554–561, 1981.[CrossRef][Medline]

Lu SM, Guido W, and Sherman SM. Effects of membrane voltage on receptive field properties of lateral geniculate neurons in the cat: contributions of the low-threshold Ca2+ conductance. J Neurophysiol 68: 2185–2198, 1992.[Abstract/Free Full Text]

Lu SM, Guido W, and Sherman SM. The brain-stem parabrachial region controls mode of response to visual stimulation of neurons in the cat's lateral geniculate nucleus. Vis Neurosci 10: 631–642, 1993.[Web of Science][Medline]

Macmillan NA and Creelman CD. Detection Theory: A User's Guide. New York: Cambridge, 1991.

Massaux A, Dutrieux G, Cotillon-Williams N, Manunta Y, and Edeline JM. Auditory thalamus bursts in anesthetized and non-anesthetized states: contribution to functional properties. J Neurophysiol 91: 2117–2134, 2004.[Abstract/Free Full Text]

McCarley RW, Benoit O, and Barrionuevo G. Lateral geniculate nucleus unitary discharge in sleep and waking: state- and rate-specific aspects. J Neurophysiol 50: 798–818, 1983.[Abstract/Free Full Text]

McCormick DA and Bal T. Sleep and arousal: thalamocortical mechanisms. Annu Rev Neurosci 20: 185–215, 1997.[CrossRef][Web of Science][Medline]

McCormick DA and Feeser HR. Functional implications of burst firing and single spike activity in lateral geniculate relay neurons. Neuroscience 39: 103–113, 1990.[CrossRef][Web of Science][Medline]

Mukherjee P and Kaplan E. Dynamics of neurons in the cat lateral geniculate nucleus: in vivo electrophysiology and computational modeling. J Neurophysiol 74: 1222–1243, 1995.[Abstract/Free Full Text]

Paxinos G and Franklin K. The Mouse Brain in Stereotaxic Coordinates: New York: Academic, 2001.

Radhakrishnan V, Tsoukatos J, Davis KD, Tasker RR, Lozano AM, and Dostrovsky JO. A comparison of the burst activity of lateral thalamic neurons in chronic pain and non-pain patients. Pain 80: 567–575, 1999.[CrossRef][Web of Science][Medline]

Ramcharan EJ, Cox CL, Zhan XJ, Sherman SM, and Gnadt JW. Cellular mechanisms underlying activity patterns in the monkey thalamus during visual behavior. J Neurophysiol 84: 1982–1987, 2000a.[Abstract/Free Full Text]

Ramcharan EJ, Gnadt JW, and Sherman SM. Burst and tonic firing in thalamic cells of unanesthetized, behaving monkeys. Vis Neurosci 17: 55–62, 2000b.[CrossRef][Web of Science][Medline]

Reinagel P, Godwin D, Sherman SM, and Koch C. Encoding of visual information by LGN bursts. J Neurophysiol 81: 2558–2569, 1999.[Abstract/Free Full Text]

Rivadulla C, Martinez L, Grieve KL, and Cudeiro J. Receptive field structure of burst and tonic firing in feline lateral geniculate nucleus. J Physiol 553: 601–610, 2003.[Abstract/Free Full Text]

Rodieck RW. Quantitative analysis of cat retinal ganglion cell response to visual stimuli. Vision Res 5: 583–601, 1965.[CrossRef][Medline]

Rowe MH and Fischer Q. Dynamic properties of retino-geniculate synapses in the cat. Vis Neurosci 18: 219–231, 2001.[CrossRef][Web of Science][Medline]

Schwarte LA, Zuurbier CJ, and Ince C. Mechanical ventilation of mice. Basic Res Cardiol 95: 510–520, 2000.[CrossRef][Web of Science][Medline]

Shapley R and Lennie P. Spatial frequency analysis in the visual system. Annu Rev Neurosci 8: 547–583, 1985.[CrossRef][Web of Science][Medline]

Sherman SM. Tonic and burst firing: dual modes of thalamocortical relay. Trends Neurosci 24: 122–126, 2001.[CrossRef][Web of Science][Medline]

Sherman SM and Guillery RW. Functional organization of thalamocortical relays. J Neurophysiol 76: 1367–1395, 1996.[Abstract/Free Full Text]

Sherman SM and Guillery RW. The role of the thalamus in the flow of information to the cortex. Philos Trans R Soc Lond B Biol Sci 357: 1695–1708, 2002.[Abstract/Free Full Text]

Sillito AM, Jones HE, Gerstein GL, and West DC. Feature-linked synchronisation of thalamic relay cell firing induced by feedback from the visual cortex. Nature 369: 479–482, 1994.[CrossRef][Medline]

Smith GD, Cox CL, Sherman SM, and Rinzel J. Fourier analysis of sinusoidally driven thalamocortical relay neurons and a minimal integrate-and-fire-or-burst model. J Neurophysiol 83: 588–610, 2000.[Abstract/Free Full Text]

So YT and Shapley R. Spatial properties of X and Y cells in the lateral geniculate nucleus of the cat and conduction velocities of their inputs. Exp Brain Res 36: 533–550, 1979.[Web of Science][Medline]

So YT and Shapley R. Spatial tuning of cells in and around lateral geniculate nucleus of the cat: X and Y relay cells and perigeniculate interneurons. J Neurophysiol 45: 107–120, 1981.[Free Full Text]

Swadlow HA and Gusev AG. The impact of "bursting" thalamic impulses at a neocortical synapse. Nat Neurosci 4: 402–408, 2001.[CrossRef][Web of Science][Medline]

Swadlow HA, Gusev AG, and Bezdudnaya T. Activation of a cortical column by a thalamocortical impulse. J Neurosci 22: 7766–7773, 2002.[Abstract/Free Full Text]

Uhlrich DJ, Manning KA, and Xue JT. Effects of activation of the histaminergic tuberomammillary nucleus on visual responses of neurons in the dorsal lateral geniculate nucleus. J Neurosci 22: 1098–1107, 2002.[Abstract/Free Full Text]

Uhlrich DJ, Tamamaki N, Murphy PC, and Sherman SM. Effects of brain stem parabrachial activation on receptive field properties of cells in the cat's lateral geniculate nucleus. J Neurophysiol 73: 2428–2447, 1995.[Abstract/Free Full Text]

Usrey WM, Reppas JB, and Reid RC. Specificity and strength of retinogeniculate connections. J Neurophysiol 82: 3527–3540, 1999.[Abstract/Free Full Text]

Wagor E, Mangini NJ, and Pearlman AL. Retinotopic organization of striate and extrastriate visual cortex in the mouse. J Comp Neurol 193: 187–202, 1980.[CrossRef][Web of Science][Medline]

Weyand TG, Boudreaux M, and Guido W. Burst and tonic response modes in thalamic neurons during sleep and wakefulness. J Neurophysiol 85: 1107–1118, 2001.[Abstract/Free Full Text]

Yu Y and Bradley A. Engineering chromosomal rearrangements in mice. Nature Rev Genet 2: 780–790, 2001.[CrossRef][Web of Science][Medline]

Zhan XJ, Cox CL, Rinzel J, and Sherman SM. Current clamp and modeling studies of low-threshold calcium spikes in cells of the cat's lateral geniculate nucleus. J Neurophysiol 81: 2360–2373, 1999.[Abstract/Free Full Text]




This article has been cited by other articles:


Home page
J. Neurophysiol.Home page
R. R. Llinas and M. Steriade
Bursting of Thalamic Neurons and States of Vigilance
J Neurophysiol, June 1, 2006; 95(6): 3297 - 3308.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
93/6/3224    most recent
00445.2004v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Grubb, M. S.
Right arrow Articles by Thompson, I. D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Grubb, M. S.
Right arrow Articles by Thompson, I. D.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2005 by the The American Physiological Society.