High-throughput morphometric and transcriptomic profiling uncovers composition of naïve and sensory-deprived cortical cholinergic VIP/CHAT neurons.

Cortical neuronal networks control cognitive output, but their composition and modulation remain elusive. Here, we studied the morphological and transcriptional diversity of cortical cholinergic VIP/ChAT interneurons (VChIs), a sparse population with a largely unknown function. We focused on VChIs from the whole barrel cortex and developed a high-throughput automated reconstruction framework, termed PopRec, to characterize hundreds of VChIs from each mouse in an unbiased manner, while preserving 3D cortical coordinates in multiple cleared mouse brains, accumulating thousands of cells. We identified two fundamentally distinct morphological types of VChIs, bipolar and multipolar that differ in their cortical distribution and general morphological features. Following mild unilateral whisker deprivation on postnatal day seven, we found after three weeks both ipsi- and contralateral dendritic arborization differences and modified cortical depth and distribution patterns in the barrel fields alone. To seek the transcriptomic drivers, we developed NuNeX, a method for isolating nuclei from fixed tissues, to explore sorted VChIs. This highlighted differentially expressed neuronal structural transcripts, altered exitatory innervation pathways and established Elmo1 as a key regulator of morphology following deprivation.

Strong and reliable synaptic communication between pyramidal neurons in adult human cerebral cortex.

Synaptic transmission constitutes the primary mode of communication between neurons. It is extensively studied in rodent but not human neocortex. We characterized synaptic transmission between pyramidal neurons in layers 2 and 3 using neurosurgically resected human middle temporal gyrus (MTG, Brodmann area 21), which is part of the distributed language circuitry. We find that local connectivity is comparable with mouse layer 2/3 connections in the anatomical homologue (temporal association area), but synaptic connections in human are 3-fold stronger and more reliable (0% vs 25% failure rates, respectively). We developed a theoretical approach to quantify properties of spinous synapses showing that synaptic conductance and voltage change in human dendritic spines are 3-4-folds larger compared with mouse, leading to significant NMDA receptor activation in human unitary connections. This model prediction was validated experimentally by showing that NMDA receptor activation increases the amplitude and prolongs decay of unitary excitatory postsynaptic potentials in human but not in mouse connections. Since NMDA-dependent recurrent excitation facilitates persistent activity (supporting working memory), our data uncovers cortical microcircuit properties in human that may contribute to language processing in MTG.

Synaptic Input and ACh Modulation Regulate Dendritic Ca2+ Spike Duration in Pyramidal Neurons, Directly Affecting Their Somatic Output.

Nonlinear synaptic integration in dendrites is a fundamental aspect of neural computation. One such key mechanism is the Ca2+ spike at the apical tuft of pyramidal neurons. Characterized by a plateau potential sustained for tens of milliseconds, the Ca2+ spike amplifies excitatory input, facilitates somatic action potentials (APs), and promotes synaptic plasticity. Despite its essential role, the mechanisms regulating it are largely unknown. Using a compartmental model of a layer 5 pyramidal cell (L5PC), we explored the plateau and termination phases of the Ca2+ spike under input current perturbations, long-step current-injections, and variations in the dendritic high-voltage-activated Ca2+ conductance (that occur during cholinergic modulation). We found that, surprisingly, timed excitatory input can shorten the Ca2+ spike duration while inhibitory input can either elongate or terminate it. A significant elongation also occurs when the high-voltage-activated Ca2+ channels (CaHVA) conductance is increased. To mechanistically understand these phenomena, we analyzed the currents involved in the spike. The plateau and termination phases are almost exclusively controlled by the CaHVA inward current and the Im outward K+ current. We reduced the full model to a single-compartment model that faithfully preserved the responses of the Ca2+ spike to interventions and consisted of two dynamic variables: the membrane potential and the K+-channel activation level. A phase-plane analysis of the reduced model provides testable predictions for modulating the Ca2+ spike and reveals various dynamical regimes that explain the robust nature of the spike. Regulating the duration of the Ca2+ spike significantly impacts the cell synaptic-plasticity window and, as we show, its input-output relationship.SIGNIFICANCE STATEMENT Pyramidal neurons are the cortex's principal projection neurons. In their apical tuft, dendritic Ca2+ spikes significantly impact information processing, synaptic plasticity, and the cell's input-output relationship. Therefore, it is essential to understand the mechanisms regulating them. Using a compartmental model of a layer 5 pyramidal cell (L5PC), we explored the Ca2+ spike responses to synaptic perturbations and cholinergic modulation. We showed a counterintuitive phenomenon: early excitatory input shortens the spike, whereas weak inhibition elongates it. Also, we demonstrated that acetylcholine (ACh) extends the spike. Through a reduced model containing only the membrane potential and the K+-channel activation level, we explained these phenomena using a phase-plane analysis. Our work provides new information about the robustness of the Ca2+ spike and its controlling mechanisms.

