arxiv: v1 [q-bio.nc] 29 Oct 2017

Transcription

1 State Dependence of Stimulus-Induced Variability Tuning in Macaque MT Joseph Lombardo,, Matthew Macellaio,, Bing Liu, Stephanie E. Palmer, *, with Leslie C. Osborne, * arxiv:7.8v [q-bio.nc] 9 Oct 7 Computational Neuroscience Graduate Program, University of Chicago, Chicago, IL USA Department of Organismal Biology and Anatomy, University of Chicago, Chicago, IL USA Neurobiology Graduate Program, University of Chicago, Chicago, IL USA Department of Neurobiology, University of Chicago, Chicago, IL USA Department of Physics, University of Chicago, Chicago, IL USA These authors contributed equally to this work. * sepalmer@uchicago.edu, osborne@uchicago.edu Abstract Average response properties of cortical sensory neurons such as firing rate and feature tuning are known to change with behavioral state. Here, we show anesthetic state also modifies the tuning of the variance of spiking responses in a fashion that is neither predicted by a Poisson model of spike count variance nor the neuron s tuning curve. Neurons recorded in the visual motion-sensitive middle temporal cortex (area MT) in sufentanil-anesthetized non-human primates tend to have lower average spike counts and greater trial-to-trial variability in response to constant direction stimuli compared to the alert state. Anesthesia-induced increases in response variation create supra-poisson variance-to-mean ratios (s) on time scales of ms. We find that alert-state MT neurons have higher mean firing rates, lower spike count variance, and largely sub-poisson variability that is stimulus-dependent and cannot be explained by firing rate differences alone. The absence of such stimulus-induced variability tuning in the anesthetized state suggests different sources of variability between states. A simple model explains qualitative shifts in the distribution of observed s in the two states via a suppression in the variance of gain fluctuations in the alert state. Author Summary Variability in neuronal firing is hypothesized to be important for read-out of the stimulus by downstream areas that control behavior, yet little is known about how the stimulus dependence of variability is modulated by behavioral context. Here, we measure the stability of variability under different anesthetic states. We show that variability can have qualitatively different stimulus tuning in different contexts. We develop a simple model for firing rate gain fluctuations that can explain changes in variability with behavioral state. Decreased variability during wakefulness may be an active process that suppresses synchronization and makes information transmission more reliable. /7

2 Introduction Many neurons modulate their firing rates smoothly in response to parametric changes in environmental stimuli, giving rise to the familiar tuning functions observed in many sensory areas, including motion-sensitive visual cortex in primates (Albright, 98; Maunsell and van Essen, 98). Single neuron contributions to stimulus decoding accuracy depend critically on this stimulus tuning, and their reliability determines their impact on decoding precision (Montijn et al., ). Much is made of how a population code might be either robust or sensitive to neuronal variability, but an important aspect of variability is often left out of the discussion: how does neuronal variability, itself, depend on the stimulus? This question becomes particularly important in the context of how response reliability impacts information transmission (Brunel and Nadal, 998; Butts and Goldman, ; Josić et al., 9) and decoding (Ma et al., ). The impacts of variability on population decoding can be analytically derived for the case of an independent or correlated population of Poisson neurons (see for example Brunel and Nadal (998)), where the mean and variance of the response are the same, yielding a =. The Poisson model of spike generation replicates many of the features of cortical spike trains recorded under some conditions (Shadlen and Newsome, 998), but the prediction is often violated in real sensory neurons (Arabzadeh et al., ; Carandini, ; Kara et al., ; Montgomery and Wehr, ; Softky and Koch, 99). MT variability changes under anesthesia (Goris et al., ), and is usually > in the anesthetized state (Buracas et al., 998). Decreases in are observed at the onset of stimulation, alongside decreases in neuronal correlation, in MT and throughout cortex (Churchland et al., ). s < have been observed in alert subjects (Barbour and Wang, ; DeWeese et al., ). To assess how robust or sensitive the sensory population code may be to state dependent changes in spiking variability, we measure the stimulus tuning of the. In this paper, we show that scaled variability in neuronal response is tuned to stimulus features, as previously observed for drifting grating stimuli in V and MT (Ponce-Alvarez et al., ; Priebe and Ferster, ; Sadagopan and Ferster, ), and can be qualitatively modulated by the anesthetic state of the subject. Awake responses display sub-poisson variability that is inversely tuned to the stimulus, decreasing at the preferred direction of the cell. Anesthetized responses show flat, supra-poisson tuning. A simple model shows how changes in a single parameter, the size of gain fluctuations in the response, can give rise to both qualitative regimes of tuning. We also show how these changes in the tuning of variability enhance stimulus discriminability in the awake state. Results Our results show that variability in MT neurons as well as the mean rate can be tuned to the direction of motion of the visual stimulus (Fig ), beyond a simple scaling with the mean as predicted by a Poisson model. Neurons in MT show increased firing rates in the awake state (Fig C), while the variance in spike count from cells in anesthetized and awake subjects show significant overlap (Fig D). The variance-to-mean ratio, or, is qualitatively different in these two states (Fig ). s from neurons recorded in MT display tuning to the direction of motion (Fig ). This is particularly pronounced in recordings from awake subjects, but can also be seen in some cells recorded from anesthetized subjects. We show that this (FF) tuning is not explained solely by increases in firing rate in the awake state (Fig ), can be accounted for by simple gain fluctuation model (Fig ), and supports lower thresholds for direction discrimination from populations of MT neurons (Fig ). /7

