High-frequency burst spiking in layer 5 thick-tufted pyramids of rat primary somatosensory cortex encodes exploratory touch

Probing multiplexed basal dendritic computations using two-photon 3D holographic uncaging

Basal dendrites of layer 5 cortical pyramidal neurons exhibit Na+ and N-methyl-D-aspartate receptor (NMDAR) regenerative spikes and are uniquely poised to influence somatic output. Nevertheless, due to technical limitations, how multibranch basal dendritic integration shapes and enables multiplexed barcoding of synaptic streams remains poorly mapped. Here, we combine 3D two-photon holographic transmitter uncaging, whole-cell dynamic clamp, and biophysical modeling to reveal how synchronously activated synapses (distributed and clustered) across multiple basal dendritic branches are multiplexed under quiescent and in vivo-like conditions. While dendritic regenerative Na+ spikes promote millisecond somatic spike precision, distributed synaptic inputs and NMDAR spikes regulate gain. These concomitantly occurring dendritic nonlinearities enable multiplexed information transfer amid an ongoing noisy background, including under back-propagating voltage resets, by barcoding the axo-somatic spike structure. Our results unveil a multibranch dendritic integration framework in which dendritic nonlinearities are critical for multiplexing different spatial-temporal synaptic input patterns, enabling optimal feature binding.

Functional Dynamics and Selectivity of Two Parallel Corticocortical Pathways from Motor Cortex to Layer 5 Circuits in Somatosensory Cortex

In the rodent whisker system, active sensing and sensorimotor integration are mediated in part by the dynamic interactions between the motor cortex (M1) and somatosensory cortex (S1). However, understanding these dynamic interactions requires knowledge about the synapses and how specific neurons respond to their input. Here, we combined optogenetics, retrograde labeling, and electrophysiology to characterize the synaptic connections between M1 and layer 5 (L5) intratelencephalic (IT) and pyramidal tract (PT) neurons in S1 of mice (both sexes). We found that M1 synapses onto IT cells displayed modest short-term depression, whereas synapses onto PT neurons showed robust short-term facilitation. Despite M1 inputs to IT cells depressing, their slower kinetics resulted in summation and a response that increased during short trains. In contrast, summation was minimal in PT neurons due to the fast time course of their M1 responses. The functional consequences of this reduced summation, however, were outweighed by the strong facilitation at these M1 synapses, resulting in larger response amplitudes in PT neurons than IT cells during repetitive stimulation. To understand the impact of facilitating M1 inputs on PT output, we paired trains of inputs with single backpropagating action potentials, finding that repetitive M1 activation increased the probability of bursts in PT cells without impacting the time dependence of this coupling. Thus, there are two parallel but dynamically distinct systems of M1 synaptic excitation in L5 of S1, each defined by the short-term dynamics of its synapses, the class of postsynaptic neurons, and how the neurons respond to those inputs.

Functional dynamics and selectivity of two parallel corticocortical pathways from motor cortex to layer 5 circuits in somatosensory cortex

In the rodent whisker system, active sensing and sensorimotor integration are mediated in part by the dynamic interactions between the motor cortex (M1) and somatosensory cortex (S1). However, understanding these dynamic interactions requires knowledge about the synapses and how specific neurons respond to their input. Here, we combined optogenetics, retrograde labeling, and electrophysiology to characterize the synaptic connections between M1 and layer 5 (L5) intratelencephalic (IT) and pyramidal tract (PT) neurons in S1 of mice (both sexes). We found that M1 synapses onto IT cells displayed modest short-term depression, whereas synapses onto PT neurons showed robust short-term facilitation. Despite M1 inputs to IT cells depressing, their slower kinetics resulted in summation and a response that increased during short trains. In contrast, summation was minimal in PT neurons due to the fast time course of their M1 responses. The functional consequences of this reduced summation, however, were outweighed by the strong facilitation at these M1 synapses, resulting in larger response amplitudes in PT neurons than IT cells during repetitive stimulation. To understand the impact of facilitating M1 inputs on PT output, we paired trains of inputs with single backpropagating action potentials, finding that repetitive M1 activation increased the probability of bursts in PT cells without impacting the time-dependence of this coupling. Thus, there are two parallel but dynamically distinct systems of M1 synaptic excitation in L5 of S1, each defined by the short-term dynamics of its synapses, the class of postsynaptic neurons, and how the neurons respond to those inputs.

