Firing Rate Modulation Research Articles

Progressively shifting patterns of co-modulation among premotor cortex neurons carry dynamically similar signals during action execution and observation.

Many neurons in the premotor cortex show firing rate modulation whether the subject performs an action or observes another individual performing a similar action. Although such "mirror neurons" have been thought to have highly congruent discharge during execution and observation, many if not most actually show non-congruent activity. Studies of neuronal populations active during both execution and observation have shown that the most prevalent patterns of co-modulation-captured as neural trajectories-pass through subspaces which are shared in part, but in part are visited exclusively during either execution or observation. These studies focused on reaching movements for which low-dimensional neural trajectories exhibit comparatively simple dynamical motifs. But the neural dynamics of hand movements are more complex. We developed a novel approach to examine prevalent patterns of co-modulation during execution and observation of a task that involved reaching, grasping, and manipulation. Rather than following neural trajectories in subspaces that contain their entire time course, we identified time series of instantaneous subspaces, calculated principal angles among them, sampled trajectory segments at the times of selected behavioral events, and projected those segments into the time series of instantaneous subspaces. We found that instantaneous neural subspaces most often remained distinct during execution versus observation. Nevertheless, latent dynamics during execution and observation could be partially aligned with canonical correlation, indicating some similarity of the relationships among neural representations of different movements relative to one another during execution and observation. We also found that during action execution, mirror neurons showed consistent patterns of co-modulation both within and between sessions, but other non-mirror neurons that were modulated only during action execution and not during observation showed considerable variability of co-modulation.

Impairment of Neuronal Activity in the Dorsolateral Prefrontal Cortex Occurs Early in Parkinsonism.

Parkinson's disease (PD) is often characterized by altered rates and patterns of neuronal activity in the sensorimotor regions of the basal ganglia thalamocortical network. Little is known, however, regarding how neuronal activity in the executive control network of the brain changes in the parkinsonian condition. Investigate the impact of parkinsonism on neuronal activity in the dorsolateral prefrontal cortex (DLPFC), a key region in executive control, during a go/nogo reaching task. Using a within-subject design, single and multi-unit neuronal activity was recorded in the DLPFC of a nonhuman primate before and after the induction of mild parkinsonism using the neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP). Coincident with development of mild parkinsonian motor signs, there was a marked reduction in the percentage of DLPFC cells with significant task-related firing rate modulation during go and nogo conditions. These results suggest that DLPFC dysfunction may occur early in parkinsonism and contribute to cognitive impairments and disrupted executive function often observed in PD patients.

Hypoxia evokes a sequence of raphe-pontomedullary network operations for inspiratory drive amplification and gasping.

Hypoxia can trigger a sequence of breathing-related behaviors, from augmentation to apneusis to apnea and gasping. Gasping is an autoresuscitative behavior that, via large tidal volumes and altered intrathoracic pressure, can enhance coronary perfusion, carotid blood flow, and sympathetic activity, and thereby coordinate cardiac and respiratory functions. We tested the hypotheses that hypoxia-evoked gasps are amplified through a disinhibitory microcircuit within the inspiratory neuron chain and that this drive is distributed via an efference copy mechanism. This generates coordinated gasplike discharges concurrently in other circuits of the raphe-pontomedullary respiratory network. Data were obtained from six decerebrate, vagotomized, neuromuscularly blocked, and artificially ventilated adult cats. Arterial blood pressure, phrenic nerve activity, end-tidal CO2, and other parameters were monitored. Hypoxia was produced by ventilation with a gas mixture of 5% O2 in nitrogen. Neuron spike trains were recorded at multiple pontomedullary sites simultaneously and evaluated for firing rate modulations and short-timescale correlations indicative of functional connectivity. Experimental perturbations evoked reconfiguration of raphe-pontomedullary circuits during initial augmentation, apneusis and augmented bursts, apnea, and gasping. Functional connectivity, altered firing rates, efference copy of gasp drive, and coordinated incremental blood pressure increases support a distributed brain stem network model for amplification and broadcasting of inspiratory drive during autoresuscitative gasping. Gasping begins with a reduction in inhibition by expiratory neurons and an initial loss of inspiratory drive during hypoxic apnea and culminates in autoresuscitative efforts. NEW & NOTEWORTHY Severe hypoxia evokes a sequence of breathing-related behaviors culminating in gasping. We report firing rate modulations and short-timescale correlations in spike trains recorded simultaneously in the raphe-pontomedullary respiratory network during hypoxia. Our findings support a disinhibitory microcircuit and a distributed efference copy mechanism for amplification of gasping. Coordinated increments in blood pressure lead to a model for autoresuscitative bootstrapping of peripheral chemoreceptor reflexes, breathing, and sympathetic activity, complementing and extending prior work.

