Generation of dense statistical connectomes from sparse morphological data

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.

Translating single-neuron axonal reconstructions into meso-scale connectivity statistics in the mouse somatosensory thalamus

Characterizing the connectomic and morphological diversity of thalamic neurons is key for better understanding how the thalamus relays sensory inputs to the cortex. The recent public release of complete single-neuron morphological reconstructions enables the analysis of previously inaccessible connectivity patterns from individual neurons. Here we focus on the Ventral Posteromedial (VPM) nucleus and characterize the full diversity of 257 VPM neurons, obtained by combining data from the MouseLight and Braintell projects. Neurons were clustered according to their most dominantly targeted cortical area and further subdivided by their jointly targeted areas. We obtained a 2D embedding of morphological diversity using the dissimilarity between all pairs of axonal trees. The curved shape of the embedding allowed us to characterize neurons by a 1-dimensional coordinate. The coordinate values were aligned both with the progression of soma position along the dorsal-ventral and lateral-medial axes and with that of axonal terminals along the posterior-anterior and medial-lateral axes, as well as with an increase in the number of branching points, distance from soma and branching width. Taken together, we have developed a novel workflow for linking three challenging aspects of connectomics, namely the topography, higher order connectivity patterns and morphological diversity, with VPM as a test-case. The workflow is linked to a unified access portal that contains the morphologies and integrated with 2D cortical flatmap and subcortical visualization tools. The workflow and resulting processed data have been made available in Python, and can thus be used for modeling and experimentally validating new hypotheses on thalamocortical connectivity.

Scaling up cell-counting efforts in neuroscience through semi-automated methods

Quantifying how the cellular composition of brain regions vary across development, aging, sex, and disease, is crucial in experimental neuroscience, and the accuracy of different counting methods is continuously debated. Due to the tedious nature of most counting procedures, studies are often restricted to one or a few brain regions. Recently, there have been considerable methodological advances in combining semi-automated feature extraction with brain atlases for cell quantification. Such methods hold great promise for scaling up cell-counting efforts. However, little focus has been paid to how these methods should be implemented and reported to support reproducibility. Here, we provide an overview of practices for conducting and reporting cell counting in mouse and rat brains, showing that critical details for interpretation are typically lacking. We go on to discuss how novel methods may increase efficiency and reproducibility of cell counting studies. Lastly, we provide practical recommendations for researchers planning cell counting.

Simulation-based inference for efficient identification of generative models in computational connectomics

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.

Circuit-selective cell-autonomous regulation of inhibition in pyramidal neurons by Ste20-like kinase

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.

Structured sampling of olfactory input by the fly mushroom body

Associative memory formation and recall in the fruit fly Drosophila melanogaster is subserved by the mushroom body (MB). Upon arrival in the MB, sensory information undergoes a profound transformation from broadly tuned and stereotyped odorant responses in the olfactory projection neuron (PN) layer to narrowly tuned and nonstereotyped responses in the Kenyon cells (KCs). Theory and experiment suggest that this transformation is implemented by random connectivity between KCs and PNs. However, this hypothesis has been challenging to test, given the difficulty of mapping synaptic connections between large numbers of brain-spanning neurons. Here, we used a recent whole-brain electron microscopy volume of the adult fruit fly to map PN-to-KC connectivity at synaptic resolution. The PN-KC connectome revealed unexpected structure, with preponderantly food-responsive PN types converging at above-chance levels on downstream KCs. Axons of the overconvergent PN types tended to arborize near one another in the MB main calyx, making local KC dendrites more likely to receive input from those types. Overconvergent PN types preferentially co-arborize and connect with dendrites of αβ and α'β' KC subtypes. Computational simulation of the observed network showed degraded discrimination performance compared with a random network, except when all signal flowed through the overconvergent, primarily food-responsive PN types. Additional theory and experiment will be needed to fully characterize the impact of the observed non-random network structure on associative memory formation and recall.

The impact of neuron morphology on cortical network architecture

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.

•

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

Computational synthesis of cortical dendritic morphologies

Neuronal morphologies provide the foundation for the electrical behavior of neurons, the connectomes they form, and the dynamical properties of the brain. Comprehensive neuron models are essential for defining cell types, discerning their functional roles, and investigating brain-disease-related dendritic alterations. However, a lack of understanding of the principles underlying neuron morphologies has hindered attempts to computationally synthesize morphologies for decades. We introduce a synthesis algorithm based on a topological descriptor of neurons, which enables the rapid digital reconstruction of entire brain regions from few reference cells. This topology-guided synthesis generates dendrites that are statistically similar to biological reconstructions in terms of morpho-electrical and connectivity properties and offers a significant opportunity to investigate the links between neuronal morphology and brain function across different spatiotemporal scales. Synthesized cortical networks based on structurally altered dendrites associated with diverse brain pathologies revealed principles linking branching properties to the structure of large-scale networks.

