Spatial arrangement of the whiskers of harbor seals ( Phoca vitulina ) compared to whisker arrangements of mice ( Mus musculus ) and rats ( Rattus norvegicus )

Most mammals have specialized facial hairs known as vibrissae (whiskers), sensitive tactile structures that subserve both touch and flow sensing. Different animals have different numbers and geometric arrangements of whiskers, and it seems nearly self-evident that these differences would correlate with functional and behavioral use. To date, however, cross-species comparisons of three-dimensional (3D) whisker array geometry have been limited because standard morphometric techniques cannot be applied. Our laboratory recently developed a novel approach to enable quantitative, cross-species vibrissal array comparisons. Here we quantify the 3D morphology of the vibrissal array of the harbor seal ( Phoca vitulina ), construct a CAD model of the array, and compare array morphologies of harbor seals, mice ( Mus musculus ) and rats ( Rattus norvegicus ). In all three species whisker arclength decreases from caudal to rostral, whisker curvature increases from caudal to rostral, and whiskers emerge from the face in smooth orientation gradients. Two aspects of whisker orientation are strikingly consistent across species: the elevation angle is constant within a row, and the twist of the whisker about its own axis varies smoothly in a diagonal gradient across the array. We suggest that invariant whisker elevation within a row may aid localization behaviors, while variable twist-orientation may help the animal sense stimulus direction. We anticipate this work will serve as a starting point for quantitative comparisons of vibrissal arrays across species, help clarify the mechanical basis by which seal vibrissae enable efficient wake detection and following, and enable the creation of whole-body biomechanical models for neuroscience and robotics.

Comparative morphology of the whiskers and faces of mice (Mus musculus) and rats (Rattus norvegicus)

Understanding neural function requires quantification of the sensory signals that an animal's brain evolved to interpret. These signals in turn depend on the morphology and mechanics of the animal's sensory structures. Although the house mouse (Mus musculus) is one of the most common model species used in neuroscience, the spatial arrangement of its facial sensors has not yet been quantified. To address this gap, the present study quantifies the facial morphology of the mouse, with a particular focus on the geometry of its vibrissae (whiskers). The study develops equations that establish relationships between the three-dimensional (3D) locations of whisker basepoints, whisker geometry (arclength, curvature) and the 3D angles at which the whiskers emerge from the face. Additionally, the positions of facial sensory organs are quantified relative to bregma-lambda. Comparisons with the Norway rat (Rattus norvegicus) indicate that when normalized for head size, the whiskers of these two species have similar spacing density. The rostral-caudal distances between facial landmarks of the rat are a factor of ∼2.0 greater than the mouse, while the scale of bilateral distances is larger and more variable. We interpret these data to suggest that the larger size of rats compared with mice is a derived (apomorphic) trait. As rodents are increasingly important models in behavioral neuroscience, the morphological model developed here will help researchers generate naturalistic, multimodal patterns of stimulation for neurophysiological experiments and allow the generation of synthetic datasets and simulations to close the loop between brain, body and environment.

Impaired trigeminal control of ingestive behavior in the Prrxl1-/- mouse is associated with a lemniscal-biased orosensory deafferentation

Although peripheral deafferentation studies have demonstrated a critical role for trigeminal afference in modulating the orosensorimotor control of eating and drinking, the central trigeminal pathways mediating that control, as well as the timescale of control, remain to be elucidated. In rodents, three ascending somatosensory pathways process and relay orofacial mechanosensory input: the lemniscal, paralemniscal, and extralemniscal. Two of these pathways (the lemniscal and extralemniscal) exhibit highly structured topographic representations of the orofacial sensory surface, as exemplified by the one-to-one somatotopic mapping between vibrissae on the animals’ face and barrelettes in brainstem, barreloids in thalamus, and barrels in cortex. Here we use the Prrxl1 knockout mouse model (also known as the DRG11 knockout) to investigate ingestive behavior deficits that may be associated with disruption of the lemniscal pathway. The Prrxl1 deletion disrupts somatotopic patterning and axonal projections throughout the lemniscal pathway but spares patterning in the extralemniscal nucleus. Our data reveal an imprecise and inefficient ingestive phenotype. Drinking behavior exhibits deficits on the timescales of milliseconds to seconds. Eating behavior shows deficits over an even broader range of timescales. An analysis of food acquisition and consummatory rate showed deficits on the timescale of seconds, and analysis of body weight suggested deficits on the scale of long term appetitive control. We suggest that ordered assembly of trigeminal sensory information along the lemniscal pathway is critical for the rapid and precise modulation of motor circuits driving eating and drinking action sequences.

Continuous, multidimensional coding of 3D complex tactile stimuli by primary sensory neurons of the vibrissal system

Across all sensory modalities, first-stage sensory neurons are an information bottleneck: they must convey all information available for an animal to perceive and act in its environment. Our understanding of coding properties of primary sensory neurons in the auditory and visual systems has been aided by the use of increasingly complex, naturalistic stimulus sets. By comparison, encoding properties of primary somatosensory afferents are poorly understood. Here, we use the rodent whisker system to examine how tactile information is represented in primary sensory neurons of the trigeminal ganglion (Vg). Vg neurons have long been thought to segregate into functional classes associated with separate streams of information processing. However, this view is based on Vg responses to restricted stimulus sets which potentially underreport the coding capabilities of these neurons. In contrast, the current study records Vg responses to complex three-dimensional (3D) stimulation while quantifying the complete 3D whisker shape and mechanics, thereby beginning to reveal their full representational capabilities. The results show that individual Vg neurons simultaneously represent multiple mechanical features of a stimulus, do not preferentially encode principal components of the stimuli, and represent continuous and tiled variations of all available mechanical information. These results directly contrast with proposed codes in which subpopulations of Vg neurons encode select stimulus features. Instead, individual Vg neurons likely overcome the information bottleneck by encoding large regions of a complex sensory space. This proposed tiled and multidimensional representation at the Vg directly constrains the computations performed by more central neurons of the vibrissotrigeminal pathway.

A dynamical model for generating synthetic data to quantify active tactile sensing behavior in the rat

As it becomes possible to simulate increasingly complex neural networks, it becomes correspondingly important to model the sensory information that animals actively acquire: the biomechanics of sensory acquisition directly determines the sensory input and therefore neural processing. Here, we exploit the tractable mechanics of the well-studied rodent vibrissal ("whisker") system to present a model that can simulate the signals acquired by a full sensor array actively sampling the environment. Rodents actively "whisk" ∼60 vibrissae (whiskers) to obtain tactile information, and this system is therefore ideal to study closed-loop sensorimotor processing. The simulation framework presented here, WHISKiT Physics, incorporates realistic morphology of the rat whisker array to predict the time-varying mechanical signals generated at each whisker base during sensory acquisition. Single-whisker dynamics were optimized based on experimental data and then validated against free tip oscillations and dynamic responses to collisions. The model is then extrapolated to include all whiskers in the array, incorporating each whisker's individual geometry. Simulation examples in laboratory and natural environments demonstrate that WHISKiT Physics can predict input signals during various behaviors, currently impossible in the biological animal. In one exemplary use of the model, the results suggest that active whisking increases in-plane whisker bending compared to passive stimulation and that principal component analysis can reveal the relative contributions of whisker identity and mechanics at each whisker base to the vibrissotactile response. These results highlight how interactions between array morphology and individual whisker geometry and dynamics shape the signals that the brain must process.

Constraints on the deformation of the vibrissa within the follicle

Nearly all mammals have a vibrissal system specialized for tactile sensation, composed of whiskers growing from sensor-rich follicles in the skin. When a whisker deflects against an object, it deforms within the follicle and exerts forces on the mechanoreceptors inside. In addition, during active whisking behavior, muscle contractions around the follicle and increases in blood pressure in the ring sinus will affect the whisker deformation profile. To date, however, it is not yet possible to experimentally measure how the whisker deforms in an intact follicle or its effects on different groups of mechanoreceptors. The present study develops a novel model to predict vibrissal deformation within the follicle sinus complex. The model is based on experimental results from a previous ex vivo study on whisker deformation within the follicle, and on a new histological analysis of follicle tissue. It is then used to simulate whisker deformation within the follicle during passive touch and active whisking. Results suggest that the most likely whisker deformation profile is “S-shaped,” crossing the midline of the follicle right below the ring sinus. Simulations of active whisking indicate that an increase in overall muscle stiffness, an increase in the ratio between deep and superficial intrinsic muscle stiffness, and an increase in sinus blood pressure will all enhance tactile sensitivity. Finally, we discuss how the deformation profiles might map to the responses of primary afferents of each mechanoreceptor type. The mechanical model presented in this study is an important first step in simulating mechanical interactions within whisker follicles.

Many mammals rely on whiskers as a mode of tactile sensation, especially when exploring in darkness. Active, rhythmic protraction and retraction of the whiskers, commonly referred to as “whisking,” is observed among many whisker specialist animals. During whisker-based sensing, forces and moments generated by external stimuli are transmitted to the base of the whisker shaft inside the follicle. Within the follicle, the interaction between the whisker’s deformation and the surrounding tissue determines how different groups of mechanoreceptors will deform, thereby transducing the mechanical signals into electrical signals. However, it is not yet possible to experimentally measure this interaction in vivo. We therefore created a mechanical model of the follicle sinus complex to simulate whisker deformation within the follicle resulting from external whisker deflection. Our results provide the first estimate of whisker shape as it deforms in the follicle, during both passive touch and active whisking. In turn, these shape estimates allow us to predict how the whisker will deform against different types of mechanoreceptors at different locations within the follicle. In addition, we find that both intrinsic muscle contraction and an increase in blood pressure will improve the tactile sensitivity of the whisker system.

Defining "active sensing" through an analysis of sensing energetics: homeoactive and alloactive sensing

The term "active sensing" has been defined in multiple ways. Most strictly, the term refers to sensing that uses self-generated energy to sample the environment (e.g., echolocation). More broadly, the definition includes all sensing that occurs when the sensor is moving (e.g., tactile stimuli obtained by an immobile versus moving fingertip) and, broader still, includes all sensing guided by attention or intent (e.g., purposeful eye movements). The present work offers a framework to help disambiguate aspects of the "active sensing" terminology and reveals properties of tactile sensing unique among all modalities. The framework begins with the well-described "sensorimotor loop," which expresses the perceptual process as a cycle involving four subsystems: environment, sensor, nervous system, and actuator. Using system dynamics, we examine how information flows through the loop. This "sensory-energetic loop" reveals two distinct sensing mechanisms that subdivide active sensing into homeoactive and alloactive sensing. In homeoactive sensing, the animal can change the state of the environment, while in alloactive sensing the animal can alter only the sensor's configurational parameters and thus the mapping between input and output. Given these new definitions, examination of the sensory-energetic loop helps identify two unique characteristics of tactile sensing: 1) in tactile systems, alloactive and homeoactive sensing merge to a mutually controlled sensing mechanism, and 2) tactile sensing may require fundamentally different predictions to anticipate reafferent input. We expect this framework may help resolve ambiguities in the active sensing community and form a basis for future theoretical and experimental work regarding alloactive and homeoactive sensing.