Dynamic spatial representation of self and others’ actions in the macaque frontal cortex

Modulation of neuronal firing rates by the spatial locations of physical objects is a widespread phenomenon in the brain. However, little is known about how neuronal responses to the actions of biological entities are spatially tuned and whether such spatially tuned responses are affected by social contexts. These issues are of key importance for understanding the neural basis of embodied social cognition, such as imitation and perspective-taking. Here, we show that spatial representation of actions can be dynamically changed depending on others' social relevance and agents of action. Monkeys performed a turn-taking choice task with a real monkey partner sitting face-to-face or a filmed partner in prerecorded videos. Three rectangular buttons (left, center, and right) were positioned in front of the subject and partner as their choice targets. We recorded from single neurons in two frontal nodes in the social brain, the ventral premotor cortex (PMv) and the medial prefrontal cortex (MPFC). When the partner was filmed rather than real, spatial preference for partner-actions was markedly diminished in MPFC, but not PMv, neurons. This social context-dependent modulation in the MPFC was also evident for self-actions. Strikingly, a subset of neurons in both areas switched their spatial preference between self-actions and partner-actions in a diametrically opposite manner. This observation suggests that these cortical areas are associated with coordinate transformation in ways consistent with an actor-centered perspective-taking coding scheme. The PMv may subserve such functions in context-independent manners, whereas the MPFC may do so primarily in social contexts.

Dissecting attention: Rate modulation vs. phase locking

In this issue of Neuron, Spyropoulos etal.1 report direct top-down modulation of neuronal firing rates underlying selective visual attentional modulation.

A large-scale neuronal network modelling study: Stimulus size modulates gamma oscillations in the primary visual cortex by long-range connections.

Stimulus size modulation of neuronal firing activity is a fundamental property of the primary visual cortex. Numerous biological experiments have shown that stimulus size modulation is affected by multiple factors at different spatiotemporal scales, but the exact pathways and mechanisms remain incompletely understood. In this paper, we establish a large-scale neuronal network model of primary visual cortex withlayer 2/3 to study how gamma oscillation properties are modulated by stimulus size and especially how long-range connections affect the modulation as realistic neuronal properties and spatial distributions of synaptic connections are considered. It is shown that long-range horizontal synaptic connections are sufficient to produce dimensional modulation of firing rates and gamma oscillations. In particular, with increasing grating stimulus size, the firing rate increases and then decreases, the peak frequency of gamma oscillations decreases and the spectral power increases. These are consistent with biological experimental observations. Furthermore, we explain in detail how the number and spatial distribution of long-range connections affect the size modulation of gamma oscillations by using the analysis of neuronal firing activity and synaptic current fluctuations. Our results provide a mechanism explanation for size modulation of gamma oscillations in the primary visual cortex and reveal the important and unique role played by long-range connections, which contributes to a deeper understanding of the cognitive function of gamma oscillations in visual cortex.

Theta Phase Entrainment of Single-Cell Spiking in Rat Somatosensory Barrel Cortex and Secondary Visual Cortex Is Enhanced during Multisensory Discrimination Behavior.

