Stimulus Dependence of Barrel Cortex Directional Selectivity

Neurons throughout the rat vibrissa somatosensory pathway are sensitive to the angular direction of whisker movement. Could this sensitivity help rats discriminate stimuli? Here we use a simple computational model of cortical neurons to analyze the robustness of directional selectivity. In the model, directional preference emerges from tuning of synaptic conductance amplitude and latency, as in recent experimental findings. We find that directional selectivity during stimulation with random deflection sequences is strongly dependent on the mean deflection frequency: Selectivity is weakened at high frequencies even when each individual deflection evokes strong directional tuning. This variability of directional selectivity is due to generic properties of synaptic integration by the neuronal membrane, and is therefore likely to hold under very general physiological conditions. Our results suggest that directional selectivity depends on stimulus context. It may participate in tasks involving brief whisker contact, such as detection of object position, but is likely to be weakened in tasks involving sustained whisker exploration (e.g., texture discrimination)

Dynamics of Action Potential Initiation in the GABAergic Thalamic Reticular Nucleus In Vivo

Understanding the neural mechanisms of action potential generation is critical to establish the way neural circuits generate and coordinate activity. Accordingly, we investigated the dynamics of action potential initiation in the GABAergic thalamic reticular nucleus (TRN) using in vivo intracellular recordings in cats in order to preserve anatomically-intact axo-dendritic distributions and naturally-occurring spatiotemporal patterns of synaptic activity in this structure that regulates the thalamic relay to neocortex. We found a wide operational range of voltage thresholds for action potentials, mostly due to intrinsic voltage-gated conductances and not synaptic activity driven by network oscillations. Varying levels of synchronous synaptic inputs produced fast rates of membrane potential depolarization preceding the action potential onset that were associated with lower thresholds and increased excitability, consistent with TRN neurons performing as coincidence detectors. On the other hand the presence of action potentials preceding any given spike was associated with more depolarized thresholds. The phase-plane trajectory of the action potential showed somato-dendritic propagation, but no obvious axon initial segment component, prominent in other neuronal classes and allegedly responsible for the high onset speed. Overall, our results suggest that TRN neurons could flexibly integrate synaptic inputs to discharge action potentials over wide voltage ranges, and perform as coincidence detectors and temporal integrators, supported by a dynamic action potential threshold

Millisecond-Timescale Local Network Coding in the Rat Primary Somatosensory Cortex

Correlation among neocortical neurons is thought to play an indispensable role in mediating sensory processing of external stimuli. The role of temporal precision in this correlation has been hypothesized to enhance information flow along sensory pathways. Its role in mediating the integration of information at the output of these pathways, however, remains poorly understood. Here, we examined spike timing correlation between simultaneously recorded layer V neurons within and across columns of the primary somatosensory cortex of anesthetized rats during unilateral whisker stimulation. We used Bayesian statistics and information theory to quantify the causal influence between the recorded cells with millisecond precision. For each stimulated whisker, we inferred stable, whisker-specific, dynamic Bayesian networks over many repeated trials, with network similarity of 83.3±6% within whisker, compared to only 50.3±18% across whiskers. These networks further provided information about whisker identity that was approximately 6 times higher than what was provided by the latency to first spike and 13 times higher than what was provided by the spike count of individual neurons examined separately. Furthermore, prediction of individual neurons' precise firing conditioned on knowledge of putative pre-synaptic cell firing was 3 times higher than predictions conditioned on stimulus onset alone. Taken together, these results suggest the presence of a temporally precise network coding mechanism that integrates information across neighboring columns within layer V about vibrissa position and whisking kinetics to mediate whisker movement by motor areas innervated by layer V

Neural Computation via Neural Geometry: A Place Code for Inter-whisker Timing in the Barrel Cortex?

The place theory proposed by Jeffress (1948) is still the dominant model of how the brain represents the movement of sensory stimuli between sensory receptors. According to the place theory, delays in signalling between neurons, dependent on the distances between them, compensate for time differences in the stimulation of sensory receptors. Hence the location of neurons, activated by the coincident arrival of multiple signals, reports the stimulus movement velocity. Despite its generality, most evidence for the place theory has been provided by studies of the auditory system of auditory specialists like the barn owl, but in the study of mammalian auditory systems the evidence is inconclusive. We ask to what extent the somatosensory systems of tactile specialists like rats and mice use distance dependent delays between neurons to compute the motion of tactile stimuli between the facial whiskers (or ‘vibrissae’). We present a model in which synaptic inputs evoked by whisker deflections arrive at neurons in layer 2/3 (L2/3) somatosensory ‘barrel’ cortex at different times. The timing of synaptic inputs to each neuron depends on its location relative to sources of input in layer 4 (L4) that represent stimulation of each whisker. Constrained by the geometry and timing of projections from L4 to L2/3, the model can account for a range of experimentally measured responses to two-whisker stimuli. Consistent with that data, responses of model neurons located between the barrels to paired stimulation of two whiskers are greater than the sum of the responses to either whisker input alone. The model predicts that for neurons located closer to either barrel these supralinear responses are tuned for longer inter-whisker stimulation intervals, yielding a topographic map for the inter-whisker deflection interval across the surface of L2/3. This map constitutes a neural place code for the relative timing of sensory stimuli

A Threshold Equation for Action Potential Initiation

In central neurons, the threshold for spike initiation can depend on the stimulus and varies between cells and between recording sites in a given cell, but it is unclear what mechanisms underlie this variability. Properties of ionic channels are likely to play a role in threshold modulation. We examined in models the influence of Na channel activation, inactivation, slow voltage-gated channels and synaptic conductances on spike threshold. We propose a threshold equation which quantifies the contribution of all these mechanisms. It provides an instantaneous time-varying value of the threshold, which applies to neurons with fluctuating inputs. We deduce a differential equation for the threshold, similar to the equations of gating variables in the Hodgkin-Huxley formalism, which describes how the spike threshold varies with the membrane potential, depending on channel properties. We find that spike threshold depends logarithmically on Na channel density, and that Na channel inactivation and K channels can dynamically modulate it in an adaptive way: the threshold increases with membrane potential and after every action potential. Our equation was validated with simulations of a previously published multicompartemental model of spike initiation. Finally, we observed that threshold variability in models depends crucially on the shape of the Na activation function near spike initiation (about −55 mV), while its parameters are adjusted near half-activation voltage (about −30 mV), which might explain why many models exhibit little threshold variability, contrary to experimental observations. We conclude that ionic channels can account for large variations in spike threshold