Burst control: Synaptic conditions for burst generation in cortical layer 5 pyramidal neurons.

The output of neocortical layer 5 pyramidal cells (L5PCs) is expressed by a train of single spikes with intermittent bursts of multiple spikes at high frequencies. The bursts are the result of nonlinear dendritic properties, including Na+, Ca2+, and NMDA spikes, that interact with the ~10,000 synapses impinging on the neuron's dendrites. Output spike bursts are thought to implement key dendritic computations, such as coincidence detection of bottom-up inputs (arriving mostly at the basal tree) and top-down inputs (arriving mostly at the apical tree). In this study we used a detailed nonlinear model of L5PC receiving excitatory and inhibitory synaptic inputs to explore the conditions for generating bursts and for modulating their properties. We established the excitatory input conditions on the basal versus the apical tree that favor burst and show that there are two distinct types of bursts. Bursts consisting of 3 or more spikes firing at < 200 Hz, which are generated by stronger excitatory input to the basal versus the apical tree, and bursts of ~2-spikes at ~250 Hz, generated by prominent apical tuft excitation. Localized and well-timed dendritic inhibition on the apical tree differentially modulates Na+, Ca2+, and NMDA spikes and, consequently, finely controls the burst output. Finally, we explored the implications of different burst classes and respective dendritic inhibition for regulating synaptic plasticity.

The gradient clusteron: A model neuron that learns to solve classification tasks via dendritic nonlinearities, structural plasticity, and gradient descent.

Synaptic clustering on neuronal dendrites has been hypothesized to play an important role in implementing pattern recognition. Neighboring synapses on a dendritic branch can interact in a synergistic, cooperative manner via nonlinear voltage-dependent mechanisms, such as NMDA receptors. Inspired by the NMDA receptor, the single-branch clusteron learning algorithm takes advantage of location-dependent multiplicative nonlinearities to solve classification tasks by randomly shuffling the locations of "under-performing" synapses on a model dendrite during learning ("structural plasticity"), eventually resulting in synapses with correlated activity being placed next to each other on the dendrite. We propose an alternative model, the gradient clusteron, or G-clusteron, which uses an analytically-derived gradient descent rule where synapses are "attracted to" or "repelled from" each other in an input- and location-dependent manner. We demonstrate the classification ability of this algorithm by testing it on the MNIST handwritten digit dataset and show that, when using a softmax activation function, the accuracy of the G-clusteron on the all-versus-all MNIST task (~85%) approaches that of logistic regression (~93%). In addition to the location update rule, we also derive a learning rule for the synaptic weights of the G-clusteron ("functional plasticity") and show that a G-clusteron that utilizes the weight update rule can achieve ~89% accuracy on the MNIST task. We also show that a G-clusteron with both the weight and location update rules can learn to solve the XOR problem from arbitrary initial conditions.

Maximally efficient prediction in the early fly visual system may support evasive flight maneuvers.

