Abstract
Temporal limits on perceptual decisions set strict boundaries on the possible underlying neural computations. How odor information is encoded in the olfactory system is still poorly understood. Here, we sought to define the limit on the speed of olfactory processing. To achieve this, we trained mice to discriminate different odor concentrations in a novel behavioral setup with precise odor delivery synchronized to the sniffing cycle. Mice reported their choice by moving a horizontal treadmill with their front limbs. We found that mice reported discriminations of 75% accuracy in 70–90 ms after odor inhalation. For a low concentration and nontrigeminal odorant, this time was 90–140 ms, showing that mice process odor information rapidly even in the absence of trigeminal stimulation. These response times establish, after accounting for odor transduction and motor delays, that olfactory processing can take tens of milliseconds. This study puts a strong limit on the underlying neural computations and suggests that the action potentials forming the neural basis for these decisions are fired in a few tens of milliseconds.
SIGNIFICANCE STATEMENT Understanding how sensory information is processed requires different approaches that span multiple levels of investigation from genes to neurons to behavior. Limits on behavioral performance constrain the possible neural mechanisms responsible for specific computations. Using a novel behavioral paradigm, we established that mice can make decisions about odor intensity surprisingly fast. After accounting for sensory and motor delays, the limit on some olfactory neural computations can be as low as a few tens of milliseconds, which suggests that only the first action potentials across a population of neurons contribute to these computations.
Introduction
Recently, the mouse has become a very powerful model organism for studying neural computations because it allows the combination of a variety of approaches at different levels from molecular and genetic to computational and behavioral. Despite great progress in mouse neurobiology, little is known about the temporal limits on mouse behavioral performance.
Olfaction, being presumably the most genetically tractable and anatomically compact sensory system, presents a great opportunity for testing the basic principles of sensory processing and decision making. In olfaction, a single sniff is considered a snapshot of the olfactory world and a unit of perception (Uchida et al., 2006). In fact, rats only need to sample an odor with a single sniff, 125–300 ms, to achieve maximum discrimination accuracy (Uchida and Mainen, 2003; Zariwala et al., 2013). These short sampling times suggest that processing of olfactory information is rapid. Wesson et al. (2008) showed that processing is indeed rapid by demonstrating that the minimal response time with above-chance accuracy from odor inhalation to an action for an innate behavior is 140 ms, which corresponded to the duration of an entire sniff cycle.
In the current work, we sought to define the limits on the speed of olfactory processing. To achieve this, we developed a novel, two-alternative choice, head-fixed odor concentration discrimination paradigm based on temporally precise stimulus delivery.
Materials and Methods
Animals.
Data were collected from six male C57BL/6 mice (8 weeks old at the start of training). Mice were on a 12 h light/dark cycle (lights on at 8:00 pm). All experimental procedures were approved by the Howard Hughes Medical Institute Institutional Animal Care and Use Committee.
Sniff recording.
We implanted a cannula (6 mm long, 25 gauge hypodermic tubing) in the nasal cavity. During experiments, the cannula was connected to a pressure sensor (MPX5050DP; Freescale Semiconductor) with capillary tubing (2.4 mm inner diameter, 7 cm long) and the pressure signal was sent to an amplifier (Cygnus Technology) with gain set at 50× and no filters selected. We verified the timing of the pressure signal with a mass flow sensor (AWM 3300; Honeywell), with a typical response time of 1 ms, as described previously (Shusterman et al., 2011). The time delay between the air flow measurement and the pressure measurement did not exceed 4 ms. The sniffing signal was sent to an external microprocessor (Mega1280; Arduino), which waited for the sniffing pressure to cross a threshold value before sending out a digital trigger to the odor valve driver (CoolDrive; NResearch). We set the threshold value to the value of the smallest exhalation peak that we observed at the beginning of each behavioral session. This was sufficient to ensure that, on >85% of trials, the odor came on at least 40 ms before inhalation (see Fig. 2 G, inset).
Surgery.
Each mouse was implanted with a head bar for head fixation and a sniffing cannula, as described previously (Smear et al., 2011), and allowed to recover for 3 d before restriction to 1 ml of water per day.
Behavioral setup.
See Figure 1 A for the behavioral study setup. The odor vial had either 2 ml of amyl acetate (Sigma-Aldrich) or 2 ml of vanillin dissolved in water (1 g/100 ml). The blank vial was either empty or contained water or mineral oil.
We assembled the treadmill using plastic Lego parts: a belt (part #779938), two wheels (part #970625), and two axles (part #970617), and a custom 3D printed main frame. To monitor the treadmill movements, we coupled a DC motor (FK-260SA; Jameco Electronics) to one of the axles of the treadmill. The motor signal was acquired continuously as an analog signal at 1000 Hz with a National Instruments acquisition board (NI USB-6259). Custom-written MATLAB software (The MathWorks) extracted the data from the acquisition board every ∼60 ms, filtered it (fourth-order Butterworth LP filter with a cutoff frequency of 50 Hz, applied bidirectionally, see Fig. 4 A for example traces), and then compared it with a threshold value to determine whether a movement had occurred. We used the filter because we occasionally saw voltage spikes in the DC motor signal. We specifically chose the filter to not introduce time shifts or amplitude changes. To verify that the filter did not induce time shifts, we estimated the response time using the unfiltered trace and calculated the difference for each trial between the two response times, which was 0.7 ± 1.8 ms (mean ± absolute mean deviation, all data pooled together). We chose the threshold value of 10 mV to be high enough to disregard the occasional rise in the voltage signal when mice grabbed the belt and lifted it up, which produced a signal of ∼5 mV. In our offline analysis, we identified movement onset using the continuous motor signal and a lower threshold of 2 mV (see Fig. 4 A for the range of peak values; scale bar, 10 cm/s = 55 mV). We occasionally saw changes of direction of movement; however, the choice (and response time) was always scored based on the direction (and onset) of the earliest movement.
To deliver the water reward (1–2 μl), we used gravitation flow under the control of a solenoid valve (NResearch). The valve was connected to a lick spout (5 cm hypodermic tubing, gauge 21) via Tygon tubing.
Behavioral training.
Training began at least 7 d after the start of water deprivation. In the first session, mice received water every time they licked. In the second session, mice had to move the treadmill in any direction to receive water. In the third and fourth sessions, we alternated which direction of the treadmill movement was rewarded in blocks of 100 trials to reduce bias for a specific direction. In the next three to five sessions, mice had to move the treadmill leftward/rightward for the low/high odor concentration. We alternated which concentration was presented in blocks of 100 trials.
In the next eight to nine sessions, the stimuli were in pseudorandom order: in a block of 20 trials, there was an equal number of low and high concentration stimuli in a random order. At this stage, mice typically made errors because either they had a bias for a particular movement direction or they made premature movements before odor inhalation. If in a block of 40 trials mice performed at <75% accuracy, we introduced in the next block of 40 trials both of these measures: (1) to reduce the bias, a 90% probability of repeating a trial if it was an error, and (2) to reduce premature movements, a “grace period” of 200 ms after the time of odor valve opening, during which any choice made was not considered. The grace period was gradually decreased to 150 ms and then to 100 ms. The training continued until mice could perform consistently with accuracy >75% with no grace period.
For data collection, we used a corrective algorithm: if overall accuracy in the last 55–70 trials dropped below 70–72% correct, the next set of 55–70 trials were “corrective trials” and included a grace period of 50–100 ms and a 90% probability of repeating an error trial. These corrective trials were not included in the analysis. We included in our analysis all sessions that had at least 100 consecutive trials regardless of their accuracy (no nonconsecutive trials were included). The number of trials per session in our analysis was 315 ± 140 (mean ± SD) and the total number of trials per mouse was 478–1226.
Data analysis.
To estimate for each mouse the response time corresponding to 75% accuracy, we binned the response times using a bin size of 20 ms. We then fitted a logistic function to the choice accuracy as a function of response time using maximum likelihood estimation (MLE) procedures, assuming a binomial distribution of correct and incorrect choices at each time bin (Dobson and Barnett, 2008). In these fits, the probability of being correct is given by the following: P correct = (1 + exp (b 0 + b 1 RT))−1, where RT is the response time. We calculated the SEs for the parameters of the logistic function using the inverse of the information matrix from the MLE fits (Dobson and Barnett, 2008). We then used these SEs to calculate a normal distribution for each parameter. We estimated 95% confidence intervals for the response time that corresponded to 75% accuracy by randomly drawing parameter values from their respective distributions 10,000 times.
To create the distribution of peak velocities shown in Figure 4 A, we excluded values as outliers if they were larger than q3 + 1.5(q3 − q1) or smaller than q1 − 1.5(q3 − q1), where q1 and q3 are the 25th and 75th percentiles, respectively. This range corresponds to approximately ±2.7σ and 99.3% coverage for normally distributed data.
To compute autocorrelation coefficients, we used the following equation: where Ck is the coefficient at a particular lag k, N is the total number of trials, yi is the value of peak velocity (or odor concentration or direction of movement) for trial i, and ȳ is the mean value across trials.
Statistical analysis.
To evaluate whether a longer time interval from the odor onset to the inhalation onset, t prior, leads to improved accuracy, we calculated the mean, μ, and SE of the distribution shown in Figure 2 H and evaluated the null hypothesis: {H 0 : μ ≥ 0} using the t test.
To evaluate whether higher peak velocity correlates with higher probability of being correct, we fitted the logistic function: P correct = (1 + exp (b 0 + b 1 v peak))−1, where v peak is peak velocity, and then evaluated the null hypothesis: {H 0 : b 1 ≤ 0}. The stated p-value is the probability that a single observation from the normal distribution with mean b 1 and its associated SE will fall in the interval [−∞, 0].
To investigate whether peak velocity carries additional information about choice accuracy that is not explained by response time, we divided the data into three equal groups according to peak velocity (high, medium, and low peak velocity; see Fig. 4 D) and fitted the logistic function: P correct = (1 + exp (b 0 + b 1 RT + b 2 I + b 3 I RT + b 4 J + b 5 J RT))−1, where RT is the response time and I and J are indicator variables (I = 1 if the trial is in the medium peak velocity group and zero otherwise and J = 1 if the trial is in the high peak velocity group and zero otherwise). We evaluated the null hypothesis: {H 0: b 4 ≤ 0}.
Results
We trained head-fixed mice to discriminate two different concentrations of an odor by moving a horizontal treadmill leftward or rightward with their front limbs (Fig. 1 A) in a task similar to that of Sanders and Kepecs (2012). The key feature of our setup was precise triggering of the odor stimulus delivery during exhalation so that the odor stimulus stabilized before the first stimulus inhalation, thus increasing the trial to trial reproducibility (Fig. 1 B). To ensure that the behavioral task was not limited by stimulus noise, we chose two stimuli that were clearly distinguishable from each other (0.02% and 0.8% of saturated vapor pressure of amyl acetate). The temporal precision of the stimulus delivery system was higher than the temporal resolution of the olfactory system, ∼10 ms, as reported in prior work (Smear et al., 2011).
Under these conditions, mice responded very rapidly. Their median response time, defined as the time from inhalation onset to movement onset, varied from 60 to 130 ms (Figs. 1 C, 2 A). Their choice accuracy increased for longer response times and this functional relationship was well fit by a logistic function (R 2 = 0.89–0.98 across mice; Figs. 1 D, 2 B). Based on these fits, mice achieved 75% accuracy for movements initiated 70–90 ms after the onset of inhalation (Fig. 2 B), which was before a single sniff was over (Fig. 2 C,D), and achieved 95% accuracy for movements initiated 100–150 ms after the onset of inhalation. Eliminating odor cues by having the same odor concentration in both vials decreased the choice accuracy to chance levels (p > 0.23, binomial test; Fig. 2 E), thus showing that mice relied on odor information to solve the task. These experiments show that mice can make discriminations of different odor concentrations in tens of milliseconds.
Our estimation of response times is based on the assumption that mice need to inhale the odor to make a decision. If mice can acquire odor information during exhalation, then our response time estimates would be extended by 50–90 ms (range of median times across mice), which is the interval from the onset of odor delivery to the onset of odor inhalation, t prior (Fig. 2 G, top). To exclude this possibility, we triggered odor delivery 10–25 ms before inhalation (median t prior across mice). We found that responses initiated 100 ms after odor onset had an accuracy of ∼75% (Fig. 2 F). Given that the odor reaches the odor port ∼25 ms after final valve onset (Fig. 1 B), mice can use up to 75 ms to reach 75% accuracy. This is in agreement with response times in the main experiment, in which mice reached 75% accuracy 70–90 ms after the onset of odor inhalation. Furthermore, we divided the data in the main experiment into three groups depending on t prior. We found that choices with similar response times after the onset of inhalation were not significantly more accurate if they were preceded by a longer odor interval before inhalation (p = 0.99, see Materials and Methods; Fig. 2 G–H). Finally, to assess whether a bout of rapid sniffs resulted in shorter response times because there may have been some odor inhaled during t prior, we divided the data in the main experiment into two groups: the 25% of trials with the shortest prior intersniff interval and the remainder of the trials. We found that shorter prior intersniff intervals were not followed by significantly shorter response times in 3/4 of mice (p = 0.06–0.93, Wilcoxon rank-sum test; we excluded one mouse because we did not acquire its full prior sniff). Therefore, we found no evidence that mice used odor information before inhalation.
To address whether the trigeminal pathway is involved in these rapid discriminations, we trained two mice to discriminate two concentrations of vanillin, an odorant that does not activate the trigeminal (Doty et al., 1978; Frasnelli et al., 2011). Some of the trigeminal neurons in the nasal epithelium have direct projections to the brainstem (Schaefer et al., 2002), so they create a direct and potentially rapid pathway for processing of odor information. For this odorant, mice reached 75% discrimination accuracy for responses initiated 90–140 ms after inhalation onset (Fig. 3 A), which is about 50 ms longer than the response times for amyl acetate (Fig. 3 A); however, this time is still shorter than the duration of the sniff cycle (Fig. 3 B). This shift may also be attributed to a significantly lower absolute concentration of vanillin, which was 105 lower than that of amyl acetate. To further compare our results for amyl acetate and vanillin, we computed the response time in sniff phase coordinates or warped time coordinates using procedures proposed previously by Shusterman et al. (2011). Briefly, we stretched or compressed individual sniffs so that each inhalation/exhalation was the length of the average inhalation/exhalation computed from all mice; the response time was stretched/compressed with its corresponding interval. In this warped sniff time, we found that the shift between amyl acetate and vanillin was only about 20 ms (Fig. 3 C), consistent with the observation that mice take longer sniffs for vanillin (p < 0.05, Wilcoxon rank-sum test). Overall, this experiment shows that mice can process odor information rapidly, even in the absence of trigeminal stimulation.
To further understand the processing of odor information at these short timescales, we quantified the peak movement velocity (Fig. 4 A), which can be thought of as the vigor of the limb movement. Does the behavioral action simply relay the outcome of these rapid decisions or does it also carry information about the likelihood that the decision is correct? We found that the higher the peak velocity, the more likely it was that the choice was correct (p < 0.01; Fig. 4 B). For this analysis, we excluded all the trials in which peak velocity occurred after the mice had received feedback on their choice (8–16% of the trials across mice; similar results were obtained with all trials included). The correlation between peak velocity and accuracy may be explained by intertrial correlation in performance. A mouse may perform with higher accuracy and with stronger limb movements during one block of trials and then decrease its accuracy and peak velocity in another block, for example, due to changes in motivation. However, this possibility cannot explain the results for the four of five mice for which the different values for peak velocity were largely intermixed across trials, as shown by the autocorrelation coefficients for the trial-to-trial sequence of peak velocities (Fig. 4 C). Finally, we have shown that accuracy covaries with both response time and peak movement velocity. Peak movement velocity, however, does not simply reflect the same accuracy information as response time, because movements initiated at the same time are more likely to be correct if their peak velocity is higher (p < 0.05; Fig. 4 D; we excluded two mice that had too few trials because this analysis requires a large number of trials).
Discussion
Together, this study shows that mice can make olfactory discriminations of 75% accuracy in 70–90 ms after odor inhalation. For a low concentration and nontrigeminal odorant, this time was 90–140 ms, showing that processing of olfactory information is rapid even in the absence of trigeminal stimulation.
Further, this study shows that the triggered behavioral action does not simply relay the outcome of the decision, but it also carries information about the likelihood that the decision is correct, consistent with previous work in humans in which arm movements were used as a motor response (Selen et al., 2012). What might peak movement velocity indicate? We speculate that high peak movement velocity is indicative of a “prepared” motor system to act and is therefore indicative of shorter motor delays, perhaps similar to the expression of saccades in primates (Fischer and Boch, 1983). It is known that the response time is made up of both a decision time, which improves accuracy, and a nondecision time, which reflects sensory and motor delays (Resulaj et al., 2009). We suspect that movements with higher peak velocity may have shorter nondecision times and thus may be more likely to reflect longer decision times and hence higher accuracies.
Our study sought to define the limit on the speed of olfactory processing. Ultimately, the distribution of response times chosen by the animal depends on weighing many factors, such as the cost of a wrong choice (for example punishment), the value of a reward now versus a surer reward later, the cost of waiting for another trial (i.e., the intertrial period), etc. Although the response time associated with 75% accuracy in our task was 70–90 ms for amyl acetate, the median response time across mice was 60–130 ms.
Previous work has reported that innate (nonlearned) responses about the novelty of a stimulus can be made by rodents in as little as 140 ms (for amyl acetate among other odors; Wesson et al., 2008). This is the minimal response time, which is the earliest response time with above chance accuracy reported for odor detection or discrimination by previous studies. In our current work, the minimal response time for concentration discrimination was 60 ms (Fig. 2 I). One explanation for this difference is that we used limb movements to indicate choice, whereas the study by Wesson et al. (2008) used the increase in sniffing rate in response to novel odors to measure choice. Specifically, rats initiate a second sniff faster when the odor is novel compared with when the odor is a learned one, so the minimum duration of the first sniff for novel odors can be used as a measure of response time. However, this measure of response time depends on how fast rats sniff. In contrast, limb movements are decoupled from sniffing and can be initiated at any time.
Response times of 70–90 ms (mean 76 ms) are very short considering that the first odor responses in sensory areas are detected at ∼30 ms after odor inhalation (Cury and Uchida, 2010). Further, motor delays for limb movements in mice after optogenetic stimulation of motor cortex are ∼25 ms (Hira et al., 2009). Therefore, sensory and motor delays could amount to ∼55 ms, leaving as little as ∼20 ms for processing. This window puts a strong limit on the underlying neural computations and it implies that only the first action potentials across a population of neurons encode this sensory information.
This estimate of processing time is also short for the nontrigeminal odorant: response times are 90–140 ms (mean 115 ms), leaving ∼60 ms for sensory processing. We expect that this is an upper bound for the processing time of nontrigeminal odorants because vanillin has a poor solubility and thus a low concentration (105 lower than that of amyl acetate), which presumably results in a longer delay from odor inhalation to activation of the olfactory receptors.
Previous studies have reported a processing time of as little as 50 ms for mixed trigeminal and olfactory odorants (amyl acetate among other odorants; Wesson et al., 2008). This time is estimated by subtracting the time of the earliest responding glomerulus, which is detected at 120 ± 45 ms after odor inhalation (mean ± SD), from the response time in the same trials. However, a key distinction is that mice in our task frequently responded before the earliest responding glomeruli reported by Wesson et al. (2008): the median response time in our task was 60–130 ms. It remains unclear what contributes to this difference, but the explanation will likely necessitate neural recordings.
Future studies that combine this behavior with neural recordings will elucidate the neural basis of these decisions. In sensory areas, olfactory information can be read out in the sequence of action potentials that tile the entire sniff cycle (Cury and Uchida, 2010; Shusterman et al., 2011). This temporal arrangement may allow discrimination of two different stimuli early in the sniff cycle—that is, within a few action potentials. Indeed, our study suggests that the relevant action potentials forming the neural basis of these decisions are fired in just a few tens of milliseconds.
Footnotes
-
This work was supported by the Howard Hughes Medical Institute. We thank J. Dudman, M. Smear, R. Shusterman, and A. Resulaj for helpful discussions; T. Tabachnik for help with designing the behavioral rig; and T. Bozza, R. Kiani, M. Smear, and C. Wilson for helpful comments on the manuscript.
-
The authors declare no competing financial interests.
- Correspondence should be addressed to Dmitry Rinberg, New York University Neuroscience Institute, New York University Langone Medical Center, 450 East 29th Street, Room 935, New York, NY 10016. rinberg{at}nyu.edu