|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Group in Vision Science, School of Optometry, and Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, California
Submitted 8 August 2007; accepted in final form 20 October 2007
|
|
ABSTRACT |
|---|
|
|
|
INTRODUCTION |
|---|
|
In this suppressive field model, as stimulus contrast increases, the local contrast within the suppressive field also increases, which decreases the response gain. This decrease in gain causes a saturation of response with increasing stimulus contrast. An assumption of this model is that the static nonlinear component is simply a linear rectifying function. This is of primary significance because it implies an approximately linear response at low contrast levels where the suppressive field is minimally activated. If in fact there is a clear nonlinearity at low contrasts, this addition can improve the predictive power of the suppressive field model, and it also could have implications for the interpretation of visual function at low contrasts.
We have examined this issue by obtaining contrast–response functions from extracellular recordings of neurons in the cat's LGN. We fit the data with a Naka–Rushton function (Naka and Rushton 1966
). Results show significant expansive nonlinearities at low contrasts. Additionally, we have measured directly the static nonlinear function by comparing a linear prediction with actual spike responses. The linear predictions are generated from spatiotemporal receptive fields obtained by m-sequence stimulation. Static nonlinear functions for the majority of neurons of our sample exhibit power-law nonlinearities with a mean exponent of 1.58. These results show a clear expansive nonlinear component in LGN neurons for low-contrast visual stimuli. It is likely that this has important consequences for basic response properties of cortical neurons.
|
|
METHODS |
|---|
|
All procedures complied with the National Institutes of Health Guide for the Care and Use of Laboratory Animals. Extracellular recordings were made using epoxy-coated tungsten microelectrodes in the LGN of anesthetized and paralyzed mature cats. Cats were initially anesthetized with isoflurane (1–4%). After catheterization, a continuous infusion was given of a combination of fentanyl citrate (10 µg·kg–1·h–1) and thiopental sodium (6 mg·kg–1·h–1). Bolus injections of thiopental sodium were given as required during surgery. After a tracheal cannula was positioned, isoflurane was discontinued and the animal was artificially ventilated with a mixture of 25% O2 and 75% N2O. Respiration rate was manually adjusted to maintain an end-tidal CO2 of 34–38 mmHg. Body temperature was maintained at 38°C with a closed-loop controlled heating pad. A craniotomy was performed over the LGN and the dura was resected and covered with agar and wax to form a closed chamber. After completion of all surgical procedures, continuous injection of fentanyl citrate was discontinued, and thiopental sodium concentration was lowered gradually to a level at which the cat was stabilized for
1 h. The level of anesthetic used was determined individually for each cat. The range used was 1.0–3.0 mg·kg–1·h–1 and a typical level was 1.5 mg·kg–1·h–1. Once a stabilized anesthetic level was reached, it was kept constant throughout the experiment. To minimize eye movements during visual tests, animals were immobilized with pancuronium bromide (0.2 mg·kg–1·h–1). EEG, ECG, heart rate, temperature, end-tidal CO2, and intratracheal pressure were monitored for the entire duration of the experiment. Contact lenses were used that were opaque except for a central 4-mm-diameter window to create an artificial pupil. To focus the eyes on the stimulation screen, ophthalmoscopic refraction was used to determine appropriate lens power.
For cortical data, electrode penetrations were made at Horsley–Clarke coordinates P4L2 angled 20° anterior and 10° medial. Electrodes were advanced until visually responsive cells were found. LGN electrode penetrations were made perpendicular to the cortical surface at Horsley–Clarke coordinates of approximately A6L9. Electrodes were then advanced until visually responsive cells with LGN response characteristics were found (typically around 12 mm below the cortical surface). Recordings were made from all layers of the LGN.
Extracellular recording
Single units were isolated in real time by the shape of their spike waveforms using custom software. An initial estimate of the tuning parameters was made qualitatively by computer-controlled manipulation of drifting sinusoidal gratings. The spatial extent of visual stimulation was kept larger than the receptive field size. Temporal frequency tuning curves were measured with drifting sinusoidal gratings at 50% contrast. Spatial frequency and contrast tuning curves were measured at optimal temporal frequencies determined for each cell, typically between 4 and 15 cycles/s.
Visual stimulation
Visual patterns consisting of sinusoidal gratings or noise patterns were presented on a large CRT at a frame rate of 75 Hz. The 47.8-cm-diameter CRT was positioned at an optical distance of 41.8 cm in front of the cat's eyes, and split so that each half of the display stimulated the left or right eye. Luminance from the CRT was calibrated for a linear range with maximum and minimum values of 90 and 0.1 cd/m2, respectively.
Data analysis
For contrast tuning data, the first harmonic (F1 component) was used for analysis. All data fitting was done by minimizing the sum squared error using fminsearch in Matlab (The MathWorks, Natick, MA), which implements the Nelder–Mead nonlinear minimization algorithm. For m-sequence analysis, spikes were first binned with a window of one m-sequence frame (typically
13 ms). Reverse correlation was done using the fast m-transform method (Reid et al. 1997
; Sutter 1991
), and the first-order kernel was extracted from m-transformed data. Estimation of the static nonlinearity was done by comparing the linear prediction to the actual spike response. Linear predictions were generated by convolving the linear spatiotemporal receptive fields with the m-sequence stimulus. Details of this analysis are given elsewhere (Anzai et al. 1999a
; Chichilnisky 2001
).
|
|
RESULTS |
|---|
|
Two visual stimulation protocols were used in this study. In the first protocol, responses to drifting gratings were measured with approximately optimal spatial and temporal frequencies at different contrast levels, as illustrated in Fig. 1A. Contrast in this experiment is defined in the standard way as
![]() | (1) |
|
![]() | (2) |
![]() | (3) |
![]() | (4) |
Recent studies and models assume that contrast–response data for LGN neurons follow the Naka–Rushton function with n = 1 (Bonin et al. 2005
; Li et al. 2006
; Priebe and Ferster 2006
). This assumption means that the contrast–response function of LGN neurons does not exhibit an expansive component. It implies that LGN neurons undergo saturation at all contrast levels above firing threshold. However, our measurements, as subsequently described, are at odds with this assumption. They demonstrate that LGN neurons exhibit expansive power-law nonlinearities when stimulated with low-contrast gratings. To quantify this observation, we measured the neural responses to gratings at different contrasts for our population of LGN neurons and calculated the best-fit Naka–Rushton function to these data. Figure 2, A, B, and C shows the best fits for n, c50, and ci, respectively, for 168 LGN cells. The mean and SE values for n, c50, and ci are 2.47 ± 0.162, 68.89 ± 5.45, and 27.42 ± 2.01, respectively. The medians for n, c50, and ci are 2.03, 38.87, and 17.72, respectively, and the modes are 1.758, 27.51, and 13.91, respectively. Mean values are indicated by unfilled arrows above the histograms. Note that these distributions are non-normal.
|
![]() | (5) |
|
|
|
Finally, for comparison to data from visual cortex, we show the exponent distribution for an expansive static nonlinearity for 160 cortical neurons in Fig. 5B. This distribution has a mean (open arrow) and SE of 2.4 ± 0.2, which is higher than that for the LGN population. Note that many neurons in the cortical population have large exponents (>4). This reflects more pronounced expansive nonlinearities in visual cortex as compared with LGN.
|
|
DISCUSSION |
|---|
|
Models of visual processing in the LGN and visual cortex have largely ignored any thalamic expansive nonlinearity (Bonin et al. 2005
; Li et al. 2006
; Priebe and Ferster 2006
). In the current study, we show that neurons in the LGN also exhibit a power-law expansive nonlinearity when activated by low-contrast visual stimuli. This expansive nonlinearity is likely to be the origin of expansive nonlinearity in membrane potentials of cortical neurons when stimulated at low contrasts (Ahmed et al. 1997
; Contreras and Palmer 2003
). We suggest, therefore, that expansive nonlinearity in the visual pathway originates early and is enhanced at various stages.
Although our measurements were made in LGN, the origin of expansive nonlinearities is probably in the retina. However, we should point out that direct neural input is not the only way to produce an expansive nonlinearity. It can also be produced as a by-product of neural noise. Random fluctuations in the membrane potential can make subthreshold responses "visible" in the presence of a threshold spiking mechanism, which can cause spike responses near threshold to simulate an expansive nonlinearity (Miller and Troyer 2002
). Given these factors, it is possible that an expansive nonlinearity originates in the retina and gradually increases in magnitude as transmission progresses in a feedforward manner along the visual pathway. Our data for LGN and visual cortex are consistent with this hypothesis.
The existence of a low-contrast expansive nonlinearity is not consistent with recent feedforward models of cross-orientation suppression in the primary visual cortex. In these models, response saturation with increasing contrast in LGN is thought to underlie cross-orientation suppression in primary visual cortex (Li et al. 2006
; Priebe and Ferster 2006
). For this to be true, the level of contrast saturation must match that of cross-orientation suppression at all contrast levels. We show here that, on average, the contrast–response function is expansive for contrasts <27% (see Fig. 2C). This would cause cross-orientation facilitation, not suppression, in the visual cortex. However, cross-orientation suppression is present in the visual cortex even with gratings at contrasts <27% (DeAngelis et al. 1992
; Freeman et al. 2002
; Li et al. 2006
). Therefore at low stimulus contrast, another mechanism must be involved in cross-orientation suppression (Li et al. 2006
).
Finally, it is relevant to consider possible consequences of a low-contrast expansive nonlinearity. In general, an expansive nonlinearity should contribute to a sharpening of tuning curves for different stimulus dimensions. This would occur by an increase in the slope of response functions so that small changes in the stimulus would generate relatively large changes in spike rates of neurons. This could apply to spatial frequency selectivity. It could also be relevant to orientation because tuning properties of orientation and spatial frequency of neurons in the visual cortex are related (Webster et al. 1990
). In a similar fashion, an expansive nonlinearity at low-contrast levels could increase contrast sensitivity by steeper slopes in contrast tuning functions. This would also yield low thresholds or high-contrast sensitivities. This accentuation of sensitivity could be highly significant in a practical sense because most visual performance occurs in a relatively low contrast environment (Mante et al. 2005
).
|
|
GRANTS |
|---|
|
|
|
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: R. D. Freeman, 360 Minor Hall, University of California, Berkeley, Berkeley, CA 94720-2020 (E-mail: freeman{at}neurovision.berkeley.edu)
|
|
REFERENCES |
|---|
|
Albrecht DG, Geisler WS. Motion selectivity and the contrast-response function of simple cells in the visual cortex. Vis Neurosci 7: 531–546, 1991.[Web of Science][Medline]
Albrecht DG, Hamilton DB. Striate cortex of monkey and cat: contrast response function. J Neurophysiol 48: 217–237, 1982.
Anzai A, Ohzawa I, Freeman RD. Neural mechanisms for processing binocular information. I. Simple cells. J Neurophysiol 82: 891–908, 1999a.
Anzai A, Ohzawa I, Freeman RD. Neural mechanisms for processing binocular information. II. Complex cells. J Neurophysiol 82: 909–924, 1999b.
Bonin V, Mante V, Carandini M. The suppressive field of neurons in lateral geniculate nucleus. J Neurosci 25: 10844–10856, 2005.
Bonin V, Mante V, Carandini M. The statistical computation underlying contrast gain control. J Neurosci 26: 6346–6353, 2006.
Carandini M. Amplification of trial-to-trial response variability by neurons in visual cortex. PLoS Biol 2: E264, 2004.[CrossRef][Medline]
Carandini M, Heeger DJ. Summation and division by neurons in primate visual cortex. Science 264: 1333–1336, 1994.
Carandini M, Heeger DJ, Movshon JA. Linearity and normalization in simple cells of the macaque primary visual cortex. J Neurosci 17: 8621–8644, 1997.
Chander D, Chichilnisky EJ. Adaptation to temporal contrast in primate and salamander retina. J Neurosci 21: 9904–9916, 2001.
Chichilnisky EJ. A simple white noise analysis of neuronal light responses. Network 12: 199–213, 2001.[Web of Science][Medline]
Contreras D, Palmer L. Response to contrast of electrophysiologically defined cell classes in primary visual cortex. J Neurosci 23: 6936–6945, 2003.
DeAngelis GC, Ohzawa I, Freeman RD. Spatiotemporal organization of simple-cell receptive fields in the cat's striate cortex. II. Linearity of temporal and spatial summation. J Neurophysiol 69: 1118–1135, 1993.
DeAngelis GC, Robson JG, Ohzawa I, Freeman RD. Organization of suppression in receptive fields of neurons in cat visual cortex. J Neurophysiol 68: 144–163, 1992.
Freeman TC, Durand S, Kiper DC, Carandini M. Suppression without inhibition in visual cortex. Neuron 35: 759–771, 2002.[CrossRef][Web of Science][Medline]
Gardner JL, Anzai A, Ohzawa I, Freeman RD. Linear and nonlinear contributions to orientation tuning of simple cells in the cat's striate cortex. Vis Neurosci 16: 1115–1121, 1999.[CrossRef][Web of Science][Medline]
Li B, Thompson JK, Duong T, Peterson MR, Freeman RD. Origins of cross-orientation suppression in the visual cortex. J Neurophysiol 96: 1755–1764, 2006.
Mante V, Frazor RA, Bonin V, Geisler WS, Carandini M. Independence of luminance and contrast in natural scenes and in the early visual system. Nat Neurosci 8: 1690–1697, 2005.[CrossRef][Web of Science][Medline]
Miller KD, Troyer TW. Neural noise can explain expansive, power-law nonlinearities in neural response functions. J Neurophysiol 87: 653–659, 2002.
Naka KI, Rushton WA. S-potentials from luminosity units in the retina of fish (Cyprinidae). J Physiol 185: 587–599, 1966.
Nykamp DQ, Ringach DL. Full identification of a linear-nonlinear system via cross-correlation analysis. J Vis 2: 1–11, 2002.[CrossRef][Medline]
Priebe NJ, Ferster D. Mechanisms underlying cross-orientation suppression in cat visual cortex. Nat Neurosci 9: 552–561, 2006.[CrossRef][Web of Science][Medline]
Reid RC, Victor JD, Shapley RM. The use of m-sequences in the analysis of visual neurons: linear receptive field properties. Vis Neurosci 14: 1015–1027, 1997.[Web of Science][Medline]
Sclar G. Expression of "retinal" contrast gain control by neurons of the cat's lateral geniculate nucleus. Exp Brain Res 66: 589–596, 1987.[CrossRef][Web of Science][Medline]
Solomon SG, Peirce JW, Dhruv NT, Lennie P. Profound contrast adaptation early in the visual pathway. Neuron 42: 155–162, 2004.[CrossRef][Web of Science][Medline]
Sutter EE. The fast m-transform: a fast computation of cross-correlations with binary m-sequences. SIAM J Comput 20: 686–694, 1991.[CrossRef]
Victor JD. The dynamics of the cat retinal X cell centre. J Physiol 386: 219–246, 1987.
Victor JD, Shapley RM. Receptive field mechanisms of cat X and Y retinal ganglion cells. J Gen Physiol 74: 275–298, 1979.
Webster MA, De Valois KK, Switkes E. Orientation and spatial-frequency discrimination for luminance and chromatic gratings. J Opt Soc Am A 7: 1034–1049, 1990.[Web of Science][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |