Search

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. Author summary: Revealing how synaptic inputs drive action potential output is one of the major challenges in neuroscience research. An increasing number of approaches therefore seek to combine detailed measurements at synaptic, cellular and network scales into biologically realistic brain models. Indeed, such models have started to make empirically testable predictions about the inputs that underlie in vivo observed activity patterns. However, the enormous complexity of these models generally prevents the derivation of interpretable descriptions that explain how neurons transform synaptic input into action potential output, and how these input-output computations depend on synaptic, cellular and network properties. Here we introduce an approach to reveal input-output computations that neurons in the cerebral cortex perform upon sensory stimulation. For this purpose, we reduce a realistic multi-scale model of the rat barrel cortex to the minimal description that accounts for in vivo observed responses to whisker stimuli. Thereby, we identify the input-output computation that these cortical neurons perform under this in vivo condition, and we show that this computation is preserved across neurons despite morphological and biophysical diversity. Our approach provides analytically tractable and hence interpretable descriptions of neuronal input-output computations during specific in vivo conditions. [ABSTRACT FROM AUTHOR]

Recent advances in connectomics research enable the acquisition of increasing amounts of data about the connectivity patterns of neurons. How can we use this wealth of data to efficiently derive and test hypotheses about the principles underlying these patterns? A common approach is to simulate neuronal networks using a hypothesized wiring rule in a generative model and to compare the resulting synthetic data with empirical data. However, most wiring rules have at least some free parameters, and identifying parameters that reproduce empirical data can be challenging as it often requires manual parameter tuning. Here, we propose to use simulation-based Bayesian inference (SBI) to address this challenge. Rather than optimizing a fixed wiring rule to fit the empirical data, SBI considers many parametrizations of a rule and performs Bayesian inference to identify the parameters that are compatible with the data. It uses simulated data from multiple candidate wiring rule parameters and relies on machine learning methods to estimate a probability distribution (the 'posterior distribution over parameters conditioned on the data') that characterizes all data-compatible parameters. We demonstrate how to apply SBI in computational connectomics by inferring the parameters of wiring rules in an in silico model of the rat barrel cortex, given in vivo connectivity measurements. SBI identifies a wide range of wiring rule parameters that reproduce the measurements. We show how access to the posterior distribution over all data-compatible parameters allows us to analyze their relationship, revealing biologically plausible parameter interactions and enabling experimentally testable predictions. We further show how SBI can be applied to wiring rules at different spatial scales to quantitatively rule out invalid wiring hypotheses. Our approach is applicable to a wide range of generative models used in connectomics, providing a quantitative and efficient way to constrain model parameters with empirical connectivity data. Author summary: The brain is composed of an intricately connected network of cells—what are the principles that contribute to constructing these patterns of connectivity, and how? To answer these questions, amassing connectivity data alone is not enough. We must also be able to efficiently develop and test our ideas about the underlying connectivity principles. For example, we could simulate a hypothetical wiring rule like "neurons near each other are more likely to form connections" in a computational model and generate corresponding synthetic data. If the synthetic, simulated data resembles the real, measured data, then we have some confidence that our hypotheses might be correct. However, the proposed wiring rules usually have unknown parameters that we need to "tune" such that simulated data matches the measurements. The central challenge thus lies in finding all the potential parametrizations of a wiring rule that can reproduce the measured data, as this process is often idiosyncratic and labor-intensive. To tackle this challenge, we introduce an approach combining computational modeling in connectomics, deep learning, and Bayesian statistical inference to automatically infer a probability distribution over the model parameters likely to explain the data. We demonstrate our approach by inferring the parameters of a wiring rule in a detailed model of the rat barrel cortex and find that the inferred distribution identifies multiple data-compatible model parameters, reveals biologically plausible parameter interactions, and allows us to make experimentally testable predictions. [ABSTRACT FROM AUTHOR]

