Neuronal variability and tuning are balanced to optimize naturalistic self-motion coding in primate vestibular pathways

Abstract
Introduction
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

It is commonly assumed that the brain’s neural coding strategies are adapted to the statistics of natural stimuli. Specifically, to maximize information transmission, a sensory neuron’s tuning function should effectively oppose the decaying stimulus spectral power, such that the neural response is temporally decorrelated (i.e. ‘whitened’). However, theory predicts that the structure of neuronal variability also plays an essential role in determining how coding is optimized. Here, we provide experimental evidence supporting this view by recording from neurons in early vestibular pathways during naturalistic self-motion. We found that central vestibular neurons displayed temporally whitened responses that could not be explained by their tuning alone. Rather, computational modeling and analysis revealed that neuronal variability and tuning were matched to effectively complement natural stimulus statistics, thereby achieving temporal decorrelation and optimizing information transmission. Taken together, our findings reveal a novel strategy by which neural variability contributes to optimized processing of naturalistic stimuli.

https://doi.org/10.7554/eLife.43019.001

Introduction

A fundamental challenge in neuroscience is to understand how the brain processes sensory information in order to generate accurate perception and guide appropriate behavioral responses. During everyday life, the inputs to our sensory systems have clearly defined statistical structures characterized by a pronounced decay in spectral power with increasing frequency (visual and auditory systems as reviewed in: Simoncelli and Olshausen, 2001; vestibular system: Carriot et al., 2014). It is generally agreed that sensory pathways are adapted to the statistical properties of natural stimuli (Laughlin, 1981; Attneave, 1954; Barlow, 1961; Simoncelli and Olshausen, 2001; Olshausen and Field, 2004; Wark et al., 2007). In this context, the brain is thought to optimize information transmission by removing the redundancy of (i.e. decorrelate) time-varying input, such that the spectral power of the resulting neuronal response is constant (i.e. are ‘whitened’; Srinivasan et al., 1982; Atick and Redlich, 1990; Atick, 2011). Indeed, the fact that whitening has been reported in several early sensory pathways supports this idea (Dan et al., 1996; Huang et al., 2016).

The prevailing view is that whitening at the single neuron level is achieved as a result of a precise match between tuning and stimulus statistics. Specifically, it has been proposed that increases in the neuronal tuning function effectively compensate for the decaying spectral power of natural stimuli (Dan et al., 1996; Wang et al., 2003; Pozzorini et al., 2013; Huang et al., 2016). Such high-pass tuning is thought to be mediated by intrinsic cellular dynamics that give rise to spike frequency adaptation (Benda and Herz, 2003; Wang et al., 2003; Pozzorini et al., 2013) and/or other mechanisms such as synaptic plasticity (Jackman and Regehr, 2017). However, in general, a match between the frequency spectra of neural tuning and stimulus statistics alone will not optimize information transmission. This is because the structure of the trial-to-trial variability of the neural response can strongly influence optimal coding as predicted from theory (van Hateren, 1992; Rieke et al., 1996; Tkacik et al., 2008; Tkacik et al., 2010). Further, in most brain areas, neural variability is not independent of frequency as it displays significant temporal correlations (Jaeger and Bower, 1994; Gershon et al., 1998; Manwani and Koch, 1999; Oram et al., 1999; Goldberg, 2000; Maimon and Assad, 2009; Shinomoto et al., 2009; Massot et al., 2011). To date, however, the effects of frequency-dependent structure of trial-to-trial variability on optimal coding at the single neuron level have for the most part been ignored (Field, 1987; Körding et al., 2004; Graham et al., 2006; Huang et al., 2016), not systematically investigated (van Hateren et al., 2002), or found to have minimal effect (Pitkow and Meister, 2012). Thus, the fundamental question of whether the frequency spectrum of neuronal variability actively contributes to optimizing information transmission of natural stimuli by sensory neurons has not yet been addressed experimentally.

The vestibular system benefits from well-described organization and easily described natural stimuli (e.g. head motion as a function of time; Cullen, 2011). This essential system generates reflexes that are vital for gaze and posture stabilization, as well as for accurate spatial perception and motor control (Cullen, 2012). Vestibular afferents displaying a wide range of resting discharge variabilities innervate the receptor cells of the vestibular sensors and synapse onto neurons within the central vestibular nuclei. These project to both motor centers and higher brain areas (e.g. thalamus; Goldberg, 2000; Cullen, 2012), thereby mediating vital reflexes and self-motion perception, respectively. Previous studies have emphasized the role of the resting discharge (i.e. in the absence of stimulation) on sensory processing across systems (for review, see: Ringach, 2009). Indeed, the variability of the resting discharge contributes to determining response variability during stimulation (Ratnam and Nelson, 2000; Chacron et al., 2003; Chacron et al., 2005; Sadeghi et al., 2007a). In particular, differences in resting discharge variability for afferents and central vestibular neurons strongly influence information transmission about self-motion, as assessed by using artificial (e.g. sinusoidal or noise) stimuli (Sadeghi et al., 2007a; Massot et al., 2011). However, how resting discharge variability contributes to central vestibular neural responses to natural self-motion remains an open question.

Here, we recorded the responses of both individual afferents and their post-synaptic central neuron targets within the vestibular nuclei to naturalistic self-motion stimuli. We found that early vestibular pathways demonstrate whitening at the first stage of central processing. Specifically, the response power spectra of central vestibular neurons are independent of frequency. Importantly, whitening is not simply inherited from afferents. Additionally, the tuning properties of central neurons alone do not predict their exceptionally whitened responses to natural self-motion. By using a computational model, we predicted that the frequency spectra of both tuning and variability will both significantly influence information transmission. Our experimental data confirmed modeling predictions of how neuronal variability affects information transmission. Notably, while increasing the level of variability will decrease information transmission, the frequency spectrum of neural variability for a given level will strongly determine optimality of coding (i.e. how close is the mutual information to its maximum value). Overall, we provide experimental evidence that the frequency spectrum of neural variability is matched with that of the tuning function in order to optimize coding during natural stimulation. Indeed, we found that coding was more optimized in central neurons than in afferents as the mutual information normalized by its maximum value was significantly higher. We hypothesize that these findings generalize across systems and species.

Results

We recorded from central vestibular-only (VO) neurons within the vestibular nuclei (VN) that are responsive to rotational head motion but not eye movements (n = 27). These neurons receive direct synaptic input from vestibular afferents, and project to higher brain areas mediating self-motion perception (Figure 1A; reviewed in: Cullen, 2012). We found that central neurons displayed large and highly variable resting discharge rates (51.1 ± 2.8 sp/s; CV*=0.44 ± 0.03) in the absence of stimulation, consistent with prior characterizations (e.g. Massot et al., 2011). For comparison, recordings were also made directly from individual semicircular canal afferents that were classified as either regular (n = 14; CV*=0.06 ± 0.005) or irregular (n = 12; CV*=0.33 ± 0.05) based on a previously established bimodal distribution of resting discharge variability corresponding to differences in axon diameter, response dynamics, and response sensitivity (Goldberg, 2000; Sadeghi et al., 2007a; Eatock et al., 2008; Massot et al., 2011; Massot et al., 2012). Notably, regular and irregular afferents displayed high resting discharge rates (regular: 102.4 ± 6.7 sp/s; irregular: 93.7 ± 8.9 sp/s), consistent with previous studies (e.g. Sadeghi et al., 2007a).

Figure 1 with 1 supplement see all

Download asset Open asset

Central vestibular neurons optimally encode natural self-motion stimuli through whitening.

(A) Vestibular afferents provide input to VO neurons in the vestibular nuclei, which in turn project to the spinal cord and higher order brain areas (i.e. thalamus/cortex). Recordings were made from VO neurons. (B) During experiments, the monkey was head-fixed and comfortably seated in a chair placed on a turntable. (C) Naturalistic head velocity during stimulation. (D) Response of an example VO neuron to the head velocity stimulus shown in C. (E) Normalized power spectra of natural stimuli and neural responses for an example VO neuron. Lower left insets shows autocorrelation functions (solid lines) and exponential fits (dashed lines) used to calculate correlation times (lower left), as well as correlation times and whitening indices (upper right). (F) Normalized power spectra of natural stimuli and population averages for VO neurons. Insets show autocorrelation functions (lower left), as well as correlation times and whitening indices (upper right). Whitening indices from VO neurons were significantly higher than those computed from the stimulus (Student’s t-test, p < 0.001, t(26) = 10.9), while their correlation times were significantly lower (Student’s t-test, p < 0.001, t(26) = 16.7). Gray bands show 95% confidence interval obtained using Poisson processes with the same firing rate as the experimental data for which the power spectrum is independent of frequency by definition (see Materials and methods). Error bars show ±1 SEM.

https://doi.org/10.7554/eLife.43019.002

Central vestibular neurons display temporal whitening in response to naturalistic self-motion

To test whether neurons within early vestibular pathways optimally encode natural self-motion statistics, we recorded their responses during naturalistic vestibular stimulation (Figure 1B) whose timecourse closely mimicked that of natural self-motion signals (see Materials and methods and Figure 1—figure supplement 1A). Because the power spectra of these signals decayed by more than two orders of magnitude over the frequency range 0–20 Hz (Figure 1—figure supplement 1A), we restricted our analysis to this frequency range. We found that central vestibular neurons showed robust responses to naturalistic self-motion stimuli (Figure 1D, bottom panels).