The visual system must make predictions to compensate for inherent delays in its processing. Yet little is known, mechanistically, about how prediction aids natural behaviors. Here, we show that despite a 20-30ms intrinsic processing delay, the vertical motion sensitive (VS) network of the blowfly achieves maximally efficient prediction. This prediction enables the fly to fine-tune its complex, yet brief, evasive flight maneuvers according to its initial ego-rotation at the time of detection of the visual threat. Combining a rich database of behavioral recordings with detailed compartmental modeling of the VS network, we further show that the VS network has axonal gap junctions that are critical for optimal prediction. During evasive maneuvers, a VS subpopulation that directly innervates the neck motor center can convey predictive information about the fly's future ego-rotation, potentially crucial for ongoing flight control. These results suggest a novel sensory-motor pathway that links sensory prediction to behavior.

Realistic retinal modeling unravels the differential role of excitation and inhibition to starburst amacrine cells in direction selectivity.

Retinal direction-selectivity originates in starburst amacrine cells (SACs), which display a centrifugal preference, responding with greater depolarization to a stimulus expanding from soma to dendrites than to a collapsing stimulus. Various mechanisms were hypothesized to underlie SAC centrifugal preference, but dissociating them is experimentally challenging and the mechanisms remain debatable. To address this issue, we developed the Retinal Stimulation Modeling Environment (RSME), a multifaceted data-driven retinal model that encompasses detailed neuronal morphology and biophysical properties, retina-tailored connectivity scheme and visual input. Using a genetic algorithm, we demonstrated that spatiotemporally diverse excitatory inputs-sustained in the proximal and transient in the distal processes-are sufficient to generate experimentally validated centrifugal preference in a single SAC. Reversing these input kinetics did not produce any centrifugal-preferring SAC. We then explored the contribution of SAC-SAC inhibitory connections in establishing the centrifugal preference. SAC inhibitory network enhanced the centrifugal preference, but failed to generate it in its absence. Embedding a direction selective ganglion cell (DSGC) in a SAC network showed that the known SAC-DSGC asymmetric connectivity by itself produces direction selectivity. Still, this selectivity is sharpened in a centrifugal-preferring SAC network. Finally, we use RSME to demonstrate the contribution of SAC-SAC inhibitory connections in mediating direction selectivity and recapitulate recent experimental findings. Thus, using RSME, we obtained a mechanistic understanding of SACs' centrifugal preference and its contribution to direction selectivity.

Differential Structure of Hippocampal CA1 Pyramidal Neurons in the Human and Mouse.

Pyramidal neurons are the most common cell type and are considered the main output neuron in most mammalian forebrain structures. In terms of function, differences in the structure of the dendrites of these neurons appear to be crucial in determining how neurons integrate information. To further shed light on the structure of the human pyramidal neurons we investigated the geometry of pyramidal cells in the human and mouse CA1 region-one of the most evolutionary conserved archicortical regions, which is critically involved in the formation, consolidation, and retrieval of memory. We aimed to assess to what extent neurons corresponding to a homologous region in different species have parallel morphologies. Over 100 intracellularly injected and 3D-reconstructed cells across both species revealed that dendritic and axonal morphologies of human cells are not only larger but also have structural differences, when compared to mouse. The results show that human CA1 pyramidal cells are not a stretched version of mouse CA1 cells. These results indicate that there are some morphological parameters of the pyramidal cells that are conserved, whereas others are species-specific.

Efficient encoding of motion is mediated by gap junctions in the fly visual system.

