Reconstruction and Simulation of Neocortical Microcircuitry.

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

In-silico EEG biomarkers of reduced inhibition in human cortical microcircuits in depression.

Reduced cortical inhibition by somatostatin-expressing (SST) interneurons has been strongly associated with treatment-resistant depression. However, due to technical limitations it is impossible to establish experimentally in humans whether the effects of reduced SST interneuron inhibition on microcircuit activity have signatures detectable in clinically-relevant brain signals such as electroencephalography (EEG). To overcome these limitations, we simulated resting-state activity and EEG using detailed models of human cortical microcircuits with normal (healthy) or reduced SST interneuron inhibition (depression), and found that depression microcircuits exhibited increased theta, alpha and low beta power (4-16 Hz). The changes in depression involved a combination of an aperiodic broadband and periodic theta components. We then demonstrated the specificity of the EEG signatures of reduced SST interneuron inhibition by showing they were distinct from those corresponding to reduced parvalbumin-expressing (PV) interneuron inhibition. Our study thus links SST interneuron inhibition level to distinct features in EEG simulated from detailed human microcircuits, which can serve to better identify mechanistic subtypes of depression using EEG, and non-invasively monitor modulation of cortical inhibition.

Age-dependent increased sag amplitude in human pyramidal neurons dampens baseline cortical activity.

Aging involves various neurobiological changes, although their effect on brain function in humans remains poorly understood. The growing availability of human neuronal and circuit data provides opportunities for uncovering age-dependent changes of brain networks and for constraining models to predict consequences on brain activity. Here we found increased sag voltage amplitude in human middle temporal gyrus layer 5 pyramidal neurons from older subjects and captured this effect in biophysical models of younger and older pyramidal neurons. We used these models to simulate detailed layer 5 microcircuits and found lower baseline firing in older pyramidal neuron microcircuits, with minimal effect on response. We then validated the predicted reduced baseline firing using extracellular multielectrode recordings from human brain slices of different ages. Our results thus report changes in human pyramidal neuron input integration properties and provide fundamental insights into the neuronal mechanisms of altered cortical excitability and resting-state activity in human aging.

Orientation processing by synaptic integration across first-order tactile neurons.

Our ability to manipulate objects relies on tactile inputs from first-order tactile neurons that innervate the glabrous skin of the hand. The distal axon of these neurons branches in the skin and innervates many mechanoreceptors, yielding spatially-complex receptive fields. Here we show that synaptic integration across the complex signals from the first-order neuronal population could underlie human ability to accurately (< 3°) and rapidly process the orientation of edges moving across the fingertip. We first derive spiking models of human first-order tactile neurons that fit and predict responses to moving edges with high accuracy. We then use the model neurons in simulating the peripheral neuronal population that innervates a fingertip. We train classifiers performing synaptic integration across the neuronal population activity, and show that synaptic integration across first-order neurons can process edge orientations with high acuity and speed. In particular, our models suggest that integration of fast-decaying (AMPA-like) synaptic inputs within short timescales is critical for discriminating fine orientations, whereas integration of slow-decaying (NMDA-like) synaptic inputs supports discrimination of coarser orientations and maintains robustness over longer timescales. Taken together, our results provide new insight into the computations occurring in the earliest stages of the human tactile processing pathway and how they may be critical for supporting hand function.

Multiregional integration in the brain during resting-state fMRI activity.

Data-driven models of functional magnetic resonance imaging (fMRI) activity can elucidate dependencies that involve the combination of multiple brain regions. Activity in some regions during resting-state fMRI can be predicted with high accuracy from the activities of other regions. However, it remains unclear in which regions activity depends on unique integration of multiple predictor regions. To address this question, sparse (parsimonious) models could serve to better determine key interregional dependencies by reducing false positives. We used resting-state fMRI data from 46 subjects, and for each region of interest (ROI) per subject we performed whole-brain recursive feature elimination (RFE) to select the minimal set of ROIs that best predicted activity in the modeled ROI. We quantified the dependence of activity on multiple predictor ROIs, by measuring the gain in prediction accuracy of models that incorporated multiple predictor ROIs compared to models that used a single predictor ROI. We identified regions that showed considerable evidence of multiregional integration and determined the key regions that contributed to their observed activity. Our models reveal fronto-parietal integration networks, little integration in primary sensory regions, as well as redundancy between some regions. Our study demonstrates the utility of whole-brain RFE to generate data-driven models with minimal sets of ROIs that predict activity with high accuracy. By determining the extent to which activity in each ROI depended on integration of signals from multiple ROIs, we find cortical integration networks during resting-state activity.

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.

In-silico testing of new pharmacology for restoring inhibition and human cortical function in depression.