If central vestibular neurons optimally encode natural stimuli, then their response power spectra should be constant as a function of temporal frequency (i.e. white). The response power of an example central vestibular neuron is shown together with that of the stimulus in Figure 1E (green and black curves, respectively). Indeed, the response power spectrum was constant for frequencies within the physiologically relevant range (0–20 Hz; Huterer and Cullen, 2002; Carriot et al., 2017). In contrast, the stimulus power spectrum strongly decayed with increasing frequency (Figure 1E, black curve). Results were quantified by computing a whitening index (see Materials and methods) that is maximum and equal to unity for a constant spectrum. We found that this example neuron displayed a higher whitening index than the stimulus (Figure 1E, left panel of top right inset). We also evaluated to what extent whitening in the frequency domain corresponded to decorrelation in the temporal domain. Specifically, if whitening occurs in the frequency domain (i.e. the power spectrum of the response is more independent of frequency than that of the stimulus), then the width of the response autocorrelation function must decrease. As expected, the response autocorrelation function correlation time (~28 ms) was much lower than that of the stimulus (~106 ms; Figure 1E, right panel of top right inset).

Overall, analysis of our entire central neuron data gave rise to qualitatively similar results: while the stimulus power spectrum decayed rapidly for frequencies greater than 2 Hz, the response power spectrum instead remained constant as a function of frequency (Figure 1F, compare black and green curves, respectively). The response whitening index values were significantly higher than those obtained from the stimulus (Figure 1F, left panel of top right inset; p < 0.001). Correspondingly, in the temporal domain, the response autocorrelation function decayed to zero faster than that of the stimulus (Figure 1F, bottom inset, compare green and black curves, respectively), as quantified by significantly lower correlations times (Figure 1F, right panel of top right inset; p < 0.001). Thus, taken together, our results show that central vestibular neurons optimally encode natural self-motion stimuli through whitening.

Peripheral afferents are not sufficient to account for whitening by central vestibular neurons

The simplest explanation for our results shown above is that afferent responses to natural self-motion stimuli are already temporally whitened and that this is simply transmitted to central vestibular neurons. To test this, we next recorded afferent responses to naturalistic self-motion (Figure 2A). We considered the responses of both regular and irregular afferents that each innervate the canal endorgans and serve as parallel channels that project to central vestibular neurons (Highstein et al., 1987; Boyle et al., 1992). We found that, while regular and irregular afferents showed robust responses to naturalistic self-motion stimuli (Figure 2B, left and right panels, respectively), neither optimally encode naturalistic self-motion stimuli as their power spectra were not independent of frequency (Figure 2C,D,E,F). First, Figure 2C shows the stimulus (black curve) and response (blue curve) power spectra for a typical regular afferent. The response power spectrum decayed as a function of frequency, largely overlapping that of the stimulus (Figure 2C). Correspondingly, in the temporal domain, there was also strong overlap between the response and stimulus autocorrelation functions (Figure 2C, lower inset). Accordingly, both the whitening index and correlation time values computed from the response and the stimulus were comparable (Figure 2C, top inset). Overall, qualitatively similar results were seen across our regular afferent dataset (Figure 2D). The population-averaged response power spectrum and autocorrelation function also largely overlapped with those of the stimulus (Figure 2D, compare blue and black curves in the main panel and in the inset). The population-averaged response whitening index value was slightly larger than that of the stimulus (Figure 2D, left panel of top right inset; p < 0.001) but still much smaller than that computed for central vestibular neurons (compare with Figure 1E, left panel of top right inset). Further, the population-averaged response correlation time value was not significantly different than that of the stimulus (Figure 2D, right panel of top right inset; p = 0.79). Thus, our results indicate that, unlike central vestibular neurons, regular afferents do not perform whitening of natural self-motion stimuli.

Figure 2

Download asset Open asset

Vestibular afferents do not perform whitening of natural self-motion.

(A) Same diagram as in Figure 1A, except that recordings were made from afferents. (B) Response of example regular (left) and irregular (right) afferents during natural stimulation. (C) Normalized power spectra of natural stimuli (black) and neural response (blue) for the example regular afferent. The lower left inset shows the autocorrelation function of the response (blue) and of the stimulus (black). Solid lines represent autocorrelation functions, while the dashed lines are exponentials fits. The upper right inset shows the whitening index (left) and correlation time (right). (D) Normalized population-averaged power spectrum for regular afferents (blue) with that of the stimulus (black). The lower left inset shows the population-averaged response autocorrelation function (blue). The upper right inset shows the population-averaged whitening index (left) and correlation time (right) values. The whitening index was not significantly different than that of the stimulus (Student’s t-test, p = 0.79, t(13)=0.27). The correlation time was significantly lower than that of the stimulus (Student’s t-test, p < 0.001, t(13)=2.9). (E) Normalized power spectra of natural stimuli (black) and neural response (red) for the example irregular afferent. The lower left inset shows the autocorrelation function of the response (red) and of the stimulus (black). Solid lines represent autocorrelation functions, while the dashed lines are exponentials fits. The upper right inset shows the whitening index (left) and correlations time (right). (F) Normalized population-averaged power spectrum for irregular afferents (red) with that of the stimulus (black). The lower left inset shows the population-averaged response autocorrelation function (red). The upper right inset shows the population-averaged whitening index (left) and correlation time (right) values. The whitening index was significantly higher than that of the stimulus (Student’s t-test, p < 0.001, t(11)=2.25). The correlation time was significantly lower than that of the stimulus (Student’s t-test, p < 0.001, t(11)=9.64). Gray bands show 95% confidence interval obtained using Poisson processes with the same firing rate as the experimental data for which the power spectrum is independent of frequency by definition (see Materials and methods). Error bars show ±1 SEM.

https://doi.org/10.7554/eLife.43019.004

Second, a comparable analysis of irregular afferent responses revealed qualitatively similar results. Figure 2E shows the stimulus (black curve) and response (blue curve) power spectra for a typical irregular afferent, while Figure 2F shows the population-averages. Overall, the response power spectra of irregular afferents were also not constant as a function of frequency (Figure 2E and F, red curves), although the power spectra of irregular afferents decayed more slowly than those of their regular counterparts (compare Figure 2D and F). Further, the response autocorrelation functions did not rapidly decay to zero (Figure 2E and F, red curves in lower left insets). Thus, our results indicate that, like regular afferents but unlike central vestibular neurons, irregular afferents did not perform whitening of natural self-motion stimuli, as their response power spectra were not constant as a function of frequency.

Temporal whitening in early vestibular pathways occurs in stages

Thus far, our results have shown that the response power spectra of central vestibular neurons but not of afferents are constant as a function of frequency. Therefore, whitening of natural self-motion by central vestibular neurons is not simply inherited from their afferent input. To systematically compare whitening by regular afferents, irregular afferents, and central vestibular neurons, we superimposed their power spectra in Figure 3A. Comparison of the population-averaged whitening indices (Figure 3B) revealed that significantly higher values for central vestibular neurons than for either class of afferents (p < 0.01; Figure 3B). However, we found that whitening index values for irregular afferents were significantly higher than those for their regular counterparts (p = 0.002; Figure 3B). Correspondingly, to facilitate comparison in the time domain, the population-averaged autocorrelation functions of regular afferents, irregular afferents, and central vestibular neurons are superimposed (Figure 3C). Consistent with our spectral analysis above, population-averaged correlation times for central vestibular neurons were lower than for either class of afferents (regular: p < 0.001; irregular: p = 0.004; Figure 3D), while correlation times for irregular afferents were significantly lower than those of regular afferents (p = 0.006; Figure 3D). Thus, we conclude that whitening occurs in stages: some decorrelation of naturalistic self-motion input is achieved by irregular afferents and further decorrelation of this input occurs as the level of central vestibular neurons.

Figure 3

Download asset Open asset

Temporal whitening of neuronal responses occurs sequentially throughout early vestibular pathways.

(A) Normalized power spectra of natural stimuli (black), and neuronal responses of regular (blue), irregular afferents (red), and VO neurons (green). (B) Population-averaged whitening index values for regular afferents (blue, left), irregular afferents (red, middle), and VO neurons (green, right). Values for VO neurons were significantly higher than those of either regular or irregular afferents, while those of irregular afferents were significantly higher than those of regular afferents (one-way ANOVA, p < 0.001, F(2,50)=24.26). (C) Normalized autocorrelation functions of natural stimuli (black), and neuronal responses of regular (blue), irregular afferents (red), and VO neurons (green). (D) Population-averaged correlation time values for regular afferents (blue, left), irregular afferents (red, middle), and VO neurons (green, right). Values for VO neurons were significantly lower than those of either regular or irregular afferents, while those of irregular afferents were significantly lower than those of regular afferents (one-way ANOVA, p < 0.001, F(2,50)=26.68).

https://doi.org/10.7554/eLife.43019.005

The tuning properties of central vestibular neurons are not sufficient to explain whitening of naturalistic self-motion stimuli

