Controlling morpho-electrophysiological variability of neurons with detailed biophysical models

Variability, which is known to be a universal feature among biological units such as neuronal cells, holds significant importance, as, for example, it enables a robust encoding of a high volume of information in neuronal circuits and prevents hypersynchronizations. While most computational studies on electrophysiological variability in neuronal circuits were done with single-compartment neuron models, we instead focus on the variability of detailed biophysical models of neuron multi-compartmental morphologies. We leverage a Markov chain Monte Carlo method to generate populations of electrical models reproducing the variability of experimental recordings while being compatible with a set of morphologies to faithfully represent specifi morpho-electrical type. We demonstrate our approach on layer 5 pyramidal cells and study the morpho-electrical variability and in particular, find that morphological variability alone is insufficient to reproduce electrical variability. Overall, this approach provides a strong statistical basis to create detailed models of neurons with controlled variability.

A universal workflow for creation, validation, and generalization of detailed neuronal models

Detailed single-neuron modeling is widely used to study neuronal functions. While cellular and functional diversity across the mammalian cortex is vast, most of the available computational tools focus on a limited set of specific features characteristic of a single neuron. Here, we present a generalized automated workflow for the creation of robust electrical models and illustrate its performance by building cell models for the rat somatosensory cortex. Each model is based on a 3D morphological reconstruction and a set of ionic mechanisms. We use an evolutionary algorithm to optimize neuronal parameters to match the electrophysiological features extracted from experimental data. Then we validate the optimized models against additional stimuli and assess their generalizability on a population of similar morphologies. Compared to the state-of-the-art canonical models, our models show 5-fold improved generalizability. This versatile approach can be used to build robust models of any neuronal type.

Overexpression of UCP4 in astrocytic mitochondria prevents multilevel dysfunctions in a mouse model of Alzheimer's disease

Alzheimer's disease (AD) is becoming increasingly prevalent worldwide. It represents one of the greatest medical challenges as no pharmacologic treatments are available to prevent disease progression. Astrocytes play crucial functions within neuronal circuits by providing metabolic and functional support, regulating interstitial solute composition, and modulating synaptic transmission. In addition to these physiological functions, growing evidence points to an essential role of astrocytes in neurodegenerative diseases like AD. Early-stage AD is associated with hypometabolism and oxidative stress. Contrary to neurons that are vulnerable to oxidative stress, astrocytes are particularly resistant to mitochondrial dysfunction and are therefore more resilient cells. In our study, we leveraged astrocytic mitochondrial uncoupling and examined neuronal function in the 3xTg AD mouse model. We overexpressed the mitochondrial uncoupling protein 4 (UCP4), which has been shown to improve neuronal survival in vitro. We found that this treatment efficiently prevented alterations of hippocampal metabolite levels observed in AD mice, along with hippocampal atrophy and reduction of basal dendrite arborization of subicular neurons. This approach also averted aberrant neuronal excitability observed in AD subicular neurons and preserved episodic-like memory in AD mice assessed in a spatial recognition task. These findings show that targeting astrocytes and their mitochondria is an effective strategy to prevent the decline of neurons facing AD-related stress at the early stages of the disease.

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.

Cell-type-specific densities in mouse somatosensory cortex derived from scRNA-seq and in situ RNA hybridization

Cells in the mammalian cerebral cortex exhibit layer-dependent patterns in their distribution. Classical methods of determining cell type distributions typically employ a painstaking process of large-scale sampling and characterization of cellular composition. We found that by combining in situ hybridization (ISH) images with cell-type-specific transcriptomes, position-dependent cortical composition in P56 mouse could be estimated in the somatosensory cortex. The method uses ISH images from the Allen Institute for Brain Science. There are two novel aspects of the methodology. First, it is not necessary to select a subset of genes that are particular for a cell type of interest, nor is it necessary to only use ISH images with low variability among samples. Second, the method also compensated for differences in soma size and incompleteness of the transcriptomes. The soma size compensation is particularly important in order to obtain quantitative estimates since relying on bulk expression alone would overestimate the contribution of larger cells. Predicted distributions of broader classes of cell types agreed with literature distributions. The primary result is that there is a high degree of substructure in the distribution of transcriptomic types beyond the resolution of layers. Furthermore, transcriptomic cell types each exhibited characteristic soma size distributions. Results suggest that the method could also be employed to assign transcriptomic cell types to well-aligned image sets in the entire brain.

Mapping of morpho-electric features to molecular identity of cortical inhibitory neurons

Knowledge of the cell-type-specific composition of the brain is useful in order to understand the role of each cell type as part of the network. Here, we estimated the composition of the whole cortex in terms of well characterized morphological and electrophysiological inhibitory neuron types (me-types). We derived probabilistic me-type densities from an existing atlas of molecularly defined cell-type densities in the mouse cortex. We used a well-established me-type classification from rat somatosensory cortex to populate the cortex. These me-types were well characterized morphologically and electrophysiologically but they lacked molecular marker identity labels. To extrapolate this missing information, we employed an additional dataset from the Allen Institute for Brain Science containing molecular identity as well as morphological and electrophysiological data for mouse cortical neurons. We first built a latent space based on a number of comparable morphological and electrical features common to both data sources. We then identified 19 morpho-electrical clusters that merged neurons from both datasets while being molecularly homogeneous. The resulting clusters best mirror the molecular identity classification solely using available morpho-electrical features. Finally, we stochastically assigned a molecular identity to a me-type neuron based on the latent space cluster it was assigned to. The resulting mapping was used to derive inhibitory me-types densities in the cortex.

Optimum trajectory learning in musculoskeletal systems with model predictive control and deep reinforcement learning