Tapered Polymer Whiskers to Enable Three-Dimensional Tactile Feature Extraction

Many mammals use their vibrissae (whiskers) to tactually explore their surrounding environment. Vibrissae are thin tapered structures that transmit mechanical signals to a wealth of mechanical receptors (sensors) located in a follicle at each vibrissal base. A recent study has shown that-provided that the whisker is tapered-three mechanical signals at the base are sufficient to determine the three-dimensional location at which a whisker made contact with an object. However, creating biomimetic tapered whiskers has proved challenging from both materials and manufacturing standpoints. This study develops and characterizes an artificial whisker for use as part of a sensory input device that is a biomimic of the biological rat whisker neurosensory system. A novel manufacturing process termed surface conforming fiber drawing (SCFD) is developed to produce artificial whiskers that meet the requirements to be a successful mechanical and geometric mimic of the biological rat vibrissae. Testing the sensory capabilities of the artificial whisker shows improved performance over previous nontapered filaments. SCFD-manufactured tapered whiskers demonstrate the ability to predict contact point locations with a median distance error of 0.47 cm.

The Cellular and Mechanical Basis for Response Characteristics of Identified Primary Afferents in the Rat Vibrissal System

Compared to our understanding of the response properties of receptors in the auditory and visual systems, we have only a limited understanding of the mechanoreceptor responses that underlie tactile sensation. Here, we exploit the stereotyped morphology of the rat vibrissal (whisker) array to investigate coding and transduction properties of identified primary tactile afferents. We performed in vivo intra-axonal recording and labeling experiments to quantify response characteristics of four different types of identified mechanoreceptors in the vibrissal follicle: ring-sinus Merkel; lanceolate; clublike; and rete-ridge collar Merkel. Of these types, only ring-sinus Merkel endings exhibited slowly adapting properties. A weak inverse relationship between response magnitude and onset response latency was found across all types. All afferents exhibited strong "angular tuning," i.e., their response magnitude and latency depended on the whisker's deflection angle. Although previous studies suggested that this tuning should be aligned with the angular location of the mechanoreceptor in the follicle, such alignment was observed only for Merkel afferents; angular tuning of the other afferent types showed no clear alignment with mechanoreceptor location. Biomechanical modeling suggested that this tuning difference might be explained by mechanoreceptors' differential sensitivity to the force directed along the whisker length. Electron microscopic investigations of Merkel endings and lanceolate endings at the level of the ring sinus revealed unique anatomical features that may promote these differential sensitivities. The present study systematically integrates biomechanical principles with the anatomical and morphological characterization of primary afferent endings to describe the physical and cellular processing that shapes the neural representation of touch.

Quantification of vibrissal mechanical properties across the rat mystacial pad

Recent work has quantified the geometric parameters of individual rat vibrissae (whiskers) and developed equations that describe how these parameters vary as a function of row and column position across the array. This characterization included a detailed quantification of whisker base diameter and arc length as well as the geometry of the whisker medulla. The present study now uses these equations for whisker geometry to quantify several properties of the whisker that govern its mechanical behavior. We first show that the average density of a whisker is lower in its proximal region than in its distal region. This density variation appears to be largely attributable to the presence of the whisker cuticle rather than the medulla. The density variation has very little effect on the center of mass of the whisker. We next show that the presence of the medulla decreases the deflection of the whisker under its own weight and also decreases its mass moment of inertia while sacrificing <1% stiffness at the whisker base compared with a solid whisker. Finally, we quantify two dimensionless parameters across the array. First, the deflection-to-length ratio decreases from caudal to rostral: caudal whiskers are longer but deflect more under their own weight. Second, the nondimensionalized radius of gyration is approximately constant across the array, which may simplify control of whisking by the intrinsic muscles. We anticipate that future work will exploit the mechanical properties computed in the present study to improve simulations of the mechanosensory signals associated with vibrissotactile exploratory behavior. NEW & NOTEWORTHY The mechanical signals transmitted by a whisker depend critically on its geometry. We used measurements of whisker geometry and mass to quantify the center of mass, mass moment of inertia, radius of gyration, and deflection under gravity of the whisker. We describe how variations in these quantities across the array could enhance sensing behaviors while reducing energy costs and simplifying whisking control. Most importantly, we provide derivations for these quantities for use in future simulation work.

Tactile Sensing with Whiskers of Various Shapes: Determining the Three-Dimensional Location of Object Contact Based on Mechanical Signals at the Whisker Base

Almost all mammals use their mystacial vibrissae (whiskers) as important tactile sensors. There are no sensors along the length of a whisker: all sensing is performed by mechanoreceptors at the whisker base. To use artificial whiskers as a sensing tool in robotics, it is essential to be able to determine the three-dimensional (3D) location at which a whisker has made contact with an object. With the assumption of quasistatic, frictionless, single-point contact, previous work demonstrated that the 3D contact point can be uniquely determined if all six components of force and moment are measured at the whisker base, but these measurements require a six-axis load cell. Here, we perform simulations to investigate the extent to which each of the 20 possible "triplet" combinations of the six mechanical signals at the whisker base uniquely determine 3D contact point location. We perform this analysis for four different whisker profiles (shapes): tapered with and without intrinsic curvature, and cylindrical with and without intrinsic curvature. We show that whisker profile has a strong influence on the particular triplet(s) of signals that uniquely map to the 3D contact point. The triplet of bending moment, bending moment direction, and axial force produces unique mappings for tapered whiskers. Four different mappings are unique for a cylindrical whisker without intrinsic curvature, but only when large deflections are excluded. These results inform the neuroscience of vibrissotactile sensing and represent an important step toward the development of artificial whiskers for robotic applications.

