Nature Neuroscience: doi: /nn.2283

Similar documents
Lateral organization & computation

Consider the following spike trains from two different neurons N1 and N2:

Exercises. Chapter 1. of τ approx that produces the most accurate estimate for this firing pattern.

Inventory of Supplemental Information

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

The homogeneous Poisson process

Frequency Adaptation and Bursting

Single-Trial Neural Correlates. of Arm Movement Preparation. Neuron, Volume 71. Supplemental Information

Temporal context calibrates interval timing

The Spike Response Model: A Framework to Predict Neuronal Spike Trains

Neural Spike Train Analysis 1: Introduction to Point Processes

Geophysical Research Letters. Supporting Information for

Internally generated preactivation of single neurons in human medial frontal cortex predicts volition

Collective Dynamics in Human and Monkey Sensorimotor Cortex: Predicting Single Neuron Spikes

Comparing integrate-and-fire models estimated using intracellular and extracellular data 1

Transformation of stimulus correlations by the retina

Neural Encoding: Firing Rates and Spike Statistics

CORRELATION TRANSFER FROM BASAL GANGLIA TO THALAMUS IN PARKINSON S DISEASE. by Pamela Reitsma. B.S., University of Maine, 2007

Evolution of the Average Synaptic Update Rule

Nature Methods: doi: /nmeth Supplementary Figure 1. In vitro screening of recombinant R-CaMP2 variants.

Phenomenological Models of Neurons!! Lecture 5!

+ + ( + ) = Linear recurrent networks. Simpler, much more amenable to analytic treatment E.g. by choosing

Mid Year Project Report: Statistical models of visual neurons

AT2 Neuromodeling: Problem set #3 SPIKE TRAINS

Lecture 7 Random Signal Analysis

Temporal whitening by power-law adaptation in neocortical neurons

Hypothesis testing (cont d)

System Identification of Adapting Neurons

Neuronal Dynamics: Computational Neuroscience of Single Neurons

Chapter 3 Functional architecture of posterior parietal cortex during tactile sensorimotor behavior (in preparation)

The Bayesian Brain. Robert Jacobs Department of Brain & Cognitive Sciences University of Rochester. May 11, 2017

Supporting Information

SUPPLEMENTAL INFORMATION. Real-time Measurement of Biomagnetic Vector Fields in Functional Syncytium Using Amorphous Metal

An Introductory Course in Computational Neuroscience

Describing Spike-Trains

The Multivariate Gaussian Distribution

The idiosyncratic nature of confidence

CHARACTERIZATION OF NONLINEAR NEURON RESPONSES

2.004 Dynamics and Control II Spring 2008

Fig. S5: Same as Fig. S3 for the combination Sample

Single-Compartment Neural Models

Adaptation in the Neural Code of the Retina

Comparison of receptive fields to polar and Cartesian stimuli computed with two kinds of models

Introduction to Biomedical Engineering

Neural Coding: Integrate-and-Fire Models of Single and Multi-Neuron Responses

Single Molecule Electrochemistry on a Porous Silica-Coated Electrode

Fitting a Stochastic Neural Network Model to Real Data

Computing with Inter-spike Interval Codes in Networks of Integrate and Fire Neurons

Auto-correlation of retinal ganglion cell mosaics shows hexagonal structure

Rice University Physics 332 LIFETIME OF THE MUON I. INTRODUCTION...2! II. MEASUREMENT PROCEDURES...3! III. ANALYSIS PROCEDURES...7!

When is an Integrate-and-fire Neuron like a Poisson Neuron?

Supplementary material: Orbitofrontal and striatal circuits dynamically encode the shift. between goal-directed and habitual actions

Physics 178/278 - David Kleinfeld - Winter 2017

Finding a Basis for the Neural State

Supporting Online Material for

Mathematical Tools for Neuroscience (NEU 314) Princeton University, Spring 2016 Jonathan Pillow. Homework 8: Logistic Regression & Information Theory

Response-Field Dynamics in the Auditory Pathway

Adaptation to a 'spatial-frequency doubled' stimulus

Event-related fmri. Christian Ruff. Laboratory for Social and Neural Systems Research Department of Economics University of Zurich

EXPERIMENT 5A RC Circuits

CS 5630/6630 Scientific Visualization. Elementary Plotting Techniques II

Introduction to Neural Networks U. Minn. Psy 5038 Spring, 1999 Daniel Kersten. Lecture 2a. The Neuron - overview of structure. From Anderson (1995)

Synchrony in Neural Systems: a very brief, biased, basic view

High-dimensional geometry of cortical population activity. Marius Pachitariu University College London

Fast neural network simulations with population density methods

Supplementary Figure 1 Change of the Tunnelling Transmission Coefficient from the Bulk to the Surface as a result of dopant ionization Colour-map of

Model neurons!!poisson neurons!

10/8/2014. The Multivariate Gaussian Distribution. Time-Series Plot of Tetrode Recordings. Recordings. Tetrode