Phase entrainment of cells by theta oscillations is thought to globally coordinate the activity of cell assemblies across different structures, such as the hippocampus and neocortex. This coordination is likely required for optimal processing of sensory input during recognition and decision-making processes. In quadruple-area ensemble recordings from male rats engaged in a multisensory discrimination task, we investigated phase entrainment of cells by theta oscillations in areas along the corticohippocampal hierarchy: somatosensory barrel cortex (S1BF), secondary visual cortex (V2L), perirhinal cortex (PER), and dorsal hippocampus (dHC). Rats discriminated between two 3D objects presented in tactile-only, visual-only, or both tactile and visual modalities. During task engagement, S1BF, V2L, PER, and dHC LFP signals showed coherent theta-band activity. We found phase entrainment of single-cell spiking activity to locally recorded as well as hippocampal theta activity in S1BF, V2L, PER, and dHC. While phase entrainment of hippocampal spikes to local theta oscillations occurred during sustained epochs of task trials and was nonselective for behavior and modality, somatosensory and visual cortical cells were only phase entrained during stimulus presentation, mainly in their preferred modality (S1BF, tactile; V2L, visual), with subsets of cells selectively phase-entrained during cross-modal stimulus presentation (S1BF: visual; V2L: tactile). This effect could not be explained by modulations of firing rate or theta amplitude. Thus, hippocampal cells are phase entrained during prolonged epochs, while sensory and perirhinal neurons are selectively entrained during sensory stimulus presentation, providing a brief time window for coordination of activity.

An interpretable mechanism for grating-induced cross-inhibition and gamma oscillation based on a visual cortical neuronal network model

<abstract><p>Biological experiments targeting the mammalian primary visual cortex have shown that neuronal response to a preferred orientation grating is cross-inhibited by an orthogonal orientation mask grating. The plaid formed by the overlap of the two gratings not only causes a decrease in the neuronal firing rate but also shifts the gamma oscillation to a weaker oscillation at a higher frequency. The mechanism for the above phenomena is unclarified. In this paper, a large-scale cortical neuronal network model with biological details is constructed. In this model, two modes of connectivity that may contribute to cross-inhibition are considered: the thalamo-cortical feedforward pathway and the push-pull organization of cortical layer 4. Based on this model, the modulation of firing rate and gamma oscillation by a plaid stimulation are successfully reproduced, which is consistent with biological experiments and suggests that it is the thalamo-cortical feedforward pathway that leads to cross-inhibition. Furthermore, our analysis of the neuronal spike clusters and current fluctuations suggests that the push-pull organization leads to an increase in gamma frequency during the transition of visual stimuli from grating to plaid by modulating the source of synaptic inhibition in local neuronal populations. Such results will help to understand the visual processing under multi-input integration.</p></abstract>

A multivariate physiological model of vagus nerve signalling during metabolic challenges in anaesthetised rats for diabetes treatment

Objective. This study aims to develop a comprehensive decoding framework to create a multivariate physiological model of vagus nerve transmission that reveals the complex interactions between the nervous and metabolic systems. Approach. Vagus nerve activity was recorded in female Sprague-Dawley rats using gold hook microwires implanted around the left cervical vagus nerve. The rats were divided into three experimental cohorts (intact nerve, ligation nerve for recording afferent activation, and ligation for recording efferent activation) and metabolic challenges were administered to change glucose levels while recording the nerve activity. The decoding methodology involved various techniques, including continuous wavelet transformation, extraction of breathing rate (BR), and correlation of neural metrics with physiological signals. Main results. Decrease in glucose level was consistently negatively correlated with an increase in the firing activity of the intact vagus nerve that was found to be conveyed by both afferent and efferent pathways, with the afferent response being more similar to the one on the intact nerve. A larger variability was observed in the sensory and motor responses to hyperglycaemia. A novel strategy to extract the BR over time based on inter-burst-interval is also presented. The vagus afferent was found to encode breathing information through amplitude and firing rate modulation. Modulations of the signal amplitude were also observed due to changes in heart rate in the intact and efferent recordings, highlighting the parasympathetic control of the heart. Significance. The analytical framework presented in this study provides an integrative understanding that considers the relationship between metabolic, cardiac, and breathing signals and contributes to the development of a multivariable physiological model for the transmission of vagus nerve signals. This work progresses toward the development of closed-loop neuro-metabolic therapeutic systems for diabetes.