Diversity of cell-types that collectively shape the cortical microcircuit ensures the necessary computational richness to orchestrate a wide variety of behaviors. The information content embedded in spiking activity of identified cell-types remain unclear to a large extent. Here, we recorded spike responses upon whisker touch of anatomically identified excitatory cell-types in primary somatosensory cortex in naive, untrained rats. We find major differences across layers and cell-types. The temporal structure of spontaneous spiking contains high-frequency bursts (≥100 Hz) in all morphological cell-types but a significant increase upon whisker touch is restricted to layer L5 thick-tufted pyramids (L5tts) and thus provides a distinct neurophysiological signature. We find that whisker touch can also be decoded from L5tt bursting, but not from other cell-types. We observed high-frequency bursts in L5tts projecting to different subcortical regions, including thalamus, midbrain and brainstem. We conclude that bursts in L5tts allow accurate coding and decoding of exploratory whisker touch. In order to investigate the information encoded by spiking activity in different neuronal cell types in the primary somatosensory cortex, de Kock et al performed electrophysiological recordings in untrained rats. They demonstrated that an increase in high-frequency burst spiking in thick tufted pyramids in layer 5 of the cortex allow accurate encoding of exploratory whisker touch. [ABSTRACT FROM AUTHOR]

Abstract: Singing behavior in the adult male zebra finch is dependent upon the activity of a cortical region known as HVC (proper name). The vast majority of HVC projection neurons send primary axons to either the downstream premotor nucleus RA (robust nucleus of the arcopallium, or primary motor cortex) or Area X (basal ganglia), which play important roles in song production or song learning, respectively. In addition to these long‐range outputs, HVC neurons also send local axon collaterals throughout that nucleus. Despite their implications for a range of circuit models, these local processes have never been completely reconstructed. Here, we use in vivo single‐neuron Neurobiotin fills to examine 40 projection neurons across 31 birds with somatic positions distributed across HVC. We show that HVC(RA) and HVC(X) neurons have categorically distinct dendritic fields. Additionally, these cell classes send axon collaterals that are either restricted to a small portion of HVC (“local neurons”) or broadly distributed throughout the entire nucleus (“broadcast neurons”). Overall, these processes within HVC offer a structural basis for significant local processing underlying behaviorally relevant population activity. [ABSTRACT FROM AUTHOR]

The cytoarchitectonic subdivision of the neocortex into six layers is often used to describe the organization of the cortical circuitry, sensory-evoked signal flow or cortical functions. However, each layer comprises neuronal cell types that have different genetic, functional and/or structural properties. Here, we reanalyze structural data from some of our recent work in the posterior-medial barrel-subfield of the vibrissal part of rat primary somatosensory cortex (vS1). We quantify the degree to which somata, dendrites and axons of the 10 major excitatory cell types of the cortex are distributed with respect to the cytoarchitectonic organization of vS1.We show that within each layer, somata of multiple cell types intermingle, but that each cell type displays dendrite and axon distributions that are aligned to specific cytoarchitectonic landmarks. The resultant quantification of the structural composition of each layer in terms of the cell type-specific number of somata, dendritic and axonal path lengths will aid future studies to bridge between layer- and cell type-specific analyses. [ABSTRACT FROM AUTHOR]

Pyramidal tract neurons (PTs) represent the major output cell type of the neocortex. To investigate principles of how the results of cortical processing are broadcasted to different downstream targets thus requires experimental approaches, which provide access to the in vivo electrophysiology of PTs, whose subcortical target regions are identified. On the example of rat barrel cortex (vS1), we illustrate that retrograde tracer injections into multiple subcortical structures allow identifying the long-range axonal targets of individual in vivo recorded PTs. Here we report that soma depth and dendritic path lengths within each cortical layer of vS1, as well as spiking patterns during both periods of ongoing activity and during sensory stimulation, reflect the respective subcortical target regions of PTs. We show that these cellular properties result in a structure–function parameter space that allows predicting a PT’s subcortical target region, without the need to inject multiple retrograde tracers. [ABSTRACT FROM AUTHOR]