Supplemental Information. Noise and Correlations. in Parallel Perceptual Decision Making. Thomas U. Otto and Pascal Mamassian

Spike Count Correlation Increases with Length of Time Interval in the Presence of Trial-to-Trial Variation

!) + log(t) # n i. The last two terms on the right hand side (RHS) are clearly independent of θ and can be

Motion Perception 1. PSY305 Lecture 12 JV Stone

Spike-Frequency Adaptation: Phenomenological Model and Experimental Tests

Comparison Figures from the New 22-Year Daily Eddy Dataset (January April 2015)

Real and Modeled Spike Trains: Where Do They Meet?

LIST of SUPPLEMENTARY MATERIALS

Nature Neuroscience: doi: /nn Supplementary Figure 1. Expression of GCaMP6s in LGN and their axons in V1.

When do Correlations Increase with Firing Rates? Abstract. Author Summary. Andrea K. Barreiro 1* and Cheng Ly 2

arxiv: v1 [q-bio.nc] 13 Feb 2018

What happens if one pulls on the spring? The spring exerts a restoring force which is proportional to the distance it is stretched, F = - k x (1)

What can spike train distances tell us about the neural code?

Discrete Time Rescaling Theorem: Determining Goodness of Fit for Statistical Models of Neural Spiking

Neural Data Analysis

SPIKE TRIGGERED APPROACHES. Odelia Schwartz Computational Neuroscience Course 2017

Spike-Frequency Adaptation of a Generalized Leaky Integrate-and-Fire Model Neuron

Surround effects on the shape of the temporal contrast-sensitivity function

Acceleration Feedback

The Mössbauer Effect

Introduction to Statistics and Error Analysis II

Recipes for the Linear Analysis of EEG and applications

Limulus. The Neural Code. Response of Visual Neurons 9/21/2011

Oscillatory Motion and Wave Motion

Mathematical Analysis of Bursting Electrical Activity in Nerve and Endocrine Cells

Those data on the time required for the pitch to traverse various distances were measured in

EXPERIMENT 07 TO STUDY DC RC CIRCUIT AND TRANSIENT PHENOMENA

Dam Monitoring With Tiltmeters

ECE 5670/6670 Lab 8. Torque Curves of Induction Motors. Objectives

Phase Locking. 1 of of 10. The PRC's amplitude determines which frequencies a neuron locks to. The PRC's slope determines if locking is stable

Robust regression and non-linear kernel methods for characterization of neuronal response functions from limited data

Data analysis methods in the neuroscience Zoltán Somogyvári

Transcription:

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S. Overall distribution of whisk cycle phases of touch events for all units with rapid response to touch. Plot of the relative frequency of touch events at different phases of the whisk cycle. The blue horizontal line indicates the mean probability of touch over all phase intervals. The shaded region around this line indicates 95 % confidence limits for a uniform distribution. Although touch events near retraction are significantly less frequent, a sufficient number of these events occur in each session to produce stable touch response estimates, e.g., the mean number of touch events per session for the interval at the start of protraction (the lowest bar on the plot) is 2 touches per session. In general, each interval had at least 8 touches per individual session in the data sets that we analyzed. Nature Neuroscience: doi:.38/nn.2283

.2 Protraction Protracted Retraction Probability of Touch. /2 /2 Phase in Whisk Cycle (radians) Nature Neuroscience: doi:.38/nn.2283 Figure S - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S2. Example of raw data for a unit with a touch response modulated by whisking. Plots of high speed video frames, videographic position signal, rectified touch sensor signal, rectified (and low pass filtered) EMG signal, and S cortex recording channels surrounding two touch events during a recording session. Dots indicate times of contact events and sorted spike times for two units (waveforms and autocorrelations are shown at the bottom). The left plots show the rat touching the sensor during retraction and touch in the right plots is during protraction. Unit is rapidly excited by touch and Unit 2 does not respond to touch. The spike activity of Unit is modulated by free whisking and its touch response is modulated by whisking. The preferred phases of free whisking and touch for Unit are near retraction, i.e., φ whsik =.84 for whisking and φ touch =.82 for touch. The touch event to the left is within the preferred phase of Unit as can be seen by increased spiking briefly after touch is initiated. The touch event to the right occurs nearly opposite to the preferred phase of Unit and spiking briefly after touch remains similar to background activity. Nature Neuroscience: doi:.38/nn.2283

Epoch with Contact Near Retraction Epoch with Contact Near Protraction θ = 6 o θ = 6 o θ = 44 o θ = 34 o θ = 37 o θ = 4 o θ = 54 o θ = 59 o θ = 5 o θ = 37 o Video frames Video frames 6 o 3 o Touch Angle, θ 6 o 3 o EMG S wire 2 µv 5 ms S wire Multi-units Unit 2 Unit Unit Unit 2 Stereo waveform µv 5 µs Spikes / Bin 5-2 2 Time from spike (ms) Autocorrelation 2-2 2 Time from spike (ms) Nature Neuroscience: doi:.38/nn.2283 Figure S2 - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S3. Distribution of touch response widths. Histogram of the width of the response to touch (full width at half maximum amplitude) for all excited or inhibited touch responses, i.e., rapid excited, slow excited, and slow inhibited. Units with no response to touch are not included. Note that all rapid excited touch response widths fall within the first bin of the plot. Nature Neuroscience: doi:.38/nn.2283