5-HT2AR and NMDAR psychedelics induce similar hyper-synchronous states in the rat cognitive-limbic cortex-basal ganglia system

The profound changes in perception and cognition induced by psychedelic drugs are thought to act on several levels, including increased glutamatergic activity, altered functional connectivity and an aberrant increase in high-frequency oscillations. To bridge these different levels of observation, we have here performed large-scale multi-structure recordings in freely behaving rats treated with 5-HT2AR psychedelics (LSD, DOI) and NMDAR psychedelics (ketamine, PCP). While interneurons and principal cells showed disparate firing rate modulations for the two classes of psychedelics, the local field potentials revealed a shared pattern of synchronized high-frequency oscillations in the ventral striatum and several cortical areas. Remarkably, the phase differences between structures were close to zero, corresponding to <1 ms delays. Likely, this hypersynchrony has major effects on the integration of information across neuronal systems and we propose that it is a key contributor to changes in perception and cognition during psychedelic drug use. Potentially, similar mechanisms could induce hallucinations and delusions in psychotic disorders and would constitute promising targets for new antipsychotic treatments.

Stable Neural Population Dynamics in the Regression Subspace for Continuous and Categorical Task Parameters in Monkeys.

Neural population dynamics provide a key computational framework for understanding information processing in the sensory, cognitive, and motor functions of the brain. They systematically depict complex neural population activity, dominated by strong temporal dynamics as trajectory geometry in a low-dimensional neural space. However, neural population dynamics are poorly related to the conventional analytical framework of single-neuron activity, the rate-coding regime that analyzes firing rate modulations using task parameters. To link the rate-coding and dynamic models, we developed a variant of state-space analysis in the regression subspace, which describes the temporal structures of neural modulations using continuous and categorical task parameters. In macaque monkeys, using two neural population datasets containing either of two standard task parameters, continuous and categorical, we revealed that neural modulation structures are reliably captured by these task parameters in the regression subspace as trajectory geometry in a lower dimension. Furthermore, we combined the classical optimal-stimulus response analysis (usually used in rate-coding analysis) with the dynamic model and found that the most prominent modulation dynamics in the lower dimension were derived from these optimal responses. Using those analyses, we successfully extracted geometries for both task parameters that formed a straight geometry, suggesting that their functional relevance is characterized as a unidimensional feature in their neural modulation dynamics. Collectively, our approach bridges neural modulation in the rate-coding model and the dynamic system, and provides researchers with a significant advantage in exploring the temporal structure of neural modulations for pre-existing datasets.

Spatial Attention Modulates Neuronal Interactions between Simple and Complex Cells in V1.

Visual perception is profoundly modulated by spatial attention, which can selectively prioritize goal-related information. Previous studies found spatial attention facilitated the efficacy of neuronal communication between visual cortices with hierarchical organizations. In the primary visual cortex (V1), there is also a hierarchical connection between simple (S) and complex (C) cells. We wonder whether and how spatial attention modulates neuronal communication within V1, especially for neuronal pairs with heterogeneous visual input. We simultaneously recorded the pairs' activity from macaque monkeys when they performed a spatial-attention-involved task, then applied likelihood-based Granger causality analysis to explore attentional modulation of neuronal interactions. First, a significant attention-related decrease in Granger causality was found in S-C pairs, which primarily displayed in the S-to-C feedforward connection. Second, the interaction strength of the feedforward connection was significantly higher than that of the feedback under attend toward (AT) conditions. Although information flow did not alter as the attentional focus shifted, the strength of communications between target- and distractor-stimuli-covered neurons differed only when attending to complex cells' receptive fields (RFs). Furthermore, pairs' communications depended on the attentional modulation of neurons' firing rates. Our findings demonstrate spatial attention does not induce specific information flow but rather amplifies directed communication within V1.