Maintaining an appropriate balance between excitation and inhibition is critical for neuronal information processing. Cortical neurons can cell-autonomously adjust the inhibition they receive to individual levels of excitatory input, but the underlying mechanisms are unclear. We describe that Ste20-like kinase (SLK) mediates cell-autonomous regulation of excitation-inhibition balance in the thalamocortical feedforward circuit, but not in the feedback circuit. This effect is due to regulation of inhibition originating from parvalbumin-expressing interneurons, while inhibition via somatostatin-expressing interneurons is unaffected. Computational modeling shows that this mechanism promotes stable excitatory-inhibitory ratios across pyramidal cells and ensures robust and sparse coding. Patch-clamp RNA sequencing yields genes differentially regulated by SLK knockdown, as well as genes associated with excitation-inhibition balance participating in transsynaptic communication and cytoskeletal dynamics. These data identify a mechanism for cell-autonomous regulation of a specific inhibitory circuit that is critical to ensure that a majority of cortical pyramidal cells participate in information coding. [Display omitted] • SLK regulates excitation-inhibition balance cell-autonomously • SLK in cortical neurons regulates feedforward but not feedback inhibition • SLK selectively regulates inhibition by parvalbumin-expressing interneurons Royero et al. identify a mechanism relying on Ste20-like kinase that allows single cortical neurons to cell-autonomously adjust feedforward inhibition they receive to the cell-specific levels of excitatory input. They propose that this mechanism is critical to ensure that a majority of cortical pyramidal cells participate in information coding. [ABSTRACT FROM AUTHOR]

Cortical inhibitory interneurons (INs) are subdivided into a variety of morphologically and functionally specialized cell types. How the respective specific properties translate into mechanisms that regulate sensory-evoked responses of pyramidal neurons (PNs) remains unknown. Here, we investigated how INs located in cortical layer 1 (L1) of rat barrel cortex affect whisker-evoked responses of L2 PNs. To do so we combined in vivo electrophysiology and morphological reconstructions with computational modeling. We show that whisker-evoked membrane depolarization in L2 PNs arises from highly specialized spatiotemporal synaptic input patterns. Temporally L1 INs and L2--5 PNs provide near synchronous synaptic input. Spatially synaptic contacts from L1 INs target distal apical tuft dendrites, whereas PNs primarily innervate basal and proximal apical dendrites. Simulations of such constrained synaptic input patterns predicted that inactivation of L1 INs increases trial-to-trial variability of whisker-evoked responses in L2 PNs. The in silico predictions were confirmed in vivo by L1-specific pharmacological manipulations. We present a mechanism--consistent with the theory of distal dendritic shunting--that can regulate the robustness of sensory-evoked responses in PNs without affecting response amplitude or latency. [ABSTRACT FROM AUTHOR]

Glycoprotein-deleted rabies virus (RABV ∆G) is a powerful tool for the analysis of neural circuits. Here, we demonstrate the utility of an anterograde RABV ∆G variant for novel neuroanatomical approaches involving either bulk or sparse neuronal populations. This technology exploits the unique features of RABV ∆G vectors, namely autonomous, rapid high-level expression of transgenes, and limited cytotoxicity. Our vector permits the unambiguous long-range and fine-scale tracing of the entire axonal arbor of individual neurons throughout the brain. Notably, this level of labeling can be achieved following infection with a single viral particle. The vector is effective over a range of ages (>14 months) aiding the studies of neurodegenerative disorders or aging, and infects numerous cell types in all brain regions tested. Lastly, it can also be readily combined with retrograde RABV ∆G variants. Together with other modern technologies, this tool provides new possibilities for the investigation of the anatomy and physiology of neural circuits. [ABSTRACT FROM AUTHOR]

Sensory-evoked signal flow, at cellular and network levels, is primarily determined by the synaptic wiring of the underlying neuronal circuitry. Measurements of synaptic innervation, connection probabilities and subcellular organization of synaptic inputs are thus among the most active fields of research in contemporary neuroscience. Methods to measure these quantities range from electrophysiological recordings over reconstructions of dendrite-axon overlap at light-microscopic levels to dense circuit reconstructions of small volumes at electron-microscopic resolution. However, quantitative and complete measurements at subcellular resolution and mesoscopic scales to obtain all local and long-range synaptic in/outputs for any neuron within an entire brain region are beyond present methodological limits. Here, we present a novel concept, implemented within an interactive software environment called NeuroNet, which allows (i) integration of sparsely sampled (sub)cellular morphological data into an accurate anatomical reference frame of the brain region(s) of interest, (ii) up-scaling to generate an average dense model of the neuronal circuitry within the respective brain region(s) and (iii) statistical measurements of synaptic innervation between all neurons within the model. We illustrate our approach by generating a dense average model of the entire rat vibrissal cortex, providing the required anatomical data, and illustrate how to measure synaptic innervation statistically. Comparing our results with data from paired recordings in vitro and in vivo, as well as with reconstructions of synaptic contact sites at light- and electron-microscopic levels, we find that our in silico measurements are in line with previous results. [ABSTRACT FROM AUTHOR]