4 Number of Single Units 35 3 25 2 5 Rapidly Excited Slowly Excited and Inhibited 5..2.4.6.8..2.4.6.8 2. Touch Response Width at Half Maximum Amplitude (s) Nature Neuroscience: doi:.38/nn.2283 Figure S3 - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S4. Free whisking modulation depths and spike waveform widths for all units according to touch response type. (Top) The distribution of different classes of single units according to the width, measured from peak to trough, of the waveform for the unit. The inserts show typical waveforms. (Lower) Free whisking modulation depth, measured from minimum to maximum divided by mean spike rates, plotted as functions of spike waveform width for the three difference responses, i.e., rapid excited, slow excited, and slow inhibited, as well as unresponsive cells. Nature Neuroscience: doi:.38/nn.2283

Count 2 Rapidly Excited Spike Rate Mean Min Max - Phase (radians) Free Whisk Modulation Depth Slowly Inhibited Slowly Excited No Touch Response 5 Spike Waveform Width (μs) Nature Neuroscience: doi:.38/nn.2283 Figure S4 - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S5. Probability of burst generation by single units in S vibrissa cortex. To quantify bursting activity in neurons excited by touch as opposed to free whisking in air, we produced separate interspike interval (ISI) histograms (bin width <.5 ms) for spikes that occurred during touch versus free whisking trials. To model normal spike activity, we used a generic form of the Inverse Gaussian equation to capture the refractory period spiking and exponential decay in number of ISI s at higher values. Bursting activity was modeled as a Gaussian function where the center and width were chosen manually such that the function was confined to the relative refractory period. If a particular ISI histogram contained a distinct peak at small values, which is indicative of bursting, the center and width of the Gaussian function were adjusted to best fit the peak. Functions were fit to the ISI histograms using Poisson distributed maximum likelihood estimates. Since particular ISI distributions were modeled as a linear combination of the Gaussian and Inverse Gaussian functions, the number of ISI values contributed by each could be calculated separately. Thus the number of ISI values that resulted from bursting was calculated by taking the difference in number of ISI values between the Gaussian and Inverse Gaussian fitted curves (bursting and normal spiking) and the Inverse Gaussian fitted curve (normal spiking only). The number of ISI values due to bursting was then divided by the total number of ISI values in order to compare the relative amount of bursting among different ISI distributions, i.e., between free whisking and touch trials of the same neuron as well as among different neurons. Nature Neuroscience: doi:.38/nn.2283

.6 Fraction of Spikes that Arrive in Bursts During Touch.2.8.4.4.8.2.6 Fraction of Spikes that Arrive in Bursts During Free Whisking Count per Bin 4 2 Burst Spikes Poisson Spikes 2 4 Interspike Interval (ms) Nature Neuroscience: doi:.38/nn.2283 Figure S5 - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S6. Laminar distribution and spike widths of different classes of single units. (A) Representative local field potential (LFP) measurements along with the current source density (CSD) calculated from these measurements. The CSD, shown also in false color to emphasize sinks and sources of current, was used to classify the cortical layer and depth for each recording. The plot on the right confirms that at least the short-time measurements concur with the past results of Swadlow, Gusev and Bezdudnaya (Activation of a cortical column by a thalamocortical impulse. 22. J Neurosci 22, 7766-7779). (B) Compendium of responses across all single-units. (C) The laminar distribution of single-units that responded to active touch, independent of their whisking related response. Nature Neuroscience: doi:.38/nn.2283

A LFP CSD CSD (Swadlow et al) Electrode Depth (µm) B 2 4 6 8 2 4 6 8 5 µv - - - Time from Vibrissa Stimulus Onset (ms) Touch Response Whisking Response Yes No Rapidly Excited 2 % 3 % Slowly Inhibited Slowly Excited None 9 % 5 % 9 % 4 % 3 % 7 % 63 % 37 % C Fraction of Units..5 mv / mm 2 Source Sink - -2 Charge Density ( e / nm 3 ) Supragranular Granular Infragranular 2 4 Time (ms) 4 6 8 2 4 6 Electrode Depth (µm) Nature Neuroscience: doi:.38/nn.2283 Figure S6 - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S7. Examples of whisking activity that surrounds touch events. Each plot shows whisking activity that surrounds touch events that are sorted by where touch occurs in the whisk cycle. The data applies to the three units included in figure 4. Single whisking trajectories are thin lines overlaid in grey and mean trajectories are thick black lines. The lower plots show the relative frequencies of touch at different whisk cycle phases for each unit s recording session. Nature Neuroscience: doi:.38/nn.2283