Beta band oscillations in the motor thalamus are modulated by visuomotor coordination in essential tremor patients.

Beta oscillations in sensorimotor structures contribute to the planning, sequencing, and stopping of movements, functions that are typically associated with the role of the basal ganglia. The presence of beta oscillations (13-30 Hz) in the cerebellar zone of the thalamus (the ventral intermediate nucleus - Vim) indicates that this rhythm may also be involved in cerebellar functions such as motor learning and visuomotor adaptation. To investigate the possible role of Vim beta oscillations in visuomotor coordination, we recorded local field potential (LFP) and multiunit activity from the Vim of essential tremor (ET) patients during neurosurgery for the implantation of deep brain stimulation (DBS) electrodes. Using a computer, patients performed a visuomotor adaptation task that required coordinating center-out movements with incongruent visual feedback imposed by inversion of the computer display. The results show that, in ET, Vim beta oscillations of the LFP were lower during the incongruent center-out task than during the congruent orientation. Vim firing rates increased significantly during periods of low beta power, particularly on approach to the peripheral target. In contrast, beta power in the subthalamic nucleus of Parkinson's disease (PD) patients did not differ significantly between the incongruent and the congruent orientation of the center-out task. The findings support the hypothesis that beta oscillations of the Vim are modulated by novel visuomotor tasks. The inverse relationship between the power of Vim-LFP beta oscillations and Vim firing rates suggest that the suppression of beta oscillations may facilitate information throughput to the thalamocortical circuit by modulation of Vim firing rates.

Shared Neural Activity But Distinct Neural Dynamics for Cognitive Control in Monkey Prefrontal and Parietal Cortex.

To better understand how prefrontal networks mediate forms of cognitive control disrupted in schizophrenia, we translated a variant of the AX continuous performance task that measures specific deficits in the human disease to 2 male monkeys and recorded neurons in PFC and parietal cortex during task performance. In the task, contextual information instructed by cue stimuli determines the response required to a subsequent probe stimulus. We found parietal neurons encoding the behavioral context instructed by cues that exhibited nearly identical activity to their prefrontal counterparts (Blackman et al., 2016). This neural population switched their preference for stimuli over the course of the trial depending on whether the stimuli signaled the need to engage cognitive control to override a prepotent response. Cues evoked visual responses that appeared in parietal neurons first, whereas population activity encoding contextual information instructed by cues was stronger and more persistent in PFC. Increasing cognitive control demand biased the representation of contextual information toward the PFC and augmented the temporal correlation of task-defined information encoded by neurons in the two areas. Oscillatory dynamics in local field potentials differed between cortical areas and carried as much information about task conditions as spike rates. We found that, at the single-neuron level, patterns of activity evoked by the task were nearly identical between the two cortical areas. Nonetheless, distinct population dynamics in PFC and parietal cortex were evident. suggesting differential contributions to cognitive control.SIGNIFICANCE STATEMENT We recorded neural activity in PFC and parietal cortex of monkeys performing a task that measures cognitive control deficits in schizophrenia. This allowed us to characterize computations performed by neurons in the two areas to support forms of cognitive control disrupted in the disease. Subpopulations of neurons in the two areas exhibited parallel modulations in firing rate; and as a result, all patterns of task-evoked activity were distributed between PFC and parietal cortex. This included the presence in both cortical areas of neurons reflecting proactive and reactive cognitive control dissociated from stimuli or responses in the task. However, differences in the timing, strength, synchrony, and correlation of information encoded by neural activity were evident, indicating differential contributions to cognitive control.

The neural representation of time distributed across multiple brain regions differs between implicit and explicit time demands.