A Translaminar Spacetime Code Supports Touch-Evoked Traveling Waves

Linking sensory-evoked traveling waves to underlying circuit patterns is critical to understanding the neural basis of sensory perception. To form this link, we performed simultaneous electrophysiology and two-photon calcium imaging through transparent NeuroGrids and mapped touch-evoked cortical traveling waves and their underlying microcircuit dynamics. In awake mice, both passive and active whisker touch elicited traveling waves within and across barrels, with a fast early component followed by a variable late wave that lasted hundreds of milliseconds post-stimulus. Strikingly, late-wave dynamics were modulated by stimulus value and correlated with task performance. Mechanistically, the late wave component was i) modulated by motor feedback, ii) complemented by a sparse ensemble pattern across layer 2/3, which a balanced-state network model reconciled via inhibitory stabilization, and iii) aligned to regenerative Layer-5 apical dendritic Ca 2+ events. Our results reveal a translaminar spacetime pattern organized by cortical feedback in the sensory cortex that supports touch-evoked traveling waves.

GRAPHICAL ABSTRACT AND HIGHLIGHTS

Whisker touch evokes both early- and late-traveling waves in the barrel cortex over 100's of millisecondsReward reinforcement modulates wave dynamics Late wave emergence coincides with network sparsity in L23 and time-locked L5 dendritic Ca 2+ spikes Experimental and computational results link motor feedback to distinct translaminar spacetime patterns.

Circuits in the motor cortex explain oscillatory responses to transcranial magnetic stimulation

Transcranial magnetic stimulation (TMS) is a popular method used to investigate brain function. Stimulation over the motor cortex evokes muscle contractions known as motor evoked potentials (MEPs) and also high-frequency volleys of electrical activity measured in the cervical spinal cord. The physiological mechanisms of these experimentally derived responses remain unclear, but it is thought that the connections between circuits of excitatory and inhibitory neurons play a vital role. Using a spiking neural network model of the motor cortex, we explained the generation of waves of activity, so called 'I-waves', following cortical stimulation. The model reproduces a number of experimentally known responses including direction of TMS, increased inhibition, and changes in strength. Using populations of thousands of neurons in a model of cortical circuitry we showed that the cortex generated transient oscillatory responses without any tuning, and that neuron parameters such as refractory period and delays influenced the pattern and timing of those oscillations. By comparing our network with simpler, previously proposed circuits, we explored the contributions of specific connections and found that recurrent inhibitory connections are vital in producing later waves that significantly impact the production of motor evoked potentials in downstream muscles (Thickbroom, 2011). This model builds on previous work to increase our understanding of how complex circuitry of the cortex is involved in the generation of I-waves.

Network-neuron interactions underlying sensory responses of layer 5 pyramidal tract neurons in barrel cortex

Neurons in the cerebral cortex receive thousands of synaptic inputs per second from thousands of presynaptic neurons. How the dendritic location of inputs, their timing, strength, and presynaptic origin, in conjunction with complex dendritic physiology, impact the transformation of synaptic input into action potential (AP) output remains generally unknown for in vivo conditions. Here, we introduce a computational approach to reveal which properties of the input causally underlie AP output, and how this neuronal input-output computation is influenced by the morphology and biophysical properties of the dendrites. We demonstrate that this approach allows dissecting of how different input populations drive in vivo observed APs. For this purpose, we focus on fast and broadly tuned responses that pyramidal tract neurons in layer 5 (L5PTs) of the rat barrel cortex elicit upon passive single whisker deflections. By reducing a multi-scale model that we reported previously, we show that three features are sufficient to predict with high accuracy the sensory responses and receptive fields of L5PTs under these specific in vivo conditions: the count of active excitatory versus inhibitory synapses preceding the response, their spatial distribution on the dendrites, and the AP history. Based on these three features, we derive an analytically tractable description of the input-output computation of L5PTs, which enabled us to dissect how synaptic input from thalamus and different cell types in barrel cortex contribute to these responses. We show that the input-output computation is preserved across L5PTs despite morphological and biophysical diversity of their dendrites. We found that trial-to-trial variability in L5PT responses, and cell-to-cell variability in their receptive fields, are sufficiently explained by variability in synaptic input from the network, whereas variability in biophysical and morphological properties have minor contributions. Our approach to derive analytically tractable models of input-output computations in L5PTs provides a roadmap to dissect network-neuron interactions underlying L5PT responses across different in vivo conditions and for other cell types.