Ergodic Exploration Using Binary Sensing for Nonparametric Shape Estimation

Current methods to estimate object shape-using either vision or touch-generally depend on high-resolution sensing. Here, we exploit ergodic exploration to demonstrate successful shape estimation when using a low-resolution binary contact sensor. The measurement model is posed as a collision-based tactile measurement, and classification methods are used to discriminate between shape boundary regions in the search space. Posterior likelihood estimates of the measurement model help the system actively seek out regions where the binary sensor is most likely to return informative measurements. Results show successful shape estimation of various objects as well as the ability to identify multiple objects in an environment. Interestingly, it is shown that ergodic exploration utilizes non-contact motion to gather significant information about shape. The algorithm is extended in three dimensions in simulation and we present two dimensional experimental results using the Rethink Baxter robot.

Variations in vibrissal geometry across the rat mystacial pad: base diameter, medulla, and taper

Many rodents tactually sense the world through active motions of their vibrissae (whiskers), which are regularly arranged in rows and columns (arcs) on the face. The present study quantifies several geometric parameters of rat whiskers that determine the tactile information acquired. Findings include the following. 1) A meta-analysis of seven studies shows that whisker base diameter varies with arc length with a surprisingly strong dependence on the whisker's row position within the array. 2) The length of the whisker medulla varies linearly with whisker length, and the medulla's base diameter varies linearly with whisker base diameter. 3) Two parameters are required to characterize whisker "taper": radius ratio (base radius divided by tip radius) and radius slope (the difference between base and tip radius, divided by arc length). A meta-analysis of five studies shows that radius ratio exhibits large variability due to variations in tip radius, while radius slope varies systematically across the array. 4) Within the resolution of the present study, radius slope does not differ between the proximal and distal segments of the whisker, where "proximal" is defined by the presence of the medulla. 5) Radius slope of the medulla is offset by a constant value from radius slope of the proximal portion of the whisker. We conclude with equations for all geometric parameters as functions of row and column position.NEW & NOTEWORTHY Rats tactually explore their world by brushing and tapping their whiskers against objects. Each whisker's geometry will have a large influence on its mechanics and thus on the tactile signals the rat obtains. We performed a meta-analysis of seven studies to generate equations that describe systematic variations in whisker geometry across the rat's face. We also quantified the geometry of the whisker medulla. A database provides access to geometric parameters of over 500 rat whiskers.

Probability distributions of whisker-surface contact: quantifying elements of the rat vibrissotactile natural scene

Analysis of natural scene statistics has been a powerful approach for understanding neural coding in the auditory and visual systems. In the field of somatosensation, it has been more challenging to quantify the natural tactile scene, in part because somatosensory signals are so tightly linked to the animal's movements. The present work takes a step towards quantifying the natural tactile scene for the rat vibrissal system by simulating rat whisking motions to systematically investigate the probabilities of whisker-object contact in naturalistic environments. The simulations permit an exhaustive search through the complete space of possible contact patterns, thereby allowing for the characterization of the patterns that would most likely occur during long sequences of natural exploratory behavior. We specifically quantified the probabilities of 'concomitant contact', that is, given that a particular whisker makes contact with a surface during a whisk, what is the probability that each of the other whiskers will also make contact with the surface during that whisk? Probabilities of concomitant contact were quantified in simulations that assumed increasingly naturalistic conditions: first, the space of all possible head poses; second, the space of behaviorally preferred head poses as measured experimentally; and third, common head poses in environments such as cages and burrows. As environments became more naturalistic, the probability distributions shifted from exhibiting a 'row-wise' structure to a more diagonal structure. Results also reveal that the rat appears to use motor strategies (e.g. head pitches) that generate contact patterns that are particularly well suited to extract information in the presence of uncertainty.

Whiskers aid anemotaxis in rats

Observation of terrestrial mammals suggests that they can follow the wind (anemotaxis), but the sensory cues underlying this ability have not been studied. We identify a significant contribution to anemotaxis mediated by whiskers (vibrissae), a modality previously studied only in the context of direct tactile contact. Five rats trained on a five-alternative forced-choice airflow localization task exhibited significant performance decrements after vibrissal removal. In contrast, vibrissal removal did not disrupt the performance of control animals trained to localize a light source. The performance decrement of individual rats was related to their airspeed threshold for successful localization: animals that found the task more challenging relied more on the vibrissae for localization cues. Following vibrissal removal, the rats deviated more from the straight-line path to the air source, choosing sources farther from the correct location. Our results indicate that rats can perform anemotaxis and that whiskers greatly facilitate this ability. Because air currents carry information about both odor content and location, these findings are discussed in terms of the adaptive significance of the interaction between sniffing and whisking in rodents.

Representation of Stimulus Speed and Direction in Vibrissal-Sensitive Regions of the Trigeminal Nuclei: A Comparison of Single Unit and Population Responses