Understanding the computational implications of specific synaptic connectivity patterns is a fundamental goal in neuroscience. In particular, the computational role of ubiquitous electrical synapses operating via gap junctions remains elusive. In the fly visual system, the cells in the vertical-system network, which play a key role in visual processing, primarily connect to each other via axonal gap junctions. This network therefore provides a unique opportunity to explore the functional role of gap junctions in sensory information processing. Our information theoretical analysis of a realistic VS network model shows that within 10 ms following the onset of the visual input, the presence of axonal gap junctions enables the VS system to efficiently encode the axis of rotation, θ, of the fly's ego motion. This encoding efficiency, measured in bits, is near-optimal with respect to the physical limits of performance determined by the statistical structure of the visual input itself. The VS network is known to be connected to downstream pathways via a subset of triplets of the vertical system cells; we found that because of the axonal gap junctions, the efficiency of this subpopulation in encoding θ is superior to that of the whole vertical system network and is robust to a wide range of signal to noise ratios. We further demonstrate that this efficient encoding of motion by this subpopulation is necessary for the fly's visually guided behavior, such as banked turns in evasive maneuvers. Because gap junctions are formed among the axons of the vertical system cells, they only impact the system's readout, while maintaining the dendritic input intact, suggesting that the computational principles implemented by neural circuitries may be much richer than previously appreciated based on point neuron models. Our study provides new insights as to how specific network connectivity leads to efficient encoding of sensory stimuli.

Comprehensive Morpho-Electrotonic Analysis Shows 2 Distinct Classes of L2 and L3 Pyramidal Neurons in Human Temporal Cortex.

There have been few quantitative characterizations of the morphological, biophysical, and cable properties of neurons in the human neocortex. We employed feature-based statistical methods on a rare data set of 60 3D reconstructed pyramidal neurons from L2 and L3 in the human temporal cortex (HL2/L3 PCs) removed after brain surgery. Of these cells, 25 neurons were also characterized physiologically. Thirty-two morphological features were analyzed (e.g., dendritic surface area, 36 333 ± 18 157 µm2; number of basal trees, 5.55 ± 1.47; dendritic diameter, 0.76 ± 0.28 µm). Eighteen features showed a significant gradual increase with depth from the pia (e.g., dendritic length and soma radius). The other features showed weak or no correlation with depth (e.g., dendritic diameter). The basal dendritic terminals in HL2/L3 PCs are particularly elongated, enabling multiple nonlinear processing units in these dendrites. Unlike the morphological features, the active biophysical features (e.g., spike shapes and rates) and passive/cable features (e.g., somatic input resistance, 47.68 ± 15.26 MΩ, membrane time constant, 12.03 ± 1.79 ms, average dendritic cable length, 0.99 ± 0.24) were depth-independent. A novel descriptor for apical dendritic topology yielded 2 distinct classes, termed hereby as "slim-tufted" and "profuse-tufted" HL2/L3 PCs; the latter class tends to fire at higher rates. Thus, our morpho-electrotonic analysis shows 2 distinct classes of HL2/L3 PCs.

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.

Dendritic Excitability and Gain Control in Recurrent Cortical Microcircuits.

Layer 5 thick tufted pyramidal cells (TTCs) in the neocortex are particularly electrically complex, owing to their highly excitable dendrites. The interplay between dendritic nonlinearities and recurrent cortical microcircuit activity in shaping network response is largely unknown. We simulated detailed conductance-based models of TTCs forming recurrent microcircuits that were interconnected as found experimentally; the network was embedded in a realistic background synaptic activity. TTCs microcircuits significantly amplified brief thalamocortical inputs; this cortical gain was mediated by back-propagation activated N-methyl-D-aspartate depolarizations and dendritic back-propagation-activated Ca(2+) spike firing, ignited by the coincidence of thalamic-activated somatic spike and local dendritic synaptic inputs, originating from the cortical microcircuit. Surprisingly, dendritic nonlinearities in TTCs microcircuits linearly multiplied thalamic inputs--amplifying them while maintaining input selectivity. Our findings indicate that dendritic nonlinearities are pivotal in controlling the gain and the computational functions of TTCs microcircuits, which serve as a dominant output source for the neocortex.

Contribution of intracolumnar layer 2/3-to-layer 2/3 excitatory connections in shaping the response to whisker deflection in rat barrel cortex.