Holistic bursting cells store long-term memory in auditory cortex

The sensory neocortex has been suggested to be a substrate for long-term memory storage, yet which exact single cells could be specific candidates underlying such long-term memory storage remained neither known nor visible for over a century. Here, using a combination of day-by-day two-photon Ca2+ imaging and targeted single-cell loose-patch recording in an auditory associative learning paradigm with composite sounds in male mice, we reveal sparsely distributed neurons in layer 2/3 of auditory cortex emerged step-wise from quiescence into bursting mode, which then invariably expressed holistic information of the learned composite sounds, referred to as holistic bursting (HB) cells. Notably, it was not shuffled populations but the same sparse HB cells that embodied the behavioral relevance of the learned composite sounds, pinpointing HB cells as physiologically-defined single-cell candidates of an engram underlying long-term memory storage in auditory cortex.

In vivo spontaneous activity and coital-evoked inhibition of mouse accessory olfactory bulb output neurons

Little is known about estrous effects on brain microcircuits. We examined the accessory olfactory bulb (AOB) in vivo, in anesthetized naturally cycling females, as model microcircuit receiving coital somatosensory information. Whole-cell recordings demonstrate that output neurons are relatively hyperpolarized in estrus and unexpectedly fire high frequency bursts of action potentials. To mimic coitus, a calibrated artificial vagino-cervical stimulation (aVCS) protocol was devised. aVCS evoked stimulus-locked local field responses in the interneuron layer independent of estrous stage. The response is sensitive to α1-adrenergic receptor blockade, as expected since aVCS increases norepinephrine release in AOB. Intriguingly, only in estrus does aVCS inhibit AOB spike output. Estrus-specific output reduction coincides with prolonged aVCS activation of inhibitory interneurons. Accordingly, in estrus the AOB microcircuit sets the stage for coital stimulation to inhibit the output neurons, possibly via high frequency bursting-dependent enhancement of reciprocal synapse efficacy between inter- and output neurons.

Highly branched and complementary distributions of layer 5 and layer 6 auditory corticofugal axons in mouse

The auditory cortex exerts a powerful, yet heterogeneous, effect on subcortical targets. Auditory corticofugal projections emanate from layers 5 and 6 and have complementary physiological properties. While several studies suggested that layer 5 corticofugal projections branch widely, others suggested that multiple independent projections exist. Less is known about layer 6; no studies have examined whether the various layer 6 corticofugal projections are independent. Therefore, we examined branching patterns of layers 5 and 6 auditory corticofugal neurons, using the corticocollicular system as an index, using traditional and novel approaches. We confirmed that dual retrograde injections into the mouse inferior colliculus and auditory thalamus co-labeled subpopulations of layers 5 and 6 auditory cortex neurons. We then used an intersectional approach to relabel layer 5 or 6 corticocollicular somata and found that both layers sent extensive branches to multiple subcortical structures. Using a novel approach to separately label layers 5 and 6 axons in individual mice, we found that layers 5 and 6 terminal distributions partially spatially overlapped and that giant terminals were only found in layer 5-derived axons. Overall, the high degree of branching and complementarity in layers 5 and 6 axonal distributions suggest that corticofugal projections should be considered as 2 widespread systems, rather than collections of individual projections.

Sleep spindles in primates: Modeling the effects of distinct laminar thalamocortical connectivity in core, matrix, and reticular thalamic circuits

Sleep spindles are associated with the beginning of deep sleep and memory consolidation and are disrupted in schizophrenia and autism. In primates, distinct core and matrix thalamocortical (TC) circuits regulate sleep spindle activity through communications that are filtered by the inhibitory thalamic reticular nucleus (TRN); however, little is known about typical TC network interactions and the mechanisms that are disrupted in brain disorders. We developed a primate-specific, circuit-based TC computational model with distinct core and matrix loops that can simulate sleep spindles. We implemented novel multilevel cortical and thalamic mixing, and included local thalamic inhibitory interneurons, and direct layer 5 projections of variable density to TRN and thalamus to investigate the functional consequences of different ratios of core and matrix node connectivity contribution to spindle dynamics. Our simulations showed that spindle power in primates can be modulated based on the level of cortical feedback, thalamic inhibition, and engagement of model core versus matrix, with the latter having a greater role in spindle dynamics. The study of the distinct spatial and temporal dynamics of core-, matrix-, and mix-generated sleep spindles establishes a framework to study disruption of TC circuit balance underlying deficits in sleep and attentional gating seen in autism and schizophrenia.