Statistical Emulation of Neural Simulators: Application to Neocortical L2/3 Large Basket Cells

Many scientific systems are studied using computer codes that simulate the phenomena of interest. Computer simulation enables scientists to study a broad range of possible conditions, generating large quantities of data at a faster rate than the laboratory. Computer models are widespread in neuroscience, where they are used to mimic brain function at different levels. These models offer a variety of new possibilities for the neuroscientist, but also numerous challenges, such as: where to sample the input space for the simulator, how to make sense of the data that is generated, and how to estimate unknown parameters in the model. Statistical emulation can be a valuable complement to simulator-based research. Emulators are able to mimic the simulator, often with a much smaller computational burden and they are especially valuable for parameter estimation, which may require many simulator evaluations. This work compares different statistical models that address these challenges, and applies them to simulations of neocortical L2/3 large basket cells, created and run with the NEURON simulator in the context of the European Human Brain Project. The novelty of our approach is the use of fast empirical emulators, which have the ability to accelerate the optimization process for the simulator and to identify which inputs (in this case, different membrane ion channels) are most influential in affecting simulated features. These contributions are complementary, as knowledge of the important features can further improve the optimization process. Subsequent research, conducted after the process is completed, will gain efficiency by focusing on these inputs.

Single cortical neurons as deep artificial neural networks

Utilizing recent advances in machine learning, we introduce a systematic approach to characterize neurons' input/output (I/O) mapping complexity. Deep neural networks (DNNs) were trained to faithfully replicate the I/O function of various biophysical models of cortical neurons at millisecond (spiking) resolution. A temporally convolutional DNN with five to eight layers was required to capture the I/O mapping of a realistic model of a layer 5 cortical pyramidal cell (L5PC). This DNN generalized well when presented with inputs widely outside the training distribution. When NMDA receptors were removed, a much simpler network (fully connected neural network with one hidden layer) was sufficient to fit the model. Analysis of the DNNs' weight matrices revealed that synaptic integration in dendritic branches could be conceptualized as pattern matching from a set of spatiotemporal templates. This study provides a unified characterization of the computational complexity of single neurons and suggests that cortical networks therefore have a unique architecture, potentially supporting their computational power.

Neural anatomy and optical microscopy (NAOMi) simulation for evaluating calcium imaging methods

BACKGROUND

The past decade has seen a multitude of new in vivo functional imaging methodologies. However, the lack of ground-truth comparisons or evaluation metrics makes the large-scale, systematic validation vital to the continued development and use of optical microscopy impossible.

NEW-METHOD

We provide a new framework for evaluating two-photon microscopy methods via in silico Neural Anatomy and Optical Microscopy (NAOMi) simulation. Our computationally efficient model generates large anatomical volumes of mouse cortex, simulates neural activity, and incorporates optical propagation and scanning to create realistic calcium imaging datasets.

RESULTS

We verify NAOMi simulations against in vivo two-photon recordings from mouse cortex. We leverage this in silico ground truth to directly compare different segmentation algorithms and optical designs. We find modern segmentation algorithms extract strong neural time-courses comparable to estimation using oracle spatial information, but with an increase in the false positive rate. Comparison between optical setups demonstrate improved resilience to motion artifacts in sparsely labeled samples using Bessel beams, increased signal-to-noise ratio and cell-count using low numerical aperture Gaussian beams and nuclear GCaMP, and more uniform spatial sampling with temporal focusing versus multi-plane imaging.

COMPARISON WITH EXISTING METHODS

NAOMi is a first-of-its kind framework for assessing optical imaging modalities. Existing methods are either anatomical simulations or do not address functional imaging. Thus there is no competing method for simulating realistic functional optical microscopy data.

CONCLUSIONS

By leveraging the rich accumulated knowledge of neural anatomy and optical physics, we provide a powerful new tool to assess and develop important methods in neural imaging.

Cellular connectomes as arbiters of local circuit models in the cerebral cortex

With the availability of cellular-resolution connectivity maps, connectomes, from the mammalian nervous system, it is in question how informative such massive connectomic data can be for the distinction of local circuit models in the mammalian cerebral cortex. Here, we investigated whether cellular-resolution connectomic data can in principle allow model discrimination for local circuit modules in layer 4 of mouse primary somatosensory cortex. We used approximate Bayesian model selection based on a set of simple connectome statistics to compute the posterior probability over proposed models given a to-be-measured connectome. We find that the distinction of the investigated local cortical models is faithfully possible based on purely structural connectomic data with an accuracy of more than 90%, and that such distinction is stable against substantial errors in the connectome measurement. Furthermore, mapping a fraction of only 10% of the local connectome is sufficient for connectome-based model distinction under realistic experimental constraints. Together, these results show for a concrete local circuit example that connectomic data allows model selection in the cerebral cortex and define the experimental strategy for obtaining such connectomic data.

