Coding of stimulus frequency by latency in thalamic networks through the interplay of GABAB-mediated feedback and stimulus shape

Theory of hierarchically organized neuronal oscillator dynamics that mediate rodent rhythmic whisking

Rodents explore their environment through coordinated orofacial motor actions, including whisking. Whisking can free-run via an oscillator of inhibitory neurons in the medulla and can be paced by breathing. Yet, the mechanics of the whisking oscillator and its interaction with breathing remain to be understood. We formulate and solve a hierarchical model of the whisking circuit. The first whisk within a breathing cycle is generated by inhalation, which resets a vibrissa oscillator circuit, while subsequent whisks are derived from the oscillator circuit. Our model posits, consistent with experiment, that there are two subpopulations of oscillator neurons. Stronger connections between the subpopulations support rhythmicity, while connections within each subpopulation induce variable spike timing that enhances the dynamic range of rhythm generation. Calculated cycle-to-cycle changes in whisking are consistent with experiment. Our model provides a computational framework to support longstanding observations of concurrent autonomous and driven rhythmic motor actions that comprise behaviors.

A vibrissa pathway that activates the limbic system

Vibrissa sensory inputs play a central role in driving rodent behavior. These inputs transit through the sensory trigeminal nuclei, which give rise to the ascending lemniscal and paralemniscal pathways. While lemniscal projections are somatotopically mapped from brainstem to cortex, those of the paralemniscal pathway are more widely distributed. Yet the extent and topography of paralemniscal projections are unknown, along with the potential role of these projections in controlling behavior. Here, we used viral tracers to map paralemniscal projections. We find that this pathway broadcasts vibrissa-based sensory signals to brainstem regions that are involved in the regulation of autonomic functions and to forebrain regions that are involved in the expression of emotional reactions. We further provide evidence that GABAergic cells of the Kölliker-Fuse nucleus gate trigeminal sensory input in the paralemniscal pathway via a mechanism of presynaptic or extrasynaptic inhibition.

Corticothalamic Pathways From Layer 5: Emerging Roles in Computation and Pathology

Large portions of the thalamus receive strong driving input from cortical layer 5 (L5) neurons but the role of this important pathway in cortical and thalamic computations is not well understood. L5-recipient "higher-order" thalamic regions participate in cortico-thalamo-cortical (CTC) circuits that are increasingly recognized to be (1) anatomically and functionally distinct from better-studied "first-order" CTC networks, and (2) integral to cortical activity related to learning and perception. Additionally, studies are beginning to elucidate the clinical relevance of these networks, as dysfunction across these pathways have been implicated in several pathological states. In this review, we highlight recent advances in understanding L5 CTC networks across sensory modalities and brain regions, particularly studies leveraging cell-type-specific tools that allow precise experimental access to L5 CTC circuits. We aim to provide a focused and accessible summary of the anatomical, physiological, and computational properties of L5-originating CTC networks, and outline their underappreciated contribution in pathology. We particularly seek to connect single-neuron and synaptic properties to network (dys)function and emerging theories of cortical computation, and highlight information processing in L5 CTC networks as a promising focus for computational studies.

Simulation of electromyographic recordings following transcranial magnetic stimulation

Transcranial magnetic stimulation (TMS) is a technique that enables noninvasive manipulation of neural activity and holds promise in both clinical and basic research settings. The effect of TMS on the motor cortex is often measured by electromyography (EMG) recordings from a small hand muscle. However, the details of how TMS generates responses measured with EMG are not completely understood. We aim to develop a biophysically detailed computational model to study the potential mechanisms underlying the generation of EMG signals following TMS. Our model comprises a feed-forward network of cortical layer 2/3 cells, which drive morphologically detailed layer 5 corticomotoneuronal cells, which in turn project to a pool of motoneurons. EMG signals are modeled as the sum of motor unit action potentials. EMG recordings from the first dorsal interosseous muscle were performed in four subjects and compared with simulated EMG signals. Our model successfully reproduces several characteristics of the experimental data. The simulated EMG signals match experimental EMG recordings in shape and size, and change with stimulus intensity and contraction level as in experimental recordings. They exhibit cortical silent periods that are close to the biological values and reveal an interesting dependence on inhibitory synaptic transmission properties. Our model predicts several characteristics of the firing patterns of neurons along the entire pathway from cortical layer 2/3 cells down to spinal motoneurons and should be considered as a viable tool for explaining and analyzing EMG signals following TMS. NEW & NOTEWORTHY A biophysically detailed model of EMG signal generation following transcranial magnetic stimulation (TMS) is proposed. Simulated EMG signals match experimental EMG recordings in shape and amplitude. Motor-evoked potential and cortical silent period properties match experimental data. The model is a viable tool to analyze, explain, and predict EMG signals following TMS.

