Reconstruction and Simulation of Neocortical Microcircuitry.

We present a first-draft digital reconstruction of the microcircuitry of somatosensory cortex of juvenile rat. The reconstruction uses cellular and synaptic organizing principles to algorithmically reconstruct detailed anatomy and physiology from sparse experimental data. An objective anatomical method defines a neocortical volume of 0.29 ± 0.01 mm(3) containing ~31,000 neurons, and patch-clamp studies identify 55 layer-specific morphological and 207 morpho-electrical neuron subtypes. When digitally reconstructed neurons are positioned in the volume and synapse formation is restricted to biological bouton densities and numbers of synapses per connection, their overlapping arbors form ~8 million connections with ~37 million synapses. Simulations reproduce an array of in vitro and in vivo experiments without parameter tuning. Additionally, we find a spectrum of network states with a sharp transition from synchronous to asynchronous activity, modulated by physiological mechanisms. The spectrum of network states, dynamically reconfigured around this transition, supports diverse information processing strategies. PAPERCLIP: VIDEO ABSTRACT.

The cellular coding of temperature in the mammalian cortex.

Temperature is a fundamental sensory modality separate from touch, with dedicated receptor channels and primary afferent neurons for cool and warm1-3. Unlike for other modalities, however, the cortical encoding of temperature remains unknown, with very few cortical neurons reported that respond to non-painful temperature, and the presence of a 'thermal cortex' is debated4-8. Here, using widefield and two-photon calcium imaging in the mouse forepaw system, we identify cortical neurons that respond to cooling and/or warming with distinct spatial and temporal response properties. We observed a representation of cool, but not warm, in the primary somatosensory cortex, but cool and warm in the posterior insular cortex (pIC). The representation of thermal information in pIC is robust and somatotopically arranged, and reversible manipulations show a profound impact on thermal perception. Despite being positioned along the same one-dimensional sensory axis, the encoding of cool and that of warm are distinct, both in highly and broadly tuned neurons. Together, our results show that pIC contains the primary cortical representation of skin temperature and may help explain how the thermal system generates sensations of cool and warm.

Neuronal Circuits in Barrel Cortex for Whisker Sensory Perception.

The array of whiskers on the snout provides rodents with tactile sensory information relating to the size, shape and texture of objects in their immediate environment. Rodents can use their whiskers to detect stimuli, distinguish textures, locate objects and navigate. Important aspects of whisker sensation are thought to result from neuronal computations in the whisker somatosensory cortex (wS1). Each whisker is individually represented in the somatotopic map of wS1 by an anatomical unit named a 'barrel' (hence also called barrel cortex). This allows precise investigation of sensory processing in the context of a well-defined map. Here, we first review the signaling pathways from the whiskers to wS1, and then discuss current understanding of the various types of excitatory and inhibitory neurons present within wS1. Different classes of cells can be defined according to anatomical, electrophysiological and molecular features. The synaptic connectivity of neurons within local wS1 microcircuits, as well as their long-range interactions and the impact of neuromodulators, are beginning to be understood. Recent technological progress has allowed cell-type-specific connectivity to be related to cell-type-specific activity during whisker-related behaviors. An important goal for future research is to obtain a causal and mechanistic understanding of how selected aspects of tactile sensory information are processed by specific types of neurons in the synaptically connected neuronal networks of wS1 and signaled to downstream brain areas, thus contributing to sensory-guided decision-making.

Cortical processing of flexible and context-dependent sensorimotor sequences.

The brain generates complex sequences of movements that can be flexibly configured based on behavioural context or real-time sensory feedback1, but how this occurs is not fully understood. Here we developed a 'sequence licking' task in which mice directed their tongue to a target that moved through a series of locations. Mice could rapidly branch the sequence online based on tactile feedback. Closed-loop optogenetics and electrophysiology revealed that the tongue and jaw regions of the primary somatosensory (S1TJ) and motor (M1TJ) cortices2 encoded and controlled tongue kinematics at the level of individual licks. By contrast, the tongue 'premotor' (anterolateral motor) cortex3-10 encoded latent variables including intended lick angle, sequence identity and progress towards the reward that marked successful sequence execution. Movement-nonspecific sequence branching signals occurred in the anterolateral motor cortex and M1TJ. Our results reveal a set of key cortical areas for flexible and context-informed sequence generation.

Pyramidal neuron subtype diversity governs microglia states in the neocortex.