As described above, the prevailing view is that sensory pathways achieve whitening by matching neural tuning to stimulus statistics. Specifically, increases in the neuronal tuning function should effectively compensate for the decaying stimulus power spectrum, such that the response power spectrum is constant as a function of frequency (Figure 4A). To test whether whitening by central vestibular neurons can be explained based on their tuning properties, we first characterized their sensitivity to naturalistic self-motion as a function of frequency (Figure 4B, inset). We then used these tuning curves (i.e. sensitivity as a function of frequency) to predict the response power spectrum to naturalistic self-motion (see Figure 4A and Materials and methods). Comparison of the population-averaged predicted and actual response power spectra revealed a poor match (Figure 4B, compare dashed green and solid green curves). Quantification of our results showed that predicted whitening index values were consistently lower than actual values (p < 0.001; Figure 4C). Further, whitening by central vestibular neurons could not be accounted for by static nonlinearities (Figure 4—figure supplement 1 and Figure 4—figure supplement 2; see Materials and methods), or by their tuning to artificial sinusoidal stimulation at discrete frequencies (Figure 5, see Materials and methods). In contrast, the tuning functions of both regular and irregular afferents were sufficient to predict their response power spectra to naturalistic self-motion (Figure 4—figure supplement 3), thereby showing that greater temporal whitening of naturalistic self-motion stimuli by irregular afferents is due to their greater high-pass tuning properties. Thus, we conclude that the tuning function of central vestibular neurons does not fully compensate for the decaying power spectrum of natural self-motion stimuli and thus cannot account for their whitened responses.

Figure 4 with 4 supplements see all

Download asset Open asset

Neural tuning to naturalistic stimulation does not account for the temporally whitened responses of central vestibular neurons to natural stimuli.

(A) Schematic showing how the response to a natural stimulus (right) is assumed to be determined from the stimulus spectrum (left) and the neural tuning function (middle). (B) Normalized population-averaged actual (solid green) and predicted (dashed green) response power spectra for central vestibular neurons. Inset: Population-averaged tuning curve showing gain as a function of frequency for central vestibular neurons. (C) Predicted versus actual whitening indices for central vestibular neurons. Most data points were below the identity line (dashed black).

https://doi.org/10.7554/eLife.43019.006

Figure 5

Download asset Open asset

Neural tuning to artificial sinusoidal stimulation does not account for responses of central vestibular neurons to natural stimuli.

(A₁) Response of example central vestibular neuron to sinusoidal yaw rotations with frequencies 0.5, 8 and 16 Hz. (A₂) Post-stimulus time histogram (PSTH) of this example central vestibular neuron to sinusoidal yaw rotations with frequencies 0.5, 8 and 16 Hz. (B) Schematic showing how the response to a natural stimulus (right) is assumed to be determined from the stimulus spectrum (left) and the neural tuning function (middle). (C) Gain as a function of frequency for an example central vestibular neuron with best fit (solid line). (D) Population-averaged gain as a function of frequency. The solid line is the best fit with bands showing ±1 SEM. (E) Normalized power spectra of neural responses (solid green) and predictions (dashed green) for central vestibular neurons. Inset: Population-averaged tuning curves obtained using naturalistic (gray) and sinusoidal (green) stimulation. Note that the tuning obtained using naturalistic stimulation was lower than that obtained using sinusoidal stimulation at low frequencies, consistent with the fact that central vestibular neurons display a boosting nonlinearity (Massot et al., 2012). (F) Predicted versus actual whitening indices.

https://doi.org/10.7554/eLife.43019.011

Effects of resting discharge variability and tuning on optimized information transmission through whitening

Theory predicts that the response power spectrum to natural stimuli is not solely determined by the tuning function (Shannon, 1948; Rieke et al., 1996). Specifically, the response power spectrum is equal to the sum of the noise power spectrum (i.e. the power spectrum of the neural variability) and the power spectrum of the response that would be predicted by tuning alone (i.e. transmitted power; Figure 4—figure supplement 4A). Thus, in order to maximize information through whitening, it is the response power spectrum that should be constant as a function of frequency, rather than the transmitted power spectrum.

To better understand how the noise and the transmitted power spectra interact in order to influence information transmission, we simulated a simple model (see Materials and methods). We initially considered a given noise spectrum (Figure 4—figure supplement 4B, black curve) and systematically varied the rate at which the transmitted power spectrum increased while keeping the area under the curve constant (Figure 4—figure supplement 4B, green curves). When the increase in the transmitted power fully compensated for the decrease in noise power (Figure 4—figure supplement 4B, condition b, solid green), the response power was constant as a function of frequency (Figure 4—figure supplement 4C, condition b, solid light green) and the mutual information was maximum (Figure 4—figure supplement 4D, point b). In contrast, when the increase in the transmitted power either undercompensated (Fig. Figure 4—figure supplement 4B, condition a, thin dashed green) or overcompensated (Fig. Figure 4—figure supplement 4B, condition c, thick dashed green) the decrease in the noise power, the response power was not constant as a function of frequency and the mutual information was less than the maximum value (Fig. Figure 4—figure supplement 4D, points a and c, respectively). Thus, our results show that, for a given noise intensity and spectrum, varying the transmitted power spectrum can have significant influence on mutual information.

We then considered the situation for which the transmitted power spectrum is fixed and increased the noise intensity (Figure 4—figure supplement 4E). We found that, with increasing noise intensity, noise as well as response power increased uniformly for all frequencies (Figure 4—figure supplement 4E,F), but that mutual information decreased (Figure 4—figure supplement 4G). This is expected given that mutual information is related to the signal-to-noise ratio, which decreases with increasing noise intensity. Figure 4—figure supplement 4H illustrates the effects of systematically varying the rate at which the transmitted power increases as well as the noise intensity. Note that the conditions illustrated in Figure 4—figure supplement 4B–D are denoted as a horizontal white line, while those illustrated in Figure 4—figure supplement 4E–G are denoted as a vertical white line, respectively. Thus, taken together, our modeling results show how both the noise and the transmitted power spectra interact in order to affect information transmission. Specifically, to maximize information, the spectral frequency content of the neural variability should perfectly compensate the transmitted power, such that the sum of the two is constant as a function of frequency.

Whitening of natural self-motion stimuli by central vestibular neurons can be explained by taking into account both their tuning properties and their resting discharge variability

We next tested our modeling predictions in the vestibular system. Theory predicts that the response power should then be equal to the sum of the noise power and of the transmitted power spectra (Figure 6A; Risken, 1996). Since central vestibular neurons display a substantial and highly variable resting discharge (e.g. Massot et al., 2011), we first used the power spectrum of the resting discharge in order to estimate that of trial-to-trial variability during stimulation. We found that the resting discharge power spectrum increases at higher frequencies (Figure 6B, inset, black curve) and thus, in conjunction with the neural tuning function (Figure 4B, inset), effectively compensated for decaying stimulus power. Indeed, the predicted and actual response power spectra were in excellent agreement (Figure 6B, compare green and light green curves). Further, predicted and actual whitening index values were not significantly different from one another (p = 0.39; Figure 6C). To confirm that the resting discharge is an accurate measurement of trial-to-trial variability during stimulation, we quantified the trial-to-trial variability of the neural response to naturalistic self-motion for a subset of neurons (see Materials and methods). We found that the power spectrum of the resting discharge and the residual were quite similar (Figure 6—figure supplement 1A), thus validating our assumption that the power spectrum of the resting discharge can be used to estimate that of the trial-to-trial variability during stimulation. Consequently, using either the resting discharge or trial-to-trial variability accurately predicted the temporally whitened neural responses to naturalistic stimuli (Figure 6—figure supplement 1B–C). We further found that the response power spectra during naturalistic self-motion were well predicted by summing the power spectrum of the trial-to-trial variability and the trial-averaged response spectrum (Figure 6—figure supplement 1).

Figure 6 with 2 supplements see all

Download asset Open asset

Neural variability and tuning determine whitening of natural self-motion stimuli by central vestibular neurons.

(A) Schematic showing how the response power spectrum (right) can be predicted from the resting activity power spectrum (left) as well as the stimulus power spectrum (middle left) and the neural tuning function (middle right). Specifically, the resting activity power spectrum should compensate the transmitted power spectrum (i.e. that obtained by multiplying the stimulus power by the tuning function) such that the response power spectrum (i.e., the sum of the two) is constant as a function of frequency. (B) Population-averaged actual response power spectrum (dark green) together with the prediction from the resting discharge and tuning (light green). The transmitted power spectrum (grey) is also shown. The inset shows the population-averaged resting discharge (black) and actual response (dark green) power spectra. (C) Predicted versus actual whitening index values for central vestibular neurons. Data points were scattered across the identity line. (D) *Left*: Population-averaged actual (left) and maximum (right) information rate values for central vestibular neurons. *Right:* The population-averaged mutual information normalized by the optimal value for VO neurons (green) was significantly greater than that of regular (blue) and irregular (red) afferents (one-way ANOVA, p < 0.001, F(2,50)=14.79). Error bars represent ±1 SEM.

https://doi.org/10.7554/eLife.43019.012

Finally, we quantified optimal coding by central vestibular neurons by comparing their rates of information transmission to optimal values (see Materials and methods). We found that actual information rates transmitted by central vestibular neurons were in good agreement with optimal values (Figure 6D, left panel), and comparison with afferents revealed that the actual information rate of central vestibular neurons was significantly closer to its maximum possible value (Figure 6D, right panel and Figure 6—figure supplement 2). Moreover, our modeling correctly predicted that, due to their higher variability, information rates computed from central vestibular neurons were lower than those of afferents (Figure 6—figure supplement 2). Thus, we conclude that the experimentally observed whitening by central vestibular neurons is achieved because both the power spectra of their variability and their tuning function effectively compensate for the decaying power of natural self-motion. Taken together, our findings show that the power spectrum of neural variability can have important consequences on determining whether sensory neurons optimally encode natural stimuli; thereby overturning the prevailing view that whitening is achieved by neural tuning alone.