EMG Activity Surrounding Touch Events at Different Phases of the Whisk Cycle - Unit Unit 2 Unit 3 3/4 Phase in Whisk Cycle (radians) /2 /4 /4 /2 3/4 - - - Time from Touch (ms) Distributions of Whisk Cycle Phases for All Touch Events Number of Touch Events 2 2 4 2 - -/2 /2 - -/2 /2 - -/2 /2 Phase in Whisk Cycle (radians) Nature Neuroscience: doi:.38/nn.2283 Figure S7 - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S8. Free whisking and touch responses for all units with rapid response to touch that are significantly modulated by phase in the whisk cycle. Autocorrelations and waveforms for each neuron are plotted at the right end of each row. (A A4). Average spike response histograms in response to free whisking centered on EMG peaks with sinusoid fits (green lines) and resultant parameters, i.e., DC = average spike rate across histograms, Amp = amplitude of sinusoid fit (DC to peak), φ = preferred phase, mod = response modulation depth. (B B4) Average spike responses to touch with BARS fits (green lines). Peak and baseline rates calculated from smooth curve are listed above each plot. (C C4) Left: Average touch response histograms parsed by phase in the whisk cycle using EMG data to determine the phase. Right: Plots of the peak values of the touch response from fits to each of the eight intervals of the smooth curve touch responses in the left hand panels. The parameters listed above each plot were estimated by fitting a sinusoid to each curve and have the same designation as in panel A. Dotted lines indicate 95 % confidence intervals for mean spike rate estimates. (D D4) Left: Average touch response histograms parsed by phase in the whisk cycle using videographic data to determine the phase. Right: Plots of the peak values of the touch response from fits to each of the eight intervals of the smooth curve touch responses in the left hand panels. (E E4) Left: Average touch response histograms parsed by speed in the whisk cycle using videographic data to determine the speed. Right: Plots of the peak values of the touch response from GLM fits to each of the eight intervals of the smoothed touch responses in the left hand panels. (F F4) Left: Average touch response histograms parsed by angular position in the whisk cycle using videographic data to determine the position. Right: Plots of the peak values of the touch response from GLM fits to each of the eight intervals of the smoothed touch responses in the left hand panels. Nature Neuroscience: doi:.38/nn.2283

DC = 68, P2T =, Ang = 26, Mod* =.8 5 5 µs µv DC =, Amp = 2.2, f =.85, Mod =.4 A Trial Averaged Free BTrial Averaged Touch C Touch Response versus phase in D Touch Response Versus Phase in E Touch Response Versus Velocity of Whisking Response Response Whisk cycle (from EMG) Whisk Cycle (from video) Whisking Motion 5 Peak = 48, Baseline = 2 - DC = 5, Amp = 3, f =.76, Mod =.2 5 DC = 54, Amp = 35, f =.87, Mod =.3 5 DC = 56, P 2T = 33, P kvel =.53, Mod =. 59 5 DC = 52, P2T = 3 Ang = 9, Mod* =.57 5 8 FTouch Response Versus Vibrissal Position - 2 2 DC = 6, Amp = 2.8, f =.85, Mod =.35 2 5 Peak = 72, Baseline = 7 - DC = 75, Amp = 27, f =.98, Mod =.74 DC = 73, Amp = 32, f =.97, Mod =.88 DC = 74, P 2T = 34, P kvel = -.38, Mod* =. 46 DC = 72, P2T = 2, Ang = 9, Mod* =.33 8. 5 -. 5 2 Spike Rate (Hz) DC = 2, Amp =., f =.35, Mod =.2 4 DC = 2, Amp =., f =.48, Mod =.2 4 DC = 9.6, Amp =.2, f =.54, Mod =.25 5 Peak = 6, Baseline = 4 Peak = 4, Baseline = 6 Peak = 79, Baseline = 9.2 Phase in Whisk Cycle (radians) - - - DC = 56, Amp = 37, f =.27, Mod =.3 DC = 4, Amp = 24, f =.44, Mod =.2 DC = 83, Amp = 48, f =.36, Mod =.2 DC = 54, Amp = 3, f =.3, Mod =. DC = 38, Amp = 22, f =.54, Mod =.2 DC = 85, Amp = 48, f =.4, Mod =. Velocity of Vibrissa Motion (deg / ms). 5 -. 5. 5 -. 5 DC = 56, P 2T = 2, P kvel = -.98, Mod* =. 36 DC = 38, P 2T = 29, P kvel = -.48, Mod =.76 DC = 82, P 2T = 63, P kvel = -., Mod =.77 Angle Relative to Head Centerline (deg) 2 4 2 4 2 4 DC = 6, P2T = 2, Ang = 25, Mod* =.35 DC = 4, P2T = 8, Ang = 42, Mod* =.45 DC = 85, P2T = 5, Ang = 56, Mod* =.6 5 2 5 2-6 DC = 9, Amp =.8, f =.84, Mod =.8 Peak = 69, Baseline = 8 - DC = 68, Amp = 26, f =.82, Mod =.76 DC = 69, Amp = 29, f =. 9, Mod =. 84 DC = 66, P 2T = 3, P kvel = -.98, Mod* =. 47 2 5 2-6 DC = 8.6, Amp =.7, f =.36, Mod =.6 Peak = 78, Baseline = 6 DC = 82, Amp = 3, f =.38, Mod =.73-8 DC = 77, Amp = 37, f =.32, Mod =.96.5 DC = 8, P2T = 3, Vel = -.2, Mod* =.38 2 3 DC = 8, P2T = 43, Ang = 38, Mod* =.53 Time from Contact (ms).5 4 - - - Spike Rate (Hz) - Spike Rate (Hz) - Spike Rate (Hz) - Spike Rate (Hz) Time from Peak of Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) 2 2 Time from Spikes (ms) EMG (ms) Figure S8 (-7) - Curtis and Kleinfeld Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Nature Neuroscience: doi:.38/nn.2283