3 Anesthetic state modulates mean firing and in MT neurons Our recordings from MT neurons in awake, behaving primates show a pronounced increase in their mean firing rate compared to neurons recorded in MT in anesthetized primates (Fig C). The mean and variability of responses were computed by counting action potentials within a ms window starting from motion onset in the visual stimulus, for each trial. An increase in neural responsiveness in LGN neurons was previously observed in alert versus anesthetized primates (Alitto et al., ), but in our study, we also characterized the variance in the response. The variance is tuned to the direction of motion in recordings from both awake and anesthetized subjects, and displays a largely overlapping distribution of magnitude (Fig D). The maintenance of roughly equal magnitude spike count variance with increasing mean firing rate implies a substantial change in the in these two states. This is illustrated in Fig. Very few neurons in either state display a of (dashed unity line in Fig A), and most awake neurons fall below a FF of, while most anesthetized neurons display FF s above. The distribution of measured FF s taken across all cells in our recordings shows a significant rightward shift in the data from anesthetized subjects (Fig B). While FF increases with spike count in anesthetized data, the FF shows a slight decrease with increasing spike count in the awake data (Fig C). Overall, this indicates a marked shift in the scaling of variability with responsiveness in these two states. While firing rates are lower in area MT in anesthetized subjects, the relative variability of this response is higher, potentially making readout of direction of motion much more ambiguous. In awake subjects, is tuned to motion direction All of the neurons showed directional tuning of both the mean and trial-to-trial variance (Fig C and D). However, in a simple linear regression, firing rate fails to explain much of the spike count variance in the awake experiments, with R =.7. We observed a dependence of the of the spike count on stimulus direction. This tuning was similar in shape and amplitude to that observed by Ponce-Alvarez et al. () using drifting grating stimuli. Neurons recorded in anesthetized subjects showed no such stimulus dependence of the, and a linear regression of spike count variance on firing rate explained more variance, with R =. (see Materials and Methods for statistics). Figs C-D show the difference in tuning between the populations of neurons in the awake and anesthetized states. The tuning shown here is averaged over the entire population. Plotting the median through the populations at each stimulus direction reveals a clear tuning in the awake state, whereas in the anesthetized state, the median population is flat across all directions. We next quantify this tuning of the. tuning index We quantified the degree of tuning in indiviual neurons by defining a Fano factor tuning index (FFTI). FFTI = FF orth FF pref FF orth + FF pref () For each neuron recorded, FF pref is the for the preferred stimulus direction of the neuron, and FF orth is the for the orthogonal stimulus directions. This tuning index is analogous to standard methods for quantifying direction /7

