Marcel Oberlaender

Cell Type–Specific Three-Dimensional Structure of Thalamocortical Circuits in a Column of Rat Vibrissal Cortex

Soma location, dendrite morphology, and synaptic innervation may represent key determinants of functional responses of individual neurons, such as sensory-evoked spiking. Here, we reconstruct the 3D circuits formed by thalamocortical afferents from the lemniscal pathway and excitatory neurons of an anatomically defined cortical column in rat vibrissal cortex. We objectively classify 9 cortical cell types and estimate the number and distribution of their somata, dendrites, and thalamocortical synapses. Somata and dendrites of most cell types intermingle, while thalamocortical connectivity depends strongly upon the cell type and the 3D soma location of the postsynaptic neuron. Correlating dendrite morphology and thalamocortical connectivity to functional responses revealed that the lemniscal afferents can account for some of the cell type- and location-specific subthreshold and spiking responses after passive whisker touch (e.g., in layer 4, but not for other cell types, e.g., in layer 5). Our data provides a quantitative 3D prediction of the cell type–specific lemniscal synaptic wiring diagram and elucidates structure–function relationships of this physiologically relevant pathway at single-cell resolution.

Generation of dense statistical connectomes from sparse morphological data.

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.

Number and Laminar Distribution of Neurons in a Thalamocortical Projection Column of Rat Vibrissal Cortex

This is the second article in a series of three studies that investigate the anatomical determinants of thalamocortical (TC) input to excitatory neurons in a cortical column of rat primary somatosensory cortex (S1). Here, we report the number and distribution of NeuN-positive neurons within the C2, D2, and D3 TC projection columns in P27 rat somatosensory barrel cortex based on an exhaustive identification of 89 834 somata in a 1.15 mm3 volume of cortex. A single column contained 19 109 ± 444 neurons (17 560 ± 399 when normalized to a standard-size projection column). Neuron density differences along the vertical column axis delineated “cytoarchitectonic” layers. The resulting neuron numbers per layer in the average column were 63 ± 10 (L1), 2039 ± 524 (L2), 3735 ± 905 (L3), 4447 ± 439 (L4), 1737 ± 251 (L5A), 2235 ± 99 (L5B), 3786 ± 168 (L6A), and 1066 ± 170 (L6B). These data were then used to derive the layer-specific action potential (AP) output of a projection column. The estimates confirmed previous reports suggesting that the ensembles of spiny L4 and thick-tufted pyramidal neurons emit the major fraction of APs of a column. The number of APs evoked in a column by a sensory stimulus (principal whisker deflection) was estimated as 4441 within 100 ms post-stimulus.

Large-Scale Automated Histology in the Pursuit of Connectomes

How does the brain compute? Answering this question necessitates neuronal connectomes, annotated graphs of all synaptic connections within defined brain areas. Further, understanding the energetics of the brains computations requires vascular graphs. The assembly of a connectome requires sensitive hardware tools to measure neuronal and neurovascular features in all three dimensions, as well 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.

Three-dimensional axon morphologies of individual layer 5 neurons indicate cell type-specific intracortical pathways for whisker motion and touch

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.

Sensory Experience Restructures Thalamocortical Axons during Adulthood

The brains capacity to rewire is thought to diminish with age. It is widely believed that development stabilizes the synapses from thalamus to cortex and that adult experience alters only synaptic connections between cortical neurons. Here we show that thalamocortical (TC) inputs themselves undergo massive plasticity in adults. We combined whole-cell recording from individual thalamocortical neurons in adult rats with a recently developed automatic tracing technique to reconstruct individual axonal trees. Whisker trimming substantially reduced thalamocortical axon length in barrel cortex but not the density of TC synapses along a fiber. Thus, sensory experience alters the total number of TC synapses. After trimming, sensory stimulation evoked more tightly time-locked responses among thalamorecipient layer 4 cortical neurons. These findings indicate that thalamocortical input itself remains plastic in adulthood, raising the possibility that the axons of other subcortical structures might also remain in flux throughout life.