Spikes / bin 5 5 5 µs µv A2 Trial Averaged Free B2 Trial Averaged Touch C2 Touch Response versus phase in D2 Touch Response Versus Phase in E2 Whisking Response Response Whisk cycle (from EMG) Whisk Cycle (from video) DC = 5, Amp =, Φ =.7, Mod =.3 2 DC = 34, Amp = 3.8, f =.33, Mod =.23 4 4 Peak = 6, Baseline = 4 Peak = 9, Baseline = 28 - DC = 7, Amp = 4, f =.3, Mod =.77 - DC = 64, Amp = 43, Φ =.6, Mod =.3 2 DC = 62, Amp = 44, ph =. 3, Mod =. 4 DC =, Amp = 57, f =.29, Mod =. 2 - Touch Response Versus Velocity of Whisking Motion DC = 6, P 2T = 25 P kvel =.7, Mod* =. 4 DC = 8, P2T = 68, PkVel = -.6, Mod* =.63 2 2 4 6 F2 Touch Response Versus Vibrissal Position 6 DC = 62, P 2T = 36, Ang = 55, Mod* =.58 DC = 8, P2T = 26, Ang = 6, Mod* =.24 2 Spike Rate (Hz) DC = 4, Amp =.7, f =.3, Mod =.24 Peak = 94, Baseline = 3 DC = 94, Amp = 58, f =.2, Mod =.2 2-2 - Time from Peak of EMG (ms) DC = 5, Amp =.7, f = -.3, Mod =.22 2 DC = 3, Amp = 3., f = -.38, Mod =.47 2 Peak = 32, Baseline = 3 Peak = 99, Baseline = DC = 35, Amp = 6, f = -.28, Mod =.89 2 - DC = 4, Amp = 59, f = -.35, Mod =. 2 - DC = 5, Amp = 2.2, f =.32, Mod =.29 Peak = 48, Baseline = 5 DC = 49, Amp = 23, f =.7, Mod =.94-2 DC =, Amp = 2.3, f = -.38, Mod =.42 2 4 4 Peak = 49, Baseline = Phase in Whisk Cycle (radians) DC = 47, Amp = 7, f = -.57, Mod =.72 5 - DC = 96, Amp = 6, f =.32, Mod =.3 2 DC = 38, Amp = 67, f = -. 7, Mod =. 97 2 DC = 3, Amp = 54, f = -.58, Mod =. 2 DC = 48, Amp = 22, f =., Mod =. 92 DC = 45, Amp = 24, f = -.32, Mod =. 5 Velocity of Vibrissa Motion (deg / ms) DC = 96, P2T = 9, PkVel =.56, Mod =.7 2.5 DC = 37, P2T = 6, PkVel =.3, Mod* =. 45 2 DC = 4, P2T = 64, PkVel =.72, Mod =. 62 2 -. 5 4 - - Spike Rate (hz) - Spike Rate (hz) - Spike Rate (hz) - Spike Rate (hz) 2 2 Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) Time from Spikes -.5 4. 5 2 -. 5 4. 5 -. 5 DC = 47, P 2T = 26, P kvel = -.22, Mod* =. 55 DC = 46, P 2T = 43, P kvel = -.68, Mod =.93 5 Angle Relative to Head Centerline (deg) 6 8 2 2 4. 5 2 DC = 94, Amp = 32, Ang = 98, Mod* =.34 2 DC = 36, P2T = 58, Ang = 45, Mod* =. 43 2 DC =, P2T = 28, Ang = 98, Mod* =. 28 2 DC = 48, P 2T = 5, Ang = 38, Mod* =. 3 DC = 45, P 2T =, Ang = 32, Mod* =. 24 5 Figure S8 (8-4) - Curtis and Kleinfeld Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin 2 2 5 Nature Neuroscience: doi:.38/nn.2283