Microglia are specialized macrophages in the brain parenchyma that exist in multiple transcriptional states and reside within a wide range of neuronal environments1-4. However, how and where these states are generated remains poorly understood. Here, using the mouse somatosensory cortex, we demonstrate that microglia density and molecular state acquisition are determined by the local composition of pyramidal neuron classes. Using single-cell and spatial transcriptomic profiling, we unveil the molecular signatures and spatial distributions of diverse microglia populations and show that certain states are enriched in specific cortical layers, whereas others are broadly distributed throughout the cortex. Notably, conversion of deep-layer pyramidal neurons to an alternate class identity reconfigures the distribution of local, layer-enriched homeostatic microglia to match the new neuronal niche. Leveraging the transcriptional diversity of pyramidal neurons in the neocortex, we construct a ligand-receptor atlas describing interactions between individual pyramidal neuron subtypes and microglia states, revealing rules of neuron-microglia communication. Our findings uncover a fundamental role for neuronal diversity in instructing the acquisition of microglia states as a potential mechanism for fine-tuning neuroimmune interactions within the cortical local circuitry.

Maturation and circuit integration of transplanted human cortical organoids.

Self-organizing neural organoids represent a promising in vitro platform with which to model human development and disease1-5. However, organoids lack the connectivity that exists in vivo, which limits maturation and makes integration with other circuits that control behaviour impossible. Here we show that human stem cell-derived cortical organoids transplanted into the somatosensory cortex of newborn athymic rats develop mature cell types that integrate into sensory and motivation-related circuits. MRI reveals post-transplantation organoid growth across multiple stem cell lines and animals, whereas single-nucleus profiling shows progression of corticogenesis and the emergence of activity-dependent transcriptional programs. Indeed, transplanted cortical neurons display more complex morphological, synaptic and intrinsic membrane properties than their in vitro counterparts, which enables the discovery of defects in neurons derived from individuals with Timothy syndrome. Anatomical and functional tracings show that transplanted organoids receive thalamocortical and corticocortical inputs, and in vivo recordings of neural activity demonstrate that these inputs can produce sensory responses in human cells. Finally, cortical organoids extend axons throughout the rat brain and their optogenetic activation can drive reward-seeking behaviour. Thus, transplanted human cortical neurons mature and engage host circuits that control behaviour. We anticipate that this approach will be useful for detecting circuit-level phenotypes in patient-derived cells that cannot otherwise be uncovered.

Temporal controls over inter-areal cortical projection neuron fate diversity.

Interconnectivity between neocortical areas is critical for sensory integration and sensorimotor transformations1-6. These functions are mediated by heterogeneous inter-areal cortical projection neurons (ICPN), which send axon branches across cortical areas as well as to subcortical targets7-9. Although ICPN are anatomically diverse10-14, they are molecularly homogeneous15, and how the diversity of their anatomical and functional features emerge during development remains largely unknown. Here we address this question by linking the connectome and transcriptome in developing single ICPN of the mouse neocortex using a combination of multiplexed analysis of projections by sequencing16,17 (MAPseq, to identify single-neuron axonal projections) and single-cell RNA sequencing (to identify corresponding gene expression). Focusing on neurons of the primary somatosensory cortex (S1), we reveal a protracted unfolding of the molecular and functional differentiation of motor cortex-projecting ([Formula: see text]) ICPN compared with secondary somatosensory cortex-projecting ([Formula: see text]) ICPN. We identify SOX11 as a temporally differentially expressed transcription factor in [Formula: see text] versus [Formula: see text] ICPN. Postnatal manipulation of SOX11 expression in S1 impaired sensorimotor connectivity and disrupted selective exploratory behaviours in mice. Together, our results reveal that within a single cortical area, different subtypes of ICPN have distinct postnatal paces of molecular differentiation, which are subsequently reflected in distinct circuit connectivities and functions. Dynamic differences in the expression levels of a largely generic set of genes, rather than fundamental differences in the identity of developmental genetic programs, may thus account for the emergence of intra-type diversity in cortical neurons.

The human somatosensory cortex contributes to the encoding of newly learned movements.