From the computational point of view, musculoskeletal control is the problem of controlling high degrees of freedom and dynamic multi-body system that is driven by redundant muscle units. A critical challenge in the control perspective of skeletal joints with antagonistic muscle pairs is finding methods robust to address this ill-posed nonlinear problem. To address this computational problem, we implemented a twofold optimization and learning framework to be specialized in addressing the redundancies in the muscle control . In the first part, we used model predictive control to obtain energy efficient skeletal trajectories to mimick human movements. The second part is to use deep reinforcement learning to obtain a sequence of stimulus to be given to muscles in order to obtain the skeletal trajectories with muscle control. We observed that the desired stimulus to muscles is only efficiently constructed by integrating the state and control input in a closed-loop setting as it resembles the proprioceptive integration in the spinal cord circuits. In this work, we showed how a variety of different reference trajectories can be obtained with optimal control and how these reference trajectories are mapped to the musculoskeletal control with deep reinforcement learning. Starting from the characteristics of human arm movement to obstacle avoidance experiment, our simulation results confirm the capabilities of our optimization and learning framework for a variety of dynamic movement trajectories. In summary, the proposed framework is offering a pipeline to complement the lack of experiments to record human motion-capture data as well as study the activation range of muscles to replicate the specific trajectory of interest. Using the trajectories from optimal control as a reference signal for reinforcement learning implementation has allowed us to acquire optimum and human-like behaviour of the musculoskeletal system which provides a framework to study human movement in-silico experiments. The present framework can also allow studying upper-arm rehabilitation with assistive robots given that one can use healthy subject movement recordings as reference to work on the control architecture of assistive robotics in order to compensate behavioural deficiencies. Hence, the framework opens to possibility of replicating or complementing labour-intensive, time-consuming and costly experiments with human subjects in the field of movement studies and digital twin of rehabilitation.

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser).

SIGNIFICANCE

There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

A calcium-based plasticity model for predicting long-term potentiation and depression in the neocortex

Pyramidal cells (PCs) form the backbone of the layered structure of the neocortex, and plasticity of their synapses is thought to underlie learning in the brain. However, such long-term synaptic changes have been experimentally characterized between only a few types of PCs, posing a significant barrier for studying neocortical learning mechanisms. Here we introduce a model of synaptic plasticity based on data-constrained postsynaptic calcium dynamics, and show in a neocortical microcircuit model that a single parameter set is sufficient to unify the available experimental findings on long-term potentiation (LTP) and long-term depression (LTD) of PC connections. In particular, we find that the diverse plasticity outcomes across the different PC types can be explained by cell-type-specific synaptic physiology, cell morphology and innervation patterns, without requiring type-specific plasticity. Generalizing the model to in vivo extracellular calcium concentrations, we predict qualitatively different plasticity dynamics from those observed in vitro. This work provides a first comprehensive null model for LTP/LTD between neocortical PC types in vivo, and an open framework for further developing models of cortical synaptic plasticity.

Representing Stimulus Information in an Energy Metabolism Pathway

We explored a computational model of astrocytic energy metabolism and demonstrated the theoretical plausibility that this type of pathway might be capable of coding information about stimuli in addition to its known functions in cellular energy and carbon budgets. Simulation results indicate that glycogenolytic glycolysis triggered by activation of adrenergic receptors can capture the intensity and duration features of a neuromodulator waveform and can respond in a dose-dependent manner, including non-linear state changes that are analogous to action potentials. We show how this metabolic pathway can translate information about external stimuli to production profiles of energy-carrying molecules such as lactate with a precision beyond simple signal transduction or non-linear amplification. The results suggest the operation of a metabolic state-machine from the spatially discontiguous yet interdependent metabolite elements. Such metabolic pathways might be well-positioned to code an additional level of salient information about a cell's environmental demands to impact its function. Our hypothesis has implications for the computational power and energy efficiency of the brain.

A Standardized Brain Molecular Atlas: A Resource for Systems Modeling and Simulation

Accurate molecular concentrations are essential for reliable analyses of biochemical networks and the creation of predictive models for molecular and systems biology, yet protein and metabolite concentrations used in such models are often poorly constrained or irreproducible. Challenges of using data from different sources include conflicts in nomenclature and units, as well as discrepancies in experimental procedures, data processing and implementation of the model. To obtain a consistent estimate of protein and metabolite levels, we integrated and normalized data from a large variety of sources to calculate Adjusted Molecular Concentrations. We found a high degree of reproducibility and consistency of many molecular species across brain regions and cell types, consistent with tight homeostatic regulation. We demonstrated the value of this normalization with differential protein expression analyses related to neurodegenerative diseases, brain regions and cell types. We also used the results in proof-of-concept simulations of brain energy metabolism. The standardized Brain Molecular Atlas overcomes the obstacles of missing or inconsistent data to support systems biology research and is provided as a resource for biomolecular modeling.

Digital Reconstruction of the Neuro-Glia-Vascular Architecture

Astrocytes connect the vasculature to neurons mediating the supply of nutrients and biochemicals. They are involved in a growing number of physiological and pathophysiological processes that result from biophysical, physiological, and molecular interactions in this neuro-glia-vascular ensemble (NGV). The lack of a detailed cytoarchitecture severely restricts the understanding of how they support brain function. To address this problem, we used data from multiple sources to create a data-driven digital reconstruction of the NGV at micrometer anatomical resolution. We reconstructed 0.2 mm3 of the rat somatosensory cortex with 16 000 morphologically detailed neurons, 2500 protoplasmic astrocytes, and its microvasculature. The consistency of the reconstruction with a wide array of experimental measurements allows novel predictions of the NGV organization, allowing the anatomical reconstruction of overlapping astrocytic microdomains and the quantification of endfeet connecting each astrocyte to the vasculature, as well as the extent to which they cover the latter. Structural analysis showed that astrocytes optimize their positions to provide uniform vascular coverage for trophic support and signaling. However, this optimal organization rapidly declines as their density increases. The NGV digital reconstruction is a resource that will enable a better understanding of the anatomical principles and geometric constraints, which govern how astrocytes support brain function.

The Role of Hub Neurons in Modulating Cortical Dynamics