Vibrissa Self-Motion and Touch Are Reliably Encoded along the Same Somatosensory Pathway from Brainstem through Thalamus

Active sensing involves the fusion of internally generated motor events with external sensation. For rodents, active somatosensation includes scanning the immediate environment with the mystacial vibrissae. In doing so, the vibrissae may touch an object at any angle in the whisk cycle. The representation of touch and vibrissa self-motion may in principle be encoded along separate pathways, or share a single pathway, from the periphery to cortex. Past studies established that the spike rates in neurons along the lemniscal pathway from receptors to cortex, which includes the principal trigeminal and ventral-posterior-medial thalamic nuclei, are substantially modulated by touch. In contrast, spike rates along the paralemniscal pathway, which includes the rostral spinal trigeminal interpolaris, posteromedial thalamic, and ventral zona incerta nuclei, are only weakly modulated by touch. Here we find that neurons along the lemniscal pathway robustly encode rhythmic whisking on a cycle-by-cycle basis, while encoding along the paralemniscal pathway is relatively poor. Thus, the representations of both touch and self-motion share one pathway. In fact, some individual neurons carry both signals, so that upstream neurons with a supralinear gain function could, in principle, demodulate these signals to recover the known decoding of touch as a function of vibrissa position in the whisk cycle.

Coding of object location in the vibrissal thalamocortical system

In whisking rodents, object location is encoded at the receptor level by a combination of motor and sensory related signals. Recoding of the encoded signals can result in various forms of internal representations. Here, we examined the coding schemes occurring at the first forebrain level that receives inputs necessary for generating such internal representations--the thalamocortical network. Single units were recorded in 8 thalamic and cortical stations in artificially whisking anesthetized rats. Neuronal representations of object location generated across these stations and expressed in response latency and magnitude were classified based on graded and binary coding schemes. Both graded and binary coding schemes occurred across the entire thalamocortical network, with a general tendency of graded-to-binary transformation from thalamus to cortex. Overall, 63% of the neurons of the thalamocortical network coded object position in their firing. Thalamocortical responses exhibited a slow dynamics during which the amount of coded information increased across 4-5 whisking cycles and then stabilized. Taken together, the results indicate that the thalamocortical network contains dynamic mechanisms that can converge over time on multiple coding schemes of object location, schemes which essentially transform temporal coding to rate coding and gradual to labeled-line coding.

Reinforcement active learning in the vibrissae system: optimal object localization

Rats move their whiskers to acquire information about their environment. It has been observed that they palpate novel objects and objects they are required to localize in space. We analyze whisker-based object localization using two complementary paradigms, namely, active learning and intrinsic-reward reinforcement learning. Active learning algorithms select the next training samples according to the hypothesized solution in order to better discriminate between correct and incorrect labels. Intrinsic-reward reinforcement learning uses prediction errors as the reward to an actor-critic design, such that behavior converges to the one that optimizes the learning process. We show that in the context of object localization, the two paradigms result in palpation whisking as their respective optimal solution. These results suggest that rats may employ principles of active learning and/or intrinsic reward in tactile exploration and can guide future research to seek the underlying neuronal mechanisms that implement them. Furthermore, these paradigms are easily transferable to biomimetic whisker-based artificial sensors and can improve the active exploration of their environment.

LTS and FS inhibitory interneurons, short-term synaptic plasticity, and cortical circuit dynamics