Discussion

Summary of results

Here, we show that neurons at the first central stage of vestibular processing optimally encode natural self-motion stimuli through whitening due to a match between their frequency tuning and response variability. Specifically, we found that the spectral power of neuronal responses to naturalistic self-motion was constant as a function of temporal frequency. Given that we did not observe such whitening in their afferent input, we hypothesized that whitening occurs as a result of response properties, which are inherent to the central neurons. We thus first tested the prevailing view that whitening is a result of neural tuning, by establishing whether the frequency-dependent sensitivity of central neurons could account for the observed whitening. However, we found that the responses of vestibular neurons to naturalistic self-motion could not be explained based on their tuning alone. We hypothesized that a precise match between stimulus statistics, neural tuning, and response variability accounts for this unexpected result. To test this, we formulated a model that explicitly considered the dependence of both neuronal response variability and tuning on temporal frequency, and then demonstrated the generality of this model by predicting the spectra of neuronal responses to naturalistic stimulation. We then confirmed that this model could explain our neurophysiological findings, by showing that indeed response power predictions made using both resting activity and tuning were in excellent agreement with our experimental data. Taken together, these results both demonstrate that the responses of central vestibular neurons are optimal, and provide experimental evidence that, in early vestibular pathways, whitening (i.e. efficient coding) requires a precise match between the input distribution, neural tuning, and neural variability.

Impact of resting discharge variability on coding of self-motion in the vestibular system

The vestibular system is well-suited for understanding how neural variability affects the coding of information. Theory predicts that the trial-to-trial variability during stimulation can be well approximated by that of the resting activity in the absence of stimulation (Risken, 1996) and our results show that this is the case for central vestibular neurons under naturalistic stimulation. Further, it has long been known that the resting rates of central vestibular neurons, despite being lower than those of afferents (~50 sp/s vs. ~100 sp/s, respectively), are substantially more variable (Goldberg, 2000; Massot et al., 2011). Previous theoretical studies have emphasized the role of the resting discharge on information transmission (Chacron et al., 2004; Chacron et al., 2005). Notably, sensory stimulation must perturb the resting discharge in order to be detected. In particular, the power spectrum of the resting discharge can be considered as a ‘noise spectrum’ as it represents the amount of noise present at each frequency, and thus can approximate neural variability under stimulation (Ratnam and Nelson, 2000; Chacron et al., 2003; Chacron et al., 2005; Massot et al., 2011; Jamali et al., 2013). In the vestibular system, the fact that central vestibular neurons display higher levels of resting discharge than afferents thus provides an explanation as to why they display higher detection thresholds, and transmit less information about artificial self-motion stimuli than single afferents (Massot et al., 2011). Consequently, the increased resting discharge variability observed at higher stages of vestibular processing is detrimental to information transmission. This leads to the question of how central vestibular neurons encode self-motion stimuli in order to mediate self-motion perception during natural behaviors.

Our results obtained using naturalistic self-motion stimuli are consistent with prior studies showing that central vestibular neurons transmit less information than single afferents. Our modeling predicts that this is due to their greater levels of resting discharge variability. Importantly, our results show that the spectral content of neural variability strongly influences information transmission. Thus, in order to optimize neural coding of natural self-motion stimuli through whitening, such that response power is constant as a function of frequency, it is essential that the spectral frequency content of the resting discharge perfectly compensate the power transmitted from the stimulus (i.e. the transmitted power, Figure 4—figure supplement 4). This theoretical concept has been previously referred to as the ‘water filling analogy’ where the transmitted power is distributed to ‘fill in’ a vessel shaped as the variability power (de Ruyter van Steveninck and Laughlin, 1996; Rieke et al., 1996). Our results provide experimental evidence for this concept in the context of optimizing information transmission of natural sensory input. Thus, while the higher resting discharge variability of central vestibular neurons actually decreases information transmission relative to that of afferents (Figure 6—figure supplement 2), the spectral content of this variability is set such as to optimize the information being transmitted since the mutual information for central vestibular neurons is closer to its maximum value than for afferents (Figure 6D).

This leads to the question: what are the implications of increased neural variability in central vestibular pathways for naturalistic self-motion coding? To answer this question, it is important to consider that, to be useful to an organism, transmitted information must be decoded by downstream brain areas. We hypothesize that a match between stimulus statistics, neural tuning, and trial-to-trial variability further optimizes information transmission at the population level as proposed in other systems (Doi et al., 2012; Kastner et al., 2015). Further studies involving multi-unit recording from central vestibular neurons however will be needed to understand how optimized coding of natural self-motion by single central vestibular neurons observed here affects population coding as well as subsequent decoding of the population response by downstream neurons.

Role of variability in the coding of natural stimuli in other systems

Neural variability is seen ubiquitously in the central nervous system (for review see: Stein et al., 2005). Previous studies performed in multiple systems have shown that neural variability actually increases at the level of central neurons relative to that observed in more peripheral areas (visual: Kara et al., 2000; auditory: Wang et al., 2008; electrosensory: Maler, 2009), which is thought to be due to intense synaptic bombardment (Destexhe et al., 2001; Destexhe et al., 2003). However, the functional role of variability in neural coding continues to be the focus of much debate (Stein et al., 2005; McDonnell and Ward, 2011). Here, we showed that there is a match between the tuning properties and variability of central vestibular neurons which enables optimized encoding of natural self-motion through temporal whitening. This result constitutes an experimental demonstration that the spectral frequency content of neural variability is matched to that of the tuning function to optimize encoding of natural stimuli, as predicted by theoretical studies (Rieke et al., 1996). Indeed, most previous studies have focused on tuning and either assumed that neural variability was negligible or did not display temporal correlations (Dan et al., 1996; Wang et al., 2003; Pitkow and Meister, 2012). However, neurons displaying large and variable spontaneous activity (i.e. resting discharge activity) with temporal correlations are found ubiquitously within the sensory periphery and central brain regions (Hubel, 1959; Evarts, 1964; Pfeiffer and Kiang, 1965; Steriade et al., 1978; Jaeger and Bower, 1994; Köppl, 1997; Aizenman and Linden, 1999; Schneidman et al., 2006; Greschner et al., 2011). Accordingly, we speculate that our results showing that neural variability, together with tuning, plays a fundamental role towards optimizing information transmission about natural stimuli, will be generally applicable across sensory systems and species.

Future directions

In this study, we considered the coding of naturalistic self-motion stimuli that were passively applied along the yaw rotation axis. However, it is important to note that, under natural conditions, vestibular stimulation results from both active and passive self-motion. Previous studies have established that afferents respond similarly to both classes of stimuli (Cullen and Minor, 2002; Sadeghi et al., 2007b; Jamali et al., 2009). Thus, the results presented here for afferents are expected to be applicable to conditions where self-motion is actively generated. However, the central vestibular neurons which were the focus of the present study (i.e. VO neurons) display markedly attenuated responses to active self-motion consisting of head and/or body orienting movements (reviewed in: Cullen, 2012) due to integration of vestibular and extra-vestibular (e.g., proprioceptive and motor) signals. Thus, further studies will be needed to establish how these encode natural active self-motion signals (e.g. those experienced during locomotion). We further note that natural self-motion stimuli are generally not restricted to one axis of motion, but rather comprise of both three-dimensional rotational and translational components (Carriot et al., 2014; Carriot et al., 2017). Thus, an important direction for future studies will be to consider the effects of potential motor synergies in the encoding of self-motion by central vestibular neurons.

Materials and methods

Surgical procedures and data acquisition

Request a detailed protocol

All experimental protocols were approved by the McGill University Animal Care Committee and were in compliance with the guidelines of the Canadian Council on Animal Care. Three male macaque monkeys (2 Macaca mulatta and 1 Macaca fascicularis) were prepared for chronic extracellular recording using aseptic surgical techniques as previously described (Massot et al., 2011). Animals (aged 7, 8, and 8 years old) were housed in pairs on a 12 hr light/dark cycle.

Head movement recording. In the first stage of this study, head movements of freely moving animals were recorded as described previously (Carriot et al., 2017). Briefly, head movement recordings were made using a microelectromechanical systems (MEMS) module (iNEMO platform, STEVAL-MKI062V2; STMicroelectronics) that was fixed to the animal’s head. The MEMS module recorded linear acceleration in all three axes of motion (i.e. fore-aft, lateral, and vertical) as well as angular velocity for rotations along each of these. Each monkey was then released into a large expansive space (240 m³) where it was able to freely move (e.g. forage, groom, walk, run) and interact with another monkey from our colony. Enrichment materials (scattered treats, toys, etc.) as well as multiple structures that the monkeys could climb, were distributed throughout this space (for details see: Carriot et al., 2017).