4 selectivity in first-order response statistics (Priebe et al., ), and takes values from - to. A positive FFTI indicates a decrease in for the preferred stimulus direction relative to the off-preferred stimuli (U-shaped or M-shaped tuning), and a negative FFTI indicates an increase in for the preferred stimulus direction (Gaussian-like tuning). The distribution of FFTI in the awake state is broader and more positive than in the anesthetized state (Fig E). The population of MT neurons recorded in the awake state ( FFTI =.7) showed a significantly greater tuning on average compared to the anesthetized state recordings ( FFTI =.) (two-tailed t-test, p =.). We chose FFTI. as a cutoff for a positively tuned. In the awake state, % of neurons had positively tuned s (7/). In the anesthetized state, only % of neurons had positively tuned s (/). One neuron from the anesthetized group was excluded from the FFTI calculations as it did not emit any spikes to the orthogonal stimulus direction, and thus FF orth was not defined. Firing rate changes do not explain differences in tuning The firing rates in our awake MT recordings are much higher on average than those in the anesthetized recordings (Fig A). We ask whether this difference in the response distribution is sufficient to explain the observed difference in tuning. It is notable that the relationship between and firing rate differs in the awake and anesthetized recordings. increases with spike count in the anesthetized state, and decreases with spike count in the awake state, over the same range of firing rates (Fig C). This suggests that changes in mean spike count alone can explain the observed differences in and its tuning. Mean-matched analysis: Ideally, we would like to be able to exclude all effects changes of mean rate might have on observed s. One way to achieve this is through analyzing data from mean-matched s, as in (Churchland et al., ). This method excludes data from the groups being compared so as to match their distributions of mean spike count. Mean matching does not maintain cell identity, and thus can only describe population-wide rather than cell-specific tuning. We generated histograms of the spike counts for both states at the preferred stimulus direction and at the orthogonal directions. We then randomly excluded points from each dataset until the histograms of spike counts in the awake and anesthetized states matched (Fig A). With this reduced data set, we calculated the mean for both stimulus conditions and both states. Population FFTI was then calculated on the mean values. This process was repeated one million times. The difference in tuning was maintained after mean matching, with FFTI awake =. (sd =.7) and FFTI anesth =. (sd =.7) (Fig B). The mean-matched population s were fit to cosine tunings (Fig C). The best cosine fits corresponded to FFTI awake =. and FFTI anesth =.9 for the mean-matched populations. Part of this difference is due to the high average in the anesthetized state, which leads to a greater denominator in the normalization of FFTI. With this in mind, we can examine the raw change in, where FF = FF orth FF pref. In awake recordings we see FF =. (sd =.) and for anesthetized, FF =.9 (sd =.). Thus, the higher spike counts recorded in the awake state are not the source of the difference in tuning. Effect of tuning bandwidth on tuning While differences in mean spike count fail to explain the differences in observed Fano factor tuning between states, this does not rule out the possibility tuning /7

5 arises from first-order statistics. It is possible that the second-order statistics are comparable between states, but differences in mean spike count tuning bandwidths determine the tuning. If tuning curves are narrower or steeper in the awake state, this could result in different tuning curves when normalized by the same variance tuning curves. Conversely, if the mean spike count tuning curves are the same in both states, tuning could arise from relatively broader or flatter count variance tuning curves. An ANOVA was used to compare the widths of mean spike count tuning curves and spike count variance tuning curves in both states to the FFTI. There was no effect of tuning curve width or variance tuning width on FFTI, and no interaction with conscious state, suggesting that arises from systematic differences in the structure of the variance rather than differences in the first-order tuning properties. A simple model accounts for changes in tuning Several models of tuning have been developed to explain the stimulus-dependent effects of neuronal response variability (Franke et al., ; Ponce-Alvarez et al., ; Zylberberg et al., ), while other models have modeled variability as rate-dependent, or due to stimulus-independent gain fluctuations (Goris et al., ). Our model combines contributions from rate dependent variability as well as rate-independent gain-fluctuations to reproduce the diversity of tunings observed in the data. A simple and common model of spike counts in cortical neurons is given by a Poisson process, or a Poisson mixture model (Goris et al., ). A Poisson process is characterized by a of. Including a multiplicative gain to describe the underlying rate of the Poisson process results in a super-poisson spike count distribution, with a that increases with firing rate. In the anesthetized state, we found increasing with mean spike count, and a greater than one (Fig ). These observations are consistent with a Poisson mixture model, and agree with previous experiments (Goris et al., ; Shadlen and Newsome, 998). Our recordings in the awake state, however, show a sub-poisson spike count, with a Fano factor that decreases with increased firing rate. The inverse rate-dependence of Fano factor was previously documented by (Ecker et al., ) in V. Our model captures the difference in rate-dependence of between states as well as the difference in tuning. We model the mean mean spike count response of neurons by a tuning curve f(θ), where θ is the stimulus direction. On any given trial, the mean rate is scaled by a multiplicative gain, g, yielding the underlying rate for the neuron on that trial µ = f(θ) g. () The gain g is taken to be a gamma-distributed variable with mean of. The spike count for a given trial sampled from a Gaussian distribution with a mean value µ, and a variance µ λ, where λ is an intrinsic property of the cell that determines how variance scales with firing rate. For a fixed value of g, a value of λ less than corresponds to a less than, and a value of λ greater than corresponds to a greater than. If we allow g to fluctuate, the variance of the spike count x is given by var(x θ) = f(θ) λ g λ + f(θ) var(g). () Because g = and the variance of g is relatively small, we can approximate g λ. Thus, we can approximate the spike count variance var(x θ) f(θ) λ + f(θ) var(g). () /7