Recent studies have indicated somatosensory cortex involvement in motor learning and retention. However, the nature of its contribution is unknown. One possibility is that the somatosensory cortex is transiently engaged during movement. Alternatively, there may be durable learning-related changes which would indicate sensory participation in the encoding of learned movements. These possibilities are dissociated by disrupting the somatosensory cortex following learning, thus targeting learning-related changes which may have occurred. If changes to the somatosensory cortex contribute to retention, which, in effect, means aspects of newly learned movements are encoded there, disruption of this area once learning is complete should lead to an impairment. Participants were trained to make movements while receiving rotated visual feedback. The primary motor cortex (M1) and the primary somatosensory cortex (S1) were targeted for continuous theta-burst stimulation, while stimulation over the occipital cortex served as a control. Retention was assessed using active movement reproduction, or recognition testing, which involved passive movements produced by a robot. Disruption of the somatosensory cortex resulted in impaired motor memory in both tests. Suppression of the motor cortex had no impact on retention as indicated by comparable retention levels in control and motor cortex conditions. The effects were learning specific. When stimulation was applied to S1 following training with unrotated feedback, movement direction, the main dependent variable, was unaltered. Thus, the somatosensory cortex is part of a circuit that contributes to retention, consistent with the idea that aspects of newly learned movements, possibly learning-updated sensory states (new sensory targets) which serve to guide movement, may be encoded there.

Evolutionarily conserved neural responses to affective touch in monkeys transcend consciousness and change with age.

Affective touch-a slow, gentle, and pleasant form of touch-activates a different neural network than which is activated during discriminative touch in humans. Affective touch perception is enabled by specialized low-threshold mechanoreceptors in the skin with unmyelinated fibers called C tactile (CT) afferents. These CT afferents are conserved across mammalian species, including macaque monkeys. However, it is unknown whether the neural representation of affective touch is the same across species and whether affective touch's capacity to activate the hubs of the brain that compute socioaffective information requires conscious perception. Here, we used functional MRI to assess the preferential activation of neural hubs by slow (affective) vs. fast (discriminative) touch in anesthetized rhesus monkeys (Macaca mulatta). The insula, anterior cingulate cortex (ACC), amygdala, and secondary somatosensory cortex were all significantly more active during slow touch relative to fast touch, suggesting homologous activation of the interoceptive-allostatic network across primate species during affective touch. Further, we found that neural responses to affective vs. discriminative touch in the insula and ACC (the primary cortical hubs for interoceptive processing) changed significantly with age. Insula and ACC in younger animals differentiated between slow and fast touch, while activity was comparable between conditions for aged monkeys (equivalent to >70 y in humans). These results, together with prior studies establishing conserved peripheral nervous system mechanisms of affective touch transduction, suggest that neural responses to affective touch are evolutionarily conserved in monkeys, significantly impacted in old age, and do not necessitate conscious experience of touch.

Supra-orbital whiskers act as wind-sensing antennae in rats.

We know little about mammalian anemotaxis or wind sensing. Recently, however, Hartmann and colleagues showed whisker-based anemotaxis in rats. To investigate how whiskers sense airflow, we first tracked whisker tips in anesthetized rats under low (0.5 m/s) and high (1.5 m/s) airflow. Whisker tips showed increasing movement from low to high airflow conditions, with all whisker tips moving during high airflow. Low airflow conditions-most similar to naturally occurring wind stimuli-engaged whisker tips differentially. Most whiskers moved little, but the long supra-orbital (lSO) whisker showed maximal displacement, followed by the α, ß, and A1 whiskers. The lSO whisker differs from other whiskers in its exposed dorsal position, upward bending, length and thin diameter. Ex vivo extracted lSO whiskers also showed exceptional airflow displacement, suggesting whisker-intrinsic biomechanics mediate the unique airflow-sensitivity. Micro computed tomography (micro-CT) revealed that the ring-wulst-the follicle structure receiving the most sensitive afferents-was more complete/closed in the lSO, and other wind-sensitive whiskers, than in non-wind-sensitive whiskers, suggesting specialization of the supra-orbital for omni-directional sensing. We localized and targeted the cortical supra-orbital whisker representation in simultaneous Neuropixels recordings with D/E-row whisker barrels. Responses to wind-stimuli were stronger in the supra-orbital whisker representation than in D/E-row barrel cortex. We assessed the behavioral significance of whiskers in an airflow-sensing paradigm. We observed that rats spontaneously turn towards airflow stimuli in complete darkness. Selective trimming of wind-responsive whiskers diminished airflow turning responses more than trimming of non-wind-responsive whiskers. Lidocaine injections targeted to supra-orbital whisker follicles also diminished airflow turning responses compared to control injections. We conclude that supra-orbital whiskers act as wind antennae.