Neuroanatomical analysis, such as classification of cell types, depends on reliable reconstruction of large numbers of complete 3D dendrite and axon morphologies. At present, the majority of neuron reconstructions are obtained from preparations in a single tissue slice in vitro, thus suffering from cut off dendrites and, more dramatically, cut off axons. In general, axons can innervate volumes of several cubic millimeters and may reach path lengths of tens of centimeters. Thus, their complete reconstruction requires in vivo labeling, histological sectioning and imaging of large fields of view. Unfortunately, anisotropic background conditions across such large tissue volumes, as well as faintly labeled thin neurites, result in incomplete or erroneous automated tracings and even lead experts to make annotation errors during manual reconstructions. Consequently, tracing reliability renders the major bottleneck for reconstructing complete 3D neuron morphologies. Here, we present a novel set of tools, integrated into a software environment named 'Filament Editor', for creating reliable neuron tracings from sparsely labeled in vivo datasets. The Filament Editor allows for simultaneous visualization of complex neuronal tracings and image data in a 3D viewer, proof-editing of neuronal tracings, alignment and interconnection across sections, and morphometric analysis in relation to 3D anatomical reference structures. We illustrate the functionality of the Filament Editor on the example of in vivo labeled axons and demonstrate that for the exemplary dataset the final tracing results after proof-editing are independent of the expertise of the human operator. [ABSTRACT FROM AUTHOR]

The neurons in the cerebral cortex are not randomly interconnected. This specificity in wiring can result from synapse formation mechanisms that connect neurons, depending on their electrical activity and genetically defined identity. Here, we report that the morphological properties of the neurons provide an additional prominent source by which wiring specificity emerges in cortical networks. This morphologically determined wiring specificity reflects similarities between the neurons' axo-dendritic projections patterns, the packing density, and the cellular diversity of the neuropil. The higher these three factors are, the more recurrent is the topology of the network. Conversely, the lower these factors are, the more feedforward is the network's topology. These principles predict the empirically observed occurrences of clusters of synapses, cell type-specific connectivity patterns, and nonrandom network motifs. Thus, we demonstrate that wiring specificity emerges in the cerebral cortex at subcellular, cellular, and network scales from the specific morphological properties of its neuronal constituents. [Display omitted] • Neuronal network architectures reflect the morphologies of their constituents • Morphology predicts nonrandom connectivity from subcellular to network scales • Morphology predicts connectivity patterns consistent with those observed empirically • Neuron morphology is a major source for wiring specificity in the cerebral cortex Udvary et al. reveal four basic principles by which the morphological properties of the neurons shape the specific architecture of the networks they form. These principles can account for nonrandom connectivity patterns that are observed empirically between the neurons in the cerebral cortex at subcellular, cellular, and network scales. [ABSTRACT FROM AUTHOR]

The three-dimensional (3D) structure of neural circuits is commonly studied by reconstructing individual or small groups of neurons in separate preparations. Investigation of structural organization principles or quantification of dendritic and axonal innervation thus requires integration of many reconstructed morphologies into a common reference frame. Here we present a standardized 3D model of the rat vibrissal cortex and introduce an automated registration tool that allows for precise placement of single neuron reconstructions. We (1) developed an automated image processing pipeline to reconstruct 3D anatomical landmarks, i.e., the barrels in Layer 4, the pia and white matter surfaces and the blood vessel pattern from high-resolution images, (2) quantified these landmarks in 12 different rats, (3) generated an average 3D model of the vibrissal cortex and (4) used rigid transformations and stepwise linear scaling to register 94 neuron morphologies, reconstructed from in vivo stainings, to the standardized cortex model. We find that anatomical landmarks vary substantially across the vibrissal cortex within an individual rat. In contrast, the 3D layout of the entire vibrissal cortex remains remarkably preserved across animals. This allows for precise registration of individual neuron reconstructions with approximately 30 μm accuracy. Our approach could be used to reconstruct and standardize other anatomically defined brain areas and may ultimately lead to a precise digital reference atlas of the rat brain. [ABSTRACT FROM AUTHOR]