6 For small values of var(g), i.e., small gain fluctuations, the variance is dominated by the first term. This results in a U-shaped tuning for λ <. As the gain fluctuations increase, the second term dominates the spike count variance, resulting in a flat or Gaussian tuning for the same value of λ. This model suggests that a difference in the amplitude of gain fluctuations is sufficient to reverse or eliminate the observed tuning. The gain fluctuations also explain the respective increase and decrease of with spike count observed in the anesthetized and awake states. For small gain fluctuations (var(g) ) and λ <, the for a given firing rate, r, is r λ, which decreases as r increases. For stronger gain fluctuations, the scales linearly with r, increasing with increasing rate. We used this model to fit the distribution of observed FFTI for the given distributions of spike counts in the two experiments. We found the best-fit values of α and var(g) for each population that generated distributions of s that matched the observed FFTI distributions (Fig A,B). The distributions of s generated by the model (Fig C,D) were similar to those observed in the data (Fig C,D). The fitted values of the intrinsic variability parameter were found to be α =. in the awake condition and α =.7 in the anesthetized condition. While these values are different, they both result in similarly U-shaped tuning curves in the absence of gain fluctuations. The model predicts a U-shaped tuning for small gain fluctuations in both states, and the reverse tuning for larger gain fluctuations. The variance of the gain parameters were var(g) =.9 in the awake condition and var(g) =.7 in the anesthetized condition, nearly an order of magnitude difference. It is this difference gain variance that results in the qualitative difference in the Fano factor tuning between the two conditions. To confirm that these results were not affected by the difference in α, we fit the same model using the same optimized value of α =. for both populations of tuning curves. The results were comparable to those shown, and the difference in gain variance between states was even larger. We also fit the model for individual neurons to match their values of FFTI. For our tuning curves f(θ), we fit the observed tuning curves to Gaussians. We then averaged the parameters within state, resulting in average awake and anesthetized tuning curves (Fig E, inset). The model parameters were then fit, such that all neurons had the same value of α and var(g) was fit separately for each neuron to reproduce the distribution of FFTI. This model predicts a distribution of gain fluctuations that is much larger for the anesthetized condition than for the awake condition (Fig F). We similarly found the gain variance distribution of the model using the observed tuning curves (Fig G,H). In both cases, the model predicts the proportion of neurons with high values of gain variance is much greater in the anesthetized condition. Impact on information transmission It seems reasonable to assume that higher in the anesthetized state will result in lower rates of information transmission and stimulus discriminability from these neural populations. Conversely, the reduction in relative variability at the preferred direction for most cells recorded from awake, behaving subjects seems to imply enhanced information transmission and decoding. To test this intuition, we simulated responses of MT neuron populations and varied the stimulus-dependence of the noise. The performance of the populations was quantified by the Fisher information, which determines the bound on the performance of an unbiased estimator reading out the population code. We first tested a homogeneous population with identically shaped tuning curves (Fig A). The first-order response properties of the simulated populations were matched to the mean statistics of MT neuron responses measured in awake macaques in the first /7

7 ms following the onset of stimulus motion. The mean tuning curve was fit to a vonmises function as in Ecker et al. (); Zylberberg et al. () and rotated such that the preferred directions evenly tiled all directions. In order to separate the effects of the magnitude of the spike-count variability from the stimulus-dependent part of the variability, we imposed various s tunings on the populations while keeping the average across all directions constant (FF=). The tunings were modeled as von Mises functions and imposed directly, rather than using the spike count model described above, so that the average could be held constant. The positive and negative tunings are symmetric by reflections over the line FF=. Short-range noise correlations were also imposed on the population such that neurons with similar preferred directions were more strongly correlated, and correlation strength decreased monotonically with the difference in preferred direction. For the results shown in Figure, the maximum pairwise correlation c max was chosen to be., with an average pairwise correlation of c =.8. In the homogeneous population, the population with the positive FFTI performed better than a population with either a flat tuning or a negative FFTI (Fig B). Performance was quantified by the Cramer-Rao bound, given by J, where J is the Fisher information. Because MT coding performance is seen as a constraint on the precision of downstream smooth pursuit behavior, we can compare the decoding performance of the population to the experimental data. Therefore, we estimated the size of the population needed to reach a Cramer-Rao bound of. The population with positive FFTI crossed this threshold with neurons (Fig B). The population in which the variability was structured to give a constant (FFTI = ) required 7 neurons to reach same level of discrimination, and the FFTI < population appeared to asymptote before reaching Cramer-Rao bound of. We also looked at the time course of decoding performance by finding the average tuning curves after cumulative intervals of time after motion onset (Fig C). We assumed a population size neurons and c max =.. We estimated the time after motion onset when the Cramer-Rao bound for these popuations would cross the threshold. The population with FFTI > crossed the threshold after 7 ms, while the FFTI = population crossed the threshold in ms. This shows that a population with stimulus-dependent variability tuning can reach the same level of stimulus discrimination in less time. The FFTI < population did not reach threshold within the ms analysis window. We then performed the same analyses for population of heterogeneous neurons drawn from the distribution of measured tuning curves (Fig D-F). Adding heterogeneity would be expected to improve overall coding performance relative to the homogeneous case in a correlated population of neurons (Shamir and Sompolinsky ), and we found this to be the case in our simulation. The tuning curves were again fit to von Mises functions that matched the cumulative spike counts ms following stimulus motion onset and the preferred directions were spaced evenly to tile the space of stimulus directions. Short-range correlations and stimulus dependent variability were imposed in the same manner as with the homogeneous population. We found that the effect of stimulus-dependent variability was smaller in the heterogeneous population relative to the homogeneous population, but the same trend was observed. The population with FFTI > required only neurons to reach the Cramer-Rao bound of, while the population with constant Fano Factor required neurons and the population with FFTI < required 7 neurons to reach the same level of stimulus discriminability. Again, estimating the time course of stimulus discriminability in populations of neurons, the FFTI > population reached the discriminability threshold in ms 7/7