The rat vibrissal (whisker) system is one of the oldest and most important models for the study of active tactile sensing and sensorimotor integration. It is well established that primary sensory neurons in the trigeminal ganglion respond to deflections of one and only one whisker, and that these neurons are strongly tuned for both the speed and direction of individual whisker deflections. During active whisking behavior, however, multiple whiskers will be deflected simultaneously. Very little is known about how neurons at central levels of the trigeminal pathway integrate direction and speed information across multiple whiskers. In the present work, we investigated speed and direction coding in the trigeminal brainstem nuclei, the first stage of neural processing that exhibits multi-whisker receptive fields. Specifically, we recorded both single-unit spikes and local field potentials from fifteen sites in spinal trigeminal nucleus interpolaris and oralis while systematically varying the speed and direction of coherent whisker deflections delivered across the whisker array. For 12/15 neurons, spike rate was higher when the whisker array was stimulated from caudal to rostral rather than rostral to caudal. In addition, 10/15 neurons exhibited higher firing rates for slower stimulus speeds. Interestingly, using a simple decoding strategy for the local field potentials and spike trains, classification of speed and direction was higher for field potentials than for single unit spike trains, suggesting that the field potential is a robust reflection of population activity. Taken together, these results point to the idea that population responses in these brainstem regions in the awake animal will be strongest during behaviors that stimulate a population of whiskers with a directionally coherent motion.

Whisking Kinematics Enables Object Localization in Head-Centered Coordinates Based on Tactile Information from a Single Vibrissa

During active tactile exploration with their whiskers (vibrissae), rodents can rapidly orient to an object even though there are very few proprioceptors in the whisker muscles. Thus a long-standing question in the study of the vibrissal system is how the rat can localize an object in head-centered coordinates without muscle-based proprioception. We used a three-dimensional model of whisker bending to simulate whisking motions against a peg to investigate the possibility that the 3D mechanics of contact from a single whisker are sufficient for localization in head-centered coordinates. Results show that for nearly all whiskers in the array, purely tactile signals at the whisker base - as would be measured by mechanoreceptors, in whisker-centered coordinates - could be used to determine the location of a vertical peg in head-centered coordinates. Both the "roll" and the "elevation" components of whisking kinematics contribute to the uniqueness and resolution of the localization. These results offer an explanation for a behavioral study showing that rats can more accurately determine the horizontal angle of an object if one column, rather than one row, of whiskers is spared.

Mechanical responses of rat vibrissae to airflow

The survival of many animals depends in part on their ability to sense the flow of the surrounding fluid medium. To date, however, little is known about how terrestrial mammals sense airflow direction or speed. The present work analyzes the mechanical response of isolated rat macrovibrissae (whiskers) to airflow to assess their viability as flow sensors. Results show that the whisker bends primarily in the direction of airflow and vibrates around a new average position at frequencies related to its resonant modes. The bending direction is not affected by airflow speed or by geometric properties of the whisker. In contrast, the bending magnitude increases strongly with airflow speed and with the ratio of the whisker's arc length to base diameter. To a much smaller degree, the bending magnitude also varies with the orientation of the whisker's intrinsic curvature relative to the direction of airflow. These results are used to predict the mechanical responses of vibrissae to airflow across the entire array, and to show that the rat could actively adjust the airflow data that the vibrissae acquire by changing the orientation of its whiskers. We suggest that, like the whiskers of pinnipeds, the macrovibrissae of terrestrial mammals are multimodal sensors - able to sense both airflow and touch - and that they may play a particularly important role in anemotaxis.

Simulations of a Vibrissa Slipping along a Straight Edge and an Analysis of Frictional Effects during Whisking

During tactile exploration, rats sweep their whiskers against objects in a motion called whisking. Here, we investigate how a whisker slips along an object's edge and how friction affects the resulting tactile signals. First, a frictionless model is developed to simulate whisker slip along a straight edge and compared with a previous model that incorporates friction but cannot simulate slip. Results of both models are compared to behavioral data obtained as a rat whisked against a smooth, stainless steel peg. As expected, the frictionless model predicts larger magnitudes of vertical slip than observed experimentally. The frictionless model also predicts forces and moments at the whisker base that are smaller and have a different direction than those predicted by the model with friction. Estimates for the friction coefficient yielded values near 0.48 (whisker/stainless steel). The present work provides the first assessments of the effects of friction on the mechanical signals received by the follicle during active whisking. It also demonstrates a proof-of-principle approach for reducing whisker tracking requirements during experiments and demonstrates the feasibility of simulating a full array of vibrissae whisking against a peg.

Evidence for Functional Groupings of Vibrissae across the Rodent Mystacial Pad

During natural exploration, rats exhibit two particularly conspicuous vibrissal-mediated behaviors: they follow along walls, and “dab” their snouts on the ground at frequencies related to the whisking cycle. In general, the walls and ground may be located at any distance from, and at any orientation relative to, the rat’s head, which raises the question of how the rat might determine the position and orientation of these surfaces. Previous studies have compellingly demonstrated that rats can accurately determine the horizontal angle at which a vibrissa first touches an object, and we therefore asked whether this parameter could provide the rat with information about the pitch, distance, and yaw of a surface relative to its head. We used a three-dimensional model of the whisker array to construct mappings between the horizontal angle of contact of each vibrissa and every possible (pitch, distance, and yaw) configuration of the head relative to a flat surface. The mappings revealed striking differences in the patterns of contact for vibrissae in different regions of the array. The exterior (A, D, E) rows provide information about the relative pitch of the surface regardless of distance. The interior (B, C) rows provide distance cues regardless of head pitch. Yaw is linearly correlated with the difference between the number of right and left whiskers touching the surface. Compared to the long reaches that whiskers can make to the side and below the rat, the reachable distance in front of the rat’s nose is relatively small. We confirmed key predictions of these functional groupings in a behavioral experiment that monitored the contact patterns that the vibrissae made with a flat vertical surface. These results suggest that vibrissae in different regions of the array are not interchangeable sensors, but rather functionally grouped to acquire particular types of information about the environment.