3D Reconstruction and Standardization of the Rat Vibrissal Cortex for Precise Registration of Single Neuron Morphology

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.

Cellular organization of cortical barrel columns is whisker-specific

Significance Cortical columns are thought to be the elementary functional building blocks of sensory cortices. Here we show that the cellular architecture of cortical “barrel” columns in rodent somatosensory cortex is not stereotypic, but specific for each whisker on the animals’ snout. Our findings challenge the concepts underlying contemporary simulation efforts that build up large-scale network models of repeatedly occurring identical cortical circuits. The cellular organization of the cortex is of fundamental importance for elucidating the structural principles that underlie its functions. It has been suggested that reconstructing the structure and synaptic wiring of the elementary functional building block of mammalian cortices, the cortical column, might suffice to reverse engineer and simulate the functions of entire cortices. In the vibrissal area of rodent somatosensory cortex, whisker-related “barrel” columns have been referred to as potential cytoarchitectonic equivalents of functional cortical columns. Here, we investigated the structural stereotypy of cortical barrel columns by measuring the 3D neuronal composition of the entire vibrissal area in rat somatosensory cortex and thalamus. We found that the number of neurons per cortical barrel column and thalamic “barreloid” varied substantially within individual animals, increasing by ∼2.5-fold from dorsal to ventral whiskers. As a result, the ratio between whisker-specific thalamic and cortical neurons was remarkably constant. Thus, we hypothesize that the cellular architecture of sensory cortices reflects the degree of similarity in sensory input and not columnar and/or cortical uniformity principles.

Automated three-dimensional detection and counting of neuron somata

We present a novel approach for automated detection of neuron somata. A three-step processing pipeline is described on the example of confocal image stacks of NeuN-stained neurons from rat somato-sensory cortex. It results in a set of position landmarks, representing the midpoints of all neuron somata. In the first step, foreground and background pixels are identified, resulting in a binary image. It is based on local thresholding and compensates for imaging and staining artifacts. Once this pre-processing guarantees a standard image quality, clusters of touching neurons are separated in the second step, using a marker-based watershed approach. A model-based algorithm completes the pipeline. It assumes a dominant neuron population with Gaussian distributed volumes within one microscopic field of view. Remaining larger objects are hence split or treated as a second neuron type. A variation of the processing pipeline is presented, showing that our method can also be used for co-localization of neurons in multi-channel images. As an example, we process 2-channel stacks of NeuN-stained somata, labeling all neurons, counterstained with GAD67, labeling GABAergic interneurons, using an adapted pre-processing step for the second channel. The automatically generated landmark sets are compared to manually placed counterparts. A comparison yields that the deviation in landmark position is negligible and that the difference between the numbers of manually and automatically counted neurons is less than 4%. In consequence, this novel approach for neuron counting is a reliable and objective alternative to manual detection.

Transmitted light brightfield mosaic microscopy for three-dimensional tracing of single neuron morphology

A fundamental challenge in neuroscience is the determination of the three-dimensional (3D) morphology of neurons in the cortex. Here we describe a semiautomated method to trace single biocytin-filled neurons using a transmitted light brightfield microscope. The method includes 3D tracing of dendritic trees and axonal arbors from image stacks of serial 100-microm-thick tangential brain sections. Key functionalities include mosaic scanning and optical sectioning, high-resolution image restoration, and fast, parallel computing for neuron tracing. The mosaic technique compensates for the limited field of view at high magnification, allowing the acquisition of high-resolution image stacks on a scale of millimeters. The image restoration by deconvolution is based on experimentally verified assumptions about the optical system. Restoration yields a significant improvement of signal-to-noise ratio and resolution of neuronal structures in the image stack. Application of local threshold and thinning filters result in a 3D graph representation of dendrites and axons in a section. The reconstructed branches are then manually edited and aligned. Branches from adjacent sections are spliced, resulting in a complete 3D reconstruction of a neuron. A comparison with 3D reconstructions from manually traced neurons shows that the semiautomated system is a fast and reliable alternative to the manual tracing systems currently available.