8 while the FFTI = population reached the same threshold in ms and the FFTI < population reached the threshold in ms. Discussion This work quantifies how variance in spiking responses to moving dot stimuli are tuned to the direction of motion of the stimulus, and how this might underlie enhanced stimulus readout. Under usual assumptions of Poisson firing statistics, the variance in spiking should scale with the mean rate in a one-to-one fashion. While variations on the Poisson assumption can give rise to other slopes in the variance-to-mean ratio, these do not produce tuning in the. We observe significant tuning of the during wakefulness, which may enhance stimulus readout. Lower variance at and around the preferred direction of each neuron in a population leads, unsurprisingly, to a finer discrimination threshold, as estimated via the Fisher information. Fewer cells are needed to achieve the same level of direction discriminability, and stimulus information is available earlier in populations that have this kind of tuning. Comparing results from anesthetized subjects to those that were awake and behaving, we were able to show that a model in which a single parameter, the gain variance, could broadly explain the shifts in tuning observed across the population of MT neurons sampled. This suggests that a higher gain state underlies the observed flat tuning of the FF during anesthesia. This aligns with results that suggest that anesthesia corresponds to a more synchronous mode of brain coupling. The active suppression of gain fluctuations during wakefulness may be directed specifically at those neurons encoding stimulus variables in an active task (here, fixation). This suggests a simple knob that attentional modulation can turn to drive more reliable decoding of the stimulus. Reducing gain fluctuations may require processes, such as activation of inhibitory networks, that are metabolically costly and are only engaged when needed to drive behavior. Materials and Methods Experimental methods We made extracellular single-unit microelectrode recordings of MT neurons in both alert, behaving and anesthetized monkeys. All procedures were performed in compliance with Institutional Animal Care and Use Committee guidelines. The anesthetized motion-step experiments consisted of recordings from four adult male macaques (Macaca fasicularis). Animals were implanted with a head-restraint and a craniotomy was performed under isoflurane anesthesia using sterile technique. During the experiments, anesthesia was maintained with the opioid sufentanil, and paralysis was induced with vecuronium bromide to minimize eye movements. Unit recordings were made using tungsten-in-glass microelectrodes. Visual stimuli comprised randomly drawn patterns of white dots that moved coherently against the dark screen of analog oscilloscopes (models A and B, P Phosphor; Hewlett-Packard, Palo Alto, CA). The size and position of the stimulus aperture was chosen to maximally excite each isolated unit, as was motion speed. The direction of stimulus motion was pseudo-randomly chosen from a set of at least directions that spanned ±9 around the preferred direction including increments. A stationary random dot stimulus was interleaved with the motion stimuli in of recorded neurons. For these neurons, stimuli were repeated - times. The remaining neurons had - stimulus repetitions. In all, dot textures appeared and remained stationary for ms, translated for ms, and were again stationary for ms. When the motion of a dot 8/7