5 2 2 4 5 5 µs µv DC = 22, Amp = 2., f =.2, Mod =.9 Peak = 85, Baseline = 24 DC = 87, Amp = 43, f =.8, Mod =.99 2-2 A3 Trial Averaged Free Whisking Response B3 Trial Averaged Touch C3 Touch Response versus phase in D3 Touch Response Versus Phase in E3 Touch Response Versus Velocity of Response Whisk cycle (from EMG) Whisk Cycle (from video) Whisking Motion 5 DC =, Amp =.8, f = -.6, Mod =.32 Peak = 65, Baseline = DC = 68, Amp = 27, f = -.2, Mod =.8-2 DC = 85, Amp = 38, f =. 3, Mod =. 89 2 DC = 69, Amp = 33, f = -. 5, Mod =. 96. 5 -. 5. 5 DC = 83, P 2T = 45, P kvel = -., Mod* =. 54 2 DC = 69, P 2T = 2, P kvel = -.6, Mod* =. 3 8 2 F3 Touch Response Versus Vibrissal Position DC = 84, P 2T = 27, Ang =, Mod* =. 32 2 DC = 7, P 2T = 22, Ang = 32, Mod* =. 3 Spike Rate (Hz) 4 DC = 7, Amp =.9, f =.7, Mod =.54 Peak = 72, Baseline = 7 DC = 75, Amp = 43, f =.23, Mod =. - 8 8 4 DC =, Amp =.7, f = -.33, Mod =.34 Peak = 6, Baseline = DC = 58, Amp = 45, f = -.29, Mod =.6-4 DC = 9, Amp =.6, f =.48, Mod =.36 Peak = 72, Baseline = DC = 75, Amp = 49, f =.43, Mod =.3-5 DC = 2, Amp = 2.7, f =.4, Mod =.45 Peak = 47, Baseline = 7 4 Phase in Whisk Cycle (radians) DC = 48, Amp = 43, f =.2, Mod =.8 - DC = 76, Amp = 44, f =. 22, Mod =. 2 DC = 59, Amp = 38, f = -. 36, Mod =. 3 DC = 73, Amp = 48, f =. 45, Mod =. 3 DC = 5, Amp = 49, f =. 2, Mod =. 9 Velocity of Vibrissa Motion (deg / ms) -. 5 - -. 5 -. 5 - DC = 76, P 2T = 5, P kvel = -.79, Mod =. 66 DC = 6, P 2T = 37, P kvel =.96, Mod =.62 DC = 72, P 2T = 69, P kvel = -., Mod =.95 DC = 49, P 2T = 35, P kvel = -.27, Mod =.7 Angle Relative to Head Centerline (deg) 3 2 2 4 3 6 2 4 DC = 74, P 2T = 2, Ang = 23, Mod* =. 27 DC = 6, P 2T = 2, Ang = 44, Mod* =. 33 DC = 75, P 2T = 23, Ang = 46, Mod* =. 3 DC = 5, P 2T = 28, Ang = 58, Mod =.56-6 DC = 5, Amp = 2., f =.2, Mod =.27 Peak = 77, Baseline = 5 DC = 78, Amp = 42, f =.7, Mod =. - 2 8 DC = 76, Amp = 43, f =. 8, Mod =.. 5 DC = 78, P 2T = 52, P kvel = -.49, Mod =. 66 2 DC = 75, P 2T = 23, Ang = 8, Mod* =. 3 - Spike Rate (Hz) Time from Contact (ms) 4 4 -. 5 - - - Spike Rate (Hz) - Spike Rate (Hz) - Spike Rate (Hz) Time from Peak of Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) 2 2 Time from Spikes (mx) EMG (ms) Figure S8 (5-2) - Curtis and Kleinfeld Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Nature Neuroscience: doi:.38/nn.2283