Monosynaptic trans-collicular pathways link mouse whisker circuits to integrate somatosensory and motor cortical signals.

The superior colliculus (SC), a conserved midbrain node with extensive long-range connectivity throughout the brain, is a key structure for innate behaviors. Descending cortical pathways are increasingly recognized as central control points for SC-mediated behaviors, but how cortico-collicular pathways coordinate SC activity at the cellular level is poorly understood. Moreover, despite the known role of the SC as a multisensory integrator, the involvement of the SC in the somatosensory system is largely unexplored in comparison to its involvement in the visual and auditory systems. Here, we mapped the connectivity of the whisker-sensitive region of the SC in mice with trans-synaptic and intersectional tracing tools and in vivo electrophysiology. The results reveal a novel trans-collicular connectivity motif in which neurons in motor- and somatosensory cortices impinge onto the brainstem-SC-brainstem sensory-motor arc and onto SC-midbrain output pathways via only one synapse in the SC. Intersectional approaches and optogenetically assisted connectivity quantifications in vivo reveal convergence of motor and somatosensory cortical input on individual SC neurons, providing a new framework for sensory-motor integration in the SC. More than a third of the cortical recipient neurons in the whisker SC are GABAergic neurons, which include a hitherto unknown population of GABAergic projection neurons targeting thalamic nuclei and the zona incerta. These results pinpoint a whisker region in the SC of mice as a node for the integration of somatosensory and motor cortical signals via parallel excitatory and inhibitory trans-collicular pathways, which link cortical and subcortical whisker circuits for somato-motor integration.

Recurrent interactions in local cortical circuits.

Most cortical synapses are local and excitatory. Local recurrent circuits could implement amplification, allowing pattern completion and other computations1-4. Cortical circuits contain subnetworks that consist of neurons with similar receptive fields and increased connectivity relative to the network average5,6. Cortical neurons that encode different types of information are spatially intermingled and distributed over large brain volumes5-7, and this complexity has hindered attempts to probe the function of these subnetworks by perturbing them individually8. Here we use computational modelling, optical recordings and manipulations to probe the function of recurrent coupling in layer 2/3 of the mouse vibrissal somatosensory cortex during active tactile discrimination. A neural circuit model of layer 2/3 revealed that recurrent excitation enhances sensory signals by amplification, but only for subnetworks with increased connectivity. Model networks with high amplification were sensitive to damage: loss of a few members of the subnetwork degraded stimulus encoding. We tested this prediction by mapping neuronal selectivity7 and photoablating9,10 neurons with specific selectivity. Ablation of a small proportion of layer 2/3 neurons (10-20, less than 5% of the total) representing touch markedly reduced responses in the spared touch representation, but not in other representations. Ablations most strongly affected neurons with stimulus responses that were similar to those of the ablated population, which is also consistent with network models. Recurrence among cortical neurons with similar selectivity therefore drives input-specific amplification during behaviour.

Value-guided remapping of sensory cortex by lateral orbitofrontal cortex.

Adaptive behaviour crucially depends on flexible decision-making, which in mammals relies on the frontal cortex, specifically the orbitofrontal cortex (OFC)1-9. How OFC encodes decision variables and instructs sensory areas to guide adaptive behaviour are key open questions. Here we developed a reversal learning task for head-fixed mice, monitored the activity of neurons of the lateral OFC using two-photon calcium imaging and investigated how OFC dynamically interacts with primary somatosensory cortex (S1). Mice learned to discriminate 'go' from 'no-go' tactile stimuli10,11 and adapt their behaviour upon reversal of stimulus-reward contingency ('rule switch'). Imaging individual neurons longitudinally across all behavioural phases revealed a distinct engagement of S1 and lateral OFC, with S1 neural activity reflecting initial task learning, whereas lateral OFC neurons responded saliently and transiently to the rule switch. We identified direct long-range projections from lateral OFC to S1 that can feed this activity back to S1 as value prediction error. This top-down signal updated sensory representations in S1 by functionally remapping responses in a subpopulation of neurons that was sensitive to reward history. Functional remapping crucially depended on top-down feedback as chemogenetic silencing of lateral OFC neurons disrupted reversal learning, as well as plasticity in S1. The dynamic interaction of lateral OFC with sensory cortex thus implements computations critical for value prediction that are history dependent and error based, providing plasticity essential for flexible decision-making.