Many neurodegenerative diseases are associated with the death of specific neuron types in particular brain regions. What makes the death of specific neuron types particularly harmful for the integrity and dynamics of the respective network is not well understood. To start addressing this question we used the most up-to-date biologically realistic dense neocortical microcircuit (NMC) of the rodent, which has reconstructed a volume of 0.3 mm³ and containing 31,000 neurons, ∼37 million synapses, and 55 morphological cell types arranged in six cortical layers. Using modern network science tools, we identified hub neurons in the NMC, that are connected synaptically to a large number of their neighbors and systematically examined the impact of abolishing these cells. In general, the structural integrity of the network is robust to cells’ attack; yet, attacking hub neurons strongly impacted the small-world topology of the network, whereas similar attacks on random neurons have a negligible effect. Such hub-specific attacks are also impactful on the network dynamics, both when the network is at its spontaneous synchronous state and when it was presented with synchronized thalamo-cortical visual-like input. We found that attacking layer 5 hub neurons is most harmful to the structural and functional integrity of the NMC. The significance of our results for understanding the role of specific neuron types and cortical layers for disease manifestation is discussed.

Morphology, physiology and synaptic connectivity of local interneurons in the mouse somatosensory thalamus

Abstract

The thalamic reticular nucleus (TRN) neurons, projecting across the external medullary lamina, have long been considered to be the only significant source of inhibition of the somatosensory ventral posterior (VP) nuclei of the thalamus. Here we report for the first time effective local inhibition and disinhibition in the VP. Inhibitory interneurons were found in GAD67–GFP‐expressing mice and studied using in vitro multiple patch clamp. Inhibitory interneurons have expansive bipolar or tripolar morphologies, reach across most of the VP nucleus and display low threshold bursting behaviour. They form triadic and non‐triadic synaptic connections onto thalamocortical relay neurons and other interneurons, mediating feedforward inhibition and disinhibition. Synaptic inputs arrive before those expected from the TRN neurons, suggesting that local inhibition plays an early and significant role in the functioning of the somatosensory thalamus.

Key points

The physiology and structure of local interneurons in the mouse somatosensory thalamus is described for the first time.

Inhibitory interneurons have extensive dendritic arborization providing significant local dendro‐dendritic inhibition in the somatosensory thalamus.

Triadic and non‐triadic synaptic connectivity onto thalamic relay neurons and other interneurons provides both local feedforward inhibition and disinhibition.

Interneurons of the somatosensory thalamus provide inhibition before the thalamic reticular nucleus, suggesting they play an important role in sensory perception.

Supervised Learning With Perceptual Similarity for Multimodal Gene Expression Registration of a Mouse Brain Atlas

The acquisition of high quality maps of gene expression in the rodent brain is of fundamental importance to the neuroscience community. The generation of such datasets relies on registering individual gene expression images to a reference volume, a task encumbered by the diversity of staining techniques employed, and by deformations and artifacts in the soft tissue. Recently, deep learning models have garnered particular interest as a viable alternative to traditional intensity-based algorithms for image registration. In this work, we propose a supervised learning model for general multimodal 2D registration tasks, trained with a perceptual similarity loss on a dataset labeled by a human expert and augmented by synthetic local deformations. We demonstrate the results of our approach on the Allen Mouse Brain Atlas (AMBA), comprising whole brain Nissl and gene expression stains. We show that our framework and design of the loss function result in accurate and smooth predictions. Our model is able to generalize to unseen gene expressions and coronal sections, outperforming traditional intensity-based approaches in aligning complex brain structures.

A Machine-Generated View of the Role of Blood Glucose Levels in the Severity of COVID-19

SARS-CoV-2 started spreading toward the end of 2019 causing COVID-19, a disease that reached pandemic proportions among the human population within months. The reasons for the spectrum of differences in the severity of the disease across the population, and in particular why the disease affects more severely the aging population and those with specific preconditions are unclear. We developed machine learning models to mine 240,000 scientific articles openly accessible in the CORD-19 database, and constructed knowledge graphs to synthesize the extracted information and navigate the collective knowledge in an attempt to search for a potential common underlying reason for disease severity. The machine-driven framework we developed repeatedly pointed to elevated blood glucose as a key facilitator in the progression of COVID-19. Indeed, when we systematically retraced the steps of the SARS-CoV-2 infection, we found evidence linking elevated glucose to each major step of the life-cycle of the virus, progression of the disease, and presentation of symptoms. Specifically, elevations of glucose provide ideal conditions for the virus to evade and weaken the first level of the immune defense system in the lungs, gain access to deep alveolar cells, bind to the ACE2 receptor and enter the pulmonary cells, accelerate replication of the virus within cells increasing cell death and inducing an pulmonary inflammatory response, which overwhelms an already weakened innate immune system to trigger an avalanche of systemic infections, inflammation and cell damage, a cytokine storm and thrombotic events. We tested the feasibility of the hypothesis by manually reviewing the literature referenced by the machine-generated synthesis, reconstructing atomistically the virus at the surface of the pulmonary airways, and performing quantitative computational modeling of the effects of glucose levels on the infection process. We conclude that elevation in glucose levels can facilitate the progression of the disease through multiple mechanisms and can explain much of the differences in disease severity seen across the population. The study provides diagnostic considerations, new areas of research and potential treatments, and cautions on treatment strategies and critical care conditions that induce elevations in blood glucose levels.

Metaball skinning of synthetic astroglial morphologies into realistic mesh models for visual analytics and in silico simulations

Motivation

Astrocytes, the most abundant glial cells in the mammalian brain, have an instrumental role in developing neuronal circuits. They contribute to the physical structuring of the brain, modulating synaptic activity and maintaining the blood–brain barrier in addition to other significant aspects that impact brain function. Biophysically, detailed astrocytic models are key to unraveling their functional mechanisms via molecular simulations at microscopic scales. Detailed, and complete, biological reconstructions of astrocytic cells are sparse. Nonetheless, data-driven digital reconstruction of astroglial morphologies that are statistically identical to biological counterparts are becoming available. We use those synthetic morphologies to generate astrocytic meshes with realistic geometries, making it possible to perform these simulations.

Results