Single-unit recordings. In the second stage of this study, we performed electrophysiology experiments in which we recorded the extracellular activity of horizontal semicircular canal afferents and central vestibular-only (VO) neurons. Head-restrained monkeys were seated in a primate chair that was mounted on a motion platform rotating about the vertical axis (i.e. yaw rotation). The motion platform produced rotations such that the head angular velocity stimulus recorded by a gyroscope mounted on the animal’s head matched those experienced during natural activities described above (Figure 1—figure supplement 1A, compare blue and red traces). We then applied rotational stimuli whose time course and spectral frequency content matched those of the recorded head movements (Figure 1B and C, and S1; see Materials and methods). We found that that the probability distributions of these signals were well-fit by a Gaussian (Figure 1—figure supplement 1B). We also recorded the extracellular activity of the same vestibular afferents and central neurons during sinusoidal rotation stimuli delivered at frequencies of 0.5–16 Hz with peak velocities of 40°/s. Data was collected through the Cerebus Neural Signal Processor (Blackrock Microsystems). Action potentials were discriminated from extracellular recordings offline (Offline Sorter, Plexon).

To confirm that each neuron in our sample discharged in a manner consistent with previous analyses, responses were characterized during voluntary eye movements and passive whole-body rotations. Monkeys were trained to track a small visual target (HeNe Laser) projected onto a white cylindrical screen located 60 cm away from the head for a juice reward, and eye position was measured using the magnetic search-coil technique. Both afferents and VO cells in our dataset were unresponsive to saccadic eye movements, smooth pursuit, and ocular fixation made to track the target. We note that, while horizontal canal afferents only respond to the component of rotational self-motion along the yaw axis (Fernandez and Goldberg, 1971), this is not the case for central vestibular neurons. Indeed, these neurons integrate inputs from multiple canal and otoliths in a manner that is both nonlinear and frequency dependent (Dickman and Angelaki, 2002; Carriot et al., 2015). Thus, we only considered yaw rotations here in order to directly compare the response properties of horizontal canal afferents to those of their central neuron targets.

Analysis of neuronal discharges

Request a detailed protocol

Data were imported into MATLAB (MathWorks) for analysis using custom-written algorithms (Mitchell et al., 2018). Head velocity signals were sampled at 1 kHz and digitally filtered at 125 Hz. Previous studies have established that vestibular afferents display large heterogeneities in terms of their resting discharge. In particular, the resting discharge regularity of vestibular afferents displays a bimodal distribution, which defines two afferent classes, regular and irregular, further characterized by differences in axon diameter as well as response dynamics (reviewed in: Goldberg, 2000; Eatock et al., 2008). Accordingly, we quantified the regularity of each afferent’s resting discharge by computing the normalized coefficient of variation (CV*) (Massot et al., 2011). We found a bimodal distribution of CV* values for our data set (Hartigan’s Dip Test; p = 0.002) and thus classified afferents with CV*<0.1 as regular and as irregular otherwise, as done previously (Goldberg et al., 1984). Spike trains were digitized at 1 kHz. Autocorrelation functions and power spectra were computed from digitized spike trains using the MATLAB functions xcorr and pwelch, respectively. Power spectra were normalized to their value at 2 Hz. The correlation time was measured by fitting an exponential to the autocorrelation function. The whitening index was computed as the integral of the spike train or stimulus power spectrum from 0 to 20 Hz divided by the difference between the minimum and maximum power multiplied by the frequency range (i.e. 20 Hz). The stimulus power spectrum was computed from the head velocity signal using the MATLAB function pwelch.

Response Dynamics: To better visualize the firing rate responses of neurons to sinusoidal stimuli at different frequencies, spike trains were convolved by a Kaiser window whose cut-off frequency was set to double that of the stimulus frequency (Cherif et al., 2008). Applying these filters does not affect the gain values computed for our dataset (Cherif et al., 2008). We verified that the activity of each central neuron recorded from was responsive to head but not eye motion, as noted above (Roy and Cullen, 2001). To compute the tuning to sinusoidal head motion, a least-squares regression analysis was used to describe each unit’s response to head rotations:

\hat{f r} (t) = b + G a i n H H V (t + θ)

where $\hat{f r} (t)$ is the filtered firing rate, b is a bias term set to the unit’s resting discharge rate, $H H V (t)$ is the time varying horizontal head velocity, Gain is the response sensitivity, and θ is the phase shift. For each sinusoidal frequency, the values of Gain and θ were determined by maximizing the variance-accounted-for as done previously (Cullen et al., 1996), these values are shown as data points in Figure 5C,D.

We next fit a transfer function of the following form to the obtained values of Gain and θ:

H (f) = \frac{A u T_{c} ({1 + u T}_{1})}{({1 + u T}_{c}) ({1 + u T}_{2})}

where u = 2 π i f, A is a constant, T_c and T₂ are the long and short time constants of the torsion-pendulum model of canal biomechanics and T₁ is proportional to the ratio of acceleration to velocity sensitivity of neuron responses, as done previously (Fernandez and Goldberg, 1971; Hullar et al., 2005; Massot et al., 2012). The estimated gain from the transfer function is then given by $|H (f)|$ while the estimated phase shift is given by $θ_{e s t} = {t a n}^{- 1} (i m a g [H (f)] / r e a l [H (f)])$ . Here real(…) and imag (…) denote the real and imaginary parts of H(f), respectively. We fixed T_c to 5.7 s consistent with the Torsion-pendulum model. For each neuron, values of A, T₁, and T₂ were estimated by minimizing the sum of residuals between the observed gain and phase values for sinusoidal stimulation and those computed from the estimated transfer function for those frequencies. Specifically, the sum of residuals was given by:

R e s = \sum_{f} {(G a i n (f) - |H (f)|)}^{2} + \frac{1}{100} \sum_{f} {(θ (f) - θ_{e s t} (f))}^{2}

where Gain and θ are obtained using equation (1). The factor of 1/100 for phase was used to make the residuals for gain and phase approximately equal to one another. The sum of residuals was minimized by using the fminsearch function in MATLAB. The estimated gain and phase from the transfer function are shown as solid lines in Figure 5C,D.

Response dynamics for naturalistic stimuli: The time-dependent firing rate was obtained by convolving the spike train with a Kaiser window whose cut-off frequency was set to 40 Hz, which will not affect the power spectrum within the frequency range 0–20 Hz. We found that the firing rate responses of central vestibular neurons were well-fit by a Gaussian (Figure 1—figure supplement 1B). For each neuron, we assumed that the time-dependent firing rate is given by:

{\hat{f r}}_{p r e d i c t e d} (t) = (h * s) (t) + f r_{0}

where ${\hat{f r}}_{p r e d i c t e d} (t)$ is the predicted firing rate, “*” denotes a convolution with a filter h(t), ${f r}_{0}$ is a constant set to the mean firing rate of the resting discharge, and s(t) is the stimulus. The Fourier transform of h(t), H₁(f) (i.e., the transfer function), was estimated using:

H_{1} (f) = \frac{P_{r s} (f)}{P_{s s} (f)}

where $P_{r s} (f)$ is the cross-spectrum between the firing rate and the stimulus, and $P_{s s} (f)$ is the stimulus power spectrum. We also estimated a static nonlinearity for each neuron by plotting the experimentally observed time-dependent firing rate as a function of ${\hat{f r}}_{e s t i m a t e d} (t)$ and used the following sigmoidal function to fit the data (Massot et al., 2012; Schneider et al., 2015):

T (x) = \frac{c_{1}}{1 + e^{- c_{2} (x - c_{3})}}

Adding this static nonlinearity to predict the firing rate response to naturalistic stimuli did not alter our results as shown in Figure 4—figure supplement 2. Moreover, replacing the transfer function H₁(f) by that obtained in response to sinusoidal stimulation (i.e. H(f) from eq. 2) also did not alter our results as shown in Figure 5C,D. We hypothesize that this is because central vestibular neurons display large resting discharge rates (~50 sp/s) that are highly variable, which effectively reduces the effects of nonlinearities (Stemmler, 1996).

We also calculated the transmitted power spectrum using:

P_{s i g n a l} = P_{s s} \times {|H (f)|}^{2}

This prediction was accurate for afferents as shown in Figure 4—figure supplement 3 but not for VO neurons as shown in Figure 4.

To quantify the amount of information transmitted by vestibular neurons over the behaviorally relevant frequency range (0–20 Hz), we calculated the mutual information rate using the following equation (Rieke et al., 1996):

M I = \int_{0}^{20} d f {l o g}_{2} [1 + S N R (f)]

Here, SNR(f) is the signal-to-noise ratio given by:

S N R (f) = \frac{P_{s i g n a l} (f)}{P_{0} (f)}

where $P_{0} (f)$ is the power spectrum of the resting discharge activity (i.e., in the absence of stimulation). The mutual information rate was then divided by the neural firing rate during stimulation. The normalized mutual information rate was obtained by dividing the actual mutual information by its maximum value obtained by systematically varying the neuronal tuning curve while keeping the resting discharge power spectrum constant for central vestibular neurons. Error estimates for mutual information quantities were determined from the distribution of values obtained from our dataset and are shown as whisker-box plots throughout.

The response power spectrum to the natural stimulus was predicted using linear response theory (Risken, 1996):

P_{R} (f) = P_{0} (f) + {|H (f)|}^{2} P_{s s} (f)

We note that linear response theory is generally applicable provided that the system is not near a critical point (e.g., a bifurcation) and that the stimulus perturbation is ‘weak’ (Risken, 1996). Previous studies have successfully used this theory to explain experimental data across multiple systems (Chacron et al., 2005; Huang et al., 2016).