Somatostatin-expressing, low threshold-spiking (LTS) cells and fast-spiking (FS) cells are two common subtypes of inhibitory neocortical interneuron. Excitatory synapses from regular-spiking (RS) pyramidal neurons to LTS cells strongly facilitate when activated repetitively, whereas RS-to-FS synapses depress. This suggests that LTS neurons may be especially relevant at high rate regimes and protect cortical circuits against over-excitation and seizures. However, the inhibitory synapses from LTS cells usually depress, which may reduce their effectiveness at high rates. We ask: by which mechanisms and at what firing rates do LTS neurons control the activity of cortical circuits responding to thalamic input, and how is control by LTS neurons different from that of FS neurons? We study rate models of circuits that include RS cells and LTS and FS inhibitory cells with short-term synaptic plasticity. LTS neurons shift the RS firing-rate vs. current curve to the right at high rates and reduce its slope at low rates; the LTS effect is delayed and prolonged. FS neurons always shift the curve to the right and affect RS firing transiently. In an RS-LTS-FS network, FS neurons reach a quiescent state if they receive weak input, LTS neurons are quiescent if RS neurons receive weak input, and both FS and RS populations are active if they both receive large inputs. In general, FS neurons tend to follow the spiking of RS neurons much more closely than LTS neurons. A novel type of facilitation-induced slow oscillations is observed above the LTS firing threshold with a frequency determined by the time scale of recovery from facilitation. To conclude, contrary to earlier proposals, LTS neurons affect the transient and steady state responses of cortical circuits over a range of firing rates, not only during the high rate regime; LTS neurons protect against over-activation about as well as FS neurons.

The brain consists of circuits of neurons that signal to one another via synapses. There are two classes of neurons: excitatory cells, which cause other neurons to become more active, and inhibitory neurons, which cause other neurons to become less active. It is thought that the activity of excitatory neurons is kept in check largely by inhibitory neurons; when such an inhibitory “brake” fails, a seizure can result. Inhibitory neurons of the low-threshold spiking (LTS) subtype can potentially fulfill this braking, or anticonvulsant, role because the synaptic input to these neurons facilitates, i.e., those neurons are active when excitatory neurons are strongly active. Using a computational model we show that, because the synaptic output of LTS neurons onto excitatory neurons depresses (decreases with activity), the ability of LTS neurons to prevent strong cortical activity and seizures is not qualitatively larger than that of inhibitory neurons of another subtype, the fast-spiking (FS) cells. Furthermore, short-term (∼one second) changes in the strength of synapses to and from LTS interneurons allow them to shape the behavior of cortical circuits even at modest rates of activity, and an RS-LTS-FS circuit is capable of producing slow oscillations, on the time scale of these short-term changes.

Functional organization of the thalamic input to the thalamic reticular nucleus

Most axons connecting the thalamus and cortex in both directions pass through the thalamic reticular nucleus (TRN), a thin layer of GABAergic cells adjacent to the thalamus, and innervate neurons there. The TRN, therefore, is in a strategic location to regulate thalamocortical communication. We recorded neurons of the somatosensory region of the TRN in a thalamocortical slice preparation and studied the spatial organization of their thalamic input using laser scanning photostimulation. We show that the thalamoreticular pathway is organized topographically for most neurons. The somatosensory region of the TRN can be organized into three tiers. From the inner (thalamoreticular) border to the outer, in a manner roughly reciprocal to the reticulothalamic pathway, each of these tiers receives its input from one of the somatosensory relays of the thalamus--the posterior medial, ventroposterior medial, and ventroposterior lateral nuclei. What is surprising is that approximately a quarter of the recorded neurons received input from multiple thalamic regions usually located in different nuclei. These neurons distribute evenly throughout the thickness of the TRN. Our results, therefore, suggest that there exist a subpopulation of TRN neurons that receive convergent inputs from multiple thalamic sources and engage in more complex patterns of inhibition of relay cells. We propose these neurons enable the TRN to act as an externally driven "searchlight" that integrates cortical and subcortical inputs and then inhibits or disinhibits specific thalamic relay cells, so that appropriate information can get through the thalamus to the cortex.

Autism and EEG phase reset: deficient GABA mediated inhibition in thalamo-cortical circuits