9 carried it outside of the aperture, it was randomly positioned along the leading aperture edge to maintain dot number. Trials were separated by a brief pause of -s. Spike waveforms were sampled at khz and isolation was aided by a window discriminator. Spike times were determined by threshold crossings. Data from three of the four monkeys have been previously published (Osborne et al.,, 8). Analytical methods We computed the spike count in the ms following motion onset, while the dot textures were translating. We computed the of the spike count (variance divided by mean) as a function of stimulus direction using a repeated-measures ANOVA on the null hypothesis of a constant across directions. Stimulus direction was found to have an effect on in the awake experiments (p = 7. 7 ) but not in the anesthetized experiments (p =.). Because there was a narrower range of stimulus directions in the anesthetized experiments, we limited both data sets to stimulus directions within ±9 of the preferred direction. Tuning widths were obtained by fitting a Gaussian curve to the spike counts as a function of direction for each neuron. The tuning width was taken to be the standard deviation of the Gaussian of best-fit. Direction selectivity was computed as a direction index (DI), where DI = r pref r orth r pref + r orth. () Here, r pref is the mean response to the preferred stimulus direction and r orth is the mean response to the orthogonal directions. FFTI was defined analogously. FFTI = FF orth FF pref FF orth + FF pref. () The procedure for mean-matched is adapted from Churchland et al.. The distribution of mean spike counts was computed for each condition. The distribution of spike counts for each stimulus direction was approximated by binning values in evenly spaced bins that spanned the range of responses. For each bin in which one condition had more data points than the other, data points were randomly discarded from that condition until both distributions matched. The was calculated for each neuron in that stimulus condition. The mean for each stimulus condition was used to calculate FFTI for each subject condition. This resampling procedure was repeated one million times to generate a distribution of Fano factor tunings. The resampled values were fit to cosine curves by a least squares fitting procedure. Acknowledgments This study was supported by grants to L.C.O. from the Alfred P. Sloan Foundation, Whitehall Foundation, Brain Research Foundation, NEI EY7 and NSF IOS-7, grants to S.E.P. from the Alfred P. Sloan Foundation and NSF CAREER-7, and support for J.L. from the National Science Foundation under Grant No. DGE-97. References Albright, T. D. (98). Direction and orientation selectivity of neurons in visual area MT of the macaque. Journal of Neurophysiology, ():. 9/7

10 Alitto, H. J., Rathbun, D. L., Martin Usrey, W., et al. (). A comparison of visual responses in the lateral geniculate nucleus of alert and anaesthetized macaque monkeys. The Journal of physiology, 89(): Arabzadeh, E., Panzeri, S., and Diamond, M. E. (). Whisker vibration information carried by rat barrel cortex neurons. Journal of Neuroscience, ():. Barbour, D. L. and Wang, X. (). Auditory cortical responses elicited in awake primates by random spectrum stimuli. Journal of Neuroscience, (8):79 7. Brunel, N. and Nadal, J.-P. (998). Mutual information, fisher information, and population coding. Neural computation, (7):7 77. Buracas, G. T., Zador, A. M., DeWeese, M. R., and Albright, T. D. (998). Efficient discrimination of temporal patterns by motion-sensitive neurons in primate visual cortex. Neuron, (): Butts, D. A. and Goldman, M. S. (). Tuning curves, neuronal variability, and sensory coding. PLoS biology, ():e9. Carandini, M. (). Amplification of trial-to-trial response variability by neurons in visual cortex. PLoS biology, (9):e. Churchland, M. M., Yu, B. M., Cunningham, J. P., Sugrue, L. P., Cohen, M. R., Corrado, G. S., Newsome, W. T., Clark, A. M., Hosseini, P., Scott, B. B., Bradley, D. C., Smith, M. A., Kohn, A., Movshon, J. A., Armstrong, K. M., Moore, T., Chang, S. W., Snyder, L. H., Lisberger, S. G., Priebe, N. J., Finn, I. M., Ferster, D., Ryu, S. I., Santhanam, G., Sahani, M., and Shenoy, K. V. (). Stimulus onset quenches neural variability: a widespread cortical phenomenon. Nature Neuroscience, ():9 78. DeWeese, M. R., Hromádka, T., and Zador, A. M. (). Reliability and representational bandwidth in the auditory cortex. Neuron, 8(): Ecker, A. S., Berens, P., Cotton, R. J., Subramaniyan, M., Denfield, G. H., Cadwell, C. R., Smirnakis, S. M., Bethge, M., and Tolias, A. S. (). State Dependence of Noise Correlations in Macaque Primary Visual Cortex. Neuron, 8(): 8. Franke, F., Fiscella, M., Sevelev, M., Roska, B., Hierlemann, A., and Azeredo da Silveira, R. (). Structures of Neural Correlation and How They Favor Coding. Neuron, 89():9. Goris, R. L. T., Movshon, J. A., and Simoncelli, E. P. (). Partitioning neuronal variability. Nature Neuroscience, 7():88 8. Josić, K., Shea-Brown, E., Doiron, B., and de la Rocha, J. (9). Stimulus-Dependent Correlations and Population Codes. Neural Computation, ():77 8. Kara, P., Reinagel, P., and Reid, R. C. (). Low response variability in simultaneously recorded retinal, thalamic, and cortical neurons. Neuron, 7():. Ma, W. J., Beck, J. M., Latham, P. E., and Pouget, A. (). Bayesian inference with probabilistic population codes. Nature neuroscience, 9(): 8. Maunsell, J. H. and van Essen, D. C. (98). Functional properties of neurons in middle temporal visual area of the macaque monkey. I. Selectivity for stimulus direction, speed, and orientation. Journal of Neurophysiology, 9():7 7. /7