We present an unconditionally robust method capable of reconstructing high fidelity polygonal meshes of astroglial cells from algorithmically-synthesized morphologies. Our method uses implicit surfaces, or metaballs, to skin the different structural components of astrocytes and then blend them in a seamless fashion. We also provide an end-to-end pipeline to produce optimized two- and three-dimensional meshes for visual analytics and simulations, respectively. The performance of our pipeline has been assessed with a group of 5000 astroglial morphologies and the geometric metrics of the resulting meshes are evaluated. The usability of the meshes is then demonstrated with different use cases.

Availability and implementation

Our metaball skinning algorithm is implemented in Blender 2.82 relying on its Python API (Application Programming Interface). To make it accessible to computational biologists and neuroscientists, the implementation has been integrated into NeuroMorphoVis, an open source and domain specific package that is primarily designed for neuronal morphology visualization and meshing.

Supplementary information

Supplementary data are available at Bioinformatics online.

In silico voltage-sensitive dye imaging reveals the emergent dynamics of cortical populations

Voltage-sensitive dye imaging (VSDI) is a powerful technique for interrogating membrane potential dynamics in assemblies of cortical neurons, but with effective resolution limits that confound interpretation. To address this limitation, we developed an in silico model of VSDI in a biologically faithful digital reconstruction of rodent neocortical microcircuitry. Using this model, we extend previous experimental observations regarding the cellular origins of VSDI, finding that the signal is driven primarily by neurons in layers 2/3 and 5, and that VSDI measurements do not capture individual spikes. Furthermore, we test the capacity of VSD image sequences to discriminate between afferent thalamic inputs at various spatial locations to estimate a lower bound on the functional resolution of VSDI. Our approach underscores the power of a bottom-up computational approach for relating scales of cortical processing.

ARC: An Open Web-Platform for Request/Supply Matching for a Prioritized and Controlled COVID-19 Response

In 2020 the world was hit by the COVID-19 pandemic putting entire governments and civil societies in crisis mode. Around the globe unprecedented shortages of equipment and qualified personnel were reported in hospitals and diagnostic laboratories. When a crisis is global, supply chains are strained worldwide and external help may not be readily available. In Switzerland, as part of the efforts of the Swiss National COVID-19 Science Task Force, we developed a tailor-made web-based tool where needs and offers for critical laboratory equipment and expertise can be brought together, coordinated, prioritized, and validated. This Academic Resources for COVID-19 (ARC) Platform presents the specialized needs of diagnostic laboratories to academic research groups at universities, allowing the sourcing of said needs from unconventional supply channels, while keeping the entities tasked with coordination of the crisis response in control of each part of the process. An instance of the ARC Platform is operated in Switzerland (arc.epfl.ch) catering to the diagnostic efforts in Switzerland and sourcing from the Swiss academic sector. The underlying technology has been released as open source so that others can adopt the customizable web-platform for need/supply match-making in their own relief efforts, during the COVID-19 pandemic or any future disaster.

Interactive visualization and analysis of morphological skeletons of brain vasculature networks with VessMorphoVis

MOTIVATION

Accurate morphological models of brain vasculature are key to modeling and simulating cerebral blood flow in realistic vascular networks. This in silico approach is fundamental to revealing the principles of neurovascular coupling. Validating those vascular morphologies entails performing certain visual analysis tasks that cannot be accomplished with generic visualization frameworks. This limitation has a substantial impact on the accuracy of the vascular models employed in the simulation.

RESULTS

We present VessMorphoVis, an integrated suite of toolboxes for interactive visualization and analysis of vast brain vascular networks represented by morphological graphs segmented originally from imaging or microscopy stacks. Our workflow leverages the outstanding potentials of Blender, aiming to establish an integrated, extensible and domain-specific framework capable of interactive visualization, analysis, repair, high-fidelity meshing and high-quality rendering of vascular morphologies. Based on the initial feedback of the users, we anticipate that our framework will be an essential component in vascular modeling and simulation in the future, filling a gap that is at present largely unfulfilled.

AVAILABILITY AND IMPLEMENTATION

VessMorphoVis is freely available under the GNU public license on Github at https://github.com/BlueBrain/VessMorphoVis. The morphology analysis, visualization, meshing and rendering modules are implemented as an add-on for Blender 2.8 based on its Python API (application programming interface). The add-on functionality is made available to users through an intuitive graphical user interface, as well as through exhaustive configuration files calling the API via a feature-rich command line interface running Blender in background mode.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Objective Morphological Classification of Neocortical Pyramidal Cells

A consensus on the number of morphologically different types of pyramidal cells (PCs) in the neocortex has not yet been reached, despite over a century of anatomical studies, due to the lack of agreement on the subjective classifications of neuron types, which is based on expert analyses of neuronal morphologies. Even for neurons that are visually distinguishable, there is no common ground to consistently define morphological types. The objective classification of PCs can be achieved with methods from algebraic topology, and the dendritic arborization is sufficient for the reliable identification of distinct types of cortical PCs. Therefore, we objectively identify 17 types of PCs in the rat somatosensory cortex. In addition, we provide a solution to the challenging problem of whether 2 similar neurons belong to different types or to a continuum of the same type. Our topological classification does not require expert input, is stable, and helps settle the long-standing debate on whether cell-types are discrete or continuous morphological variations of each other.

Impact of higher order network structure on emergent cortical activity

Synaptic connectivity between neocortical neurons is highly structured. The network structure of synaptic connectivity includes first-order properties that can be described by pairwise statistics, such as strengths of connections between different neuron types and distance-dependent connectivity, and higher order properties, such as an abundance of cliques of all-to-all connected neurons. The relative impact of first- and higher order structure on emergent cortical network activity is unknown. Here, we compare network structure and emergent activity in two neocortical microcircuit models with different synaptic connectivity. Both models have a similar first-order structure, but only one model includes higher order structure arising from morphological diversity within neuronal types. We find that such morphological diversity leads to more heterogeneous degree distributions, increases the number of cliques, and contributes to a small-world topology. The increase in higher order network structure is accompanied by more nuanced changes in neuronal firing patterns, such as an increased dependence of pairwise correlations on the positions of neurons in cliques. Our study shows that circuit models with very similar first-order structure of synaptic connectivity can have a drastically different higher order network structure, and suggests that the higher order structure imposed by morphological diversity within neuronal types has an impact on emergent cortical activity.