To maximize information transmission (i.e., optimally encode) under the constraint of a constant noise variance (Shannon, 1948; Rieke et al., 1996), the following relationship must be satisfied:

P_{R} (f) = C

and C is a constant. We note that this assumes that the response, stimulus, and resting activity display Gaussian probability distributions. This is the case for the response and stimulus (Fig. 1 – figure supplement 1B) and we furthermore confirmed that probability distributions of the resting time-varying firing rate were well fit by Gaussian distributions for vestibular afferents (R = 0.99±0.005 and 0.98±0.01 for regular and irregular afferents, respectively), as well as for central vestibular neurons (R = 0.83±0.03). To test whether the experimentally obtained response power spectra were constant as a function of frequency, we simulated 1000 independent Poisson processes with the same number of spikes as the experimental data for each neuron. The power spectra obtained for each Poisson spike train were then used to compute the probability distribution of the power for each frequency. We found that the distributions were all Gaussian (Shapiro-Wilk test, p > 0.05 in all cases) as expected from the central limit theorem, and obtained a 95% confidence interval. These are plotted as gray bands in Figures 1 and 2.

Quantifying trial-to-trial variability to repeated stimulus presentations

Request a detailed protocol

For a subset of VO neurons (N = 8), the naturalistic stimulus was repeated multiple (i.e. 2–5) times, allowing us to quantify trial-to-trial variability. For each trial, the spiking response was convolved with a Kaiser window as described above in order to obtain the time-dependent firing rate. We then computed the residual as the difference between the firing rate response to a given trial and the mean response averaged over all trials. The power spectrum of this residual response was then computed as described above and compared to that of the resting discharge for each neuron. Our results, shown in Figure 6—figure supplement 1, show that both power spectra were quite similar. We found no significant differences in the estimated power spectrum of the trial-to-trial variability when using only a subset of trials for neurons for which the stimulus was repeated more than twice.

Modeling the effects of noise and transmitted power on information transmission: We used a simple model for which we assume that the total response power $P_{r e s p o n s e} (f)$ is given by:

P_{r e s p o n s e} (f) = P_{n o i s e} (f) + P_{t u n i n g, p r e d i c t e d} (f)

where $P_{n o i s e} (f)$ is the noise power spectrum and $P_{t u n i n g, p r e d i c t e d} (f)$ is the transmitted power spectrum. We focused on the 0–20 Hz frequency range, corresponding to that of natural self-motion stimuli as noted above. The mutual information rate was computed from (Rieke et al., 1996):

M I = \int_{0}^{20} {d f l o g}_{2} [1 + \frac{P_{t u n i n g, p r e d i c t e d} (f)}{P_{n o i s e} (f)}]

We then considered the following dependencies for the rest and transmitted power spectra:

P_{n o i s e} (f) = P_{0} + \frac{2 - \frac{f}{20}}{\int_{0}^{20} d f (2 - \frac{f}{20})}

P_{t u n i n g, p r e d i c t e d} (f) = \frac{1 + s e n s i t i v i t y \frac{f}{20}}{\int_{0}^{20} d f (1 + s e n s i t i v i t y \frac{f}{20})}

where P₀ is the background noise power which is completely determined by the noise intensity and sensitivity controls the rate at which the transmitted power increases with frequency, while keeping the area under the curve constant and equal to unity. Both P₀ and sensitivity were varied systematically and their effects on information computed.

Statistics

Our sample sizes were similar to those generally employed in the field (Massot et al., 2011). Before statistical analysis, normality of distribution was evaluated using a Shapiro-Wilk’s test. All data were tested for presence of non-stationarities using an augmented Dickey-Fuller test. We did not find any significant non-stationarities during either resting discharge or driven activities for either of afferents or central vestibular neurons (p > 0.05 in all cases). To determine if variances between groups were comparable, an F-test was used. Statistical significance (p < 0.05) was determined using parametric analysis with either two-tailed Student’s t-test or one-way ANOVA. Post-hoc pairwise comparisons were conducted using Tukey’s honestly significant difference test (Figure 3) or Dunnett’s test (Figure 6). Throughout, values are expressed as mean ±SEM.

Data availability

All data have been deposited on Figshare under the URL https://doi.org/10.6084/m9.figshare.7423724.v1.

The following data sets were generated

1. Mitchell D
2. Kwan A
3. Carriot J
4. Chacron M
5. Cullen KE
(2018) Figshare
Figure source data for "Neuronal variability and tuning are balanced to optimize coding of naturalistic self-motion in primate vestibular pathways".

https://doi.org/10.6084/m9.figshare.7423724.v1

References

1. Aizenman CD
2. Linden DJ
(1999) Regulation of the rebound depolarization and spontaneous firing patterns of deep nuclear neurons in slices of rat cerebellum
Journal of Neurophysiology 82:1697–1709.

https://doi.org/10.1152/jn.1999.82.4.1697
- PubMed
- Google Scholar
1. Atick JJ
2. Redlich AN
(1990) Towards a Theory of Early Visual Processing
Neural Computation 2:308–320.

https://doi.org/10.1162/neco.1990.2.3.308
- Google Scholar
1. Atick JJ
(2011) Could information theory provide an ecological theory of sensory processing?
Network: Computation in Neural Systems 22:4–44.

https://doi.org/10.3109/0954898X.2011.638888
- Google Scholar
1. Attneave F
(1954) Some informational aspects of visual perception
Psychological Review 61:183–193.

https://doi.org/10.1037/h0054663
- PubMed
- Google Scholar
1. Barlow HB
(1961) Possible principles underlying the transformations of sensory messages
Mitpress pp. 217–234.

https://doi.org/10.7551/mitpress/9780262518420.003.0013
- Google Scholar
1. Benda J
2. Herz AV
(2003) A universal model for spike-frequency adaptation
Neural Computation 15:2523–2564.

https://doi.org/10.1162/089976603322385063
- PubMed
- Google Scholar
(1992) Inputs from regularly and irregularly discharging vestibular nerve afferents to secondary neurons in squirrel monkey vestibular nuclei. III. Correlation with vestibulospinal and vestibuloocular output pathways
Journal of Neurophysiology 68:471–484.

https://doi.org/10.1152/jn.1992.68.2.471
- PubMed
- Google Scholar
(2014) Statistics of the vestibular input experienced during natural self-motion: implications for neural processing
Journal of Neuroscience 34:8347–8357.

https://doi.org/10.1523/JNEUROSCI.0692-14.2014
- PubMed
- Google Scholar
(2015) Integration of canal and otolith inputs by central vestibular neurons is subadditive for both active and passive self-motion: implication for perception
Journal of Neuroscience 35:3555–3565.

https://doi.org/10.1523/JNEUROSCI.3540-14.2015
- PubMed
- Google Scholar
(2017) The statistics of the vestibular input experienced during natural self-motion differ between rodents and primates
The Journal of Physiology 595:2751–2766.

https://doi.org/10.1113/JP273734
- PubMed
- Google Scholar
(2003) The effects of spontaneous activity, background noise, and the stimulus ensemble on information transfer in neurons
Network: Computation in Neural Systems 14:803–824.

https://doi.org/10.1088/0954-898X_14_4_010
- Google Scholar
(2004) Noise shaping by interval correlations increases information transfer
Physical Review Letters 92:080601.

https://doi.org/10.1103/PhysRevLett.92.080601
- PubMed
- Google Scholar
(2005) Electroreceptor neuron dynamics shape information transmission
Nature Neuroscience 8:673–678.

https://doi.org/10.1038/nn1433
- PubMed
- Google Scholar
(2008) An improved method for the estimation of firing rate dynamics using an optimal digital filter
Journal of Neuroscience Methods 173:165–181.

https://doi.org/10.1016/j.jneumeth.2008.05.021
- PubMed
- Google Scholar
1. Cullen KE
2. Rey CG
3. Guitton D
4. Galiana HL
(1996) The use of system identification techniques in the analysis of oculomotor burst neuron spike train dynamics
Journal of Computational Neuroscience 3:347–368.

https://doi.org/10.1007/BF00161093
- Google Scholar
1. Cullen KE
2. Minor LB
(2002) Semicircular canal afferents similarly encode active and passive head-on-body rotations: implications for the role of vestibular efference
The Journal of Neuroscience 22:RC226.

https://doi.org/10.1523/JNEUROSCI.22-11-j0002.2002
- PubMed
- Google Scholar
1. Cullen KE
(2011) The neural encoding of self-motion
Current Opinion in Neurobiology 21:587–595.

https://doi.org/10.1016/j.conb.2011.05.022
- PubMed
- Google Scholar
1. Cullen KE
(2012) The vestibular system: multimodal integration and encoding of self-motion for motor control
Trends in Neurosciences 35:185–196.

https://doi.org/10.1016/j.tins.2011.12.001
- PubMed
- Google Scholar
1. Dan Y
2. Atick JJ
3. Reid RC
(1996) Efficient coding of natural scenes in the lateral geniculate nucleus: experimental test of a computational theory
The Journal of Neuroscience 16:3351–3362.

https://doi.org/10.1523/JNEUROSCI.16-10-03351.1996
- PubMed
- Google Scholar
1. de Ruyter van Steveninck RR
2. Laughlin SB
(1996) The rate of information transfer at graded-potential synapses
Nature 379:642–645.

https://doi.org/10.1038/379642a0
- Google Scholar
(2001) Fluctuating synaptic conductances recreate in vivo-like activity in neocortical neurons
Neuroscience 107:13–24.