A simple method defines 3D morphology and axon projections of filled neurons in a small CNS volume: Steps toward understanding functional network circuitry

BACKGROUND

Fundamental to understanding neuronal network function is defining neuron morphology, location, properties, and synaptic connectivity in the nervous system. A significant challenge is to reconstruct individual neuron morphology and connections at a whole CNS scale and bring together functional and anatomical data to understand the whole network.

NEW METHOD

We used a PC controlled micropositioner to hold a fixed whole mount of Xenopus tadpole CNS and replace the stage on a standard microscope. This allowed direct recording in 3D coordinates of features and axon projections of one or two neurons dye-filled during whole-cell recording to study synaptic connections. Neuron reconstructions were normalised relative to the ventral longitudinal axis of the nervous system. Coordinate data were stored as simple text files.

RESULTS

Reconstructions were at 1 μm resolution, capturing axon lengths in mm. The output files were converted to SWC format and visualised in 3D reconstruction software NeuRomantic. Coordinate data are tractable, allowing correction for histological artefacts. Through normalisation across multiple specimens we could infer features of network connectivity of mapped neurons of different types.

COMPARISON WITH EXISTING METHODS

Unlike other methods using fluorescent markers and utilising large-scale imaging, our method allows direct acquisition of 3D data on neurons whose properties and synaptic connections have been studied using whole-cell recording.

CONCLUSIONS

This method can be used to reconstruct neuron 3D morphology and follow axon projections in the CNS. After normalisation to a common CNS framework, inferences on network connectivity at a whole nervous system scale contribute to network modelling to understand CNS function.

Database of literature derived cellular measurements from the murine basal ganglia