2 25 5 µs µv A4 Trial Averaged Free Whisking Response B4 Trial Averaged Touch C4 Touch Response versus phase in D4Touch Response Versus Phase in E4 Touch Response Versus Velocity of Response Whisk cycle (from EMG) Whisk Cycle (from video) Whisking Motion F4 Touch Response Versus Vibrissal Position DC = 7, Amp =., f =.45 Mod =.3 8 2 Peak = 25, Baseline = 7 DC = 23, Amp = 2 f =.48, Mod =. 4 - DC = 25, Amp = 5, f =. 65, Mod =. 2 4 DC = 25, P 2T = 4, P kvel = -.2, Mod* =.56 4 2 DC = 24, P 2T = 2, Ang = 87, Mod* =. 5 4-4 DC = 8, Amp =., f = -.66 Mod =.28 4 Peak = 47, Baseline = 9 DC = 47, Amp = 23, f = -.87, Mod =.98 - DC = 46, Amp = 2, f = -. 69, Mod =.9 DC = 49, P 2T = 35, P kvel =.5, Mod =. 7 2 DC = 48, P 2T = 22, Ang = 97, Mod =. 46-4 Spike Rate (Hz) DC = 8, Amp =., f =.29 Mod =.25 Peak = 5, Baseline = DC = 54, Amp = 39, f =.38, Mod =.7-4 DC = 2, Amp = 2., f =.2 Mod =.35 Peak = 35, Baseline = DC = 35, Amp = 2, f =.22, Mod =.2-4 DC = 8, Amp =.6, f = -.82 Mod =.5 Peak = 38, Baseline = 7 DC = 4, Amp = 27, f = -.72, Mod =.4-4 Phase in Whisk Cycle (radians) DC = 56, Amp = 35, f =. 35, Mod =. 3 DC = 37, Amp = 28, f =. 8, Mod =. 5 DC = 42, Amp = 25, f = -. 68, Mod =. 2 Whisking Velocity (deg / ms) DC = 55, P 2T = 3, P kvel = -., Mod =. 55. 5 -. 5 DC = 35, P 2T = 25, P kvel = -., Mod =.7 -. 5 -. 5 DC = 42, P 2T = 33, P kvel =.82, Mod =.79 Vibrissa Angle (deg) 3 5 2 DC = 57 P 2T = 2, Ang = 24, Mod* =. 37 DC = 34, P 2T = 8, Ang = 3, Mod =.53 4 DC = 4, P 2T = 9, Ang = 32, Mod* =. 23 3 4 5 2 2 2 DC =, Amp =.3, f =.28 Mod =.26 Peak = 45, Baseline = 9 DC = 46, Amp = 33, f =.2, Mod =.4 - DC = 45, Amp = 26, f =. 2, Mod =. 2 DC = 44, P 2T = 25, P kvel =.6, Mod* =. 57 DC = 46, P 2T = 22, Ang = 39, Mod* =. 48 4 2-4 - Spike Rate (Hz) Time from Contact (ms) DC = 7, Amp =.2, f =.7 Mod =.34 Peak = 66, Baseline = DC = 7, Amp = 5, f =., Mod =.4 DC = 67, Amp = 6, f =. 3, Mod =.8 DC = 68, P 2T = 37, P kvel = -.2, Mod =.54 DC = 69, P 2T = 53, Ang = 42, Mod =.77 -. 5 2 6 3 -. 5 4 - - - Spike Rate (Hz) - Spike Rate (Hz) - Spike Rate (Hz) Time from Peak of Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) Time from Contact (ms) 2 2 Time from Spikes (ms) EMG (ms) Figure S8 (22-28) - Curtis and Kleinfeld Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Nature Neuroscience: doi:.38/nn.2283

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S9. Free whisking and touch responses for all units with rapid response to touch that are not modulated by phase in the whisk cycle. (A). Average spike response histograms in response to free whisking centered on EMG peaks with sinusoid fits (green lines) and resultant parameters defined in caption for figure S8. An asterisk denotes insignificant fit. (B) Average spike responses to touch with BARS fits (green lines). Peak and baseline rates calculated from smooth curve are listed above each plot. (C) Average touch response histograms parsed by phase in the whisk cycle. (D) Plots of the peak values of the touch response from fits to each of the eight intervals of the smooth curve touch responses in panel C. The parameters listed above each plot were estimated by fitting a sinusoid to each curve and have the same designation as in panel A. Dotted lines indicate 95 % confidence intervals for mean spike rate estimates. (E) Autocorrelations and waveforms for each neuron. Nature Neuroscience: doi:.38/nn.2283

5 2 5 5 5 2 2 Time from Spikes (ms) 5 µs µv Nature Neuroscience: doi:.38/nn.2283 Figure S9 - Curtis and Kleinfeld A Trial Averaged Free Whisking Response DC = 8, Amp =.5, Φ* =.29, mod =.3 B Trial Averaged Touch Response C Touch Response Parsed by Phase in Whisk Cycle D Touch Response versus Phase in Whisk Cycle Peak = 32, Baseline = 8 DC = 35, Amp = 4, Φ* =.22, mod =.2 - E Unit Autocorrelations and Waveforms Spike Rate (Hz) 4 DC = 3, Amp =.8, Φ =.36, mod =.27 2 DC = 6, Amp =.2, Φ* =.8, mod =.7 4 2 DC = 2, Amp =.5, Φ =.94, mod =.25 2 DC = 2, Amp =.5, Φ* =.25, mod =.8 4 Peak = 24, Baseline = 6 DC = 28, Amp =.8, Φ* =.76, mod =.6 Peak = 5, Baseline = 29 4 2 4 4 DC = 5, Amp =.7, Φ =.48, mod =. 5 Peak = 57, Baseline = 4 Peak = 45, Baseline = Peak = 46, Baseline = 2 Peak = 5, Baseline = Phase in Whisk Cycle (radians) - - - - - - Spike Rate (Hz) 5 DC = 58, Amp = 7, Φ* =.56, Mod =.24 DC = 25, Amp = 3, Φ* = -., Mod =.24 5 DC = 48, Amp =, Φ* =.9, Mod =.46 DC = 4, Amp = 8,Φ* =.87, Mod =.35 2 DC = 45, Amp = 6, Φ* =.5, Mod =.27 5 DC = 5, Amp = 2, Φ* =.9, Mod =.2 5 Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin Spikes / bin - - - Time from Peak of Time from Contact (ms) Time from Contact (ms) EMG (ms) - Phase in Whisk Cycle (radians)

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S. Comparisons of maximum touch response and modulation depth at different preferred phases. (A) The spike rate for the maximum touch response as a function of preferred phase for units with whisking modulated rapid touch responses. (B) The modulation depth of the touch response plotted as a function of each unit s preferred phase. Mean values across all preferred phases are indicated in each plot by dark green line along with 95 % confidence band in light green. There is no significant tendency for either maximum touch response amplitude or modulation depth to be greater at any specific preferred phase relative to other preferred phases. Nature Neuroscience: doi:.38/nn.2283