Animals do not passively sense the world. They actively move their sensors to acquire the information best suited to their current behavioral needs. Thus, to study neural processing in any sensory system, it is critical to determine the spatiotemporal structure of the inputs to that system. The rat whisker system is a well-studied model of active sensing because rats rhythmically brush and tap their vibrissae against objects during tactual exploration. To date, however, the patterns of sensory input during vibrisso-tactile exploration have not been quantified. This study quantifies some of the statistics of vibrissal-surface contact, focusing specifically on the angle of contact with a flat surface. Given that during exploration the rat could encounter a surface at any position and orientation relative to its head, we simulated the spatial patterns of vibrissal contact for all possible head-surface configurations. Results reveal functional groupings of whiskers across the array, with different whiskers being better suited for different aspects of tactile sensory tasks than others. Key predictions of the simulations were validated in a behavioral experiment that monitored vibrissal contact patterns with a flat wall. These results provide some of the first quantitative insights into the “natural tactile scene” for the vibrisso-trigeminal system.

Spatiotemporal Patterns of Contact Across the Rat Vibrissal Array During Exploratory Behavior

The rat vibrissal system is an important model for the study of somatosensation, but the small size and rapid speed of the vibrissae have precluded measuring precise vibrissal-object contact sequences during behavior. We used a laser light sheet to quantify, with 1 ms resolution, the spatiotemporal structure of whisker-surface contact as five naïve rats freely explored a flat, vertical glass wall. Consistent with previous work, we show that the whisk cycle cannot be uniquely defined because different whiskers often move asynchronously, but that quasi-periodic (~8 Hz) variations in head velocity represent a distinct temporal feature on which to lock analysis. Around times of minimum head velocity, whiskers protract to make contact with the surface, and then sustain contact with the surface for extended durations (~25-60 ms) before detaching. This behavior results in discrete temporal windows in which large numbers of whiskers are in contact with the surface. These "sustained collective contact intervals" (SCCIs) were observed on 100% of whisks for all five rats. The overall spatiotemporal structure of the SCCIs can be qualitatively predicted based on information about head pose and the average whisk cycle. In contrast, precise sequences of whisker-surface contact depend on detailed head and whisker kinematics. Sequences of vibrissal contact were highly variable, equally likely to propagate in all directions across the array. Somewhat more structure was found when sequences of contacts were examined on a row-wise basis. In striking contrast to the high variability associated with contact sequences, a consistent feature of each SCCI was that the contact locations of the whiskers on the glass converged and moved more slowly on the sheet. Together, these findings lead us to propose that the rat uses a strategy of "windowed sampling" to extract an object's spatial features: specifically, the rat spatially integrates quasi-static mechanical signals across whiskers during the period of sustained contact, resembling an "enclosing" haptic procedure.

The search space of the rat during whisking behavior

Rodents move their vibrissae rhythmically to tactually explore their surroundings. We used a three-dimensional model of the vibrissal array to quantify the rat's 'search space' during whisking. Search space was quantified either as the volume encompassed by the array or as the surface formed by the vibrissal tips. At rest, the average position of the vibrissal tips lies near the rat's mouth, and the tips are all approximately equidistant from the midpoint between the rat's eyes, suggesting spatial registration with the visual system. The intrinsic curvature of the vibrissae greatly increases the volume encompassed by the array, and during a protraction, roll and elevation changes have strong effects on the trajectories of the vibrissal tips. The size of the rat's search space--as measured either by the volume of the array or by the surface area formed by the vibrissal tips--was surprisingly unaffected by protraction angle. In contrast, search space was strongly correlated with the 'spread' of the array, defined as the angle between rostral and caudal-most whiskers. We draw two conclusions: first, that with some caveats, spread can be used as a proxy for changes in search space, and second, in order to change its sensing resolution, the rat must differentially control rostral and caudal vibrissae. Finally, we show that behavioral data can be incorporated into the three-dimensional model to visualize changes in vibrissal search space and sensing resolution during natural exploratory whisking.

Tactile signals transmitted by the vibrissa during active whisking behavior

The rodent vibrissal-trigeminal system is one of the most widely used models for the study of somatosensation and tactile perception, but to date the field has been unable to quantify the complete set of mechanical input signals generated during natural whisking behavior. In this report we show that during whisking behavior of awake rats (Rattus norvegicus), the whisker will often bend out of its plane of rotation, generating sizeable mechanical (tactile) signals out of the plane. We then develop a model of whisker bending that allows us to compute the three-dimensional tactile signals at the vibrissal base during active whisking behavior. Considerable information can be lost if whisking motions are considered only in two dimensions, and we offer some suggestions for experimentalists concerned with monitoring the direction of bending. These data represent the first quantification of the physical signals transmitted to the mechanoreceptors in the follicle during active whisking behavior.

Sensory prediction on a whiskered robot: a tactile analogy to "optical flow"

When an animal moves an array of sensors (e.g., the hand, the eye) through the environment, spatial and temporal gradients of sensory data are related by the velocity of the moving sensory array. In vision, the relationship between spatial and temporal brightness gradients is quantified in the “optical flow” equation. In the present work, we suggest an analog to optical flow for the rodent vibrissal (whisker) array, in which the perceptual intensity that “flows” over the array is bending moment. Changes in bending moment are directly related to radial object distance, defined as the distance between the base of a whisker and the point of contact with the object. Using both simulations and a 1×5 array (row) of artificial whiskers, we demonstrate that local object curvature can be estimated based on differences in radial distance across the array. We then develop two algorithms, both based on tactile flow, to predict the future contact points that will be obtained as the whisker array translates along the object. The translation of the robotic whisker array represents the rat's head velocity. The first algorithm uses a calculation of the local object slope, while the second uses a calculation of the local object curvature. Both algorithms successfully predict future contact points for simple surfaces. The algorithm based on curvature was found to more accurately predict future contact points as surfaces became more irregular. We quantify the inter-related effects of whisker spacing and the object's spatial frequencies, and examine the issues that arise in the presence of real-world noise, friction, and slip.

