Abstract
During many natural behaviors the relevant sensory stimuli and motor outputs are difficult to quantify. Furthermore, the high dimensionality of the space of possible stimuli and movements compounds the problem of experimental control. Head fixation facilitates stimulus control and movement tracking, and can be combined with techniques for recording and manipulating neural activity. However, head-fixed mouse behaviors are typically trained through extensive instrumental conditioning. Here we present a whisker-based, tactile virtual reality system for head-fixed mice running on a spherical treadmill. Head-fixed mice displayed natural movements, including running and rhythmic whisking at 16 Hz. Whisking was centered on a set point that changed in concert with running so that more protracted whisking was correlated with faster running. During turning, whiskers moved in an asymmetric manner, with more retracted whisker positions in the turn direction and protracted whisker movements on the other side. Under some conditions, whisker movements were phase-coupled to strides. We simulated a virtual reality tactile corridor, consisting of two moveable walls controlled in a closed-loop by running speed and direction. Mice used their whiskers to track the walls of the winding corridor without training. Whisker curvature changes, which cause forces in the sensory follicles at the base of the whiskers, were tightly coupled to distance from the walls. Our behavioral system allows for precise control of sensorimotor variables during natural tactile navigation.
Introduction
Investigating neural circuits requires recording and manipulating defined cell types in behaving animals (O'Connor et al., 2009). Genetically encoded reporters (Chen et al., 2013) and effectors (Fenno et al., 2011) of neural activity can be combined with transgenic mice (Heintz, 2004; Madisen et al., 2012; Gerfen et al., 2013) for cell-type-specific imaging and perturbation of neural circuits (Luo et al., 2008). These advances drive the need for mouse behaviors compatible with cellular imaging and optogenetics. Such assays should (1) provide quantitative behavioral data to enable correlations with neural variables, (2) be quick for mice to learn to facilitate throughput in experimentation, and (3) be ethologically relevant to ensure study of the brain under relevant conditions.
Head-fixed behaviors allow for a high degree of stimulus control and precise monitoring of movement (Evarts, 1968; Wurtz, 1968; Boyden et al., 2006; Verhagen et al., 2007; Komiyama et al., 2010; O'Connor et al., 2010). Current head-fixed mouse behaviors involve training mice to make decisions based on sensory stimuli (Andermann et al., 2010; Komiyama et al., 2010; O'Connor et al., 2010) or allow mice to navigate visually defined virtual environments (Harvey et al., 2012). Virtual reality enables sensory stimuli to be presented coupled to the motion of the animal, in closed-loop, and independently in open-loop (Wolf and Heisenberg, 1990; Ahrens et al., 2012; Dombeck and Reiser, 2012; Keller et al., 2012; Saleem et al., 2013). Open-loop and closed-loop experiments can be used to disentangle correlated sensory and motor variables and deduce control laws underlying behavior (Poggio and Reichardt, 1973).
Mice are primarily nocturnal animals, which live in narrow underground burrows (Brecht et al., 1997; Latham and Mason, 2004). Rodents use their whiskers to recognize and locate objects in space (Krupa et al., 2001; Knutsen et al., 2006; Mehta et al., 2007; O'Connor et al., 2010; Pammer et al., 2013) and guide their locomotion (Vincent, 1912). Whiskers are thin tapered rods (Birdwell et al., 2007) with mechanosensory receptors at their base. These translate forces generated by whisker-object contact into neural excitation. Sensory stimuli and motor actions can be quantified by tracking whisker curvature and position using high-speed video (Carvell and Simons, 1990; Bermejo and Zeigler, 2000; Knutsen et al., 2006; Voigts et al., 2008; Grant et al., 2009; O'Connor et al., 2010; Clack et al., 2012).
During locomotion, rodents move their whiskers rhythmically and direct their whiskers toward objects (Mitchinson et al., 2007). Whisker movements can reveal a rodent's expectations about the environment (Grant et al., 2009). Running speed and direction is coupled with average whisker position (Towal and Hartmann, 2006; Arkley et al., 2011; Mitchinson et al., 2011); however, the finer-scale structure of this coupling has been hard to ascertain due to difficulties in measuring whisking in freely moving rodents.
We developed a tactile virtual reality system for head-fixed mice that enables quantitative analysis of whisker-guided locomotion. Here we show that whisking and running are tightly coupled, and that mice naturally guide their locomotion along walls by monitoring wall position, velocity, and acceleration.
Materials and Methods
Mice.
Mice (females or males; >3-months-old) were housed individually in cages with bedding and running wheels (Bio-Serv; K3327 and K3251) in a reverse light cycle room (Table 1). Data from the following strains were pooled: C57BL/6Crl (The Jackson Laboratory; 000664, RRID:IMSR_JAX:000664), scnn1a-tg3-cre (The Jackson Laboratory; 009613, RRID:IMSR_JAX:009613) X AI32 (The Jackson Laboratory; 012569), VGAT-ChR2-EYFP (The Jackson Laboratory; 014548, RRID:IMSR_JAX:014548; Table 1). Mice received 1.5 ml of water per day, which corresponds to ∼50% of ad libitum water consumption for C57BL/6J mice (Mouse Phenome Database from The Jackson Laboratory; http://www.jax.org/phenome). Mice obtained water either during behavioral sessions or by supplemental water given after the session. Weight gain during the behavioral session was taken as a lower bound on the amount of water consumed, and the water supplemented was determined accordingly. The weight and health (posture, condition of fur, and motor activity) of the mice were monitored daily. All procedures were in accordance with protocols approved by the Janelia Farm Institutional Animal Care and Use Committee.
Mice used
Surgical procedures.
Mice were head-fixed using a titanium head postimplant (22.5 mm length; 3.2 mm width; Guo et al., 2014). For implantation, mice were anesthetized with isoflurane, (1–2% by volume in O2; SurgiVet, Smiths Medical), and maintained at 37°C body temperature on a thermal blanket (Harvard Apparatus). Following a local injection of Marcaine (50 μl), the scalp and periosteum over the dorsal surface of the skull were removed. The skull was then coated with a thin layer of cyanoacrylate adhesive (Krazy Glue, Elmer's Products). The head post was placed on the skull with its anterior edge aligned to the lambdoid suture and fixed in place using clear dental acrylic (Jet Repair Acrylic, Lang Dental Manufacturing; P/N 1223-clear). The remaining exposed cyanoacrylate adhesive was covered with a thin layer of dental acrylic. Buprenorphine HCl (0.1 mg/kg, subcutaneous injection; Bedford Laboratories) was used for postoperative analgesia. Ketoprofen (5 mg/kg, subcutaneous injection; Fort Dodge Animal Health) was used at the time of surgery and postoperatively for two additional days to reduce inflammation. Mice were allowed at least 3 d to recover from surgery before being placed on water restriction.
Whisker trimming.
For some experiments (see Fig. 12) whiskers were trimmed progressively: full whisker field; a row of whiskers (C1–C3); a single whisker (C2); no whiskers. Whisker trimming was performed under isofluorane anesthesia and mice were given overnight to recover. At least two behavioral sessions were performed with each whisker configuration.
Data acquisition and control.
Behavioral control was implemented using a real-time system (BControl; http://brodylab.princeton.edu/bcontrol) based on real-time Linux (code.google.com/p/rt-fsm). A MATLAB (MathWorks) program running on Windows XP communicated over Ethernet with the real-time system which also interfaced with PCI-6229 and PCI-6713 data acquisition boards (National Instruments) using COMEDI drivers (comedi.org; O'Connor et al., 2010). The real-time system executed a module of C code at 6 kHz that read and wrote digital and analog voltages on the data acquisition boards for monitoring and controlling the apparatus. The system received commands from MATLAB and sent data back to MATLAB for logging at 4 Hz.
Spherical treadmill.
The spherical treadmill was a 15.56 inch diameter hollow Smoothfoam ball (Plasteel; 16 inch diameter, Ball no. 183). Balls were purchased as hollow halves with an initial wall thickness of 19.5 ± 0.2 mm and weight of 121.3 ± 0.5 g. The inside of the halves was carved with a hot wire system (Hot Wire Foam Factory) to give a wall thickness of 3.79 ± 0.66 mm. Two matched halves were glued together with expanded polystyrene foam glue (Hot Wire Foam Factory; Foam Glue, no. 028B-8). The total weight of the ball was 82.5 g. The ball was supported by 10 Ping-Pong balls (JOOLA Gold 3-Star 40 mm) in air cannons. Each air cannon consisted of a 1.577 inch diameter acrylic tube plugged at one end with an acetyl resin base plate containing a tube fitting (McMaster Carr, no. 50745K15) for air flow. The air cannons were clamped at regular intervals around the ball in a custom 19 inch diameter acrylic ribbed bowl. One cannon was located under the bottom of the ball, three cannons were located in a ring at a latitude of 60° S, and six were located in a ring at a latitude of 20° S. The airflow to the bottom cannon, the ring of cannons at 60°S, and two groups of three cannons at 20° S was controlled independently with regulators (McMaster Carr, nos. 3846K29 and 5627K511).
Ball tracking.
Rotation of the ball was tracked using two cameras containing chips that measure optic flow (Avago Technologies, ADNS-6090) with a serial interface to a microcontroller (Atmel, ATMega644p; Seelig et al., 2010). The cameras were mounted on the ribbed bowl around the equator of the ball, 45 mm from the ball surface. One camera was directly in front of the mouse, and one was on the right. The surface of the ball in front of each camera was illuminated using a 940 nm IR LED (Roithner, ELJ-940-629). Each camera imaged a 2.2 × 2.2 mm2 field of view onto a 30 × 30 pixel sensor with a lens (focal length 25 mm; Computar) and a 10 mm extension tube. The microcontroller interfaced with a PC running MATLAB via serial communication and provided shutter time (inverse of light level) and contrast metrics. These parameters were used to focus the lenses and adjust the illumination during initial setup. For tracking ball motion, every 2 ms the microcontroller received two signed integers measuring the optic flow Δx and Δy from each camera. These signals were converted to analog voltages ranging from 0 to 5 V, centered on 2.5 V, with 150 mV resolution. The microcontroller also provided a clock signal. The RTLinux was triggered by the clock signal to read and rediscretize the analog values of the optic flow. The vector of camera motion displacement signals was transformed to a vector of ball motion vball by multiplication with a calibration matrix Acalib, vball = Acalib × [Δx1, Δy1, Δx2, Δy2]′.
Ball rotation calibration.
To determine Acalib, the rotation of the ball was recorded independently of the ball tracking system using a high-speed (500 Hz) camera (Mikrotron MC1362). The field-of-view was 11 × 11 mm2 (27 pixels per mm). The camera was focused on the surface of the top of the ball, where the head of the mouse would be. Each video image was corrected for uneven illumination. For each frame, the forward and sideways rotation of the ball was computed from the x and y displacements that gave the peak cross-correlation with the next frame. Acalib was determined by least squares fitting a linear transformation from the camera-motion displacement signals to the recorded forward and sideways rotation of the ball. The predicted ball motion and actual ball motion had a correlation coefficient of 0.97. Thus, the real-time tracking system provided a high bandwidth and accurate measurement of ball motion.
Moveable walls for tactile virtual reality.
Two walls were mounted on either side of the mouse. For most experiments, each wall (height, 0.75 inches; length, 2 inches) was a 0.04-inch-thick aluminum sheet; in some cases, the walls were standard glass microscope slides (Fischerbrand, 12 -550-343). The use of glass walls allowed unobstructed video of the mouse from the side. The lower edges of the walls were 1 mm above the surface of the ball and parallel to the anterior–posterior axis of the mouse. The back edges of the walls were positioned in line with the posterior end of the whisker pad. This prevented the front paws of the mouse from touching the walls. The head of the mouse was positioned such that the distance from the surface of the ball to the center of the eyes was ∼19 mm and to the tip of the nose was ∼15 mm. Each wall was attached to an x-y translation stage moved by linear servomotors (MX80 L, Parker Hannifin). The linear servos were tuned without integral control to maximize response speed. Each axis had 50 mm of travel and could be controlled either from software or via an analog voltage provided directly to the controller (ACR9000, Parker Hannifin). The analog input was used to update the position of the walls at 500 Hz. Each axis was controlled over 40 mm of travel. The repeatability of the wall position was ±13 μm (SD), as determined by tracking wall position with the same high-speed video system that was used for calibrating the ball rotation, (37 μm pixel size). The analog input was low-pass filtered with a second order low-pass Butterworth filter (cutoff frequency, 100 Hz). To characterize the control of the servomotors, sinusoidal command signals (frequencies 0.5–7.5 Hz, amplitude 5 mm) were sent to the wall while tracking wall position with high-speed video. The average lag of the wall position was 13 ms across frequencies. The ratio of the measured motion amplitude to the command amplitude ranged from 1.00 to 1.12 across frequencies. Before the start of each experiment, the software control was used to home and position the walls in the desired place. The walls were not allowed to come closer than 4 mm to the face, ensuring that mice touched the walls only using their whiskers.
Delivery of water rewards.
Water was delivered through a 0.05 inch outer diameter stainless steel tube (Small Parts, HTX-18H-36-10). Water flow was controlled using a solenoid valve (Lee Company, LHDA0533115H) such that each droplet was ∼6.6 μl. The tube was attached to a 3D printed thermoplastic piece that could be positioned using a three-axis stepper motor system controlled via a joystick (Zaber 3X, NA11B30 and T-JOY3). Licks were detected using an electrical lick detector that uses a small current to trigger a voltage change when a lick completes a circuit and triggers a change in the output voltage of the detector (Slotnick, 2009). Lick rate was calculated by convolving the individually detected licks with a 1 s window.
Infrared illumination system.
The tactile virtual reality system was contained in a light-tight box built from 25 mm construction rails (Thorlabs, XE25 series) and aluminum sheets held in place with magnetic tape (McMaster Carr, 5759K36). For videography, the area around the mouse was illuminated with four 940 nm IR LEDs (Roithner, TO package 940-66-60) powered by LED drivers (BuckPuck 3023-D-E-1000) and focused with 40 mm focal length aspheric condenser lenses (Thorlabs, ACL5040-B). A fiber optic illuminator (Edmund Optics, Dolan-Jenner Model 180) controlled by a foot pedal (Treadlite II T-91-SC3) provided visible light for placing the animal into the head-fixation apparatus. For trials with visual illumination two 470 nm 700 mA LEDs (Luxeon Star) were used (see Fig. 13).
Video tracking of movements.
Video of whisker movements was acquired at 500 Hz, 640 × 480 pixels, using a high-speed camera (Mikrotron, MC1362) with a 0.36× telecentric lens (Edmund Optics, no. 58-257), and frame grabber (Bitflow) controlled by StreamPix5 multicamera software (Norpix). Each frame was triggered by a transistor–transistor logic (TTL) pulse from the real-time system. A mirror directed light from the whiskers to the camera (Ted Pella; 26002-G). For some sessions, a piece of laminated white paper was placed immediately on top of the surface of the ball and under the whiskers to provide a uniform background for whisker video.
Video of other mouse movements was acquired at 100 Hz from two cameras (Basler, acA640) in StreamPix5. Each frame was triggered by a TTL pulse from the real-time system. Either every trial, or every 12 s, whichever was shorter, a trigger was sent by the RTLinux to instruct StreamPix5 to create a new video file.
Locomotion.
High-resolution forward and lateral components of velocity were calculated from the ball motion displacements by filtering with a low-pass filter (four pole Butterworth, 30 Hz). Running speed was the magnitude of the velocity. A running speed threshold was defined as the minimum of the speed distribution between 0 cm/s and the peak running speed of each mouse (Fig. 1D). The lateral acceleration was the derivative of the lateral velocity component. The run angle was the angle of the velocity smoothed with a moving average (500 ms window). The smoothing was applied to average angle fluctuations produced by individual strides. Positive values of run angle correspond to running left. Power spectra of the lateral acceleration were calculated for each trial with Welch's method (1 s Hamming window, 90% overlap, 1024 fft points) and averaged together across 4 s trials where the mean running speed was above the running speed threshold (Fig. 2C). To generate a probability density function of stride frequency (Fig. 2D), the peak frequency (<6 Hz) of the power spectrum was calculated on a trial-by-trial basis. The stride signal (Fig. 2B) was the bandpass filtered lateral acceleration (four pole Butterworth, 2–6 Hz). Stride amplitude and phase were calculated as the magnitude and angle of the Hilbert transform of the stride signal. Stride frequency is the derivative of the unwrapped stride phase computed with the Savitzky Golay method (fourth order, 400 ms window). The slope of stride frequency against running speed was calculated by linear regression after discarding speeds <0.1 cm/s (Fig. 2F).
Whisker tracking and whisking variables.
Whisker angle and curvature were tracked using automated, freely available software (O'Connor et al., 2010; Clack et al., 2012; http://openwiki.janelia.org/wiki/display/MyersLab/Whisker+Tracking). All whisker tracking was done on mice with only C2 whiskers either unilaterally or bilaterally. Whisker angle was measured as the angle between the whisker and a line perpendicular to the midline of the mouse. Thus, 0° corresponds to a whisker perpendicular to the mouse, positive values correspond to protracted whiskers, and negative values correspond to retracted angles. Trials with whisker tracking errors in >10% of images (27% of trials) were rejected. For the remaining trials, whisker angle was interpolated over the mistracked frames. Whisker set point was calculated as the moving average of whisker angle (500 ms window). Whisking amplitude and phase were calculated as the magnitude and angle of the Hilbert transform of the bandpass filtered whisker angle (four pole Butterworth, 8–30 Hz; Hill et al., 2011). Whisking frequency was the derivative of the unwrapped whisker phase taken with the Savitzky Golay method (fourth order, 400 ms window). The slope of whisking variables against running speed was calculated by linear regression after discarding speeds less than the running threshold (see Fig. 4D,F, H). Power spectra of whisking were calculated from the whisking angle for each trial with Welch's method (1 s Hamming window, 90% overlap, 1024 fft points) and averaged together across trials where the mean running speed was above the running speed threshold (Fig. 3E). To generate a probability density function of whisking frequency (Fig. 3F), the peak frequency (>10 Hz) of the power spectrum was calculated on a trial-by-trial basis (Berg and Kleinfeld, 2003). Whisking frequency was also analyzed when mice were within 8 mm of a wall and when mice were running in closed-loop with the walls. A low-frequency component of the whisking signal was isolated with a bandpass filter (four pole Butterworth, 2–8 Hz). Amplitude, phase, and frequency of this component of whisking were calculated using the Hilbert transform as described above. The normalized strength of this low-frequency component was calculated as the ratio of its amplitude to the sum of its amplitude and the whisking amplitude (see Fig. 5D,E). The coherence of whisker angle and lateral acceleration was computed using Welch's method (MATLAB mscohere; 1 s Hamming window, 90% overlap, 1024 fft points) on the unfiltered signals on a trial-by-trial basis for trials where the mean running speed was above the running speed threshold (see Fig. 5B).
Whisker curvature was measured 1 mm from the base. The intrinsic curvature of the whisker in the absence of contact was subtracted off (Pammer et al., 2013). Measured whisker curvature changes during light touches are likely to underestimate the actual whisker curvature change due to resolution of the images. Whisker curvature was quantified only in mice with a single C2 whisker. In this case, most touches were retraction touches (see Fig. 7A). This resulted in negative values of whisker curvature and the absolute value of the whisker curvature was taken (see Fig. 7B). The maximum whisker curvature on each trial was calculated as the 90th percentile of the whisker curvature and averaged together across trials where the mean running speed was above the running speed threshold (see Fig. 7D). In mice with multiple whiskers both protraction and retraction touches were observed which were either of a sustained or transient nature.
Acclimatization to head-fixation and the behavioral apparatus.
Acclimatization lasted between one and three 15–30 min daily sessions (see Fig. 11). During these sessions mice were head-fixed on the ball and rewarded for forward motion. Initially the mice were rewarded for every 40 cm of movement. As mice ran more, the reward distance threshold was steadily increased up to a final value of 200 cm.
Open-loop trials.
Open-loop trials either contained no walls within reach of the mouse (see Figs. 1⇓⇓⇓⇓–6) or had the right hand wall move to a fixed distance from the mouse for 2 s before moving away (see Figs. 7, 8). Each trial lasted 4 s and was separated by a 1 s intertrial interval. The right wall moved from out of reach of the mouse to 4, 6, 8, 10, 12, 14, 16, 20, 24, 29, or 34 mm from the mouse's whisker pad during the first second of the trial. The wall returned to 34 mm from the face during the last second of the trial. Data from the middle two seconds of the trial, when the wall was in a fixed position, was analyzed. Mice were rewarded with one drop of water for every 200 cm of distance covered on the ball, regardless of trial start and end and independent of behavioral performance.
Closed-loop control of wall position.
For the tactile virtual reality system the walls moved in closed loop with the movements of the mice (see Figs. 9⇓⇓–12). Lateral motion of the walls was controlled in closed-loop using the component of the ball velocity perpendicular to a coupling vector (Table 3; Fig. 9A). The length and sign of the coupling vector, the coupling gain, γ, determines how far and in which direction the walls move for a given movement of the ball. A negative value of the coupling gain means that movement on the ball toward a wall moves that wall closer, as expected when moving toward a real wall. A positive value of the coupling gain does the opposite. A gain value of γ = −0.2 (i.e., 2 mm of wall motion for 1 cm of ball motion) was used. The choice of gamma depends on the weight of the ball and the friction of the system, as the ability of the mouse to move the ball depends on these values. γ = −0.2 was chosen to prevent large wall movements associated with the stepping of the animal during straight running. The angle of the coupling vector, referred to as the turn angle,ψ, determines an ideal movement direction in the forward-sideways plane along which motion produces no movement of the walls. When the run angle matches the turn angle, the walls do not move. ψ = 0° corresponds to a straight corridor as in this case pure forward motion results in no movement of the walls. Positive values of ψ corresponded to left turns and negative values to right turns. In the tactile virtual reality, the walls were always controlled together so that the corridor had a fixed width of 30 mm travel. The walls were 45 mm apart to allow for space for the mouse's head in the center of the corridor.
Running and whisking parameters
Closed-loop virtual reality corridor.
Mice ran through a winding corridor in tactile virtual reality (see Figs. 9⇓⇓–12). Trials were 200 cm long. Trials were either straight or contained a 100 cm long left or right 11.3° turn during the middle of the trial. Trials of different turn angles were randomly interleaved. For the majority of experiments, during the last 10 cm of each trial, the walls reset themselves such that mice started every trial in the center of the corridor. During some experiments, the walls did not reset themselves, and a new trial would not start unless the mouse was in the middle 87% of the corridor. Mice received a drop of water if they were in the middle 87% of the corridor during the 170–180th cm of the trial.
Task performance was evaluated using the absolute difference between the run angle of the mice and the turn angle over the last 50 cm of the turn. This difference is termed the “angle error”. The smallest angle error that mice can achieve without knowledge of the wall position occurs if they run perfectly straight. This behavior will result in a mean angle error of 7.53° = 2/3 × 11.3° (as 2/3 of the trials contain turns). If mice guess random turn angles, the angle error will be higher. For statistical evaluation of performance of the task (see Fig. 10C), angle errors were compared with this chance level using a one-sided t test for each animal. For statistical evaluation of performance in two different conditions (see Figs. 11A, 12A, 13A) one-sided unpaired t tests were done for each animal.
Model fitting.
Only trials where mice stayed in the middle 20 mm of the corridor were used for modeling as these trials are indicative of good performance (see Fig. 14). Only mice which had at least 20 acceptable trials on both left and right turns were included in the modeling (14 of 23 mice). For each trial, the run angle was computed as a function of forward distance in the trial by taking the derivative of the lateral run trajectory with respect to forward distance using the Savitzky–Golay method (fourth order, 12.5 cm window). The mean run angle trajectory on left turns and the negative of the mean run angle trajectory on right turns were averaged together to remove left/right biases. The run angle trajectory was further normalized by subtraction of its mean value during the 25 cm before the onset of the turn to give the run angle trajectory z used for modeling. The error in the mean run angle trajectory was calculated using bootstrapping with 1000 resamples of the included trials. Our model makes the following assumptions: (1) the mouse wants to run straight (as mean left/right biases have been subtracted), (2) if the mouse deviates from the its ideal angle it corrects its trajectory at a rate that is proportional to the deviation, (3) the mouse experiences a damping force proportional to its velocity, and (4) turns act like external perturbations. This produces dynamics of a forced damped harmonic oscillator z″ + 2ζω0z′ + ω02z = ω02F (Table 3). Where, ζ is the damping coefficient, ω0 is the undamped angular frequency of the oscillator, and F is the applied force, which should match the turn angle. The run angle trajectory z was fit to solutions of the forced damped harmonic oscillator, allowing for the damping coefficient ζ to be <1 (under-damped), equal to 1 (critically damped), or >1 (overdamped).
The three parameters, ζ, ω0, and F were fit using a nonlinear regression (MATLAB, nlinfit). Corridor trajectories from left and right trials were averaged together, after the right trials had been flipped in sign to give the wall trajectory u, which represents deviation from the baseline wall distance. The derivative of the wall trajectory was then modeled as û′ = γ(ψ − z), where ψ represents the turn angle, and γ represents the coupling gain (see Fig. 14B). This equation is a small angle approximation of the true equation coupling ball motion to wall motion. The parameters were fit using regression (MATLAB, regress). The run angle was then regressed (MATLAB, regress) against the modeled wall trajectory position, velocity, and acceleration along with a bias term, z = Aû + Bû″ + Cû″ + Dx. This equation is independent of γ and ψ, and describes how mice transform sensory stimuli into motor output.
Definition of variables
Results
Head-fixed mice exhibit natural running and whisking on a spherical treadmill
Our goal was to develop a behavioral assay for whisker-guided locomotion. To achieve precise measurement of movement and tactile input we used head-fixed mice running on a large (diameter, 15.5 inch), light (82.5 g) spherical treadmill (Fig. 1A). We first quantified free running and whisking behavior of head-fixed mice in the absence of walls. Movements of the ball were tracked in two-dimensions at 500 Hz (Seelig et al., 2010). Whisker movements were recorded using high-speed video at 500 Hz and tracked using automated software (Clack et al., 2012). For simplicity in whisker tracking, mice were trimmed to have only the C2 whisker on each side.
A whisker-based tactile virtual reality system. A, Side view of a mouse on a spherical treadmill. Ball movement was measured using two optic flow sensors (only one sensor is shown). Values of ball movement were used to update the position of two walls on either side of the mouse. B, Example data showing running speed (black), licking (magenta), and water rewards (blue). C, Higher-magnification view (corresponding to gray bar in B). D, Histograms of speed for four mice (gray line corresponds to data in B).
We motivated mice to run on the treadmill by providing water rewards for every 200 cm of forward travel. Mice thus ran in short bursts, and stopped to drink in between bursts (Fig. 1B,C). The distributions of running speed showed peaks at zero and 23 ± 4.1 cm/s (4 mice), which corresponds to the preferred running speed of the mice. The distributions also have a minimum at 8.3 ± 3.4 cm/s (4 mice), which we define as the running speed threshold (Fig. 1D,E). Mice stopped 136 ± 71 ms (4 mice) after onset of water delivery and drank for 4.8 ± 2.4 s (4 mice) before starting to run again (Fig. 1E). While drinking, mice licked at rates of 7.7 ± 0.2 Hz (4 mice; Fig. 1E).
We extracted mouse strides, the successive touchdowns of the paws, from video images (Fig. 2A). We noticed that the lateral ball acceleration was synchronized with the stride frequency, allowing us to compute the stride frequency from the ball movement (Fig. 2A,B). Stride frequency during running was 3.8 ± 0.5 Hz (4 mice) (Fig. 2C,D) and increased with running speed as expected (Fig. 2 E,F).
Stride frequency. A, Time-lapse images showing one stride. The right hindpaw was rendered red. B, Example traces of paw position along the longitudinal axis determined manually from video (top) and the filtered lateral ball acceleration (bottom). C, The power spectrum of the lateral ball acceleration (mean ± SD). The peak between 3 and 4 Hz corresponds to the stride frequency of the mouse. D, Probability density function of the peak frequency of the power spectrum, four mice, gray line corresponds to mouse in C. E, Stride frequency as a function of running speed across time points (mean ± SD; same mouse as in B, C). F, Stride frequency as a function of running speed for four mice.
Whisker movements were measured from high-speed videography (Fig. 3A). During running, mice moved their whiskers in a rhythmic fashion (Fig. 3B,C). The average whisking amplitude was 9.3 ± 1.1° (4 mice; Fig. 3C,D). Power spectra of whisking frequency showed peaks at 3.9 ± 0.6 Hz (4 mice) and 16.9 ± 0.5 Hz (4 mice; Fig. 3E), implying two timescales of whisking. The lower-frequency is close to the stride frequency, whereas the faster-frequency (Fig. 3F) is close to the typical whisking frequency seen in freely moving mice (Voigts et al., 2008).
Whisker movements of mice running on a spherical treadmill. A, Left, Definition of the whisker angle. Right, Time-lapse images of a mouse with a single whisker on the right hand side, moving over half a whisker cycle. The first 7.5 mm of the tracked whisker are highlighted in blue. B, Example trace of whisker angle. C, Blow-up of whisker trace (corresponding to gray bar in B). D, Histogram of fraction of time spent at a particular whisking amplitude when running, four mice; gray line corresponds to mouse in B. E, Power spectrum of the whisker angle for one mouse (mean ± SD). F, Probability density function of the major peak frequency (i.e., >10 Hz) of the power spectrum, four mice; gray line corresponds to mouse in C.
Running speed and whisking are tightly coupled
We investigated the coupling between whisking and running, two basic rhythmic movements. Running speed and whisker set point, the slowly changing whisker position about which the rapid whisker movements occur, were highly correlated. After initiation of running, mice protracted their whiskers by ∼30° (Fig. 4A,B). The faster mice ran, the further forward they pushed their whiskers (Fig. 4C,D). The fast component of the whisking frequency increased monotonically with running speed (Fig. 4E,F). Whisking amplitude initially increased sharply with running speed, but then monotonically decreased with speeds from 10 to 40 cm/s as the set point became more protracted (3 of 4 mice; one mouse did not run at high speeds; Fig. 4G,H).
Coupling between running speed and whisking. A, Example showing running speed (black) and whisker angle (gray). The slow (4 Hz) oscillation in the running speed corresponds to strides and is locked one-to-one to a slow oscillation in the whisker angle. B, Overlay of 300 whisker traces (600 ms) when the mouse is still (left) and during a running bout (right). C, Whisking set point as a function of running speed for all time points for one mouse (mean ± SD). D, Whisking set point as a function of running speed, four mice. E, Whisking frequency as a function of running speed for all time points (mean ± SD). F, Whisking frequency as a function of running speed; four mice. G, Whisking amplitude as a function of running speed across time points (mean ± SD). H, Same as G; four mice.
Whisking epochs either had one dominant whisking frequency ∼15 Hz (Fig. 3B), or contained an additional frequency component ∼4 Hz (Fig. 5A). The 4 Hz whisking component was coherent with the strides of the mice (Fig. 5A–C). The strength of the 4 Hz whisking component increased at higher speeds (3 of 4 mice; one mouse did not run at high speeds; Fig. 5D,E). Thus, at high running speeds whisking and running were coupled on a step-by-step basis.
Whisking and running are coupled on a step-by-step basis. A, Example trace of whisking with a 4 Hz component (black) and strides (gray dots). B, Cross-coherence between whisker angle and lateral ball acceleration across trials for one mouse (mean ± SD; same mouse as Figs. 1⇑⇑⇑–5). C, Cross-coherence, four mice. D, Normalized 4 Hz whisking strength as a function of running speed across time points (mean ± SD). E, Normalized 4 Hz whisking strength, four mice.
Running direction and whisking are tightly coupled
We recorded bilateral whisker movements to assess whisking behavior during turning. Mice moved their left and right whiskers at the same frequency (Fig. 6A). The difference in whisking set point between the left and right side was modulated by running direction. When turning to the left, mice protracted their right whiskers and retracted their left whiskers, resulting in an asymmetric arrangement of the whiskers (Fig. 6B–D). This asymmetry is similar to behaviors performed by freely moving rats (Towal and Hartmann, 2006; Mitchinson et al., 2011). This suggests that mice have not learnt unnatural actions to interact with the ball, but are instead performing natural movements.
Asymmetric whisking during turning. A, Example trace of left whisker angle (red) and right whisker angle (blue) during a transition from turning right to turning left. B, Overlay of 250 whisker traces (500 ms) when the mouse is turning right, running straight, and turning left. C, Right/left difference in whisking set point as a function of run angle for one mouse (mean ± SD). D, Right/left difference in whisking set point as a function of run angle; three mice.
Whisker curvature increased with decreasing wall distance
We next characterized the sensory variables that could underlie whisker-guided locomotion. We recorded the interactions of the C2 whisker with a wall presented at different distances from the mice (Fig. 7A). The movements of the wall were in an open-loop, independent of the movement of the mice. Mice continued to whisk at a normal frequency 16.0 ± 0.5 Hz (4 mice) when running close to the wall. During running along the wall, the C2 whisker was very protracted with a set point of 42.8 ± 2.8° (4 mice). Therefore mice mainly made retraction touches onto the wall, resulting in negative whisker curvature changes (Fig. 7A,B). We analyzed the absolute value of the curvature change. The curvature change in the whiskers during touch is proportional to the forces acting on the mechanoreceptors at the base of the whiskers (Birdwell et al., 2007; Bagdasarian et al., 2013; Pammer et al., 2013). Whisker curvature increased as mice moved their whiskers into the wall. Smaller wall distances resulted in larger slopes of whisker curvature versus whisker angle, as predicted by mechanical models (Fig. 7C,D; Birdwell et al., 2007). Whisker curvature increased monotonically with decreasing wall distance (Fig. 7E,F). The resulting distance-dependent forces likely underlie perception of wall distance.
Interactions between whiskers and the wall. A, Images of a single C2 whisker touching a wall at different distances (4, 6, and 8 mm) or without wall present. B, Cumulative distribution of curvature changes during close wall distance (≤8 mm) trials (dark purple) and trials where the wall was out of reach (pink); four mice. C, Scatter plot of absolute change in whisker curvature against whisker angle for different wall distances (dark purple, 4 mm; medium purple, 6 mm; light purple, 8 mm; pink, no wall present). D, Slope of absolute change in whisker curvature versus angle, four mice. E, Max absolute change in whisker curvature for different wall distances across trials for one mouse (mean ± SD). F, Max absolute change in whisker curvature for different wall distances; four mice.
Mice naturally keep walls at a distance
During open-loop presentations of the wall, mice turned away from the wall. The turn angle increased monotonically with decreasing distance between wall and mouse (Fig. 8A–C). Mice thus naturally kept a safe distance from the wall. Based on this observation, we reasoned that mice would try to stay in the middle of a winding corridor if they could control the position of the walls in closed-loop.
Interaction between mice and walls. A, Schematic illustrating mouse behavior in response to wall location. Mice tend to turn away from the wall, with sharper turn angles for closer wall distances. B, Change in run angle for different wall distances across trials for one mouse. C, Change in run angle; four mice.
Head-fixed mice navigate a two-dimensional tactile virtual reality environment
We next presented mice with two walls in closed-loop with their running speed and direction. This configuration allowed us to test the ability of head-fixed mice to use their whiskers to guide their locomotion through a winding corridor on the spherical treadmill (Fig. 1A). The two walls 45 mm apart, one on each side of the mouse. This width corresponds to typical sizes of mouse tunnels (Latham and Mason, 2004). The faces of mice are >5 mm wide at the whisker pad, so the minimum distance to a wall from the whisker pad was always <20 mm. When mice were in the center of the corridor, they could touch both walls with their whiskers. Experiments were conducted in the dark to remove visual cues. The walls were moved in real-time based on movements of the ball. Changing the coupling between ball motion and wall motion allowed simulation of turns at different angles (Fig. 9A).
Tactile virtual reality environment. A, Schematic illustrating coupling between ball movement and movement of the walls. Turns were defined by a turn angle ψ that was positive for left turns. The coupling gain γ determined the scaling between ball motion and wall motion. The movement of the wall, Δu, was determined based on the coupling equation. The example shown contains a leftward winding corridor (ψ > 0), and a mouse running very left (φ > ψ), resulting in the left wall moving closer to the mouse and the right wall moving farther away. If the mouse run angle matched the turn angle, (φ = ψ), the walls did not move (Δu = 0). B, Trial structure of the winding corridor task. Trials were 200 cm long and contained a single 100 cm left turn, right turn, or no turn, during their middle half. C, Overlay of one example running trajectory from a mouse, and the ideal trajectory if the mouse maintained the same position in the corridor throughout the turn. Performance was quantified by comparing the angle of the mouse's trajectory at the end of the turn and the angle of the turn.
Individual trials were 200 cm long and were either straight or contained a single left or right 11.3° turn over the middle 100 cm of the trial (Fig. 9B). Trials with different turn-directions were randomly interleaved, thus simulating a winding two-dimensional corridor. The two-dimensional trajectory of the mouse on the ball was compared with the ideal running trajectory, which occurred if mice maintained a constant position in the corridor throughout the turn (Fig. 9C). Task performance was quantified using the angle error, the absolute difference between the run angle of the mice and the turn angle during the last 50 cm of the turn. Zero angle error implies that mice track the walls at a fixed distance throughout the turn.
Mice successfully matched their run angle to the turn angle on a trial-by-trial basis (Fig. 10A). Changes in distance to the walls drove turning (Fig. 10B). During straight running, mice tended to stay in the center of the corridor, though some mice favored different sides of the corridor. On left turns, mice ran closer to the right wall, whereas on right turns mice ran closer to the left wall. Mice tracked the walls with angle errors of 2.1 ± 0.9° (23 mice), well below chance levels (Fig. 10C). The median trial duration was 14.1 ± 5.9 s (23 mice), including drinking (Fig. 10D). Mice ran 199 ± 96 m per session (Table 1). Thus, mice naturally track a winding corridor in a virtual reality environment using their whiskers.
Performance in tactile virtual reality. A, Left, Twenty randomly selected running trajectories from each turn angle for one behavioral session. Red, right turns; purple, no turns; blue, left turns; ideal running trajectories, black dashed lines. Middle, Deviation of run angles at given turn angles across trials from one mouse (mean ± SD). Right, Mean run angle versus turn angle; 23 mice. B, Corridor position trajectories for same trials as in A. Thick gray lines represent the left hand wall (LW) and right hand wall (RW). Middle, Corridor position at the end of the turn at given turn angles across trials from one mouse (mean ± SD). Right, Mean corridor position versus turn angle; 23 mice. C, Boxplots of run angle error, 23 of 23 mice showed angle errors significantly below chance levels (p < 0.001, one-sided t test). Gray shading represents range of angle errors that could be generated via chance performance. D, Histogram of trial durations (mean ± SD; 23 mice).
Corridor-tracking does not require training
We next quantified performance as a function of training. Mice were first acclimatized to head-fixation and were allowed to run on the ball in the absence of walls. During these baseline sessions, the running angle error was calculated by randomly assigning each trial a turn angle. These mice were then introduced into the tactile virtual reality environment. The running angle error dropped to steady-state performance levels on the first trial in which mice experienced the simulated winding corridor (Fig. 11A,B). These data show that following the narrow winding corridor is a natural behavior for mice.
Corridor-tracking did not require training. A, Average angle error before (green) and after (black) the walls were introduced. Seven of 7 mice showed a significant decrease in angle error following introduction of the walls (p < 0.001, one-sided unpaired t test). B, Time course of the angle error before and immediately after introduction of the walls; seven mice.
Corridor-tracking is whisker-dependent
We next tested how corridor-tracking in the dark depended on the large mystacial whiskers. Whiskers were gradually trimmed between behavioral sessions. Wall-tracking performance degraded as the number of whiskers was reduced bilaterally, with significant performance drops from all whiskers to the C rows of whiskers (4 of 5 mice, p < 0.02, one-sided unpaired t test), from the C rows to single C2 whiskers (5 of 5 mice, p < 0.01, one-sided unpaired t test), and from single C2 whiskers to no whiskers (5 of 5 mice, p < 0.001, one-sided unpaired t test; Fig. 12A). Whisker movements were also measured in a subset of these mice while they were performing the task with their C2 whiskers. These mice continued to show normal whisking while performing the closed-loop corridor-tracking task, with whisking frequency 15.6 ± 0.2 Hz (3 mice) and whisker set point 40.1 ± 6.7° (3 mice). When mice had all their whiskers, touches occurred with both positive and negative curvature on the protraction and retraction phases of whisking. Mice with fewer whiskers ran closer to the walls (Fig. 12B). Mice without whiskers are disoriented and spend most of their time at the closest allowed wall distance, farthest from the center of the corridor. Thus, running through the middle of the narrow winding corridor requires whiskers. Furthermore, multiple whiskers are normally used together for guiding locomotion.
Corridor-tracking is whisker-dependent. A, Average angle error for mice with different numbers of whiskers. Black, all whiskers; brown, C row of whiskers; dark red, single C2 whisker; bright red, no whisker); four of five mice showed a significant increase in angle error following trimming of all whiskers to a row of whiskers (p < 0.02, one-sided unpaired t test), five of five mice showed a significant increase in angle error following trimming of a row of whiskers to a single whisker (p < 0.01, one-sided unpaired t test), and trimming of a single whisker to no whiskers (p < 0.001, one-sided unpaired t test). B, Average histograms of distance to nearest wall for different whisker configurations; five mice.
Visually mediated corridor-tracking without whiskers
Mice without whiskers may have been unable to perform the corridor-tracking task due to discomfort or a lack of motivation. We tested these hypotheses by illuminating the behavioral apparatus, allowing mice to see the walls. Four of five mice without whiskers partially regained the ability to run through the simulated winding corridor in an illuminated arena (Fig. 13A). This finding shows that in the dark mice without whiskers lacked the necessary sensory information to navigate their environment.
Vision partially compensates for absence of whiskers. Average angle error for mice with all whiskers (black), no whiskers (red), and no whiskers but with illumination provided (yellow). Four of 5 mice showed a significant decrease in angle error following introduction of the light when they had no whiskers (p < 0.001, one-sided unpaired t test).
Modeling whisker-guided turning behavior
How do mice transform sensory inputs to motor control signals as they execute turns in the corridor? Because mice tried to maintain a fixed position in the corridor, u0, we modeled how mice might determine their mean run angle trajectory, z, as a function of their deviation from this ideal position, û = u − u0. This approach is inspired by the framework of control theory, where a controller, in this case the mouse, tries to limit an error, in this case deviation from the ideal corridor position, in a system with feedback, in this case the closed-loop coupling between the ball and the walls (Fig. 14A; Åström and Murray, 2008). We modeled data where mice had all their whiskers so at least one of the walls was within their contact range at all times.
Modeling of whisker-guided turning. A, Box diagram describing control theory inspired model of the behavior. B, Example showing average run angle trajectory (black) with SD across bootstrap estimates (gray) and results of model fit of transformation from sensory stimuli to turn angle (red). C, Example showing average corridor position trajectory (black) with SD across bootstrap estimates (gray) and results of model fit of coupling equation (red).
During a turn, the mean run angle trajectory, z, approached the target value of the turn angle with an oscillatory transient (Fig. 14B). Curves with such transients can be described as solutions of a forced damped harmonic oscillator z″ + 2ζω0z′ + ω02z = ω02F. We thus used a nonlinear fitting approach to determine the coefficients ζ, ω0, and F (Table 4). The distance for mice to reach the steady-state turn angle, given by (ζω0)−1, was 12.1 ± 5.2 cm (14 mice). The parameters for this equation only hold for a specific turn angle.
Parameters from model fits
To determine how mice transformed sensory stimuli into motor output, we fit an equation describing run angle as a function of wall position and its derivatives that is independent of the experimentally set turn angle and coupling gain. The open-loop experiments (Fig. 8) suggest that the mice determine their run angle as a function of wall position; however, the presence of the oscillatory transient in the run angle requires that mice also adjust their run angle as a function of the wall velocity and acceleration. Mice that only monitored wall position would show an exponential approach to the turn angle. Therefore, we fit an equation describing run angle as a function of wall distance, velocity, acceleration and an innate bias, z = Aû + Bû′ + Cû″ + Dx (Fig. 14B; Table 4). The parameters A–D no longer depend on the turn angle or coupling gain. The bias term was included to reflect that some mice naturally run more to the left and others naturally run more to the right. The run angle trajectory was related to the corridor position trajectory, u, by fitting the coupling equation û′ = γ(ψ − z). The average behavior of the mice can be well described by a function of the deviation of wall position from an ideal value, wall velocity, and wall acceleration. This suggests that that neural activity in the whisker system may code for these state variables.
Discussion
We developed a whisker-guided locomotion task for head-fixed mice. The virtual reality aspects of the system allowed for decoupling the sensory stimuli from motor behavior, which would not be possible in a study of freely moving rodents.
The running speeds of the head-fixed mice (Fig. 1) were similar to reports of running in freely moving mice (Clarke and Still, 1999; Saleem et al., 2013), but somewhat slower than some reports of head-fixed mice on smaller spherical treadmills (Harvey et al., 2012). Stride frequencies are also similar to freely moving mice (Clarke and Still, 1999). Whisking frequencies in freely moving mice range from 11.4 ± 1.0 Hz (30 mice) during locomotion with a mean speed of 13.2 ± 8.2 cm/s (30 mice; Mitchinson et al., 2011) to 15–20 Hz in mice performing a gap crossing task (Voigts et al., 2008). These values are in the same range as in our experiments (Fig. 3). Head-fixed mice running on a large spherical treadmill therefore appear to perform natural exploratory behaviors.
Few studies have examined the coupling between running and whisking. The correlation between speed and whisker set point has been reported for rats (Arkley et al., 2011). Whisking asymmetry during head turning has also been noted (Towal and Hartmann, 2006; Mitchinson et al., 2011). In these studies, during left turning the right whiskers were more protracted than the left whiskers, similar to our observed increase in the set point of the right whisker during head-fixed leftward turning. The 4 Hz whisking component and its synchronization with the stride frequency (Fig. 5) has not previously been described. This coupling could be passive, related to stretching of the skin during stepping, or might be under neural control. Breathing was not recorded in our study; however, normal breathing in mice can occur ∼4 Hz (Ranade et al., 2013) and quadrupeds tend to lock breathing and stride in a one-to-one manner (Bramble and Carrier, 1983). As the breathing oscillator can reset the whisking oscillator (Moore et al., 2013), the coherence between whisking and strides could be a consequence of the coupling between whisking and breathing.
The synchronization of whisking and stepping could provide a common time base for neural computations. It is unknown whether touch signals in barrel cortex are modulated by locomotor signals relating to either running velocity or the phase of locomotion. Such modulation would allow for self-motion information to be coupled to sensory input for information processing. During the closed-loop corridor-tracking experiments, mice continued to whisk at normal frequencies; however, they may subtly adapt their whisking strategies to maximize the amount of information about the walls available to them (Grant et al., 2009).
Whisker-guided locomotion in the simulated winding corridor did not require training (Fig. 11), demonstrating that mice experience the virtual reality as natural. Virtual reality experiments using visual cues have been mainly confined to one-dimensional tracks (Keller et al., 2012; Domnisoru et al., 2013; Ravassard et al., 2013; Saleem et al., 2013) or required extensive training to get rodents to make turns and not run through walls (Hölscher et al., 2005; Harvey et al., 2012). An advantage of a tactile virtual reality system is that naive mice correctly navigate two-dimensional environments from the first trial. In the current system, the walls only moved perpendicular to the mouse. Mice therefore lack cues related to the parallel motion of the walls that would normally be present during running along a wall. As the walls were smooth, it is unclear how significant this absence was. A more realistic system would have walls made from rotating cylinders whose speed of rotation could be coupled to the speed of the mouse. Tactile virtual reality can be combined with visual virtual reality to create a system where the boundaries of an environment are defined with tactile walls and where visual information gives the animal distal cues to it location in the environment.
Mice were only lightly water restricted to motivate more consistent running and larger total running distance per session; many mice received rewards on every single trial (Table 1). Water restriction was not necessary for performance of whisker-guided locomotion.
Corridor-tracking performance depended on the number of whiskers (Fig. 12). With multiple rows of whiskers, interactions with the walls were rich and varied, occurring both during protraction and retraction, and lasting different durations. Similar to previous object localization tasks, performance degraded with whisker trimming (Krupa et al., 2001; Knutsen et al., 2006; Mehta et al., 2007; O'Connor et al., 2010). There was a large decrement in performance from three C row whiskers to a C2 single whisker. Interactions between the wall and a C2 single whisker were mainly retraction touches. Multiple whiskers provide mice with a rich set of algorithms to estimate wall distance, some of which are formally inaccessible to mice with single whiskers. This includes for example, strategies based on labeled line coding (Brecht et al., 1997) or triangulation (Solomon and Hartmann, 2011). Because whiskers are naturally of different lengths, a simple strategy to measure distance to a wall would be to determine which whiskers are touching, and which whiskers are not touching a wall. However, mice that were able to use only this strategy would be unable to measure distances with a single whisker, in contrast to behavior that we observed during the open-loop presentations of the wall. Triangulation strategies depend on mice using knowledge of the positions of their whiskers during touch to calculate object distance (Solomon and Hartmann, 2011; Pammer et al., 2013). The large drop in performance may be because mice with multiple whiskers use multiple strategies. However, as mice were still able to track the winding corridor with a single whisker, these strategies cannot be the only ones available to mice.
Whisker curvature change is directly proportional to the forces acting on the mechanoreceptors at the base of the whisker (Birdwell et al., 2007). Because whiskers are linearly tapered, the bending stiffness decreases rapidly with distance from the follicle. The absolute and relative amplitudes of different stress variables change dramatically with the point of force application (i.e., position of contact with the wall) along the whisker length. Distance-dependent whisker mechanics could be used by mice to estimate the wall distance. Consistent with this view, induced whisker curvature increased monotonically with decreasing wall distance (Fig. 7; Birdwell et al., 2007). Whiskers can also slip along the walls. These high-frequency acceleration events could also underlie measurements of wall distance (Hires et al., 2013). Measuring wall-tracking performance for walls with different mechanical properties could further disambiguate radial distance coding strategies that depend on absolute or relative force measurements. Mice may make bilateral comparisons of the forces on their whiskers to stay in the center of the winding corridor.
The phenomenological modeling of turning behavior suggests that mice do not just compute their running angle based on an instantaneous measurement of wall position, but also take into account the velocity and acceleration of the wall (Fig. 14). Computation of these additional variables may require integration of sensory inputs over hundreds of milliseconds.
This tactile virtual reality system allows for rich natural behavior while sensory stimuli and motor output can be precisely quantified. The system can easily be extended beyond simple trials with a single turn to richer environments with multiple turns and changes in the width of the corridor that create choice points. Trials can also be presented in a repeating fashion, creating two-dimensional mazes that mice can plan routes through. This paradigm would allow for study of the role of the hippocampus in a whisker-based navigation paradigm. Tactile virtual reality can be also combined with recordings and manipulation of neural activity to examine coding of object location in the whisker system.
Footnotes
This work was supported by the Howard Hughes Medical Institute. We thank Dan O'Connor and Daniel Huber for help during the initial stages of this project, and Tansu Celikel, Zengcai Guo, Diego Gutnisky, Vivek Jayaraman, Nuo Li, Dan O'Connor, and Tony Prescott for comments on the paper.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr Karel Svoboda, Janelia Farm Research Campus, Howard Hughes Medical Institute, Ashburn, VA 20147. svobodak{at}janelia.hhmi.org