This computational study integrates anatomical and physiological data to assess the functional role of the lateral excitatory connections between layer 2/3 (L2/3) pyramidal cells (PCs) in shaping their response during early stages of intracortical processing of a whisker deflection (WD). Based on in vivo and in vitro recordings, and 3D reconstructions of connected pairs of L2/3 PCs, our model predicts that: 1) AMPAR and NMDAR conductances/synapse are 0.52 ± 0.24 and 0.40 ± 0.34 nS, respectively; 2) following WD, connection between L2/3 PCs induces a composite EPSPs of 7.6 ± 1.7 mV, well below the threshold for action potential (AP) initiation; 3) together with the excitatory feedforward L4-to-L2/3 connection, WD evoked a composite EPSP of 16.3 ± 3.5 mV and a probability of 0.01 to generate an AP. When considering the variability in L4 spiny neurons responsiveness, it increased to 17.8 ± 11.2 mV; this 3-fold increase in the SD yielded AP probability of 0.35; 4) the interaction between L4-to-L2/3 and L2/3-to-L2/3 inputs is highly nonlinear; 5) L2/3 dendritic morphology significantly affects L2/3 PCs responsiveness. We conclude that early stages of intracortical signaling of WD are dominated by a combination of feedforward L4-L2/3 and L2/3-L2/3 lateral connections.

Dendritic and Axonal Architecture of Individual Pyramidal Neurons across Layers of Adult Human Neocortex.

The size and shape of dendrites and axons are strong determinants of neuronal information processing. Our knowledge on neuronal structure and function is primarily based on brains of laboratory animals. Whether it translates to human is not known since quantitative data on "full" human neuronal morphologies are lacking. Here, we obtained human brain tissue during resection surgery and reconstructed basal and apical dendrites and axons of individual neurons across all cortical layers in temporal cortex (Brodmann area 21). Importantly, morphologies did not correlate to etiology, disease severity, or disease duration. Next, we show that human L(ayer) 2 and L3 pyramidal neurons have 3-fold larger dendritic length and increased branch complexity with longer segments compared with temporal cortex neurons from macaque and mouse. Unsupervised cluster analysis classified 88% of human L2 and L3 neurons into human-specific clusters distinct from mouse and macaque neurons. Computational modeling of passive electrical properties to assess the functional impact of large dendrites indicates stronger signal attenuation of electrical inputs compared with mouse. We thus provide a quantitative analysis of "full" human neuron morphologies and present direct evidence that human neurons are not "scaled-up" versions of rodent or macaque neurons, but have unique structural and functional properties.

Dendrites impact the encoding capabilities of the axon.

This study highlights a new and powerful direct impact of the dendritic tree (the input region of neurons) on the encoding capability of the axon (the output region). We show that the size of the dendritic arbors (its impedance load) strongly modulates the shape of the action potential (AP) onset at the axon initial segment; it is accelerated in neurons with larger dendritic surface area. AP onset rapidness is key in determining the capability of the axonal spikes to encode (phase lock to) rapid changes in synaptic inputs. Hence, our findings imply that neurons with larger dendritic arbors have improved encoding capabilities. This "dendritic size effect" was explored both analytically as well as numerically, in simplified and detailed models of 3D reconstructed layer 2/3 cortical pyramidal cells of rats and humans. The cutoff frequency of spikes phase locking to modulated inputs increased from 100 to 200 Hz in pyramidal cells of young rats to 400-600 Hz in human cells. In the latter case, phase locking reached close to 1 KHz in in vivo-like conditions. This work highlights new and functionally profound cross talk between the dendritic tree and the axon initial segment, providing new understanding of neurons as sophisticated nonlinear input/output devices.

Spatially distributed dendritic resonance selectively filters synaptic input.