Interareal Synaptic Inputs Underlying Whisking-Related Activity in the Primary Somatosensory Barrel Cortex.

Body movements influence brain-wide neuronal activities. In the sensory cortex, thalamocortical bottom-up inputs and motor-sensory top-down inputs are thought to affect the dynamics of membrane potentials (Vm ) of neurons and change their processing of sensory information during movements. However, direct perturbation of the axons projecting to the sensory cortex from other remote areas during movements has remained unassessed, and therefore the interareal circuits generating motor-related signals in sensory cortices remain unclear. Using a Gi/o -coupled opsin, eOPN3, we here inhibited interareal signals incoming to the whisker primary somatosensory barrel cortex (wS1) of awake male mice and tested their effects on whisking-related changes in neuronal activities in wS1. Spontaneous whisking in air induced the changes in spike rates of a subset of wS1 neurons, which were accompanied by depolarization and substantial reduction of slow-wave oscillatory fluctuations of Vm Despite an extensive innervation, inhibition of inputs from the whisker primary motor cortex (wM1) to wS1 did not alter the spike rates and Vm dynamics of wS1 neurons during whisking. In contrast, inhibition of axons from the whisker-related thalamus (wTLM) and the whisker secondary somatosensory cortex (wS2) to wS1 largely attenuated the whisking-related supra- and sub-threshold Vm dynamics of wS1 neurons. Notably, silencing inputs from wTLM markedly decreased the modulation depth of whisking phase-tuned neurons in wS1, while inhibiting wS2 inputs did not impact the whisking variable tuning of wS1 neurons. Thus, sensorimotor integration in wS1 during spontaneous whisking is predominantly facilitated by direct synaptic inputs from wTLM and wS2 rather than from wM1.

Thalamocortical Dynamics during Rapid Eye Movement Sleep in the Mouse Somatosensory Pathway.

Rapid eye movement (REM) sleep, also referred to as paradoxical sleep for the striking resemblance of its electroencephalogram (EEG) to the one observed in wakefulness, is characterized by the occurrence of transient events such as limb twitches or facial and rapid eye movements. Here, we investigated the local activity of the primary somatosensory or barrel cortex (S1) in naturally sleeping head-fixed male mice during REM. Through local field potential recordings, we uncovered local appearances of spindle waves in the barrel cortex during REM concomitant with strong delta power, challenging the view of a wakefulness-like activity in REM. We further performed extra- and intracellular recordings of thalamic cells in head-fixed mice. Our data show high-frequency thalamic bursts of spikes and subthreshold spindle oscillations in approximately half of the neurons of the ventral posterior medial nucleus which further confirmed the thalamic origin of local cortical spindles in S1 in REM. Cortical spindle oscillations were suppressed, while thalamus spike firing increased, associated with rapid mouse whisker movements and S1 cortical activity transitioned to an activated state. During REM, the sensory thalamus and barrel cortex therefore alternate between high (wake-like) and low (non-REM sleep-like) activation states, potentially providing a neuronal substrate for mnemonic processes occurring during this paradoxical sleep stage.

Balancing the Senses: Electrophysiological Responses Reveal the Interplay between Somatosensory and Visual Processing During Body-Related Multisensory Conflict.

In the study of bodily awareness, the predictive coding theory has revealed that our brain continuously modulates sensory experiences to integrate them into a unitary body representation. Indeed, during multisensory illusions (e.g., the rubber hand illusion, RHI), the synchronous stroking of the participant's concealed hand and a fake visible one creates a visuotactile conflict, generating a prediction error. Within the predictive coding framework, through sensory processing modulation, prediction errors are solved, inducing participants to feel as if touches originated from the fake hand, thus ascribing the fake hand to their own body. Here, we aimed to address sensory processing modulation under multisensory conflict, by disentangling somatosensory and visual stimuli processing that are intrinsically associated during the illusion induction. To this aim, we designed two EEG experiments, in which somatosensory- (SEPs; Experiment 1; N = 18; F = 10) and visual-evoked potentials (VEPs; Experiment 2; N = 18; F = 9) were recorded in human males and females following the RHI. Our results show that, in both experiments, ERP amplitude is significantly modulated in the illusion as compared with both control and baseline conditions, with a modality-dependent diametrical pattern showing decreased SEP amplitude and increased VEP amplitude. Importantly, both somatosensory and visual modulations occur in long-latency time windows previously associated with tactile and visual awareness, thus explaining the illusion of perceiving touch at the sight location. In conclusion, we describe a diametrical modulation of somatosensory and visual processing as the neural mechanism that allows maintaining a stable body representation, by restoring visuotactile congruency under the occurrence of multisensory conflicts.