Excitation states of metabolic networks predict dose-response fingerprinting and ligand pulse phase signalling

With a computational model of energy metabolism in an astrocyte, we show how a system of enzymes in a cascade can act as a functional unit of interdependent reactions, rather than merely a series of independent reactions. These systems may exist in multiple states, depending on the level of stimulation, and the effects of substrates at any point will depend on those states. Response trajectories of metabolites downstream from cAMP-stimulated glycogenolysis exhibit a host of non-linear dynamical response characteristics including hysteresis and response envelopes. Dose-dependent phase transitions predict a novel intracellular signalling mechanism and suggest a theoretical framework that could be relevant to single cell information processing, drug discovery or synthetic biology. Ligands may produce unique dose-response fingerprints depending on the state of the system, allowing selective output tuning. We conclude with the observation that state- and dose-dependent phase transitions, what we dub "ligand pulses" (LPs), may carry information and resemble action potentials (APs) generated from excitatory postsynaptic potentials. In our model, the relevant information from a cAMP-dependent glycolytic cascade in astrocytes could reflect the level of neuromodulatory input that signals an energy demand threshold. We propose that both APs and LPs represent specialized cases of molecular phase signalling with a common evolutionary root.

Estimating the Readily-Releasable Vesicle Pool Size at Synaptic Connections in the Neocortex

Previous studies based on the 'Quantal Model' for synaptic transmission suggest that neurotransmitter release is mediated by a single release site at individual synaptic contacts in the neocortex. However, recent studies seem to contradict this hypothesis and indicate that multi-vesicular release (MVR) could better explain the synaptic response variability observed in vitro. In this study we present a novel method to estimate the number of release sites per synapse, also known as the size of the readily releasable pool (N_RRP), from paired whole-cell recordings of connections between layer 5 thick tufted pyramidal cell (L5_TTPC) in the juvenile rat somatosensory cortex. Our approach extends the work of Loebel et al. (2009) by leveraging a recently published data-driven biophysical model of neocortical tissue. Using this approach, we estimated N_RRP to be between two to three for synaptic connections between L5_TTPCs. To constrain N_RRP values for other connections in the microcircuit, we developed and validated a generalization approach using published data on the coefficient of variation (CV) of the amplitudes of post-synaptic potentials (PSPs) from literature and comparing them against in silico experiments. Our study predicts that transmitter release at synaptic connections in the neocortex could be mediated by MVR and provides a data-driven approach to constrain the MVR model parameters in the microcircuit.

Caries Detection with Near-Infrared Transillumination Using Deep Learning

Dental caries is the most prevalent chronic condition worldwide. Early detection can significantly improve treatment outcomes and reduce the need for invasive procedures. Recently, near-infrared transillumination (TI) imaging has been shown to be effective for the detection of early stage lesions. In this work, we present a deep learning model for the automated detection and localization of dental lesions in TI images. Our method is based on a convolutional neural network (CNN) trained on a semantic segmentation task. We use various strategies to mitigate issues related to training data scarcity, class imbalance, and overfitting. With only 185 training samples, our model achieved an overall mean intersection-over-union (IOU) score of 72.7% on a 5-class segmentation task and specifically an IOU score of 49.5% and 49.0% for proximal and occlusal carious lesions, respectively. In addition, we constructed a simplified task, in which regions of interest were evaluated for the binary presence or absence of carious lesions. For this task, our model achieved an area under the receiver operating characteristic curve of 83.6% and 85.6% for occlusal and proximal lesions, respectively. Our work demonstrates that a deep learning approach for the analysis of dental images holds promise for increasing the speed and accuracy of caries detection, supporting the diagnoses of dental practitioners, and improving patient outcomes.

A null model of the mouse whole-neocortex micro-connectome

In connectomics, the study of the network structure of connected neurons, great advances are being made on two different scales: that of macro- and meso-scale connectomics, studying the connectivity between populations of neurons, and that of micro-scale connectomics, studying connectivity between individual neurons. We combine these two complementary views of connectomics to build a first draft statistical model of the micro-connectome of a whole mouse neocortex based on available data on region-to-region connectivity and individual whole-brain axon reconstructions. This process reveals a targeting principle that allows us to predict the innervation logic of individual axons from meso-scale data. The resulting connectome recreates biological trends of targeting on all scales and predicts that an established principle of scale invariant topological organization of connectivity can be extended down to the level of individual neurons. It can serve as a powerful null model and as a substrate for whole-brain simulations.

Cortical reliability amid noise and chaos

Typical responses of cortical neurons to identical sensory stimuli appear highly variable. It has thus been proposed that the cortex primarily uses a rate code. However, other studies have argued for spike-time coding under certain conditions. The potential role of spike-time coding is directly limited by the internally generated variability of cortical circuits, which remains largely unexplored. Here, we quantify this internally generated variability using a biophysical model of rat neocortical microcircuitry with biologically realistic noise sources. We find that stochastic neurotransmitter release is a critical component of internally generated variability, causing rapidly diverging, chaotic recurrent network dynamics. Surprisingly, the same nonlinear recurrent network dynamics can transiently overcome the chaos in response to weak feed-forward thalamocortical inputs, and support reliable spike times with millisecond precision. Our model shows that the noisy and chaotic network dynamics of recurrent cortical microcircuitry are compatible with stimulus-evoked, millisecond spike-time reliability, resolving a long-standing debate.

A Kinetic Map of the Homomeric Voltage-Gated Potassium Channel (Kv) Family