11 Montgomery, N. and Wehr, M. (). Auditory cortical neurons convey maximal stimulus-specific information at their best frequency. Journal of Neuroscience, ():. Montijn, J. S., Vinck, M., and Pennartz, C. M. A. (). Population coding in mouse visual cortex: response reliability and dissociability of stimulus tuning and noise correlation. Frontiers in Computational Neuroscience, 8:. Osborne, L. C., Bialek, W., and Lisberger, S. G. (). Time Course of Information about Motion Direction in Visual Area MT of Macaque Monkeys. The Journal of neuroscience : the official journal of the Society for Neuroscience, ():. Osborne, L. C., Palmer, S. E., Lisberger, S. G., and Bialek, W. (8). The Neural basis for combinatorial coding in a cortical population responseeural basis for combinatorial coding in a cortical population response. Journal of Neuroscience, 8():. Ponce-Alvarez, A., Thiele, A., Albright, T. D., Stoner, G. R., and Deco, G. (). Stimulus-dependent variability and noise correlations in cortical MT neurons. Proceedings of the National Academy of Sciences, (): 7. Priebe, N. J. and Ferster, D. (). Mechanisms of neuronal computation in mammalian visual cortex. Neuron, 7():9 8. Priebe, N. J., Lisberger, S. G., and Movshon, J. A. (). Tuning for Spatiotemporal Frequency and Speed in Directionally Selective Neurons of Macaque Striate Cortex. The Journal of neuroscience : the official journal of the Society for Neuroscience, ():9 9. Sadagopan, S. and Ferster, D. (). Feedforward origins of response variability underlying contrast invariant orientation tuning in cat visual cortex. Neuron, 7():9 9. Shadlen, M. N. and Newsome, W. T. (998). The Variable Discharge of Cortical Neurons: Implications for Connectivity, Computation, and Information Coding. The Journal of Neuroscience, 8(): Softky, W. R. and Koch, C. (99). The highly irregular firing of cortical cells is inconsistent with temporal integration of random epsps. Journal of Neuroscience, ():. Zylberberg, J., Cafaro, J., Turner, M. H., Shea-Brown, E., and Rieke, F. (). Direction-Selective Circuits Shape Noise to Ensure a Precise Population Code. Neuron, 89():9 8. /7

12 A B motion onset motion offset C spike count Spikes/s Time (ms) D spike count variance stim. direction stim. direction Fig. Mean and variance tuning of MT responses to % coherent random dot stimuli. (A) A cartoon of the experimental setup in awake experiments. Head-fixed macaques fixate on a central target while a random dot stimulus is presented to the receptive field of the isolated neuron. The dots translate coherently behind an aperture to probe neuronal tuning. (B) Sample peri-stimulus time histogram (PSTH), showing the time course of firing rate in a representative neuron to different stimulus directions. Dashed lines indicate motion onset and motion offset. Arrows at the right indicate motion direction. (C) Average spike count across neurons in the awake experiments (blue) and anesthetized (orange) by direction. Circles indicate mean, error bars indicate standard deviation, and the solid traces indicates a Gaussian best fit. Responses are aligned such that the preferred stimulus direction for each neuron is taken to be. (D) Average spike count variance across populations of neurons in the same manner as (C). The solid traces indicate a best-fit quadratic curve. /7

13 A B. spike count variance 8 7 Awake Anesthetized C probability density spike count spike count Fig. s differ significantly between MT responses recorded in awake versus anesthetized primates. (A) The relationship between spike count and spike count variance in awake (blue circles) and anesthetized (orange circles) experiments. Each circle indicates a single neuron under a single stimulus condition. Histogram above shows distribution of spike counts observed in both experiments. The histogram to the right shows the distribution spike count variance values in both experiments. Stimulus directions shown are within of the preferred motion direction, in increments. (B) Distribution of spike count s for awake (blue) and anesthetized (orange) experiments. Solid trace indicates the best-fit lognormal distribution for each case. (C) Distribution of s as a function of spike count for awake experiments (blue) and anesthetized experiments (orange). Black circles and triangles indicate the median binned by spike count in awake and anesthetized experiments, respectively. Error bars show the standard deviation of s in each bin. The dashed lines are the best-fit quadratic curve for each case. /7