https://doi.org/10.1016/S0306-4522(01)00344-X
- PubMed
- Google Scholar
(2003) The high-conductance state of neocortical neurons in vivo
Nature Reviews Neuroscience 4:739–751.

https://doi.org/10.1038/nrn1198
- PubMed
- Google Scholar
1. Dickman JD
2. Angelaki DE
(2002) Vestibular convergence patterns in vestibular nuclei neurons of alert primates
Journal of Neurophysiology 88:3518–3533.

https://doi.org/10.1152/jn.00518.2002
- PubMed
- Google Scholar
1. Doi E
2. Gauthier JL
3. Field GD
4. Shlens J
5. Sher A
6. Greschner M
7. Machado TA
8. Jepson LH
9. Mathieson K
10. Gunning DE
11. Litke AM
12. Paninski L
13. Chichilnisky EJ
14. Simoncelli EP
(2012) Efficient coding of spatial information in the primate retina
Journal of Neuroscience 32:16256–16264.

https://doi.org/10.1523/JNEUROSCI.4036-12.2012
- PubMed
- Google Scholar
(2008) Ion channels in mammalian vestibular afferents may set regularity of firing
Journal of Experimental Biology 211:1764–1774.

https://doi.org/10.1242/jeb.017350
- PubMed
- Google Scholar
1. Evarts E
(1964) Temporal patterns of discharge of pyramidal tract neurons during sleep and waking in the monkey
Journal of Neurophysiology 27:152–171.

https://doi.org/10.1152/jn.1964.27.2.152
- PubMed
- Google Scholar
1. Fernandez C
2. Goldberg JM
(1971) Physiology of peripheral neurons innervating semicircular canals of the squirrel monkey. II. Response to sinusoidal stimulation and dynamics of peripheral vestibular system
Journal of Neurophysiology 34:661–675.

https://doi.org/10.1152/jn.1971.34.4.661
- PubMed
- Google Scholar
1. Field DJ
(1987) Relations between the statistics of natural images and the response properties of cortical cells
Journal of the Optical Society of America A 4:2379–2394.

https://doi.org/10.1364/JOSAA.4.002379
- Google Scholar
(1998) Coding strategies in monkey V1 and inferior temporal cortices
Journal of Neurophysiology 79:1135–1144.

https://doi.org/10.1152/jn.1998.79.3.1135
- PubMed
- Google Scholar
(1984) Relation between discharge regularity and responses to externally applied galvanic currents in vestibular nerve afferents of the squirrel monkey
Journal of Neurophysiology 51:1236–1256.

https://doi.org/10.1152/jn.1984.51.6.1236
- PubMed
- Google Scholar
1. Goldberg JM
(2000) Afferent diversity and the organization of central vestibular pathways
Experimental Brain Research 130:277–297.

https://doi.org/10.1007/s002210050033
- PubMed
- Google Scholar
(2006) Can the theory of “whitening” explain the center-surround properties of retinal ganglion cell receptive fields?
Vision Research 46:2901–2913.

https://doi.org/10.1016/j.visres.2006.03.008
- Google Scholar
1. Greschner M
2. Shlens J
3. Bakolitsa C
4. Field GD
5. Gauthier JL
6. Jepson LH
7. Sher A
8. Litke AM
9. Chichilnisky EJ
(2011) Correlated firing among major ganglion cell types in primate retina
The Journal of Physiology 589:75–86.

https://doi.org/10.1113/jphysiol.2010.193888
- PubMed
- Google Scholar
(1987) Inputs from regularly and irregularly discharging vestibular nerve afferents to secondary neurons in the vestibular nuclei of the squirrel monkey. II. Correlation with output pathways of secondary neurons
Journal of Neurophysiology 58:719–738.

https://doi.org/10.1152/jn.1987.58.4.719
- PubMed
- Google Scholar
(2016) Temporal decorrelation by SK channels enables efficient neural coding and perception of natural stimuli
Nature Communications 7:11353.

https://doi.org/10.1038/ncomms11353
- PubMed
- Google Scholar
1. Hubel DH
(1959) Single unit activity in striate cortex of unrestrained cats
The Journal of Physiology 147:226–238.

https://doi.org/10.1113/jphysiol.1959.sp006238
- PubMed
- Google Scholar
(2005) Responses of irregularly discharging chinchilla semicircular canal vestibular-nerve afferents during high-frequency head rotations
Journal of Neurophysiology 93:2777–2786.

https://doi.org/10.1152/jn.01002.2004
- PubMed
- Google Scholar
1. Huterer M
2. Cullen KE
(2002) Vestibuloocular reflex dynamics during high-frequency and high-acceleration rotations of the head on body in rhesus monkey
Journal of Neurophysiology 88:13–28.

https://doi.org/10.1152/jn.2002.88.1.13
- PubMed
- Google Scholar
1. Jackman SL
2. Regehr WG
(2017) The Mechanisms and Functions of Synaptic Facilitation
Neuron 94:447–464.

https://doi.org/10.1016/j.neuron.2017.02.047
- PubMed
- Google Scholar
1. Jaeger D
2. Bower JM
(1994) Prolonged responses in rat cerebellar Purkinje cells following activation of the granule cell layer: an intracellular in vitro and in vivo investigation
Experimental Brain Research 100:200–214.

https://doi.org/10.1007/BF00227191
- PubMed
- Google Scholar
(2009) Response of vestibular nerve afferents innervating utricle and saccule during passive and active translations
Journal of Neurophysiology 101:141–149.

https://doi.org/10.1152/jn.91066.2008
- PubMed
- Google Scholar
(2013) Strong correlations between sensitivity and variability give rise to constant discrimination thresholds across the otolith afferent population
Journal of Neuroscience 33:11302–11313.

https://doi.org/10.1523/JNEUROSCI.0459-13.2013
- PubMed
- Google Scholar
(2000) Low response variability in simultaneously recorded retinal, thalamic, and cortical neurons
Neuron 27:635–646.

https://doi.org/10.1016/S0896-6273(00)00072-6
- PubMed
- Google Scholar
(2015) Critical and maximally informative encoding between neural populations in the retina
PNAS 112:2533–2538.

https://doi.org/10.1073/pnas.1418092112
- PubMed
- Google Scholar
1. Köppl C
(1997) Frequency tuning and spontaneous activity in the auditory nerve and cochlear nucleus magnocellularis of the barn owl Tyto alba
Journal of Neurophysiology 77:364–377.

https://doi.org/10.1152/jn.1997.77.1.364
- PubMed
- Google Scholar
(2004) How are complex cell properties adapted to the statistics of natural stimuli?
Journal of Neurophysiology 91:206–212.

https://doi.org/10.1152/jn.00149.2003
- PubMed
- Google Scholar
1. Laughlin S
(1981)
A simple coding procedure enhances a neuron’s information capacity

Zeitschrift für Naturforschung Sect C Biosci 36:910–912.
- Google Scholar
1. Maimon G
2. Assad JA
(2009) Beyond Poisson: increased spike-time regularity across primate parietal cortex
Neuron 62:426–440.

https://doi.org/10.1016/j.neuron.2009.03.021
- PubMed
- Google Scholar
1. Maler L
(2009) Receptive field organization across multiple electrosensory maps. II. Computational analysis of the effects of receptive field size on prey localization
The Journal of Comparative Neurology 516:394–422.

https://doi.org/10.1002/cne.22120
- PubMed
- Google Scholar
1. Manwani A
2. Koch C
(1999) Detecting and estimating signals in noisy cable structure, I: neuronal noise sources
Neural Computation 11:1797–1829.

https://doi.org/10.1162/089976699300015972
- PubMed
- Google Scholar
(2011) Information transmission and detection thresholds in the vestibular nuclei: single neurons vs. population encoding
Journal of Neurophysiology 105:1798–1814.

https://doi.org/10.1152/jn.00910.2010
- PubMed
- Google Scholar
(2012) The vestibular system implements a linear-nonlinear transformation in order to encode self-motion
PLoS Biology 10:e1001365.

https://doi.org/10.1371/journal.pbio.1001365
- PubMed
- Google Scholar
1. McDonnell MD
2. Ward LM
(2011) The benefits of noise in neural systems: bridging theory and experiment
Nature Reviews Neuroscience 12:415–426.

https://doi.org/10.1038/nrn3061
- PubMed
- Google Scholar
1. Mitchell DE
2. Kwan A
3. Carriot J
4. Chacron MJ
5. Cullen KE
(2018) Figure source data for neuronal variability and tuning are balanced to optimize coding of naturalistic self-motion in primate vestibular pathways
eLife e43019.

https://doi.org/10.7554/eLife.43019
- Google Scholar
1. Olshausen BA
2. Field DJ
(2004) Sparse coding of sensory inputs
Current Opinion in Neurobiology 14:481–487.

https://doi.org/10.1016/j.conb.2004.07.007
- PubMed
- Google Scholar
(1999) Stochastic nature of precisely timed spike patterns in visual system neuronal responses
Journal of Neurophysiology 81:3021–3033.

https://doi.org/10.1152/jn.1999.81.6.3021
- PubMed
- Google Scholar
1. Pfeiffer RR
2. Kiang NY
(1965) Spike Discharge Patterns of Spontaneous and Continuously Stimulated Activity in the Cochlear Nucleus of Anesthetized Cats
Biophysical Journal 5:301–316.