The voltage-gated potassium (Kv) channels, encoded by 40 genes, repolarize all electrically excitable cells, including plant, cardiac, and neuronal cells. Although these genes were fully sequenced decades ago, a comprehensive kinetic characterization of all Kv channels is still missing, especially near physiological temperature. Here, we present a standardized kinetic map of the 40 homomeric Kv channels systematically characterized at 15, 25, and 35°C. Importantly, the Kv kinetics at 35°C differ significantly from commonly reported kinetics, usually performed at room temperature. We observed voltage-dependent Q₁₀ for all active Kv channels and inherent heterogeneity in kinetics for some of them. Kinetic properties are consistent across different host cell lines and conserved across mouse, rat, and human. All electrophysiology data from all Kv channels are made available through a public website (Channelpedia). This dataset provides a solid foundation for exploring kinetics of heteromeric channels, roles of auxiliary subunits, kinetic modulation, and for building accurate Kv models.

A Derived Positional Mapping of Inhibitory Subtypes in the Somatosensory Cortex

Obtaining a catalog of cell types is a fundamental building block for understanding the brain. The ideal classification of cell-types is based on the profile of molecules expressed by a cell, in particular, the profile of genes expressed. One strategy is, therefore, to obtain as many single-cell transcriptomes as possible and isolate clusters of neurons with similar gene expression profiles. In this study, we explored an alternative strategy. We explored whether cell-types can be algorithmically derived by combining protein tissue stains with transcript expression profiles. We developed an algorithm that aims to distribute cell-types in the different layers of somatosensory cortex of the developing rat constrained by the tissue- and cellular level data. We found that the spatial distribution of major inhibitory cell types can be approximated using the available data. The result is a depth-wise atlas of inhibitory cell-types of the rat somatosensory cortex. In principle, any data that constrains what can occur in a particular part of the brain can also strongly constrain the derivation of cell-types. This draft inhibitory cell-type mapping is therefore dynamic and can iteratively converge towards the ground truth as further data is integrated.

NeuroMorphoVis: a collaborative framework for analysis and visualization of neuronal morphology skeletons reconstructed from microscopy stacks

Motivation

From image stacks to computational models, processing digital representations of neuronal morphologies is essential to neuroscientific research. Workflows involve various techniques and tools, leading in certain cases to convoluted and fragmented pipelines. The existence of an integrated, extensible and free framework for processing, analysis and visualization of those morphologies is a challenge that is still largely unfulfilled.

Results

We present NeuroMorphoVis, an interactive, extensible and cross-platform framework for building, visualizing and analyzing digital reconstructions of neuronal morphology skeletons extracted from microscopy stacks. Our framework is capable of detecting and repairing tracing artifacts, allowing the generation of high fidelity surface meshes and high resolution volumetric models for simulation and in silico imaging studies. The applicability of NeuroMorphoVis is demonstrated with two case studies. The first simulates the construction of three-dimensional profiles of neuronal somata and the other highlights how the framework is leveraged to create volumetric models of neuronal circuits for simulating different types of in vitro imaging experiments.

Availability and implementation

The source code and documentation are freely available on https://github.com/BlueBrain/NeuroMorphoVis under the GNU public license. The morphological analysis, visualization and surface meshing are implemented as an extensible Python API (Application Programming Interface) based on Blender, and the volume reconstruction and analysis code is written in C++ and parallelized using OpenMP. The framework features are accessible from a user-friendly GUI (Graphical User Interface) and a rich CLI (Command Line Interface).

Supplementary information

Supplementary data are available at Bioinformatics online.

A Brief History of Simulation Neuroscience

Our knowledge of the brain has evolved over millennia in philosophical, experimental and theoretical phases. We suggest that the next phase is simulation neuroscience. The main drivers of simulation neuroscience are big data generated at multiple levels of brain organization and the need to integrate these data to trace the causal chain of interactions within and across all these levels. Simulation neuroscience is currently the only methodology for systematically approaching the multiscale brain. In this review, we attempt to reconstruct the deep historical paths leading to simulation neuroscience, from the first observations of the nerve cell to modern efforts to digitally reconstruct and simulate the brain. Neuroscience began with the identification of the neuron as the fundamental unit of brain structure and function and has evolved towards understanding the role of each cell type in the brain, how brain cells are connected to each other, and how the seemingly infinite networks they form give rise to the vast diversity of brain functions. Neuronal mapping is evolving from subjective descriptions of cell types towards objective classes, subclasses and types. Connectivity mapping is evolving from loose topographic maps between brain regions towards dense anatomical and physiological maps of connections between individual genetically distinct neurons. Functional mapping is evolving from psychological and behavioral stereotypes towards a map of behaviors emerging from structural and functional connectomes. We show how industrialization of neuroscience and the resulting large disconnected datasets are generating demand for integrative neuroscience, how the scale of neuronal and connectivity maps is driving digital atlasing and digital reconstruction to piece together the multiple levels of brain organization, and how the complexity of the interactions between molecules, neurons, microcircuits and brain regions is driving brain simulation to understand the interactions in the multiscale brain.

Cellular, Synaptic and Network Effects of Acetylcholine in the Neocortex

The neocortex is densely innervated by basal forebrain (BF) cholinergic neurons. Long-range axons of cholinergic neurons regulate higher-order cognitive function and dysfunction in the neocortex by releasing acetylcholine (ACh). ACh release dynamically reconfigures neocortical microcircuitry through differential spatiotemporal actions on cell-types and their synaptic connections. At the cellular level, ACh release controls neuronal excitability and firing rate, by hyperpolarizing or depolarizing target neurons. At the synaptic level, ACh impacts transmission dynamics not only by altering the presynaptic probability of release, but also the magnitude of the postsynaptic response. Despite the crucial role of ACh release in physiology and pathophysiology, a comprehensive understanding of the way it regulates the activity of diverse neocortical cell-types and synaptic connections has remained elusive. This review aims to summarize the state-of-the-art anatomical and physiological data to develop a functional map of the cellular, synaptic and microcircuit effects of ACh in the neocortex of rodents and non-human primates, and to serve as a quantitative reference for those intending to build data-driven computational models on the role of ACh in governing brain states.

Individual differences in sensory sensitivity: Further lessons from an Autism model