Animals appear to possess an internal timer during action, based on the passage of time. However, the neural underpinnings of the perception of time, ranging from seconds to minutes, remain unclear. Herein, we considered the neural representation of time based on mounting evidence on the neural correlates of time perception. The passage of time in the brain is represented by two types of neural encoding as follows: (i) the modulation of firing rates in single neurons and (ii) the sequential activity in neural ensembles. Time-dependent neural activity reflects the relative time rather than the absolute time, similar to a clock. They emerge in multiple regions, including the hippocampus, medial and lateral entorhinal cortices, medial prefrontal cortex, and dorsal striatum. Moreover, they involve different brain regions, depending on an implicit or explicit event duration. Thus, the two types of internal timers distributed across multiple brain regions simultaneously engage in time perception, in response to implicit or explicit time demands.

Dynamic attention signalling in V4: Relation to fast-spiking/non-fast-spiking cell class and population coupling.

The computational role of a neuron during attention depends on its firing properties, neurotransmitter expression and functional connectivity. Neurons in the visual cortical area V4 are reliably engaged by selective attention but exhibit diversity in the effect of attention on firing rates and correlated variability. It remains unclear what specific neuronal properties shape these attention effects. In this study, we quantitatively characterised the distribution of attention modulation of firing rates across populations of V4 neurons. Neurons exhibited a continuum of time-varying attention effects. At one end of the continuum, neurons' spontaneous firing rates were slightly depressed with attention (compared to when unattended), whereas their stimulus responses were enhanced with attention. The other end of the continuum showed the converse pattern: attention depressed stimulus responses but increased spontaneous activity. We tested whether the particular pattern of time-varying attention effects that a neuron exhibited was related to the shape of their actions potentials (so-called 'fast-spiking' [FS] neurons have been linked to inhibition) and the strength of their coupling to the overall population. We found an interdependence among neural attention effects, neuron type and population coupling. In particular, we found neurons for which attention enhanced spontaneous activity but suppressed stimulus responses were less likely to be fast-spiking (more likely to be non-fast-spiking) and tended to have stronger population coupling, compared to neurons with other types of attention effects. These results add important information to our understanding of visual attention circuits at the cellular level.

The steady state and response to a periodic stimulation of the firing rate for a theta neuron with correlated noise

The stochastic activity of neurons is caused by various sources of correlated fluctuations and can be described in terms of simplified, yet biophysically grounded, integrate-and-fire models. One paradigmatic model is the quadratic integrate-and-fire model and its equivalent phase description by the theta neuron. Here we study the theta neuron model driven by a correlated Ornstein-Uhlenbeck noise and by periodic stimuli. We apply the matrix-continued-fraction method to the associated Fokker-Planck equation to develop an efficient numerical scheme to determine the stationary firing rate as well as the stimulus-induced modulation of the instantaneous firing rate. For the stationary case, we identify the conditions under which the firing rate decreases or increases by the effect of the colored noise and compare our results to existing analytical approximations for limit cases. For an additional periodic signal we demonstrate how the linear and nonlinear response terms can be computed and report resonant behavior for some of them. We extend the method to the case of two periodic signals, generally with incommensurable frequencies, and present a particular case for which a strong mixed response to both signals is observed, i.e. where the response to the sum of signals differs significantly from the sum of responses to the single signals. We provide Python code for our computational method: https://github.com/jannikfranzen/theta_neuron.

Selective modulation of cortical population dynamics during neuroprosthetic skill learning