The purpose of this study was to explore the relationship between electroencephalogram (EEG) phase reset in autism spectrum disorder (ASD) subjects as compared to age matched control subjects. The EEG was recorded from 19 scalp locations from 54 autistic subjects and 241 control subjects ranging in age from 2.6 years to 11 years. Complex demodulation was used to compute instantaneous phase differences between all pairs of electrodes and the 1st and 2nd derivatives were used to measure phase reset by phase shift duration and phase lock duration. In both short (6 cm) and long (21-24 cm) inter-electrode distances phase shift duration in ASD subjects was significantly shorter in all frequency bands but especially in the alpha-1 frequency band (8-10 Hz) (p < .0001). Phase lock duration was significantly longer in the alpha-2 frequency band (10-12 Hz) in ASD subjects (p < .0001). An anatomical gradient was present with the occipital-parietal regions the most significant. The findings in this study support the hypothesis that neural resource recruitment occurs in the lower frequency bands and especially the alpha-1 frequency band while neural resource allocation occurs in the alpha-2 frequency band. The results are consistent with a general GABA inhibitory neurotransmitter deficiency resulting in reduced number and/or strength of thalamo-cortical connections in autistic subjects.

Modeling the emergence of whisker direction maps in rat barrel cortex

Based on measuring responses to rat whiskers as they are mechanically stimulated, one recent study suggests that barrel-related areas in layer 2/3 rat primary somatosensory cortex (S1) contain a pinwheel map of whisker motion directions. Because this map is reminiscent of topographic organization for visual direction in primary visual cortex (V1) of higher mammals, we asked whether the S1 pinwheels could be explained by an input-driven developmental process as is often suggested for V1. We developed a computational model to capture how whisker stimuli are conveyed to supragranular S1, and simulate lateral cortical interactions using an established self-organizing algorithm. Inputs to the model each represent the deflection of a subset of 25 whiskers as they are contacted by a moving stimulus object. The subset of deflected whiskers corresponds with the shape of the stimulus, and the deflection direction corresponds with the movement direction of the stimulus. If these two features of the inputs are correlated during the training of the model, a somatotopically aligned map of direction emerges for each whisker in S1. Predictions of the model that are immediately testable include (1) that somatotopic pinwheel maps of whisker direction exist in adult layer 2/3 barrel cortex for every large whisker on the rat's face, even peripheral whiskers; and (2) in the adult, neurons with similar directional tuning are interconnected by a network of horizontal connections, spanning distances of many whisker representations. We also propose specific experiments for testing the predictions of the model by manipulating patterns of whisker inputs experienced during early development. The results suggest that similar intracortical mechanisms guide the development of primate V1 and rat S1.

Derivation and analysis of an ordinary differential equation mean-field model for studying clinically recorded epilepsy dynamics

In this paper we describe how an ordinary differential equation model of corticothalamic interactions may be obtained from a more general system of delay differential equations. We demonstrate that transitions to epileptic dynamics via changes in system parameters are qualitatively the same as in the original model with delay, as well as demonstrating that the onset of epileptic activity may arise due to regions of bistability. Hence, the model presents in one unique framework, two competing theories for the genesis of epileptiform activity. Similarities between model transitions and clinical data are presented and we argue that statistics obtained from, and a parameter estimation of this model may be a potential means of classifying and predicting the onset and offset of seizure activity.

Encoding of stimulus frequency and sensor motion in the posterior medial thalamic nucleus

In all sensory systems, information is processed along several parallel streams. In the vibrissa-to-barrel cortex system, these include the lemniscal system and the lesser-known paralemniscal system. The posterior medial nucleus (POm) is the thalamic structure associated with the latter pathway. Previous studies suggested that POm response latencies are positively correlated with stimulation frequency and negatively correlated with response duration, providing a basis for a phase locked loop-temporal decoding of stimulus frequency. We tested this hypothesis by analyzing response latencies of POm neurons, in both awake and anesthetized rats, to vibrissae deflections at frequencies between 0.3 and 11 Hz. We found no significant, systematic correlation between stimulation frequency and the latency or duration of POm responses. We obtained similar findings from recording in awake rats, in rats under different anesthetics, and in anesthetized rats in which the reticular activating system was stimulated. These findings suggest that stimulus frequency is not reliably reflected in response latency of POm neurons. We also tested the hypothesis that POm neurons respond preferentially to sensor motion, that is, they respond to whisking in air, without contacts. We recorded from awake, head-restrained rats while monitoring vibrissae movements. All POm neurons responded to passive whisker deflections, but none responded to noncontact whisking. Thus like their counterparts in the trigeminal ganglion, POm neurons may not reliably encode whisking kinematics. These observations suggest that POm neurons might not faithfully encode vibrissae inputs to provide reliable information on vibrissae movements or contacts.