Radial distance determination in the rat vibrissal system and the effects of Weber's law

Rats rhythmically tap and brush their vibrissae (whiskers) against objects to tactually explore the environment. To extract a complex feature such as the contour of an object, the rat must at least implicitly estimate radial object distance, that is, the distance from the base of the vibrissa to the point of object contact. Radial object distance cannot be directly measured, however, because there are no mechanoreceptors along the vibrissa. Instead, the mechanical signals generated by the vibrissa's interaction with the environment must be transmitted to mechanoreceptors near the vibrissa base. The first part of this paper surveys the different mechanical methods by which the rat could determine radial object distance. Two novel methods are highlighted: one based on measurement of bending moment and axial force at the vibrissa base, and a second based on measurement of how far the vibrissa rotates beyond initial contact. The second part of the paper discusses the application of Weber's law to two methods for radial distance determination. In both cases, Weber's law predicts that the rat will have greatest sensing resolution close to the vibrissa tip. These predictions could be tested with behavioural experiments that measure the perceptual acuity of the rat.

Mechanical signals at the base of a rat vibrissa: the effect of intrinsic vibrissa curvature and implications for tactile exploration

Rats actively tap and sweep their large mystacial vibrissae (whiskers) against objects to tactually explore their surroundings. When a vibrissa makes contact with an object, it bends, and this bending generates forces and bending moments at the vibrissa base. Researchers have only recently begun to quantify these mechanical variables. The present study quantifies the forces and bending moments at the vibrissa base with a quasi-static model of vibrissa deflection. The model was validated with experiments on real vibrissae. Initial simulations demonstrated that almost all vibrissa-object collisions during natural behavior will occur with the concave side of the vibrissa facing the object, and we therefore paid particular attention to the role of the vibrissa's intrinsic curvature in shaping the forces at the base. Both simulations and experiments showed that vibrissae with larger intrinsic curvatures will generate larger axial forces. Simulations also demonstrated that the range of forces and moments at the vibrissal base vary over approximately three orders of magnitude, depending on the location along the vibrissa at which object contact is made. Both simulations and experiments demonstrated that collisions in which the concave side of the vibrissa faces the object generate longer-duration contacts and larger net forces than collisions with the convex side. These results suggest that the orientation of the vibrissa's intrinsic curvature on the mystacial pad may increase forces during object contact and provide increased sensitivity to detailed surface features.

A night in the life of a rat: vibrissal mechanics and tactile exploration

The rat vibrissal (whisker) system is an increasingly important model for the study of the sense of touch. This paper describes recent results obtained from high-speed videography of rat exploratory behavior and from modeling studies of vibrissal biomechanics. We review several features of vibrissal touch, including the mechanics of contact versus noncontact whisking, the coordination between head and whisker movements, and the use of information obtained from the whiskers to detect, localize, and extract the spatial properties of objects. This work highlights the idea that mechanics are critical to an understanding of sensory systems and describes some new tools to monitor the spatiotemporal patterns of whisker-object contact during natural tactile exploratory behavior.

The morphology of the rat vibrissal array: a model for quantifying spatiotemporal patterns of whisker-object contact

In all sensory modalities, the data acquired by the nervous system is shaped by the biomechanics, material properties, and the morphology of the peripheral sensory organs. The rat vibrissal (whisker) system is one of the premier models in neuroscience to study the relationship between physical embodiment of the sensor array and the neural circuits underlying perception. To date, however, the three-dimensional morphology of the vibrissal array has not been characterized. Quantifying array morphology is important because it directly constrains the mechanosensory inputs that will be generated during behavior. These inputs in turn shape all subsequent neural processing in the vibrissal-trigeminal system, from the trigeminal ganglion to primary somatosensory ("barrel") cortex. Here we develop a set of equations for the morphology of the vibrissal array that accurately describes the location of every point on every whisker to within ±5% of the whisker length. Given only a whisker's identity (row and column location within the array), the equations establish the whisker's two-dimensional (2D) shape as well as three-dimensional (3D) position and orientation. The equations were developed via parameterization of 2D and 3D scans of six rat vibrissal arrays, and the parameters were specifically chosen to be consistent with those commonly measured in behavioral studies. The final morphological model was used to simulate the contact patterns that would be generated as a rat uses its whiskers to tactually explore objects with varying curvatures. The simulations demonstrate that altering the morphology of the array changes the relationship between the sensory signals acquired and the curvature of the object. The morphology of the vibrissal array thus directly constrains the nature of the neural computations that can be associated with extraction of a particular object feature. These results illustrate the key role that the physical embodiment of the sensor array plays in the sensing process.

Active touch, exploratory movements, and sensory prediction

The relation between somatosensory input and motor output is asymmetric. Somatosensation is associated with every movement an animal makes, but movement is not required for somatosensation. This symposium paper proposes a classification scheme for movement, in which movements are placed along a continuum that describes the role that somatosensory information plays during the movement. Fine sensorimotor control-manipulation and exploration-are found to fall to one extreme of the spectrum, and exploratory movements in particular are shown to possess characteristics that clearly distinguish them from other varieties of movement. Specifically, the exploratory process must permit animals to extract an object's features independently of the sequence of movements executed to explore the object. Based in part on our work on the rat vibrissal system, we suggest that exploration of objects may consist of two complementary levels of sensorimotor prediction operating in parallel. At the cognitive level, the animal might move so as to perform hypothesis testing about the identity or nature of the object. The particular hypothesis tests chosen by the animal might be implemented through sequences of control-level predictions that could be generated at the level of the brainstem and cerebellum.