https://doi.org/10.1016/S0006-3495(65)86718-2
- PubMed
- Google Scholar
1. Pitkow X
2. Meister M
(2012) Decorrelation and efficient coding by retinal ganglion cells
Nature Neuroscience 15:628–635.

https://doi.org/10.1038/nn.3064
- PubMed
- Google Scholar
1. Pozzorini C
2. Naud R
3. Mensi S
4. Gerstner W
(2013) Temporal whitening by power-law adaptation in neocortical neurons
Nature Neuroscience 16:942–948.

https://doi.org/10.1038/nn.3431
- PubMed
- Google Scholar
1. Ratnam R
2. Nelson ME
(2000) Nonrenewal statistics of electrosensory afferent spike trains: implications for the detection of weak sensory signals
The Journal of Neuroscience 20:6672–6683.

https://doi.org/10.1523/JNEUROSCI.20-17-06672.2000
- PubMed
- Google Scholar
Book
(1996)
Spikes: Exploring the Neural Code

MIT Press.
- Google Scholar
1. Ringach DL
(2009) Spontaneous and driven cortical activity: implications for computation
Current Opinion in Neurobiology 19:439–444.

https://doi.org/10.1016/j.conb.2009.07.005
- Google Scholar
Book
1. Risken H
(1996)
The Fokker-Planck Equation

Berlin, Heidelberg, Germany: Springer.
- Google Scholar
1. Roy JE
2. Cullen KE
(2001) Selective processing of vestibular reafference during self-generated head motion
The Journal of Neuroscience 21:2131–2142.

https://doi.org/10.1523/JNEUROSCI.21-06-02131.2001
- PubMed
- Google Scholar
(2007a) Neural variability, detection thresholds, and information transmission in the vestibular system
Journal of Neuroscience 27:771–781.

https://doi.org/10.1523/JNEUROSCI.4690-06.2007
- Google Scholar
(2007b) Response of vestibular-nerve afferents to active and passive rotations under normal conditions and after unilateral labyrinthectomy
Journal of Neurophysiology 97:1503–1514.

https://doi.org/10.1152/jn.00829.2006
- PubMed
- Google Scholar
(2015) The increased sensitivity of irregular peripheral canal and otolith vestibular afferents optimizes their encoding of natural stimuli
Journal of Neuroscience 35:5522–5536.

https://doi.org/10.1523/JNEUROSCI.3841-14.2015
- PubMed
- Google Scholar
(2006) Weak pairwise correlations imply strongly correlated network states in a neural population
Nature 440:1007–1012.

https://doi.org/10.1038/nature04701
- PubMed
- Google Scholar
1. Shannon CE
(1948) A Mathematical Theory of Communication
Bell System Technical Journal 27:379–423.

https://doi.org/10.1002/j.1538-7305.1948.tb01338.x
- Google Scholar
1. Shinomoto S
2. Kim H
3. Shimokawa T
4. Matsuno N
5. Funahashi S
6. Shima K
7. Fujita I
8. Tamura H
9. Doi T
10. Kawano K
11. Inaba N
12. Fukushima K
13. Kurkin S
14. Kurata K
15. Taira M
16. Tsutsui K
17. Komatsu H
18. Ogawa T
19. Koida K
20. Tanji J
21. Toyama K
(2009) Relating neuronal firing patterns to functional differentiation of cerebral cortex
PLoS Computational Biology 5:e1000433.

https://doi.org/10.1371/journal.pcbi.1000433
- PubMed
- Google Scholar
1. Simoncelli EP
2. Olshausen BA
(2001) Natural image statistics and neural representation
Annual Review of Neuroscience 24:1193–1216.

https://doi.org/10.1146/annurev.neuro.24.1.1193
- PubMed
- Google Scholar
(1982) Predictive Coding: A Fresh View of Inhibition in the Retina
Proceedings of the Royal Society B: Biological Sciences 216:427–459.

https://doi.org/10.1098/rspb.1982.0085
- Google Scholar
(2005) Neuronal variability: noise or part of the signal?
Nature Reviews Neuroscience 6:389–397.

https://doi.org/10.1038/nrn1668
- PubMed
- Google Scholar
1. Stemmler M
(1996) A single spike suffices: the simplest form of stochastic resonance in model neurons
Network: Computation in Neural Systems 7:687–716.

https://doi.org/10.1088/0954-898X_7_4_005
- Google Scholar
(1978) Firing rates and patterns of output and nonoutput cells in cortical areas 5 and 7 of cat during the sleep-waking cycle
Experimental Neurology 60:443–468.

https://doi.org/10.1016/0014-4886(78)90003-1
- PubMed
- Google Scholar
(2008) Information flow and optimization in transcriptional regulation
PNAS 105:12265–12270.

https://doi.org/10.1073/pnas.0806077105
- PubMed
- Google Scholar
(2010) Optimal population coding by noisy spiking neurons
PNAS 107:14419–14424.

https://doi.org/10.1073/pnas.1004906107
- PubMed
- Google Scholar
1. van Hateren JH
(1992) A theory of maximizing sensory information
Biological Cybernetics 68:23–29.

https://doi.org/10.1007/BF00203134
- PubMed
- Google Scholar
(2002) Processing of natural temporal stimuli by macaque retinal ganglion cells
The Journal of Neuroscience 22:9945–9960.

https://doi.org/10.1523/JNEUROSCI.22-22-09945.2002
- PubMed
- Google Scholar
(2003) Adaptation and temporal decorrelation by single neurons in the primary visual cortex
Journal of Neurophysiology 89:3279–3293.

https://doi.org/10.1152/jn.00242.2003
- PubMed
- Google Scholar
1. Wang X
2. Lu T
3. Bendor D
4. Bartlett E
(2008) Neural coding of temporal information in auditory thalamus and cortex
Neuroscience 154:294–303.

https://doi.org/10.1016/j.neuroscience.2008.03.065
- PubMed
- Google Scholar
(2007) Sensory adaptation
Current Opinion in Neurobiology 17:423–429.

https://doi.org/10.1016/j.conb.2007.07.001
- PubMed
- Google Scholar

Article and author information

Author details

Diana E Mitchell

Department of Physiology, McGill University, Montreal, Canada

Contribution
Data curation, Formal analysis, Methodology, Writing—original draft, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-0733-484X
Annie Kwan

Department of Physiology, McGill University, Montreal, Canada

Contribution
Data curation, Formal analysis

Competing interests
No competing interests declared
Jerome Carriot

Department of Physiology, McGill University, Montreal, Canada

Contribution
Data curation, Investigation, Methodology

Competing interests
No competing interests declared
Maurice J Chacron

Department of Physiology, McGill University, Montreal, Canada

Contribution
Conceptualization, Data curation, Formal analysis, Supervision, Methodology, Writing—original draft, Writing—review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-3032-452X
Kathleen E Cullen
1. Department of Physiology, McGill University, Montreal, Canada
2. Department of Biomedical Engineering, The Johns Hopkins University, Baltimore, United States
Present address
Department of Biomedical Engineering, The Johns Hopkins University, Baltimore, Maryland, United States

Contribution
Conceptualization, Data curation, Supervision, Funding acquisition, Methodology, Writing—original draft, Writing—review and editing

For correspondence
kathleen.cullen@jhu.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-9348-0933

Funding

Canadian Institutes of Health Research

Maurice J Chacron
Kathleen E Cullen

Canada Research Chairs

Maurice J Chacron

National Institutes of Health (DC2390)

Kathleen E Cullen

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This research study was supported by the Canadian Institutes of Health Research (MJC and KEC), Canada research chair program (MJC), and the National Institutes of Health (KEC).

Ethics

Animal experimentation: All experimental protocols were approved by the McGill University Animal Care Committee (#2001-4096) and were in compliance with the guidelines of the Canadian Council on Animal Care. Three male macaque monkeys (2 Macaca mulatta and 1 Macaca fascicularis) were prepared for chronic extracellular recording using aseptic surgical techniques as previously described (Massot et al., 2011). Animals (aged 7, 8, and 8 years old) were housed in pairs on a 12 hour light/dark cycle.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

1,490

views
230

downloads
29

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Diana E Mitchell
Annie Kwan
Jerome Carriot
Maurice J Chacron
Kathleen E Cullen

(2018)

Neuronal variability and tuning are balanced to optimize naturalistic self-motion coding in primate vestibular pathways

eLife 7:e43019.

https://doi.org/10.7554/eLife.43019

Categories and tags

Research organism

Rhesus macaque

Share this article

Cite this article

Central vestibular neurons optimally encode natural self-motion stimuli through whitening.

Vestibular afferents do not perform whitening of natural self-motion.

Temporal whitening of neuronal responses occurs sequentially throughout early vestibular pathways.

Neural tuning to naturalistic stimulation does not account for the temporally whitened responses of central vestibular neurons to natural stimuli.

Neural tuning to artificial sinusoidal stimulation does not account for responses of central vestibular neurons to natural stimuli.

Neural variability and tuning determine whitening of natural self-motion stimuli by central vestibular neurons.

Author details

Diana E Mitchell

Contribution

Competing interests

Annie Kwan

Contribution

Competing interests

Jerome Carriot

Contribution

Competing interests

Maurice J Chacron

Contribution

Competing interests

Kathleen E Cullen

Present address

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Research organism

Further reading