In this commentary, we join Ward (this issue) in the usefulness of conceptualizing neural output in terms of signal and noise relationships, to create the missing links between neural, behavioral and subjective sensory sensitivity. We draw from our work in the Intense World Theory of Autism and the valproic acid rodent model, to complement the discussion with the consideration of developmental time and function of the system, for a neural output to serve as a predictor of atypical outcome in sensory sensitivity, and guide personalized therapies.

A Topological Representation of Branching Neuronal Morphologies

Many biological systems consist of branching structures that exhibit a wide variety of shapes. Our understanding of their systematic roles is hampered from the start by the lack of a fundamental means of standardizing the description of complex branching patterns, such as those of neuronal trees. To solve this problem, we have invented the Topological Morphology Descriptor (TMD), a method for encoding the spatial structure of any tree as a "barcode", a unique topological signature. As opposed to traditional morphometrics, the TMD couples the topology of the branches with their spatial extents by tracking their topological evolution in 3-dimensional space. We prove that neuronal trees, as well as stochastically generated trees, can be accurately categorized based on their TMD profiles. The TMD retains sufficient global and local information to create an unbiased benchmark test for their categorization and is able to quantify and characterize the structural differences between distinct morphological groups. The use of this mathematically rigorous method will advance our understanding of the anatomy and diversity of branching morphologies.

A Cell Atlas for the Mouse Brain

Despite vast numbers of studies of stained cells in the mouse brain, no current brain atlas provides region-by-region neuron counts. In fact, neuron numbers are only available for about 4% of brain of regions and estimates often vary by as much as 3-fold. Here we provide a first 3D cell atlas for the whole mouse brain, showing cell positions constructed algorithmically from whole brain Nissl and gene expression stains, and compared against values from the literature. The atlas provides the densities and positions of all excitatory and inhibitory neurons, astrocytes, oligodendrocytes, and microglia in each of the 737 brain regions defined in the AMBA. The atlas is dynamic, allowing comparison with previously reported numbers, addition of cell types, and improvement of estimates as new data is integrated. The atlas also provides insights into cellular organization only possible at this whole brain scale, and is publicly available.

Cell Densities in the Mouse Brain: A Systematic Review

The mouse brain is the most extensively studied brain of all species. We performed an exhaustive review of the literature to establish our current state of knowledge on cell numbers in mouse brain regions, arguably the most fundamental property to measure when attempting to understand a brain. The synthesized information, collected in one place, can be used by both theorists and experimentalists. Although for commonly-studied regions cell densities could be obtained for principal cell types, overall we know very little about how many cells are present in most brain regions and even less about cell-type specific densities. There is also substantial variation in cell density values obtained from different sources. This suggests that we need a new approach to obtain cell density datasets for the mouse brain.

Data-Driven Modeling of Cholinergic Modulation of Neural Microcircuits: Bridging Neurons, Synapses and Network Activity

Neuromodulators, such as acetylcholine (ACh), control information processing in neural microcircuits by regulating neuronal and synaptic physiology. Computational models and simulations enable predictions on the potential role of ACh in reconfiguring network activity. As a prelude into investigating how the cellular and synaptic effects of ACh collectively influence emergent network dynamics, we developed a data-driven framework incorporating phenomenological models of the physiology of cholinergic modulation of neocortical cells and synapses. The first-draft models were integrated into a biologically detailed tissue model of neocortical microcircuitry to investigate the effects of levels of ACh on diverse neuron types and synapses, and consequently on emergent network activity. Preliminary simulations from the framework, which was not tuned to reproduce any specific ACh-induced network effects, not only corroborate the long-standing notion that ACh desynchronizes spontaneous network activity, but also predict that a dose-dependent activation of ACh gives rise to a spectrum of neocortical network activity. We show that low levels of ACh, such as during non-rapid eye movement (nREM) sleep, drive microcircuit activity into slow oscillations and network synchrony, whereas high ACh concentrations, such as during wakefulness and REM sleep, govern fast oscillations and network asynchrony. In addition, spontaneous network activity modulated by ACh levels shape spike-time cross-correlations across distinct neuronal populations in strikingly different ways. These effects are likely due to the regulation of neurons and synapses caused by increasing levels of ACh, which enhances cellular excitability and decreases the efficacy of local synaptic transmission. We conclude by discussing future directions to refine the biological accuracy of the framework, which will extend its utility and foster the development of hypotheses to investigate the role of neuromodulators in neural information processing.

A Process for Digitizing and Simulating Biologically Realistic Oligocellular Networks Demonstrated for the Neuro-Glio-Vascular Ensemble

One will not understand the brain without an integrated exploration of structure and function, these attributes being two sides of the same coin: together they form the currency of biological computation. Accordingly, biologically realistic models require the re-creation of the architecture of the cellular components in which biochemical reactions are contained. We describe here a process of reconstructing a functional oligocellular assembly that is responsible for energy supply management in the brain and creating a computational model of the associated biochemical and biophysical processes. The reactions that underwrite thought are both constrained by and take advantage of brain morphologies pertaining to neurons, astrocytes and the blood vessels that deliver oxygen, glucose and other nutrients. Each component of this neuro-glio-vasculature ensemble (NGV) carries-out delegated tasks, as the dynamics of this system provide for each cell-type its own energy requirements while including mechanisms that allow cooperative energy transfers. Our process for recreating the ultrastructure of cellular components and modeling the reactions that describe energy flow uses an amalgam of state-of the-art techniques, including digital reconstructions of electron micrographs, advanced data analysis tools, computational simulations and in silico visualization software. While we demonstrate this process with the NGV, it is equally well adapted to any cellular system for integrating multimodal cellular data in a coherent framework.

The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow

Every neuron is part of a network, exerting its function by transforming multiple spatiotemporal synaptic input patterns into a single spiking output. This function is specified by the particular shape and passive electrical properties of the neuronal membrane, and the composition and spatial distribution of ion channels across its processes. For a variety of physiological or pathological reasons, the intrinsic input/output function may change during a neuron's lifetime. This process results in high variability in the peak specific conductance of ion channels in individual neurons. The mechanisms responsible for this variability are not well understood, although there are clear indications from experiments and modeling that degeneracy and correlation among multiple channels may be involved. Here, we studied this issue in biophysical models of hippocampal CA1 pyramidal neurons and interneurons. Using a unified data-driven simulation workflow and starting from a set of experimental recordings and morphological reconstructions obtained from rats, we built and analyzed several ensembles of morphologically and biophysically accurate single cell models with intrinsic electrophysiological properties consistent with experimental findings. The results suggest that the set of conductances expressed in any given hippocampal neuron may be considered as belonging to two groups: one subset is responsible for the major characteristics of the firing behavior in each population and the other is responsible for a robust degeneracy. Analysis of the model neurons suggests several experimentally testable predictions related to the combination and relative proportion of the different conductances that should be expressed on the membrane of different types of neurons for them to fulfill their role in the hippocampus circuitry.

Norepinephrine stimulates glycogenolysis in astrocytes to fuel neurons with lactate

The mechanism of rapid energy supply to the brain, especially to accommodate the heightened metabolic activity of excited states, is not well-understood. We explored the role of glycogen as a fuel source for neuromodulation using the noradrenergic stimulation of glia in a computational model of the neural-glial-vasculature ensemble (NGV). The detection of norepinephrine (NE) by the astrocyte and the coupled cAMP signal are rapid and largely insensitive to the distance of the locus coeruleus projection release sites from the glia, implying a diminished impact for volume transmission in high affinity receptor transduction systems. Glucosyl-conjugated units liberated from glial glycogen by NE-elicited cAMP second messenger transduction winds sequentially through the glycolytic cascade, generating robust increases in NADH and ATP before pyruvate is finally transformed into lactate. This astrocytic lactate is rapidly exported by monocarboxylate transporters to the associated neuron, demonstrating that the astrocyte-to-neuron lactate shuttle activated by glycogenolysis is a likely fuel source for neuromodulation and enhanced neural activity. Altogether, the energy supply for both astrocytes and neurons can be supplied rapidly by glycogenolysis upon neuromodulatory stimulus.

Although efficient compared to computers, the human brain utilizes energy at 10-fold the rate of other organs by mass. How the brain is supplied with sufficient on-demand energy to support its activity in the absence of neuronal storage capacity remains unknown. Neurons are not capable of meeting their own energy requirements, instead energy supply in the brain is managed by an oligocellular cartel composed of neurons, glia and the local vasculature (NGV), wherein glia can provide the ergogenic metabolite lactate to the neuron in a process called the astrocyte-to-neuron shuttle (ANLS). The only means of energy storage in the brain is glycogen, a polymerized form of glucose that is localized largely to astrocytes, but its exact role and conditions of use are not clear. In this computational model we show that neuromodulatory stimulation by norepinephrine induces astrocytes to recover glucosyl subunits from glycogen for use in a glycolytic process that favors the production of lactate. The ATP and NADH produced support metabolism in the astrocyte while the lactate is exported to feed the neuron. Thus, rapid energy demands by both neurons and glia in a stimulated brain can be met by glycogen mobilization.

Timed Synaptic Inhibition Shapes NMDA Spikes, Influencing Local Dendritic Processing and Global I/O Properties of Cortical Neurons

The NMDA spike is a long-lasting nonlinear phenomenon initiated locally in the dendritic branches of a variety of cortical neurons. It plays a key role in synaptic plasticity and in single-neuron computations. Combining dynamic system theory and computational approaches, we now explore how the timing of synaptic inhibition affects the NMDA spike and its associated membrane current. When impinging on its early phase, individual inhibitory synapses strongly, but transiently, dampen the NMDA spike; later inhibition prematurely terminates it. A single inhibitory synapse reduces the NMDA-mediated Ca²⁺ current, a key player in plasticity, by up to 45%. NMDA spikes in distal dendritic branches/spines are longer-lasting and more resilient to inhibition, enhancing synaptic plasticity at these branches. We conclude that NMDA spikes are highly sensitive to dendritic inhibition; sparse weak inhibition can finely tune synaptic plasticity both locally at the dendritic branch level and globally at the level of the neuron's output.

Morphological Diversity Strongly Constrains Synaptic Connectivity and Plasticity

Synaptic connectivity between neurons is naturally constrained by the anatomical overlap of neuronal arbors, the space on the axon available for synapses, and by physiological mechanisms that form synapses at a subset of potential synapse locations. What is not known is how these constraints impact emergent connectivity in a circuit with diverse morphologies. We investigated the role of morphological diversity within and across neuronal types on emergent connectivity in a model of neocortical microcircuitry. We found that the average overlap between the dendritic and axonal arbors of different types of neurons determines neuron-type specific patterns of distance-dependent connectivity, severely constraining the space of possible connectomes. However, higher order connectivity motifs depend on the diverse branching patterns of individual arbors of neurons belonging to the same type. Morphological diversity across neuronal types, therefore, imposes a specific structure on first order connectivity, and morphological diversity within neuronal types imposes a higher order structure of connectivity. We estimate that the morphological constraints resulting from diversity within and across neuron types together lead to a 10-fold reduction of the entropy of possible connectivity configurations, revealing an upper bound on the space explored by structural plasticity.

Cliques of Neurons Bound into Cavities Provide a Missing Link between Structure and Function

The lack of a formal link between neural network structure and its emergent function has hampered our understanding of how the brain processes information. We have now come closer to describing such a link by taking the direction of synaptic transmission into account, constructing graphs of a network that reflect the direction of information flow, and analyzing these directed graphs using algebraic topology. Applying this approach to a local network of neurons in the neocortex revealed a remarkably intricate and previously unseen topology of synaptic connectivity. The synaptic network contains an abundance of cliques of neurons bound into cavities that guide the emergence of correlated activity. In response to stimuli, correlated activity binds synaptically connected neurons into functional cliques and cavities that evolve in a stereotypical sequence toward peak complexity. We propose that the brain processes stimuli by forming increasingly complex functional cliques and cavities.