Brain-machine interfaces (BMIs) provide a framework for studying how cortical population dynamics evolve over learning in a task in which the mapping between neural activity and behavior is precisely defined. Learning to control a BMI is associated with the emergence of coordinated neural dynamics in populations of neurons whose activity serves as direct input to the BMI decoder (direct subpopulation). While previous work shows differential modification of firing rate modulation in this population relative to a population whose activity was not directly input to the BMI decoder (indirect subpopulation), little is known about how learning-related changes in cortical population dynamics within these groups compare.To investigate this, we monitored both direct and indirect subpopulations as two macaque monkeys learned to control a BMI. We found that while the combined population increased coordinated neural dynamics, this increase in coordination was primarily driven by changes in the direct subpopulation. These findings suggest that motor cortex refines cortical dynamics by increasing neural variance throughout the entire population during learning, with a more pronounced coordination of firing activity in subpopulations that are causally linked to behavior.

A circuit mechanism for independent modulation of excitatory and inhibitory firing rates after sensory deprivation

Diverse interneuron subtypes shape sensory processing in mature cortical circuits. During development, sensory deprivation evokes powerful synaptic plasticity that alters circuitry, but how different inhibitory subtypes modulate circuit dynamics in response to this plasticity remains unclear. We investigate how deprivation-induced synaptic changes affect excitatory and inhibitory firing rates in a microcircuit model of the sensory cortex with multiple interneuron subtypes. We find that with a single interneuron subtype (parvalbumin-expressing [PV]), excitatory and inhibitory firing rates can only be comodulated-increased or decreased together. To explain the experimentally observed independent modulation, whereby one firing rate increases and the other decreases, requires strong feedback from a second interneuron subtype (somatostatin-expressing [SST]). Our model applies to the visual and somatosensory cortex, suggesting a general mechanism across sensory cortices. Therefore, we provide a mechanistic explanation for the differential role of interneuron subtypes in regulating firing rates, contributing to the already diverse roles they serve in the cortex.

State-Dependent Modulation of Activity in Distinct Layer 6 Corticothalamic Neurons in Barrel Cortex of Awake Mice.

Layer 6 corticothalamic (L6 CT) neurons are in a strategic position to control sensory input to the neocortex, yet we understand very little about their functions. Apart from studying their anatomic, physiological, and synaptic properties, most recent efforts have focused on the activity-dependent influences CT cells can exert on thalamic and cortical neurons through causal optogenetic manipulations. However, few studies have attempted to study them during behavior. To address this gap, we performed juxtacellular recordings from optogenetically identified CT neurons in whisker-related primary somatosensory cortex (wS1) of awake, head-fixed mice (either sex) free to rest quietly or self-initiate bouts of whisking and locomotion. We found a rich diversity of response profiles exhibited by CT cells. Their spiking patterns were either modulated by whisking-related behavior (∼28%) or not (∼72%). Whisking-responsive neurons exhibited both increases (activated-type) and decreases in firing rates (suppressed-type) that aligned with whisking onset better than locomotion. We also encountered responsive neurons with preceding modulations in firing rate before whisking onset. Overall, whisking better explained these changes in rates than overall changes in arousal. Whisking-unresponsive CT cells were generally quiet, with many having low spontaneous firing rates (sparse-type) and others being completely silent (silent-type). Remarkably, the sparse firing CT population preferentially spiked at the state transition point when pupil diameter constricted, and the mouse entered quiet wakefulness. Thus, our results demonstrate that L6 CT cells in wS1 show diverse spiking patterns, perhaps subserving distinct functional roles related to precisely timed responses during complex behaviors and transitions between discrete waking states.SIGNIFICANCE STATEMENT Layer 6 corticothalamic neurons provide a massive input to the sensory thalamus and local connectivity within cortex, but their role in thalamocortical processing remains unclear because of difficulty accessing and isolating their activity. Although several recent optogenetic studies reveal that the net influence of corticothalamic actions, suppression versus enhancement, depends critically on the rate these neurons fire, the factors that influence their spiking are poorly understood, particularly during wakefulness. Using the well-established Ntsr1-Cre line to target this elusive population in the whisker somatosensory cortex of awake mice, we found that corticothalamic neurons show diverse state-related responses and modulations in firing rate. These results suggest separate corticothalamic populations can differentially influence thalamocortical excitability during rapid state transitions in awake, behaving animals.