Figures
Abstract
When multiple stimuli appear together in the receptive field of a visual cortical neuron, the response is typically close to the average of that neuron’s response to each individual stimulus. The departure from a linear sum of each individual response is referred to as normalization. In mammals, normalization has been best characterized in the visual cortex of macaques and cats. Here we study visually evoked normalization in the visual cortex of awake mice using imaging of calcium indicators in large populations of layer 2/3 (L2/3) V1 excitatory neurons and electrophysiological recordings across layers in V1. Regardless of recording method, mouse visual cortical neurons exhibit normalization to varying degrees. The distributions of normalization strength are similar to those described in cats and macaques, albeit slightly weaker on average.
Citation: Zayyad ZA, Maunsell JHR, MacLean JN (2023) Normalization in mouse primary visual cortex. PLoS ONE 18(12): e0295140. https://doi.org/10.1371/journal.pone.0295140
Editor: Stéphane Charpier, Sorbonne Universite UFR de Biologie, FRANCE
Received: April 25, 2023; Accepted: November 16, 2023; Published: December 18, 2023
Copyright: © 2023 Zayyad et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The DOI necessary to access the repository information for our data is as follows: 10.12751/g-node.y6g8wz.
Funding: This work was supported by NIH F30 EY030733 (ZAZ), a Naomi Ragins Goldsmith Fellowship (ZAZ), NIH R01 EY022338 (JNM), NIH U19 NS107464 (JHRM), and NIH T32 GM007281 (ZAZ). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Nonlinearities in the responses of visual neurons, such as the saturation of the contrast response function and sublinear summation of responses when two visual stimuli are presented together, are referred to as normalization [1–4]. The responses of neurons in the visual cerebral cortex to multiple stimuli appearing simultaneously in their receptive field is generally sublinear in cats and macaques [2, 5]. This normalized neuronal response to multiple stimuli typically approximates a contrast-weighted average of a neuron’s responses to those stimuli presented individually and is found across sensory modalities. Normalization also describes the large modulations of responses that occur with attentional shifts between two stimuli in a neuron’s receptive field, and may contribute to pairwise spike correlations and gamma oscillations [6–10]. The prevalence of the observation of neuronal nonlinear responses regardless of the species studied or experimental details suggests that normalization is a fundamental feature of signal processing in neural systems.
Neurons can exhibit a range of different normalization strengths, sometimes depending on the specific stimulus configuration in the receptive field [8]. Variations in the strength of normalization have been shown to correlate with attention-related modulations [6–8] and the magnitude of pairwise noise correlations [9]. This latter observation suggests that normalization arises, at least in part, from local circuit mechanisms. Consequently, to better understand normalization and its role in sensory processing, it will be important to delineate the cell types and synaptic interactions that mechanistically underlie nonlinear responses. The mouse offers many advantages for such levels of analysis, including genetic and optical access [11, 12]. However, normalization is understudied in mouse sensory cortices [13, 14].
To study normalization we investigated the visual response properties of primary visual cortex (V1) neurons to cross-oriented gratings using both optophysiology and electrophysiology in awake mice [2, 3, 15]. We find that mouse V1 neurons exhibit visually-evoked normalization that varies in strength from neuron to neuron similar to that previously described in other species [5, 7]. Moreover, we find that this distribution of nonlinear responses is consistent across optophysiology and electrophysiology measurements and note that there are some differences from other species, including weaker overall strength of normalization in mice.
Materials and methods
Mice
C57BL/6J (B6) mice (Jax stock #000664) and mice expressing GCaMP6s under the Thy-1 promoter (Jax stock #025776; [16, 17]) were bred and maintained in accordance with the animal care and use regulations of the University of Chicago Institutional Animal Care and Use Committee. Heterozygous male and female mice (n = 12; n females = 6) aged P60-P200 were used across experiments.
Intrinsic imaging to identify primary visual cortex
Anesthetized adult mice were surgically implanted with a head bar and a 3.0 mm or 3.5 mm diameter chronic imaging window centered over V1, as previously described [18]. Mice were anesthetized using isoflurane (1.5–2.0% in 50% O2) with (electrophysiology experiments) or without (optophysiology experiments) ketamine (40 mg/kg, i.p.), and xylazine (2 mg/kg i.p.).
For all experiments, a retinotopic map of the visual cortex based on intrinsic imaging was produced under anesthesia (1.25% isoflurane in 50% O2) two weeks after the cranial implant surgery to guide calcium imaging or electrode placement. For optophysiology experiments, intrinsic imaging was performed as in Dechery et al. [18] using drifting bar stimuli. In electrophysiology experiments, intrinsic imaging was performed as in Cone et al. [19] using stationary Gabor stimuli.
Optophysiology
Excitatory neurons expressing GCaMP6s in Layer 2/3 of V1 were imaged using a custom-built two-photon microscope. The awake mouse was head-fixed and allowed to ambulate on a custom-made linear treadmill. Treadmill velocity was measured with a rotary encoder [18, 20]. Neurons were identified by an automated algorithm [18]. A typical field of view (FOV) contained 200 to 400 identified excitatory neurons scanned at 18.6–28.9 Hz (mean = 23.1 Hz; SD = 3.7) [18, 21]. Fluorescence data was preprocessed to calculate changes in fluorescence over baseline as previously described. Briefly, raw fluorescence time-series were thresholded to remove spike-induced fluorescence transients and smoothed with a 4th-order, 81-point Savitsky-Golay filter to obtain baseline. Fluorescence time-series were then normalized to percent change over baseline to obtain dF/F [18].
Electrophysiology
Silicone multi-channel electrodes (32 sites, model 4x8-100-200-177; NeuroNexus Technologies) were used. Electrodes sampled from depths 100–1000 μm in the lateral visual field representation of V1 of anesthetized mice (1.2–1.5% isoflurane in 50% O2) approximately normal to the cortical surface. After the electrode was inserted, anesthesia was removed. Mice were awake but head-fixed and seated in a custom plastic sled during recording sessions. Signals were amplified, filtered, and sampled at threshold crossings to identify putative spikes (BlackRock Inc.) and sorted offline (Offline Sorter, Plexon Inc.) to identify small multi-unit clusters.
Stimulus presentation
Stimuli were displayed on a calibrated video monitor positioned ~23 cm from the mouse’s right eye during data collection. Two sets of stimuli were shown: drifting gratings in one of 8 or 12 directions to measure direction tuning and preferences, and a cross-oriented stimulus to elicit normalization.
In the optophysiology experiments, drifting square-wave grating stimuli (0.04–0.06 cycles/degree, 2 Hz) in eight equally spaced directions were shown in a pseudorandom order to measure the responsivity and tuning of each neuron. The drifting grating stimuli were each presented for 4 s, separated by 3 s of a luminance matched gray screen to allow time for the calcium signal to return to baseline. Each stimulus direction was shown at least 20 times.
In the electrophysiology experiments, drifting squarewave grating stimuli (0.04 cpd, 2 Hz) in 12 equally spaced directions were shown in a pseudorandom order to determine the responsivity and tuning of each neuron. The drifting grating stimuli were shown for 1 s per direction, interleaved with 1 s of a luminance matched gray screen. Each stimulus direction was shown at least 20 times.
To assess normalization, cross-oriented stimuli (0.04 cpd, 2 Hz) composed of superimposed orthogonal gratings were presented 20 times [2, 15]. The cross-orientation stimuli consisted of orthogonal pairs of square-wave gratings of independently varying contrasts (0%, 6%, 12%, 25%, or 50%) superimposed to create a plaid. In six out of the eight optophysiology experiments, stimuli were counter-phased and in two the stimuli were drifting. Because responses to both types of stimuli were qualitatively similar, the data was pooled. In the electrophysiology experiments (n = 13), all stimuli were counter-phased.
Analysis
The data was analyzed using custom-written Matlab (MathWorks) scripts. For both optophysiology and electrophysiology data, direction and orientation tuning were assessed as described previously [18, 22, 23]. For the optophysiology data, neurons were classified as responsive if activity evoked by a stimulus was greater than activity recorded during the preceding gray epoch for at least one grating stimulus condition. Responsiveness was established using paired sample t-tests to compare the responses of neurons to the stimulus compared with the preceding gray for each direction; p values for responsiveness were Bonferroni corrected for the number of stimulus directions (p < 0.05/8). Responsive units were further tested for direction or orientation tuning using repeated measures analysis of variance (p < 0.05). For the electrophysiology data, units were tested for responsiveness using an analysis of variance comparing responses to each direction with activity during gray epochs, further evaluated with Tukey’s Honest Significant Difference test (p < 0.01). Responsive units were then tested for orientation or direction tuning using a Hotelling test or t-test of the trial vectors in orientation or direction spaces (p < 0.01) [18, 22]. Tuning properties of optophysiology and electrophysiology were further assessed by fitting an asymmetric Gaussian to visually evoked activity [18, 22].
Activity evoked by the plaid stimulus set was calculated using the mean grating response minus the mean spontaneous activity during luminance-matched gray periods. Only neurons that were responsive to one of the grating components at 50% contrast as determined by a Kruskal-Wallis test with alpha set to 0.1 and had a minimum positive response to the normalization stimulus set were included. This minimum positive response was set to 0.1 dF/F evoked response in optophysiology experiments and 0.3 spikes/s in electrophysiology experiments.
Optophysiological and electrophysiological neuronal evoked activity to cross-oriented stimuli were analyzed using a normalization index defined as:
Here, R[p] is the response to the 50% contrast in the preferred plaid component, R[n] is the response to the 50% contrast in the non-preferred plaid component, and R[p+n] is the response to the 50%/50% plaid composed of both the preferred and non-preferred components.
Results
Optophysiological and electrophysiological recordings of mouse V1 neurons
Calcium signals were recorded from a total of 2,120 excitatory neurons from layer 2/3 in eight imaging sessions from eight sites in five awake and head-fixed mice (Fig 1). Full contrast drifting grating stimuli were presented moving in eight directions, pseudo-randomly ordered, while imaging. Of the imaged neurons, 1,689 exhibited elevated activity in response to one or more gratings as compared to the response to the gray screen (Fig 1B). In a different set of experiments and a different set of mice, multi-channel silicone probes with 32 electrode sites spanning 700 μm of V1 were inserted approximately orthogonal to the cortical surface to record single and multi-unit activity from awake, head-fixed mice. Units were recorded from 13 penetrations from seven mice. Stable recordings were obtained from 114 single units and small multi-unit clusters that had visual stimulus evoked firing rates of at least 0.3 spikes/s. Of those, 72 units were more responsive to visual stimuli than to luminance matched gray.
a, Intrinsic imaging reveals retinotopic responses in the visual cortex. The white oval represents V1; the white square represents an example two-photon imaging field of view (FOV); the color map represents azimuth from approximately 20° (violet) to 160° (green). b, Representative calcium fluorescence traces from five example neurons. Different colors represent 4 s presentations of eight different direction stimuli. c, Example spikes represented as tics recorded using multichannel silicon probes. Different colors represent 1 s presentations of twelve different direction stimuli, which were presented in pseudo-random order.
To manage the impact of low amplitude fluorescence changes and its systematic effect on the normalization index metric, we only calculated a normalization index for cells with a minimum of a 10% increase in activity during plaid/preferred component presentation compared to the gray response. Of the responsive neurons recorded using optophysiology, 341 of 493 (69%) of neurons were significantly tuned to orientation or direction (Fig 2A–2D). Of the responsive units in the electrophysiology experiments, 63 of 72 (88%) were direction or orientation tuned (Fig 2E–2H). Because normalization has previously been studied in relation to feature strength and selection, only tuned neurons were included in further analyses [24].
a, An example tuning curve from a direction-selective neuron recorded optically in mouse V1. The red line is the mean dF/F of responses to drifting grating stimuli and the black line represents the tuning curve fit (r2 = 0.99). b, Mean response of the neuron shown in Fig 2A to the drifting grating stimulus in its least preferred (null) and most preferred stimuli (pref). The gray region marks the period when the stimulus was present. c-d, Same as Fig 1A and 1B, except averaged across tuned neurons (n = 341). e, An example tuning curve from a direction-selective unit recorded with a microelectrode in mouse V1. The blue line represents the mean response to drifting grating stimuli and the black line represents the tuning curve fit (r2 = 0.95). f, Mean response of the neuron shown in Fig 2E to the drifting grating stimulus in its least preferred (null) and most preferred stimuli (pref). The gray region represents the stimulus-on period. g-h, Same as Fig 2E and Fig 2F, except averaged across tuned electrophysiologically recorded units (n = 63). All error bars and bands represent ±1 SEM.
Cross-oriented stimuli elicit normalizing responses
We presented a set of counterphased cross-oriented grating stimuli [2, 3, 15] that included 25 combinations of orthogonal gratings at one of 5 contrasts ranging from 0% to 50% while recording from V1. To test whether these stimuli evoked normalizing responses, we used a normalization index ([7]; see methods). The index was calculated using the linear sum of responses to the component stimuli (gratings) minus the response to the cross-oriented (plaid) stimulus, divided by the sum of these two quantities. For this index, a value of 0.00 corresponds to a plaid response equal to linear summation and a value of 0.33 represents a plaid response that is the average of two grating responses.
When component stimuli were superimposed to create a cross-oriented stimulus neuronal responses were normalizing (Fig 3). We found a range of normalization index values. The median normalization index as measured by optophysiology was 0.15 (interquartile range (IQR) = 0.03–0.25) among tuned excitatory neurons in layer 2/3 of the mouse visual cortex (Fig 3). The bootstrapped 95% confidence interval (CI) for the median was 0.13–0.17.
a, Responses of an example neuron to plaid stimuli recorded with optophysiology. Four of the 25 contrast combinations are shown. The normalization index was calculated from the evoked activity. b, Mean responses to plaid stimuli across neurons (n = 341). All mean response SEMs were less than or equal to 0.01. c, Population histogram of normalization index values (median = 0.15, IQR = 0.03–0.25; bootstrapped median 95% CI = 0.13–0.17). The median is marked with a triangle. Dashed line, expected average of normalization index (0.33) for simple averaging. d, Visually-evoked responses to plaid stimuli in mouse V1 recorded using electrophysiology. e, Mean responses to plaid stimuli across neurons (n = 63). Data points and error bars correspond to mean rate ±1 SEM. f, Population histogram of normalization index values (median = 0.13, IQR = -0.01–0.26; bootstrapped median 95% CI = 0.06–0.19). The median is marked with a triangle.
Mice were allowed to freely ambulate in the optophysiology experiments, and did so ~5% of total stimulus presentation time. We found no difference in the medians or overall distribution of normalization indices regardless of inclusion or exclusion of running epochs (S1 Fig; with running: median = 0.15, IQR = 0.03–0.25; bootstrapped median 95% CI = 0.13–0.17; without running: 0.13, IQR = -0.02–0.25, bootstrapped median 95% CI = 0.10–0.15; Mann-Whitney U test, p = 0.09; two-sample Kolmogorov-Smirnov test, p = 0.16). Because running trials only composed ~5% of total trials, this data was not analyzed separately.
When evaluating V1 neuron responses using multi-channel silicon probes, recorded units showed normalizing responses when the stimulus components were superimposed to form plaids. The median normalization index was 0.13 (IQR = -0.01–0.26), indistinguishable from the value measured using optophysiology (Mann-Whitney U test, p = 0.68; t test p = 0.75; two-sample Kolmogorov-Smirnov test, 0.48). The bootstrapped 95% CI for the median was 0.06–0.19.
Discussion
Our findings reveal visually-evoked normalizing responses to cross-oriented plaid stimuli in awake mouse V1 as measured using optophysiology and electrophysiology. The findings from the optophysiology and electrophysiology experiments were highly consistent. This suggests that the advantages of both techniques can be deployed to dissect the circuitry of normalization and that comparable measurements can be expected across these recording modalities. The range of normalization strengths is consistent with tuned normalization, in which normalization is weighted differently for different stimulus feature values (orientations) [7, 9]. The overall normalization strength in mice reported here largely agrees with a recent report on the strength and range of responses to cross-oriented stimuli in the mouse primary visual cortex [13].
A recent study of macaque V1 suggests that normalization strength can be modulated by stimulus alignment with neurons’ preferred stimuli [25]. While direct comparison between these studies cannot be made due to stimulus set differences, subsets of our data that roughly correspond with that study’s analysis suggest that mouse V1 neurons exhibit weaker normalization than macaque V1 neurons. Additionally, the range of normalizing responses in mouse V1 is somewhat weaker than that reported in macaque area MT. Median normalization index values for electrophysiological measurements are ~0.2 in macaque MT [7], compared with 0.13 for electrophysiological measures in mouse V1.
With the myriad tools available in the mouse, this model will facilitate future efforts to establish the cell type(s) and network structure(s) that mechanistically produce normalizing responses.
Supporting information
S1 Fig. Running did not appreciably affect the median or overall distribution of normalization index.
a, Normalization index of neurons recorded during optophysiology experiments (n = 341). The triangle represents the median (median = 0.15, IQR = 0.03–0.25; bootstrapped median 95% CI = 0.13–0.17). b, Normalization index measured in the same neurons as in S1A Fig, after excluding trials in which mice ran. The triangle represents the median (0.13, IQR = -0.02–0.25, bootstrapped median 95% CI = 0.10–0.15). No difference in the medians or overall distribution of normalization indices was statistically detectable (Mann-Whitney U test, n.s.; two-sample Kolmogorov-Smirnov test, n.s.). A normalization index of 0.00 represents the expected normalization index for linear summation and the dashed line at 0.33 represents the expected normalization index for simple averaging.
https://doi.org/10.1371/journal.pone.0295140.s001
(TIF)
Acknowledgments
We thank Supriya Ghosh, Elizabeth de Laittre, Yuqing Zhu, and Chery Cherian for helpful comments on drafts of this manuscript.
References
- 1. Albrecht DG, Geisler WS. Motion selectivity and the contrast-response function of simple cells in the visual cortex. Vis Neurosci. 1991;7: 531–546. pmid:1772804
- 2. Carandini M, Heeger DJ. Summation and division by neurons in primate visual cortex. Science. 1994;264: 1333–1336. pmid:8191289
- 3. Carandini M, Heeger DJ, Movshon JA. Linearity and normalization in simple cells of the macaque primary visual cortex. J Neurosci. 1997;17: 8621–8644. pmid:9334433
- 4. Heeger DJ. Normalization of cell responses in cat striate cortex. Vis Neurosci. 1992;9: 181–197. pmid:1504027
- 5. Rust NC, Mante V, Simoncelli EP, Movshon JA. How MT cells analyze the motion of visual patterns. Nat Neurosci. 2006;9: 1421–1431. pmid:17041595
- 6. Lee J, Maunsell JHR. A normalization model of attentional modulation of single unit responses. PLoS ONE. 2009;4: e4651. pmid:19247494
- 7. Ni AM, Ray S, Maunsell JHR. Tuned normalization explains the size of attention modulations. Neuron. 2012;73: 803–813. pmid:22365552
- 8. Ni AM, Maunsell JHR. Spatially tuned normalization explains attention modulation variance within neurons. J Neurophysiol. 2017;118: 1903–1913. pmid:28701536
- 9. Verhoef B-E, Maunsell JHR. Attention-related changes in correlated neuronal activity arise from normalization mechanisms. Nat Neurosci. 2017;20: 969–977. pmid:28553943
- 10. Ray S, Ni AM, Maunsell JHR. Strength of gamma rhythm depends on normalization. PLoS Biol. 2013;11: e1001477. pmid:23393427
- 11. Luo L, Callaway EM, Svoboda K. Genetic Dissection of Neural Circuits: A Decade of Progress. Neuron. 2018;98: 256–281. pmid:29673479
- 12. Niell CM. Cell types, circuits, and receptive fields in the mouse visual cortex. Annu Rev Neurosci. 2015;38: 413–431. pmid:25938727
- 13. Barbera D, Priebe NJ, Glickfeld LL. Feedforward mechanisms of cross-orientation interactions in mouse V1. Neuron. 2022;110: 297–311.e4. pmid:34735779
- 14. Wilson NR, Runyan CA, Wang FL, Sur M. Division and subtraction by distinct cortical inhibitory networks in vivo. Nature. 2012;488: 343–348. pmid:22878717
- 15. Busse L, Wade AR, Carandini M. Representation of concurrent stimuli by population activity in visual cortex. Neuron. 2009;64: 931–942. pmid:20064398
- 16. Dana H, Chen T-W, Hu A, Shields BC, Guo C, Looger LL, et al. Thy1-GCaMP6 Transgenic Mice for Neuronal Population Imaging In Vivo. PLOS ONE. 2014;9: e108697. pmid:25250714
- 17. Chen T-W, Wardill TJ, Sun Y, Pulver SR, Renninger SL, Baohan A, et al. Ultra-sensitive fluorescent proteins for imaging neuronal activity. Nature. 2013;499: 295–300. pmid:23868258
- 18. Dechery JB, MacLean JN. Functional triplet motifs underlie accurate predictions of single-trial responses in populations of tuned and untuned V1 neurons. PLoS Comput Biol. 2018;14: e1006153. pmid:29727448
- 19. Cone JJ, Bade ML, Masse NY, Page EA, Freedman DJ, Maunsell JHR. Mice Preferentially Use Increases in Cerebral Cortex Spiking to Detect Changes in Visual Stimuli. J Neurosci. 2020;40: 7902–7920. pmid:32917791
- 20. Niell CM, Stryker MP. Modulation of Visual Responses by Behavioral State in Mouse Visual Cortex. Neuron. 2010;65: 472–479. pmid:20188652
- 21. Sadovsky AJ, Kruskal PB, Kimmel JM, Ostmeyer J, Neubauer FB, MacLean JN. Heuristically optimal path scanning for high-speed multiphoton circuit imaging. J Neurophysiol. 2011;106: 1591–1598. pmid:21715667
- 22. Mazurek M, Kager M, Van Hooser SD. Robust quantification of orientation selectivity and direction selectivity. Front Neural Circuits. 2014;8: 92. pmid:25147504
- 23. Levy M, Sporns O, MacLean JN. Network Analysis of Murine Cortical Dynamics Implicates Untuned Neurons in Visual Stimulus Coding. Cell Rep. 2020;31: 107483. pmid:32294431
- 24. Carandini M, Heeger DJ. Normalization as a canonical neural computation. Nat Rev Neurosci. 2011;13: 51–62. pmid:22108672
- 25. Guan S-C, Zhang S-H, Zhang Y-C, Tang S-M, Yu C. Plaid Detectors in Macaque V1 Revealed by Two-Photon Calcium Imaging. Curr Biol. 2020;30: 934–940.e3. pmid:32084400