Targeting operational regimes of interest in recurrent neural networks

Neural computations emerge from local recurrent neural circuits or computational units such as cortical columns that comprise hundreds to a few thousand neurons. Continuous progress in connectomics, electrophysiology, and calcium imaging require tractable spiking network models that can consistently incorporate new information about the network structure and reproduce the recorded neural activity features. However, for spiking networks, it is challenging to predict which connectivity configurations and neural properties can generate fundamental operational states and specific experimentally reported nonlinear cortical computations. Theoretical descriptions for the computational state of cortical spiking circuits are diverse, including the balanced state where excitatory and inhibitory inputs balance almost perfectly or the inhibition stabilized state (ISN) where the excitatory part of the circuit is unstable. It remains an open question whether these states can co-exist with experimentally reported nonlinear computations and whether they can be recovered in biologically realistic implementations of spiking networks. Here, we show how to identify spiking network connectivity patterns underlying diverse nonlinear computations such as XOR, bistability, inhibitory stabilization, supersaturation, and persistent activity. We establish a mapping between the stabilized supralinear network (SSN) and spiking activity which allows us to pinpoint the location in parameter space where these activity regimes occur. Notably, we find that biologically-sized spiking networks can have irregular asynchronous activity that does not require strong excitation-inhibition balance or large feedforward input and we show that the dynamic firing rate trajectories in spiking networks can be precisely targeted without error-driven training algorithms.

Malformations-related neocortical circuits in focal seizures

This review article gives an overview on the molecular, cellular and network mechanisms underlying focal seizures in neocortical networks with developmental malformations. Neocortical malformations comprise a large variety of structural abnormalities associated with epilepsy and other neurological and psychiatric disorders. Genetic or acquired disorders of neocortical cell proliferation, neuronal migration and/or programmed cell death may cause pathologies ranging from the expression of dysmorphic neurons and heterotopic cell clusters to abnormal layering and cortical misfolding. After providing a brief overview on the pathogenesis and structure of neocortical malformations in humans, animal models are discussed and how they contributed to our understanding on the mechanisms of neocortical hyperexcitability associated with developmental disorders. State-of-the-art molecular biological and electrophysiological techniques have been also used in humans and on resectioned neocortical tissue of epileptic patients and provide deep insights into the subcellular, cellular and network mechanisms contributing to focal seizures. Finally, a brief outlook is given how novel models and methods can shape translational research in the near future.

Neocortical layer 5 subclasses: From cellular properties to roles in behavior

Layer 5 (L5) serves as the main output layer of cortical structures, where long-range projecting pyramidal neurons broadcast the columnar output to other cortical and extracortical regions of the brain. L5 pyramidal neurons are grouped into two subclasses based on their projection targets; while intratelencephalic (IT) neurons project to cortical areas and the striatum, extratelencephalic (ET) neurons project to subcortical areas such as the thalamus, midbrain, and brainstem. Each L5 subclass possesses distinct morphological and electrophysiological properties and is incorporated into a unique synaptic network. Thanks to recent advances in genetic tools and methodologies, it has now become possible to distinguish between the two subclasses in the living brain. There is increasing evidence indicating that each subclass plays a unique role in sensory processing, decision-making, and learning. This review first summarizes the anatomical and physiological properties as well as the neuromodulation of IT and ET neurons in the rodent neocortex, and then reviews recent literature on their roles in sensory processing and rodent behavior. Our ultimate goal is to provide a comprehensive understanding of the role of each subclass in cortical function by examining their operational regimes based on their cellular properties.

Pareto optimality, economy-effectiveness trade-offs and ion channel degeneracy: improving population modelling for single neurons