14 A C stimulus direction awake FFTI =.9 D stimulus direction (degrees) stimulus direction (degrees) 7 B stimulus direction anesthetized FFTI = -. E probability density. FF FFTI FF 9 FF stim. dir.. Aw. An. Aw. Fit An. Fit Fig. s show significant tuning to motion direction in the awake state. across stimulus directions in an example neuron from (A) awake and (B) anesthetized recordings in MT. Black dots indicate by direction. Blue and orange traces show values smoothed by a moving window. for respective experimental conditions. Red dashed lines indicate FF =. (C) Black dots show the s by stimulus direction for all neurons neurons in the awake experiments. Stimulus directions are aligned such that the preferred stimulus direction of each neuron is at. The median (green) is significantly tuned by direction. Gray dashed line indicates median in response to null stimulus. Red dashed line is at FF = is as predicted for a Poisson process. (D) Same as (C) but for the anesthetized experiments. The median is not tuned by stimulus direction. (E) The distributions of tuning indices (FFTI) from the awake (blue) and anesthetized (orange) experiments. Arrows indicate the FFTI of the individual neurons shown in and (A) and (B). The dashed line indicates the tuning threshold of.. The cartoon at top right shows how FFTI is calculated using the median trace from (C). /7

15 A count -9 count count 9 C FFTI... spike count spike count spike count -. awake anesthetized B D 8 probability density stimulus direction FFTI Fig. Mean-matched s preserve variability tuning. (A) Histograms of spike counts by neuron for a given stimulus direction (left: -9 degrees; center: preferred direction; right: +9 degrees) for awake (blue, throughout) and anesthetized (orange, throughout) recordings. The overlapping area is shown in gray. (B) Sample mean-matched s, with spike count distributions corresponding to the gray histograms in (A). Each circle is a sample mean. Dashed traces indicate the mean Fano factor across directions. Solid traces indicate the mean across bootstrapped samples. (C) tuning index (FFTI) of bootstrapped resamples for awake (FFTI =., s.d. =.) and anesthetized (FFTI =., s.d. =.7). (D) Distribution of FFTI in for mean-matched distributions. Solid traces indicate a best-fit Gaussian curve. /7

16 -9-9 Stim. Direction -9-9 Stim. Direction A probability density... Spike Count Aw. Model An. Model Aw. Model Fit An. Model Fit Aw. Data Fit An. Data Fit B awake model C anesthetized model D probability density... Spike Count FFTI FFTI cum. prob. density stimulus direction (degrees) E F gain variance stimulus direction (degrees) 9 stimulus direction (degrees) Fig. A single parameter, gain variance, can account for the observed changes in Fano factor tuning with behavioral state. (A-D) Model fitting captures differences in Fano factor tuning through changes in gain variance. The model is fit with the optimal α and var(g) for each population. (A) The distribution of FFTI for awake (blue, throughout) and anesthetized (orange, throughout) in the observed population (dashed trace) and the model values (solid trace). (B,C) The variance model predicts direction tuning for the awake (B) and anesthetized (C) experiments. Compare model results to observed s in Fig C-D. Parameters are fit to match distribution of FFTI, but reasonably reproduce tuning as well. (D) Distribution of FFTI for alternate model fitting in which the spike count used for each neuron was replaced by the average tuning curve over all neurons recorded in each experimental condition. The gain variance was fit separately to match the observed FFTI distribution. The α parameter is identical for all neurons and was chosen to minimize mean-squared error in the fit. Inset: Mean tuning curves for the awake and anesthetized experiments. (E) The distribution of gain variance parameters from the model fit in (D), showing larger gain variance values in the awake experiment. (F) Sample tunings generated by the best-fit models shown in (A-C). /7

17 spike count spike count A B C stimulus direction stimulus direction discriminability threshold discriminability threshold 9 8 FFTI > FFTI = FFTI < 7 T= ms number of cells (N) D E F FFTI > 9 FFTI = 8 FFTI < 7 T= ms 8 number of cells (N) discriminability threshold discriminability threshold FFTI > 9 FFTI = 8 FFTI < 7 N= T (ms) FFTI > FFTI = FFTI < N= T (ms) Fig. tuning and heterogeneity both contribute to lower discriminability thresholds in MT populations. (A) Model population with homogeneous tuning curves. (B) Cramer-Rao bound in homogeneous population models of different sizes with short-range correlations (c max =., c =.8). Tuning curves are fit to the average cumulative response up to ms after motion onset. Black traces show performance of models with varying stimulus-dependent variance. The solid trace is FFTI >, the dashed trace is FFTI =, and the dotted trace is FFTI <. The red line indicates the stimulus discriminability threshold for smooth pursuit behavior in macaques ms after pursuit initiation. (C) Same as in (B) but for populations of neurons with first-order statistics matched to average response at time t after motion onset. (D) Sample population of heterogeneous tuning curves drawn from measured tuning curves in recorded neurons. (E,F) Same as in (B,C) but for heterogeneous populations. The shaded areas show the standard deviation of the Cramer-Rao bound. 7/7