How does the brain compute? Answering this question necessitates neuronal connectomes, annotated graphs of aU synaptic connections within defined brain areas. Further, understanding the energetics of the brain's computations require vascular graphs. The a embly of a connectomc requires sensitive hardware tools to measure neuronal and neurovascular features in all three dimensions, as weU as software and machine learning for data analysis and visualization. We present the state of the art on the reconstruction of circuits and vasculature that link brain anatomy and function. Analysis at the scale of tens of nanometers yields connections between identified neurons, while analysis at the micrometer scale yields probabilistic rules of connection between neurons and exact vascular connectivity. [ABSTRACT FROM AUTHOR]

The cortical output layer 5 contains two excitatory cell types, slender- and thick-tufted neurons. In rat vibrissal cortex, slender-tufted neurons carry motion and phase information during active whisking, but remain inactive after passive whisker touch. In contrast, thick-tufted neurons reliably increase spiking preferably after passive touch. By reconstructing the 3D patterns of intracortical axon projections from individual slender- and thick-tufted neurons, filled in vivo with biocytin, we were able to identify cell type-specific intracortical circuits that may encode whisker motion and touch. Individual slender-tufted neurons showed elaborate and dense innervation of supragranular layers of large portions of the vibrissal area (total length, 86.8 ± 5.5 mm). During active whisking, these long-range projections may modulate and phase-lock the membrane potential of dendrites in layers 2 and 3 to the whisking cycle. Thick-tufted neurons with soma locations intermingling with those of slender-tufted ones display less dense intracortical axon projections (total length, 31.6 ± 14.3 mm) that are primarily confined to infragranular layers. Based on anatomical reconstructions and previous measurements of spiking, we put forward the hypothesis that thick-tufted neurons in rat vibrissal cortex receive input of whisker motion from slender-tufted neurons onto their apical tuft dendrites and input of whisker touch from thalamic neurons onto their basal dendrites. During tactile-driven behavior, such as object location, near-coincident input from these two pathways may result in increased spiking activity of thick-tufted neurons and thus enhanced signaling to their subcortical targets. [ABSTRACT FROM AUTHOR]

We present a novel approach for deconvolution of 3D image stacks of cortical tissue taken by mosaic/optical-sectioning technology, using a transmitted light brightfield microscope. Mosaic/optical-sectioning offers the possibility of imaging large volumes (e.g. from cortical sections) on a millimetre scale at sub-micrometre resolution. However, a blurred contribution from out-of-focus light results in an image quality that usually prohibits 3D quantitative analysis. Such quantitative analysis is only possible after deblurring by deconvolution. The resulting image quality is strongly dependent on how accurate the point spread function used for deconvolution resembles the properties of the imaging system. Since direct measurement of the true point spread function is laborious and modelled point spread functions usually deviate from measured ones, we present a method of optimizing the microscope until it meets almost ideal imaging conditions. These conditions are validated by measuring the aberration function of the microscope and tissue using a Shack-Hartmann sensor. The analysis shows that cortical tissue from rat brains embedded in Mowiol and imaged by an oil-immersion objective can be regarded as having a homogeneous index of refraction. In addition, the amount of spherical aberration that is caused by the optics or the specimen is relatively low. Consequently the image formation is simplified to refraction between the embedding and immersion medium and to 3D diffraction at the finite entrance pupil of the objective. The resulting model point spread function is applied to the image stacks by linear or iterative deconvolution algorithms. For the presented dataset of large 3D images the linear approach proves to be superior. The linear deconvolution yields a significant improvement in signal-to-noise ratio and resolution. This novel approach allows a quantitative analysis of the cortical image stacks such as the reconstruction of biocytin-stained neuronal dendrites and axons. [ABSTRACT FROM AUTHOR]