Learning-related congruent and incongruent changes of excitation and inhibition in distinct cortical areas.

Excitatory and inhibitory neurons in diverse cortical regions are likely to contribute differentially to the transformation of sensory information into goal-directed motor plans. Here, we investigate the relative changes across mouse sensorimotor cortex in the activity of putative excitatory and inhibitory neurons-categorized as regular spiking (RS) or fast spiking (FS) according to their action potential (AP) waveform-comparing before and after learning of a whisker detection task with delayed licking as perceptual report. Surprisingly, we found that the whisker-evoked activity of RS versus FS neurons changed in opposite directions after learning in primary and secondary whisker motor cortices, while it changed similarly in primary and secondary orofacial motor cortices. Our results suggest that changes in the balance of excitation and inhibition in local circuits concurrent with changes in the long-range synaptic inputs in distinct cortical regions might contribute to performance of delayed sensory-to-motor transformation.

The tuning of tuning: How adaptation influences single cell information transfer.

Sensory neurons reconstruct the world from action potentials (spikes) impinging on them. To effectively transfer information about the stimulus to the next processing level, a neuron needs to be able to adapt its working range to the properties of the stimulus. Here, we focus on the intrinsic neural properties that influence information transfer in cortical neurons and how tightly their properties need to be tuned to the stimulus statistics for them to be effective. We start by measuring the intrinsic information encoding properties of putative excitatory and inhibitory neurons in L2/3 of the mouse barrel cortex. Excitatory neurons show high thresholds and strong adaptation, making them fire sparsely and resulting in a strong compression of information, whereas inhibitory neurons that favour fast spiking transfer more information. Next, we turn to computational modelling and ask how two properties influence information transfer: 1) spike-frequency adaptation and 2) the shape of the IV-curve. We find that a subthreshold (but not threshold) adaptation, the 'h-current', and a properly tuned leak conductance can increase the information transfer of a neuron, whereas threshold adaptation can increase its working range. Finally, we verify the effect of the IV-curve slope in our experimental recordings and show that excitatory neurons form a more heterogeneous population than inhibitory neurons. These relationships between intrinsic neural features and neural coding that had not been quantified before will aid computational, theoretical and systems neuroscientists in understanding how neuronal populations can alter their coding properties, such as through the impact of neuromodulators. Why the variability of intrinsic properties of excitatory neurons is larger than that of inhibitory ones is an exciting question, for which future research is needed.

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.

What does the mean mean? A simple test for neuroscience.

Trial-averaged metrics, e.g. tuning curves or population response vectors, are a ubiquitous way of characterizing neuronal activity. But how relevant are such trial-averaged responses to neuronal computation itself? Here we present a simple test to estimate whether average responses reflect aspects of neuronal activity that contribute to neuronal processing. The test probes two assumptions implicitly made whenever average metrics are treated as meaningful representations of neuronal activity: Reliability: Neuronal responses repeat consistently enough across trials that they convey a recognizable reflection of the average response to downstream regions.Behavioural relevance: If a single-trial response is more similar to the average template, it is more likely to evoke correct behavioural responses. We apply this test to two data sets: (1) Two-photon recordings in primary somatosensory cortices (S1 and S2) of mice trained to detect optogenetic stimulation in S1; and (2) Electrophysiological recordings from 71 brain areas in mice performing a contrast discrimination task. Under the highly controlled settings of Data set 1, both assumptions were largely fulfilled. In contrast, the less restrictive paradigm of Data set 2 met neither assumption. Simulations predict that the larger diversity of neuronal response preferences, rather than higher cross-trial reliability, drives the better performance of Data set 1. We conclude that when behaviour is less tightly restricted, average responses do not seem particularly relevant to neuronal computation, potentially because information is encoded more dynamically. Most importantly, we encourage researchers to apply this simple test of computational relevance whenever using trial-averaged neuronal metrics, in order to gauge how representative cross-trial averages are in a given context.