Reduced inhibition by somatostatin-expressing interneurons is associated with depression. Administration of positive allosteric modulators of α5 subunit-containing GABAA receptor (α5-PAM) that selectively target this lost inhibition exhibit antidepressant and pro-cognitive effects in rodent models of chronic stress. However, the functional effects of α5-PAM on the human brain in vivo are unknown, and currently cannot be assessed experimentally. We modeled the effects of α5-PAM on tonic inhibition as measured in human neurons, and tested in silico α5-PAM effects on detailed models of human cortical microcircuits in health and depression. We found that α5-PAM effectively recovered impaired cortical processing as quantified by stimulus detection metrics, and also recovered the power spectral density profile of the microcircuit EEG signals. We performed an α5-PAM dose-response and identified simulated EEG biomarker candidates. Our results serve to de-risk and facilitate α5-PAM translation and provide biomarkers in non-invasive brain signals for monitoring target engagement and drug efficacy.

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.

Models of neocortical layer 5b pyramidal cells capturing a wide range of dendritic and perisomatic active properties.

The thick-tufted layer 5b pyramidal cell extends its dendritic tree to all six layers of the mammalian neocortex and serves as a major building block for the cortical column. L5b pyramidal cells have been the subject of extensive experimental and modeling studies, yet conductance-based models of these cells that faithfully reproduce both their perisomatic Na(+)-spiking behavior as well as key dendritic active properties, including Ca(2+) spikes and back-propagating action potentials, are still lacking. Based on a large body of experimental recordings from both the soma and dendrites of L5b pyramidal cells in adult rats, we characterized key features of the somatic and dendritic firing and quantified their statistics. We used these features to constrain the density of a set of ion channels over the soma and dendritic surface via multi-objective optimization with an evolutionary algorithm, thus generating a set of detailed conductance-based models that faithfully replicate the back-propagating action potential activated Ca(2+) spike firing and the perisomatic firing response to current steps, as well as the experimental variability of the properties. Furthermore, we show a useful way to analyze model parameters with our sets of models, which enabled us to identify some of the mechanisms responsible for the dynamic properties of L5b pyramidal cells as well as mechanisms that are sensitive to morphological changes. This automated framework can be used to develop a database of faithful models for other neuron types. The models we present provide several experimentally-testable predictions and can serve as a powerful tool for theoretical investigations of the contribution of single-cell dynamics to network activity and its computational capabilities.

Reduced inhibition in depression impairs stimulus processing in human cortical microcircuits.

Cortical processing depends on finely tuned excitatory and inhibitory connections in neuronal microcircuits. Reduced inhibition by somatostatin-expressing interneurons is a key component of altered inhibition associated with treatment-resistant major depressive disorder (depression), which is implicated in cognitive deficits and rumination, but the link remains to be better established mechanistically in humans. Here we test the effect of reduced somatostatin interneuron-mediated inhibition on cortical processing in human neuronal microcircuits using a data-driven computational approach. We integrate human cellular, circuit, and gene expression data to generate detailed models of human cortical microcircuits in health and depression. We simulate microcircuit baseline and response activity and find a reduced signal-to-noise ratio and increased false/failed detection of stimuli due to a higher baseline activity in depression. We thus apply models of human cortical microcircuits to demonstrate mechanistically how reduced inhibition impairs cortical processing in depression, providing quantitative links between altered inhibition and cognitive deficits.

An in vitro whole-cell electrophysiology dataset of human cortical neurons.

BACKGROUND: Whole-cell patch-clamp electrophysiology is an essential technique for understanding how single neurons translate their diverse inputs into a functional output. The relative inaccessibility of live human cortical neurons for experimental manipulation has made it difficult to determine the unique features of how human cortical neurons differ from their counterparts in other species. FINDINGS: We present a curated repository of whole-cell patch-clamp recordings from surgically resected human cortical tissue, encompassing 118 neurons from 35 individuals (age range, 21-59 years; 17 male, 18 female). Recorded human cortical neurons derive from layers 2 and 3 (L2&3), deep layer 3 (L3c), or layer 5 (L5) and are annotated with a rich set of subject and experimental metadata. For comparison, we also provide a limited set of comparable recordings from 21-day-old mice (11 cells from 5 mice). All electrophysiological recordings are provided in the Neurodata Without Borders (NWB) format and are available for further analysis via the Distributed Archives for Neurophysiology Data Integration online repository. The associated data conversion code is made publicly available and can help others in converting electrophysiology datasets to the open NWB standard for general reuse. CONCLUSION: These data can be used for novel analyses of biophysical characteristics of human cortical neurons, including in cross-species or cross-lab comparisons or in building computational models of individual human neurons.