An important task performed by a neuron is the selection of relevant inputs from among thousands of synapses impinging on the dendritic tree. Synaptic plasticity enables this by strenghtening a subset of synapses that are, presumably, functionally relevant to the neuron. A different selection mechanism exploits the resonance of the dendritic membranes to preferentially filter synaptic inputs based on their temporal rates. A widely held view is that a neuron has one resonant frequency and thus can pass through one rate. Here we demonstrate through mathematical analyses and numerical simulations that dendritic resonance is inevitably a spatially distributed property; and therefore the resonance frequency varies along the dendrites, and thus endows neurons with a powerful spatiotemporal selection mechanism that is sensitive both to the dendritic location and the temporal structure of the incoming synaptic inputs.

A hierarchical structure of cortical interneuron electrical diversity revealed by automated statistical analysis.

Although the diversity of cortical interneuron electrical properties is well recognized, the number of distinct electrical types (e-types) is still a matter of debate. Recently, descriptions of interneuron variability were standardized by multiple laboratories on the basis of a subjective classification scheme as set out by the Petilla convention (Petilla Interneuron Nomenclature Group, PING). Here, we present a quantitative, statistical analysis of a database of nearly five hundred neurons manually annotated according to the PING nomenclature. For each cell, 38 features were extracted from responses to suprathreshold current stimuli and statistically analyzed to examine whether cortical interneurons subdivide into e-types. We showed that the partitioning into different e-types is indeed the major component of data variability. The analysis suggests refining the PING e-type classification to be hierarchical, whereby most variability is first captured within a coarse subpartition, and then subsequently divided into finer subpartitions. The coarse partition matches the well-known partitioning of interneurons into fast spiking and adapting cells. Finer subpartitions match the burst, continuous, and delayed subtypes. Additionally, our analysis enabled the ranking of features according to their ability to differentiate among e-types. We showed that our quantitative e-type assignment is more than 90% accurate and manages to catch several human errors.

Morpho-electric diversity of human hippocampal CA1 pyramidal neurons.

Hippocampal pyramidal neuron activity underlies episodic memory and spatial navigation. Although extensively studied in rodents, extremely little is known about human hippocampal pyramidal neurons, even though the human hippocampus underwent strong evolutionary reorganization and shows lower theta rhythm frequencies. To test whether biophysical properties of human Cornu Amonis subfield 1 (CA1) pyramidal neurons can explain observed rhythms, we map the morpho-electric properties of individual CA1 pyramidal neurons in human, non-pathological hippocampal slices from neurosurgery. Human CA1 pyramidal neurons have much larger dendritic trees than mouse CA1 pyramidal neurons, have a large number of oblique dendrites, and resonate at 2.9 Hz, optimally tuned to human theta frequencies. Morphological and biophysical properties suggest cellular diversity along a multidimensional gradient rather than discrete clustering. Across the population, dendritic architecture and a large number of oblique dendrites consistently boost memory capacity in human CA1 pyramidal neurons by an order of magnitude compared to mouse CA1 pyramidal neurons.

Preserving axosomatic spiking features despite diverse dendritic morphology.

Throughout the nervous system, cells belonging to a certain electrical class (e-class)-sharing high similarity in firing response properties-may nevertheless have widely variable dendritic morphologies. To quantify the effect of this morphological variability on the firing of layer 5 thick-tufted pyramidal cells (TTCs), a detailed conductance-based model was constructed for a three-dimensional reconstructed exemplar TTC. The model exhibited spike initiation in the axon and reproduced the characteristic features of individual spikes, as well as of the firing properties at the soma, as recorded in a population of TTCs in young Wistar rats. When using these model parameters over the population of 28 three-dimensional reconstructed TTCs, both axonal and somatic ion channel densities had to be scaled linearly with the conductance load imposed on each of these compartments. Otherwise, the firing of model cells deviated, sometimes very significantly, from the experimental variability of the TTC e-class. The study provides experimentally testable predictions regarding the coregulation of axosomatic membrane ion channels density for cells with different dendritic conductance load, together with a simple and systematic method for generating reliable conductance-based models for the whole population of modeled neurons belonging to a particular e-class, with variable morphology as found experimentally.