Neurons encounter unavoidable evolutionary trade-offs between multiple tasks. They must consume as little energy as possible while effectively fulfilling their functions. Cells displaying the best performance for such multi-task trade-offs are said to be Pareto optimal, with their ion channel configurations underpinning their functionality. Ion channel degeneracy, however, implies that multiple ion channel configurations can lead to functionally similar behaviour. Therefore, instead of a single model, neuroscientists often use populations of models with distinct combinations of ionic conductances. This approach is called population (database or ensemble) modelling. It remains unclear, which ion channel parameters in the vast population of functional models are more likely to be found in the brain. Here we argue that Pareto optimality can serve as a guiding principle for addressing this issue by helping to identify the subpopulations of conductance-based models that perform best for the trade-off between economy and functionality. In this way, the high-dimensional parameter space of neuronal models might be reduced to geometrically simple low-dimensional manifolds, potentially explaining experimentally observed ion channel correlations. Conversely, Pareto inference might also help deduce neuronal functions from high-dimensional Patch-seq data. In summary, Pareto optimality is a promising framework for improving population modelling of neurons and their circuits.

Relationship between high-frequency activity in the cortical sensory and the motor hand areas, and their myelin content

BACKGROUND

The human primary sensory (S1) and primary motor (M1) hand areas feature high-frequency neuronal responses. Electrical nerve stimulation evokes high-frequency oscillations (HFO) at around 650 Hz in the contralateral S1. Likewise, transcranial magnetic stimulation (TMS) of M1 can evoke a series of descending volleys in the corticospinal pathway that can be detected non-invasively with a paired-pulse TMS protocol, called short interval intracortical facilitation (SICF). SICF features several peaks of facilitation of motor evoked potentials in contralateral hand muscles, which are separated by inter-peak intervals resembling HFO rhythmicity.

HYPOTHESIS

In this study, we tested the hypothesis that the individual expressions of HFO and SICF are tightly related to each other and to the regional myelin content in the sensorimotor cortex.

METHODS

In 24 healthy volunteers, we recorded HFO and SICF, and, in a subgroup of 20 participants, we mapped the cortical myelin content using the ratio between the T1- and T2-weighted MRI signal as read-out.

RESULTS

The individual frequencies and magnitudes of HFO and SICF curves were tightly correlated: the intervals between the first and second peak of cortical HFO and SICF showed a positive linear relationship (r = 0.703, p < 0.001), while their amplitudes were inversely related (r = -0.613, p = 0.001). The rhythmicity, but not the magnitude of the high-frequency responses, was related to the cortical myelin content: the higher the cortical myelin content, the shorter the inter-peak intervals of HFO and SICF.

CONCLUSION

The results confirm a tight functional relationship between high-frequency responses in S1 (i.e., HFO) and M1 (i.e., as measured with SICF). They also establish a link between the degree of regional cortical myelination and the expression of high-frequency responses in the human sensorimotor cortex, giving further the opportunity to infer their generators.

The guide to dendritic spikes of the mammalian cortex in vitro and in vivo

Half a century since their discovery by Llinás and colleagues, dendritic spikes have been observed in various neurons in different brain regions, from the neocortex and cerebellum to the basal ganglia. Dendrites exhibit a terrifically diverse but stereotypical repertoire of spikes, sometimes specific to subregions of the dendrite. Despite their prevalence, we only have a glimpse into their role in the behaving animal. This article aims to survey the full range of dendritic spikes found in excitatory and inhibitory neurons, compare them in vivo versus in vitro, and discuss new studies describing dendritic spikes in the human cortex. We focus on dendritic spikes in neocortical and hippocampal neurons and present a roadmap to identify and understand the broader role of dendritic spikes in single-cell computation.

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.

Neural burst codes disguised as rate codes

The burst coding hypothesis posits that the occurrence of sudden high-frequency patterns of action potentials constitutes a salient syllable of the neural code. Many neurons, however, do not produce clearly demarcated bursts, an observation invoked to rule out the pervasiveness of this coding scheme across brain areas and cell types. Here we ask how detrimental ambiguous spike patterns, those that are neither clearly bursts nor isolated spikes, are for neuronal information transfer. We addressed this question using information theory and computational simulations. By quantifying how information transmission depends on firing statistics, we found that the information transmitted is not strongly influenced by the presence of clearly demarcated modes in the interspike interval distribution, a feature often used to identify the presence of burst coding. Instead, we found that neurons having unimodal interval distributions were still able to ascribe different meanings to bursts and isolated spikes. In this regime, information transmission depends on dynamical properties of the synapses as well as the length and relative frequency of bursts. Furthermore, we found that common metrics used to quantify burstiness were unable to predict the degree with which bursts could be used to carry information. Our results provide guiding principles for the implementation of coding strategies based on spike-timing patterns, and show that even unimodal firing statistics can be consistent with a bivariate neural code.