Quantitative measurements and descriptive statistics of different cellular elements in the brain are typically published in journal articles as text, tables, and example figures, and represent an important basis for the creation of biologically constrained computational models, design of intervention studies, and comparison of subject groups. Such data can be challenging to extract from publications and difficult to normalise and compare across studies, and few studies have so far attempted to integrate quantitative information available in journal articles. We here present a database of quantitative information about cellular parameters in the frequently studied murine basal ganglia. The database holds a curated and normalised selection of currently available data collected from the literature and public repositories, providing the most comprehensive collection of quantitative neuroanatomical data from the basal ganglia to date. The database is shared as a downloadable resource from the EBRAINS Knowledge Graph (https://kg.ebrains.eu), together with a workflow that allows interested researchers to update and expand the database with data from future reports.

Anatomy and Physiology of Macaque Visual Cortical Areas V1, V2, and V5/MT: Bases for Biologically Realistic Models

The cerebral cortex of primates encompasses multiple anatomically and physiologically distinct areas processing visual information. Areas V1, V2, and V5/MT are conserved across mammals and are central for visual behavior. To facilitate the generation of biologically accurate computational models of primate early visual processing, here we provide an overview of over 350 published studies of these three areas in the genus Macaca, whose visual system provides the closest model for human vision. The literature reports 14 anatomical connection types from the lateral geniculate nucleus of the thalamus to V1 having distinct layers of origin or termination, and 194 connection types between V1, V2, and V5, forming multiple parallel and interacting visual processing streams. Moreover, within V1, there are reports of 286 and 120 types of intrinsic excitatory and inhibitory connections, respectively. Physiologically, tuning of neuronal responses to 11 types of visual stimulus parameters has been consistently reported. Overall, the optimal spatial frequency (SF) of constituent neurons decreases with cortical hierarchy. Moreover, V5 neurons are distinct from neurons in other areas for their higher direction selectivity, higher contrast sensitivity, higher temporal frequency tuning, and wider SF bandwidth. We also discuss currently unavailable data that could be useful for biologically accurate models.

An efficient analytical reduction of detailed nonlinear neuron models

Detailed conductance-based nonlinear neuron models consisting of thousands of synapses are key for understanding of the computational properties of single neurons and large neuronal networks, and for interpreting experimental results. Simulations of these models are computationally expensive, considerably curtailing their utility. Neuron_Reduce is a new analytical approach to reduce the morphological complexity and computational time of nonlinear neuron models. Synapses and active membrane channels are mapped to the reduced model preserving their transfer impedance to the soma; synapses with identical transfer impedance are merged into one NEURON process still retaining their individual activation times. Neuron_Reduce accelerates the simulations by 40-250 folds for a variety of cell types and realistic number (10,000-100,000) of synapses while closely replicating voltage dynamics and specific dendritic computations. The reduced neuron-models will enable realistic simulations of neural networks at unprecedented scale, including networks emerging from micro-connectomics efforts and biologically-inspired "deep networks". Neuron_Reduce is publicly available and is straightforward to implement.

Cortical Output Is Gated by Horizontally Projecting Neurons in the Deep Layers

Pyramidal tract neurons (PTs) represent the major output cell type of the mammalian neocortex. Here, we report the origins of the PTs’ ability to respond to a broad range of stimuli with onset latencies that rival or even precede those of their intracortical input neurons. We find that neurons with extensive horizontally projecting axons cluster around the deep-layer terminal fields of primary thalamocortical axons. The strategic location of these corticocortical neurons results in high convergence of thalamocortical inputs, which drive reliable sensory-evoked responses that precede those in other excitatory cell types. The resultant fast and horizontal stream of excitation provides PTs throughout the cortical area with input that acts to amplify additional inputs from thalamocortical and other intracortical populations. The fast onsets and broadly tuned characteristics of PT responses hence reflect a gating mechanism in the deep layers, which assures that sensory-evoked input can be reliably transformed into cortical output.

•

Simulations predict in vivo responses for major output cell type of the neocortex

•

Simulations reveal strategy how to test the origins of cortical output empirically

•

Manipulations confirm that deep-layer corticocortical neurons gate cortical output

•

Gating of cortical output originates from deep-layer thalamocortical input stratum

Firing rate distributions in spiking networks with heterogeneous connectivity

Mean-field theory for networks of spiking neurons based on the so-called diffusion approximation has been used to calculate certain measures of neuronal activity which can be compared with experimental data. This includes the distribution of firing rates across the network. However, the theory in its current form applies only to networks in which there is relatively little heterogeneity in the number of incoming and outgoing connections per neuron. Here we extend this theory to include networks with arbitrary degree distributions. Furthermore, the theory takes into account correlations in the in-degree and out-degree of neurons, which would arise, e.g., in the case of networks with hublike neurons. Finally, we show that networks with broad and positively correlated degrees can generate a large-amplitude sustained response to transient stimuli which does not occur in more homogeneous networks.

Human Cortical Pyramidal Neurons: From Spines to Spikes via Models

We present detailed models of pyramidal cells from human neocortex, including models on their excitatory synapses, dendritic spines, dendritic NMDA- and somatic/axonal Na⁺ spikes that provided new insights into signal processing and computational capabilities of these principal cells. Six human layer 2 and layer 3 pyramidal cells (HL2/L3 PCs) were modeled, integrating detailed anatomical and physiological data from both fresh and postmortem tissues from human temporal cortex. The models predicted particularly large AMPA- and NMDA-conductances per synaptic contact (0.88 and 1.31 nS, respectively) and a steep dependence of the NMDA-conductance on voltage. These estimates were based on intracellular recordings from synaptically-connected HL2/L3 pairs, combined with extra-cellular current injections and use of synaptic blockers, and the assumption of five contacts per synaptic connection. A large dataset of high-resolution reconstructed HL2/L3 dendritic spines provided estimates for the EPSPs at the spine head (12.7 ± 4.6 mV), spine base (9.7 ± 5.0 mV), and soma (0.3 ± 0.1 mV), and for the spine neck resistance (50–80 MΩ). Matching the shape and firing pattern of experimental somatic Na⁺-spikes provided estimates for the density of the somatic/axonal excitable membrane ion channels, predicting that 134 ± 28 simultaneously activated HL2/L3-HL2/L3 synapses are required for generating (with 50% probability) a somatic Na⁺ spike. Dendritic NMDA spikes were triggered in the model when 20 ± 10 excitatory spinous synapses were simultaneously activated on individual dendritic branches. The particularly large number of basal dendrites in HL2/L3 PCs and the distinctive cable elongation of their terminals imply that ~25 NMDA-spikes could be generated independently and simultaneously in these cells, as compared to ~14 in L2/3 PCs from the rat somatosensory cortex. These multi-sites non-linear signals, together with the large (~30,000) excitatory synapses/cell, equip human L2/L3 PCs with enhanced computational capabilities. Our study provides the most comprehensive model of any human neuron to-date demonstrating the biophysical and computational distinctiveness of human cortical neurons.

MultiMap: A Tool to Automatically Extract and Analyse Spatial Microscopic Data From Large Stacks of Confocal Microscopy Images

The development of 3D visualization and reconstruction methods to analyse microscopic structures at different levels of resolutions is of great importance to define brain microorganization and connectivity. MultiMap is a new tool that allows the visualization, 3D segmentation and quantification of fluorescent structures selectively in the neuropil from large stacks of confocal microscopy images. The major contribution of this tool is the posibility to easily navigate and create regions of interest of any shape and size within a large brain area that will be automatically 3D segmented and quantified to determine the density of puncta in the neuropil. As a proof of concept, we focused on the analysis of glutamatergic and GABAergic presynaptic axon terminals in the mouse hippocampal region to demonstrate its use as a tool to provide putative excitatory and inhibitory synaptic maps. The segmentation and quantification method has been validated over expert labeled images of the mouse hippocampus and over two benchmark datasets, obtaining comparable results to the expert detections.

Cell Type-Specific Structural Organization of the Six Layers in Rat Barrel Cortex

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.

3D reconstruction and standardization of the rat facial nucleus for precise mapping of vibrissal motor networks

The rodent facial nucleus (FN) comprises motoneurons (MNs) that control the facial musculature. In the lateral part of the FN, populations of vibrissal motoneurons (vMNs) innervate two groups of muscles that generate movements of the whiskers. Vibrissal MNs thus represent the terminal point of the neuronal networks that generate rhythmic whisking during exploratory behaviors and that modify whisker movements based on sensory-motor feedback during tactile-based perception. Here, we combined retrograde tracer injections into whisker-specific muscles, with large-scale immunohistochemistry and digital reconstructions to generate an average model of the rat FN. The model incorporates measurements of the FN geometry, its cellular organization and a whisker row-specific map formed by vMNs. Furthermore, the model provides a digital 3D reference frame that allows registering structural data - obtained across scales and animals - into a common coordinate system with a precision of ∼60 µm. We illustrate the registration method by injecting replication competent rabies virus into the muscle of a single whisker. Retrograde transport of the virus to vMNs enabled reconstruction of their dendrites. Subsequent trans-synaptic transport enabled mapping the presynaptic neurons of the reconstructed vMNs. Registration of these data to the FN reference frame provides a first account of the morphological and synaptic input variability within a population of vMNs that innervate the same muscle.

Neuroanatomy and Global Neuroscience

Our brains are like a dense forest-a complex, seemingly impenetrable terrain of interacting cells mediating cognition and behavior. However, we should view the challenge of understanding the brain with optimism, provided that we choose appropriate strategies for the development of global neuroscience.

A Manual Segmentation Tool for Three-Dimensional Neuron Datasets

To date, automated or semi-automated software and algorithms for segmentation of neurons from three-dimensional imaging datasets have had limited success. The gold standard for neural segmentation is considered to be the manual isolation performed by an expert. To facilitate the manual isolation of complex objects from image stacks, such as neurons in their native arrangement within the brain, a new Manual Segmentation Tool (ManSegTool) has been developed. ManSegTool allows user to load an image stack, scroll down the images and to manually draw the structures of interest stack-by-stack. Users can eliminate unwanted regions or split structures (i.e., branches from different neurons that are too close each other, but, to the experienced eye, clearly belong to a unique cell), to view the object in 3D and save the results obtained. The tool can be used for testing the performance of a single-neuron segmentation algorithm or to extract complex objects, where the available automated methods still fail. Here we describe the software's main features and then show an example of how ManSegTool can be used to segment neuron images acquired using a confocal microscope. In particular, expert neuroscientists were asked to segment different neurons from which morphometric variables were subsequently extracted as a benchmark for precision. In addition, a literature-defined index for evaluating the goodness of segmentation was used as a benchmark for accuracy. Neocortical layer axons from a DIADEM challenge dataset were also segmented with ManSegTool and compared with the manual “gold-standard” generated for the competition.

From single cells and single columns to cortical networks: dendritic excitability, coincidence detection and synaptic transmission in brain slices and brains

Although patch pipettes were initially designed to record extracellularly the elementary current events from muscle and neuron membranes, the whole-cell and loose cell-attached recording configurations proved to be useful tools for examination of signalling within and between nerve cells. In this Paton Prize Lecture, I will initially summarize work on electrical signalling within single neurons, describing communication between the dendritic compartments, soma and nerve terminals via forward- and backward-propagating action potentials. The newly discovered dendritic excitability endows neurons with the capacity for coincidence detection of spatially separated subthreshold inputs. When these are occurring during a time window of tens of milliseconds, this information is broadcast to other cells by the initiation of bursts of action potentials (AP bursts). The occurrence of AP bursts critically impacts signalling between neurons that are controlled by target-cell-specific transmitter release mechanisms at downstream synapses even in different terminals of the same neuron. This can, in turn, induce mechanisms that underly synaptic plasticity when AP bursts occur within a short time window, both presynaptically in terminals and postsynaptically in dendrites. A fundamental question that arises from these findings is: 'what are the possible functions of active dendritic excitability with respect to network dynamics in the intact cortex of behaving animals?' To answer this question, I highlight in this review the functional and anatomical architectures of an average cortical column in the vibrissal (whisker) field of the somatosensory cortex (vS1), with an emphasis on the functions of layer 5 thick-tufted cells (L5tt) embedded in this structure. Sensory-evoked synaptic and action potential responses of these major cortical output neurons are compared with responses in the afferent pathway, viz. the neurons in primary somatosensory thalamus and in one of their efferent targets, the secondary somatosensory thalamus. Coincidence-detection mechanisms appear to be implemented in vivo as judged from the occurrence of AP bursts. Three-dimensional reconstructions of anatomical projections suggest that inputs of several combinations of thalamocortical projections and intra- and transcolumnar connections, specifically those from infragranular layers, could trigger active dendritic mechanisms that generate AP bursts. Finally, recordings from target cells of a column reveal the importance of AP bursts for signal transfer to these cells. The observations lead to the hypothesis that in vS1 cortex, the sensory afferent sensory code is transformed, at least in part, from a rate to an interval (burst) code that broadcasts the occurrence of whisker touch to different targets of L5tt cells. In addition, the occurrence of pre- and postsynaptic AP bursts may, in the long run, alter touch representation in cortex.

EM connectomics reveals axonal target variation in a sequence-generating network

The sequential activation of neurons has been observed in various areas of the brain, but in no case is the underlying network structure well understood. Here we examined the circuit anatomy of zebra finch HVC, a cortical region that generates sequences underlying the temporal progression of the song. We combined serial block-face electron microscopy with light microscopy to determine the cell types targeted by HVC(RA) neurons, which control song timing. Close to their soma, axons almost exclusively targeted inhibitory interneurons, consistent with what had been found with electrical recordings from pairs of cells. Conversely, far from the soma the targets were mostly other excitatory neurons, about half of these being other HVC(RA) cells. Both observations are consistent with the notion that the neural sequences that pace the song are generated by global synaptic chains in HVC embedded within local inhibitory networks.

Comparative ultrastructural features of excitatory synapses in the visual and frontal cortices of the adult mouse and monkey

The excitatory glutamatergic synapse is the principal site of communication between cortical pyramidal neurons and their targets, a key locus of action of many drugs, and highly vulnerable to dysfunction and loss in neurodegenerative disease. A detailed knowledge of the structure of these synapses in distinct cortical areas and across species is a prerequisite for understanding the anatomical underpinnings of cortical specialization and, potentially, selective vulnerability in neurological disorders. We used serial electron microscopy to assess the ultrastructural features of excitatory (asymmetric) synapses in the layers 2-3 (L2-3) neuropil of visual (V1) and frontal (FC) cortices of the adult mouse and compared findings to those in the rhesus monkey (V1 and lateral prefrontal cortex [LPFC]). Analyses of multiple ultrastructural variables revealed four organizational features. First, the density of asymmetric synapses does not differ between frontal and visual cortices in either species, but is significantly higher in mouse than in monkey. Second, the structural properties of asymmetric synapses in mouse V1 and FC are nearly identical, by stark contrast to the significant differences seen between monkey V1 and LPFC. Third, while the structural features of postsynaptic entities in mouse and monkey V1 do not differ, the size of presynaptic boutons are significantly larger in monkey V1. Fourth, both presynaptic and postsynaptic entities are significantly smaller in the mouse FC than in the monkey LPFC. The diversity of synaptic ultrastructural features demonstrated here have broad implications for the nature and efficacy of glutamatergic signaling in distinct cortical areas within and across species.

Pyramidal Neurons Are Not Generalizable Building Blocks of Cortical Networks

A key challenge in cortical neuroscience is to gain a comprehensive understanding of how pyramidal neuron heterogeneity across different areas and species underlies the functional specialization of individual neurons, networks, and areas. Comparative studies have been important in this endeavor, providing data relevant to the question of which of the many inherent properties of individual pyramidal neurons are necessary and sufficient for species-specific network and areal function. In this mini review, the importance of pyramidal neuron structural properties for signaling are outlined, followed by a summary of our recent work comparing the structural features of mouse (C57/BL6 strain) and rhesus monkey layer 3 (L3) pyramidal neurons in primary visual and frontal association cortices and their implications for neuronal and areal function. Based on these and other published data, L3 pyramidal neurons plausibly might be considered broadly “generalizable” from one area to another in the mouse neocortex due to their many similarities, but major differences in the properties of these neurons in diverse areas in the rhesus monkey neocortex rules this out in the primate. Further, fundamental differences in the dendritic topology of mouse and rhesus monkey pyramidal neurons highlight the implausibility of straightforward scaling and/or extrapolation from mouse to primate neurons and cortical networks.

The Impact of Structural Heterogeneity on Excitation-Inhibition Balance in Cortical Networks

Models of cortical dynamics often assume a homogeneous connectivity structure. However, we show that heterogeneous input connectivity can prevent the dynamic balance between excitation and inhibition, a hallmark of cortical dynamics, and yield unrealistically sparse and temporally regular firing. Anatomically based estimates of the connectivity of layer 4 (L4) rat barrel cortex and numerical simulations of this circuit indicate that the local network possesses substantial heterogeneity in input connectivity, sufficient to disrupt excitation-inhibition balance. We show that homeostatic plasticity in inhibitory synapses can align the functional connectivity to compensate for structural heterogeneity. Alternatively, spike-frequency adaptation can give rise to a novel state in which local firing rates adjust dynamically so that adaptation currents and synaptic inputs are balanced. This theory is supported by simulations of L4 barrel cortex during spontaneous and stimulus-evoked conditions. Our study shows how synaptic and cellular mechanisms yield fluctuation-driven dynamics despite structural heterogeneity in cortical circuits.

•

Structural heterogeneity threatens the dynamic balance of excitation and inhibition

•

Reconstruction of cortical networks reveals significant structural heterogeneity

•

Spike-frequency adaptation can act locally to facilitate global balance

•

Inhibitory homeostatic plasticity can compensate for structural imbalance

Morphological Characteristics of Electrophysiologically Characterized Layer Vb Pyramidal Cells in Rat Barrel Cortex

Layer Vb pyramidal cells are the major output neurons of the neocortex and transmit the outcome of cortical columnar signal processing to distant target areas. At the same time they contribute to local tactile information processing by emitting recurrent axonal collaterals into the columnar microcircuitry. It is, however, not known how exactly the two types of pyramidal cells, called slender-tufted and thick-tufted, contribute to the local circuitry. Here, we investigated in the rat barrel cortex the detailed quantitative morphology of biocytin-filled layer Vb pyramidal cells in vitro, which were characterized for their intrinsic electrophysiology with special emphasis on their action potential firing pattern. Since we stained the same slices for cytochrome oxidase, we could also perform layer- and column-related analyses. Our results suggest that in layer Vb the unambiguous action potential firing patterns "regular spiking (RS)" and "repetitive burst spiking (RB)" (previously called intrinsically burst spiking) correlate well with a distinct morphology. RS pyramidal cells are somatodendritically of the slender-tufted type and possess numerous local intralaminar and intracolumnar axonal collaterals, mostly reaching layer I. By contrast, their transcolumnar projections are less well developed. The RB pyramidal cells are somatodendritically of the thick-tufted type and show only relatively sparse local axonal collaterals, which are preferentially emitted as long horizontal or oblique infragranular collaterals. However, contrary to many previous slice studies, a substantial number of these neurons also showed axonal collaterals reaching layer I. Thus, electrophysiologically defined pyramidal cells of layer Vb show an input and output pattern which suggests RS cells to be more "locally segregating" signal processors whereas RB cells seem to act more on a "global integrative" scale.

From Neuron Biophysics to Orientation Selectivity in Electrically Coupled Networks of Neocortical L2/3 Large Basket Cells

In the neocortex, inhibitory interneurons of the same subtype are electrically coupled with each other via dendritic gap junctions (GJs). The impact of multiple GJs on the biophysical properties of interneurons and thus on their input processing is unclear. The present experimentally based theoretical study examined GJs in L2/3 large basket cells (L2/3 LBCs) with 3 goals in mind: (1) To evaluate the errors due to GJs in estimating the cable properties of individual L2/3 LBCs and suggest ways to correct these errors when modeling these cells and the networks they form; (2) to bracket the GJ conductance value (0.05-0.25 nS) and membrane resistivity (10 000-40 000 Ω cm(2)) of L2/3 LBCs; these estimates are tightly constrained by in vitro input resistance (131 ± 18.5 MΩ) and the coupling coefficient (1-3.5%) of these cells; and (3) to explore the functional implications of GJs, and show that GJs: (i) dynamically modulate the effective time window for synaptic integration; (ii) improve the axon's capability to encode rapid changes in synaptic inputs; and (iii) reduce the orientation selectivity, linearity index, and phase difference of L2/3 LBCs. Our study provides new insights into the role of GJs and calls for caution when using in vitro measurements for modeling electrically coupled neuronal networks.

High-Degree Neurons Feed Cortical Computations

Recent work has shown that functional connectivity among cortical neurons is highly varied, with a small percentage of neurons having many more connections than others. Also, recent theoretical developments now make it possible to quantify how neurons modify information from the connections they receive. Therefore, it is now possible to investigate how information modification, or computation, depends on the number of connections a neuron receives (in-degree) or sends out (out-degree). To do this, we recorded the simultaneous spiking activity of hundreds of neurons in cortico-hippocampal slice cultures using a high-density 512-electrode array. This preparation and recording method combination produced large numbers of neurons recorded at temporal and spatial resolutions that are not currently available in any in vivo recording system. We utilized transfer entropy (a well-established method for detecting linear and nonlinear interactions in time series) and the partial information decomposition (a powerful, recently developed tool for dissecting multivariate information processing into distinct parts) to quantify computation between neurons where information flows converged. We found that computations did not occur equally in all neurons throughout the networks. Surprisingly, neurons that computed large amounts of information tended to receive connections from high out-degree neurons. However, the in-degree of a neuron was not related to the amount of information it computed. To gain insight into these findings, we developed a simple feedforward network model. We found that a degree-modified Hebbian wiring rule best reproduced the pattern of computation and degree correlation results seen in the real data. Interestingly, this rule also maximized signal propagation in the presence of network-wide correlations, suggesting a mechanism by which cortex could deal with common random background input. These are the first results to show that the extent to which a neuron modifies incoming information streams depends on its topological location in the surrounding functional network.

We recorded the electrical activity of hundreds of neurons simultaneously in brain tissue from mice and we analyzed these signals using state-of-the-art tools from information theory. These tools allowed us to ascertain which neurons were transmitting information to other neurons and to characterize the computations performed by neurons using the inputs they received from two or more other neurons. We found that computations did not occur equally in all neurons throughout the networks. Surprisingly, neurons that computed large amounts of information tended to be recipients of information from neurons with a large number of outgoing connections. Interestingly, the number of incoming connections to a neuron was not related to the amount of information that neuron computed. To better understand these results, we built a network model to match the data. Unexpectedly, the model also maximized information transfer in the presence of network-wide correlations. This suggested a way that networks of cortical neurons could deal with common random background input. These results are the first to show that the amount of information computed by a neuron depends on where it is located in the surrounding network.

Robustness of sensory-evoked excitation is increased by inhibitory inputs to distal apical tuft dendrites

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.

An algorithm to predict the connectome of neural microcircuits

Experimentally mapping synaptic connections, in terms of the numbers and locations of their synapses and estimating connection probabilities, is still not a tractable task, even for small volumes of tissue. In fact, the six layers of the neocortex contain thousands of unique types of synaptic connections between the many different types of neurons, of which only a handful have been characterized experimentally. Here we present a theoretical framework and a data-driven algorithmic strategy to digitally reconstruct the complete synaptic connectivity between the different types of neurons in a small well-defined volume of tissue—the micro-scale connectome of a neural microcircuit. By enforcing a set of established principles of synaptic connectivity, and leveraging interdependencies between fundamental properties of neural microcircuits to constrain the reconstructed connectivity, the algorithm yields three parameters per connection type that predict the anatomy of all types of biologically viable synaptic connections. The predictions reproduce a spectrum of experimental data on synaptic connectivity not used by the algorithm. We conclude that an algorithmic approach to the connectome can serve as a tool to accelerate experimental mapping, indicating the minimal dataset required to make useful predictions, identifying the datasets required to improve their accuracy, testing the feasibility of experimental measurements, and making it possible to test hypotheses of synaptic connectivity.

The anatomical problem posed by brain complexity and size: a potential solution

Over the years the field of neuroanatomy has evolved considerably but unraveling the extraordinary structural and functional complexity of the brain seems to be an unattainable goal, partly due to the fact that it is only possible to obtain an imprecise connection matrix of the brain. The reasons why reaching such a goal appears almost impossible to date is discussed here, together with suggestions of how we could overcome this anatomical problem by establishing new methodologies to study the brain and by promoting interdisciplinary collaboration. Generating a realistic computational model seems to be the solution rather than attempting to fully reconstruct the whole brain or a particular brain region.

Beyond Columnar Organization: Cell Type- and Target Layer-Specific Principles of Horizontal Axon Projection Patterns in Rat Vibrissal Cortex

Vertical thalamocortical afferents give rise to the elementary functional units of sensory cortex, cortical columns. Principles that underlie communication between columns remain however unknown. Here we unravel these by reconstructing in vivo-labeled neurons from all excitatory cell types in the vibrissal part of rat primary somatosensory cortex (vS1). Integrating the morphologies into an exact 3D model of vS1 revealed that the majority of intracortical (IC) axons project far beyond the borders of the principal column. We defined the corresponding innervation volume as the IC-unit. Deconstructing this structural cortical unit into its cell type-specific components, we found asymmetric projections that innervate columns of either the same whisker row or arc, and which subdivide vS1 into 2 orthogonal [supra-]granular and infragranular strata. We show that such organization could be most effective for encoding multi whisker inputs. Communication between columns is thus organized by multiple highly specific horizontal projection patterns, rendering IC-units as the primary structural entities for processing complex sensory stimuli.