Extracting Object Contours with the Sweep of a Robotic Whisker Using Torque Information

Several recent studies have investigated the problem of object feature extraction with artificial whiskers. Many of these studies have used an approach in which the whisker is rotated against the object through a small angle. Each small-angle “tap” of the whisker provides information about the radial distance between the base of the whisker and the object. By tapping at various points on the object, a full representation of the surface can be gradually constructed in three-dimensional space. It is clear, however, that this tapping method does not exploit useful information about object contours that could be extracted by “sweeping” the whisker against the object. Rotating the whisker against the object through a large angle permits the collection of a sequence of contact points as the whisker slips along the surface. The present paper derives an algorithm based on a numerical cantilever beam model of the whisker to measure object profile shape over a single large-angle whisker rotation using only information about torque and angle at the whisker base. The algorithm is validated experimentally using three different object shapes. As the method does not require measurement of force, it is well suited for implementation on an array of robotic whiskers to accomplish quick and precise object feature extraction.

Variability in Velocity Profiles During Free-Air Whisking Behavior of Unrestrained Rats

During exploratory behaviors, the velocity of an organism's sensory surfaces can have a pronounced effect on the incoming flow of sensory information. In this study, we quantified variability in the velocity profiles of rat whisking during natural exploratory behavior that included head rotations. A wide continuum of profiles was observed, including monotonic, delayed, and reversing velocities during protractions and retractions. Three alternative hypotheses for the function of the variable velocity profiles were tested: 1) that they produce bilateral asymmetry specifically correlated with rotational head velocity, 2) that they serve to generate bilaterally asymmetric and/or asynchronous whisker movements independent of head velocity, and 3) that the different profiles—despite increasing variability in instantaneous velocity—reduce variability in the average whisking velocity. Our results favor the third hypothesis and do not support the first two. Specifically, the velocity variability within a whisk can be observed as a shift in the phase of the maximum velocity. We discuss the implications of these results for the control of whisker motion, horizontal object localization, and processing in the thalamus and cortex of the rat vibrissal system.

A two-dimensional force sensor in the millinewton range for measuring vibrissal contacts

The rat vibrissal (whisker) array is a common model system in neuroscience used to study sensorimotor integration. Recent work has suggested that during object contact, the forces and moments at the whisker base may serve as important perceptual cues to the rat. To date, however, the force/moment profile that results from a whisker sweeping against an object has yet to be characterized, because it requires the simultaneous measurement of two-dimensional forces on the order of millinewtons. Current technology for these measurements typically involves prohibitively bulky, expensive equipment with complicated fabrication techniques. We have developed a simple, yet effective two-dimensional force sensor with +/-0.02 mN resolution; it is extremely compact, has a highly linear static response with low-noise output, and is inexpensive to build. We demonstrate the advantages and limitations of the sensor in three different experimental protocols, ranging from the precise quantification of forces on isolated (plucked) whiskers, to the detection of whisker-contact times in the awake behaving animal. Given the high fidelity of the sensor, it could have utility in a broad range of applications in which measuring contact/detach occurrence and/or small magnitude forces are important.

Using hardware models to quantify sensory data acquisition across the rat vibrissal array

Our laboratory investigates how animals acquire sensory data to understand the neural computations that permit complex sensorimotor behaviors. We use the rat whisker system as a model to study active tactile sensing; our aim is to quantitatively describe the spatiotemporal structure of incoming sensory information to place constraints on subsequent neural encoding and processing. In the first part of this paper we describe the steps in the development of a hardware model (a 'sensobot') of the rat whisker array that can perform object feature extraction. We show how this model provides insights into the neurophysiology and behavior of the real animal. In the second part of this paper, we suggest that sensory data acquisition across the whisker array can be quantified using the complete derivative. We use the example of wall-following behavior to illustrate that computing the appropriate spatial gradients across a sensor array would enable an animal or mobile robot to predict the sensory data that will be acquired at the next time step.

Biomechanical models for radial distance determination by the rat vibrissal system

Rats use active, rhythmic movements of their whiskers to acquire tactile information about three-dimensional object features. There are no receptors along the length of the whisker; therefore all tactile information must be mechanically transduced back to receptors at the whisker base. This raises the question: how might the rat determine the radial contact position of an object along the whisker? We developed two complementary biomechanical models that show that the rat could determine radial object distance by monitoring the rate of change of moment (or equivalently, the rate of change of curvature) at the whisker base. The first model is used to explore the effects of taper and inherent whisker curvature on whisker deformation and used to predict the shapes of real rat whiskers during deflections at different radial distances. Predicted shapes closely matched experimental measurements. The second model describes the relationship between radial object distance and the rate of change of moment at the base of a tapered, inherently curved whisker. Together, these models can account for recent recordings showing that some trigeminal ganglion (Vg) neurons encode closer radial distances with increased firing rates. The models also suggest that four and only four physical variables at the whisker base -- angular position, angular velocity, moment, and rate of change of moment -- are needed to describe the dynamic state of a whisker. We interpret these results in the context of our evolving hypothesis that neural responses in Vg can be represented using a state-encoding scheme that includes combinations of these four variables.