Protraction Protracted Retraction Maximum Touch Response by Phase in Whisk Cycle Maximum Spike Rate of Touch Response (Hz) 2 μ+2σ μ μ 2σ.6 μ+2σ Modulation Depth of Touch Response.2.8.4 μ μ 2σ /2 /2 Preferred Phase in Whisk Cycle, φ whisk (radians) Nature Neuroscience: doi:.38/nn.2283 Figure S - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S. The preferred phase in the whisk cycle for contact estimated from EMG data versus videographic data. The slope of the best fit line is.92 (p <.). Nature Neuroscience: doi:.38/nn.2283

Count - Preferred Phase During Touch, f touch (radians), from EMG -/2 /2 - -/2 /2 Preferred Phase During Touch, Count φ touch (radians), from Videography Nature Neuroscience: doi:.38/nn.2283 Figure S - Curtis and Kleinfeld

Supplemental Material for NN-A2678-T Phase-to-rate transformations encode touch in cortical neurons of a scanning sensorimotor system by John Curtis and David Kleinfeld Figure S2. Touch responses as a function of phase in the whisk cycle for changes in whisking parameters. Shown are examples, including that for unit 3 shown in figure 6 of the main text. Nature Neuroscience: doi:.38/nn.2283

f : 5-8 Unit 26 DC = 42, Amp = 25, f = -. 68, Mod =. 2 All DC = 4, Amp = 27, f = -. 65, Mod =. 3 f : 8-2 DC = 43, Amp = 26, f = -. 5, Mod =. 2 DC = 42, Amp = 9, f = -. 64, Mod =. 9 amp: 3-5 DC = 44, Amp = 33, f = -. 6, Mod =. 5 amp: 5-28 Figure S2 - Curtis and Kleinfeld Unit 5 Unit 3* Unit 2 Unite 2 Unit 2 DC = 85, Amp = 48, f =.4, Mod =. DC = 73, Amp = 48, f =. 45, Mod =. 3 DC = 3, Amp = 54, f = -.58, Mod =. DC = 5, Amp = 49, f =. 2, Mod =. 9 DC = 76, Amp = 43, f =. 8, Mod =. All All All All All All trials - - - 2 - - - Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Top 5 % f whisk Bottom 5 % f whisk DC = 84, Amp = 4, f =. 43, Mod =. f : 4-9 DC = 83, Amp = 47, f =. 36, Mod =. f : 9-3 DC = 73, Amp = 45, ph =. 37, Mod =. 3 f : 5-9 DC = 72, Amp = 46., ph =. 4, Mod =. 3 f : 9-3 DC = 3, Amp = 4, f = -. 54, Mod =. 8 f : 4-8 DC = 2, Amp = 58, f = -. 68, Mod =. f : 8-3 DC = 53, Amp = 47, f =. 7, Mod =.8 f : 6-8 DC = 49, Amp = 46, f =., Mod =.9 f : 8-2 DC = 73, Amp = 33. 89, ph =. 22, Mod =. 4 f : 4-8 DC = 75, Amp = 5. 54, ph =. 3, Mod =. 4 f : 8-2 Top 5 % θ o Botom 5 % θ o DC = 85, Amp = 35, f =. 33, Mod =. 8 amp: - 6 DC = 83, Amp = 49, f =. 4, Mod =. 2 amp: 6-24 DC = 74, Amp = 45, ph =. 46, Mod =. 2 amp: 4-2 DC = 72, Amp = 46, ph =. 38, Mod =. 3 amp: 2-25 DC = 3, Amp = 53, f = -. 59, Mod =. amp: 7-4 DC = 4, Amp = 5, f = -. 67, Mod =. amp: 4-28 DC = 48, Amp = 42, f =. 3, Mod =. 8 amp: 9-24 DC = 52, Amp = 45, f =. 3, Mod =. 7 amp: 24-37 DC = 78, Amp = 43. 26, ph =. 8, Mod =. 8 amp: 6-4 DC = 77, Amp = 38. 53, ph =. 2, Mod =. 23 amp: 4-29 - 2 - - 2 3 - - - Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Time after contact (ms) Spike Rate (Hz) Nature Neuroscience: doi:.38/nn.2283