Thalamocortical transformations of periodic stimuli: the effect of stimulus velocity and synaptic short-term depression in the vibrissa-barrel system

Recent works on the response of barrel neurons to periodic deflections of the rat vibrissae have shown that the stimulus velocity is encoded in the corti cal spike rate (Pinto et al., Journal of Neurophysiology, 83(3), 1158-1166, 2000; Arabzadeh et al., Journal of Neuroscience, 23(27), 9146-9154, 2003). Other studies have reported that repetitive pulse stimulation produces band-pass filtering of the barrel response rate centered around 7-10 Hz (Garabedian et al., Journal of Neurophysiology, 90, 1379-1391, 2003) whereas sinusoidal stimulation gives an increasing rate up to 350 Hz (Arabzadeh et al., Journal of Neuroscience, 23(27), 9146-9154, 2003). To explore the mechanisms underlying these results we propose a simple computational model consisting in an ensemble of cells in the ventro-posterior medial thalamic nucleus (VPm) encoding the stimulus velocity in the temporal profile of their response, connected to a single barrel cell through synapses showing short-term depression. With sinusoidal stimulation, encoding the velocity in VPm facilitates the response as the stimulus frequency increases and it causes the velocity to be encoded in the cortical rate in the frequency range 20-100 Hz. Synaptic depression does not suppress the response with sinusoidal stimulation but it produces a band-pass behavior using repetitive pulses. We also found that the passive properties of the cell membrane eventually suppress the response to sinusoidal stimulation at high frequencies, something not observed experimentally. We argue that network effects not included here must be important in sustaining the response at those frequencies.

VPM and PoM nuclei of the rat somatosensory thalamus: intrinsic neuronal properties and corticothalamic feedback

Sensory information originating in individual whisker follicles ascends through focused projections to the brainstem, then to the ventral posteromedial nucleus (VPM) of the thalamus, and finally into barrels of the primary somatosensory cortex (S1). By contrast, the posteromedial complex (PoM) of the thalamus receives more diffuse sensory projections from the brainstem and projects to the interbarrel septa of S1. Both VPM and PoM receive abundant corticothalamic projections from S1. Using a thalamocortical slice preparation, we characterized differences in intrinsic neuronal properties and in responses to corticothalamic feedback in neurons of VPM and PoM. Due to the plane of the slice, the majority of our observed responses came from activation of layer VI because most or all of the layer V axons terminating in PoM are cut. We found that VPM neurons exhibit higher firing rates than PoM neurons when stimulated with injected current. Stimulation of corticothalamic fibers evoked monosynaptic excitation, disynaptic inhibition, or a combination of the two in both nuclei. A few differences in the feedback responses emerged: purely excitatory postsynaptic potentials (EPSPs) in VPM were smaller and facilitated more than those in PoM, and only the EPSPs in VPM had a strong NMDA component. For both nuclei, some of the feedback responses were purely disynaptic inhibitory postsynaptic potentials (IPSPs) from the thalamic reticular nucleus (TRN). This was due to EPSP failures within VPM and PoM combined with greater reliability of S1-originating synapses onto TRN. These findings suggest that despite the exclusively excitatory nature of corticothalamic fibers, activation of cortex can trigger excitation or inhibition in thalamic relay neurons.

Active sensation: insights from the rodent vibrissa sensorimotor system

Rats sweep their vibrissae through space to locate objects in their immediate environment. In essence, their view of the proximal world is generated through pliable hairs that tap and palpate objects. The texture and shape of those objects must be discerned for the rat to assess the value of the object. Furthermore, the location of those objects must be specified with reference to the position of the rat's head for the rat to plan its movements. Recent in vivo and in vitro electrophysiological measurements provide insight into the algorithms and mechanisms that underlie these behavioral-based computations.