Brain cells can be classified according to their transcriptomic, morphological, and electrophysiological features. Comprehensive data on the spatial density of cell types that integrate all three properties are critical for improving models of how neuronal diversity contributes to brain function, yet existing atlases lack this information. To address this gap, we created a quantitative, three-dimensional atlas of cell type density distributions in the mouse brain. We began by generating a transcriptomic cell type atlas, scaling regional density estimates from brain slices using cell counts and anatomical dimensions. For densely populated regions like the cerebellum, we further refined these estimates by applying voxel-wise corrections based on average Nissl staining intensity. To connect transcriptomic identities with functional characteristics, we leveraged patch-sequencing datasets that combine single-neuron mRNA profiles, morphological reconstructions, and electrophysiological recordings from cortical neurons. Transcriptomic types were determined from gene expression data, morphological types were assigned based on structural reconstructions, and electrophysiological types were identified using K-means clustering. The resulting whole-brain atlas (consisting of 5274 transcriptomic clusters and 458 functional morphological-electrophysiological types) and computational tools offer a high-resolution (25 μm3 voxel size), integrative resource compatible with a broad range of neuroscience applications and enable the parsing of individual cell types to reveal previously unrecognized features.

Disruptions of energy supply to the brain are associated with many neurodegenerative pathologies and are difficult to study due to numerous interlinked metabolic pathways. We explored the effects of diminished energy supply on brain metabolism using a computational model of the neuro-glia-vasculature ensemble, in the form of a neuron, an astrocyte and local blood supply. As a case study, we investigated the glucose transporter type-1 deficiency syndrome (GLUT1-DS), a childhood affliction characterized by impaired glucose utilization and associated with phenotypes including seizures. Compared to neurons, astrocytes exhibited markedly higher metabolite concentration variabilities for all but a few redox species. This effect could signal a role for astrocytes in absorbing the shock of blood nutrient fluctuations. Redox balances were disrupted in GLUT1-DS with lower levels of reducing equivalent carriers NADH and ATP. The best non-glucose nutrient or pharmacotherapies for re-establishing redox normalcy involved lactate, the keto-diet (β-hydroxybutyrate), NAD and Q10 supplementation, suggesting a possible glucose sparing mechanism. GLUT1-DS seizures resulted from after-discharge neuronal firing caused by post-stimulus ATP reductions and impaired Na+/K+-ATPase, which can be rescued by restoring either normal glucose or by relatively small increases in neuronal ATP.

The high energetic demands of the brain arise primarily from neuronal activity. Neurons consume substantial energy to transmit information as electrical signals and maintain their resting membrane potential. These energetic requirements are met by the neuro-glial-vascular (NGV) ensemble, which generates energy in a coupled metabolic process. In ageing, metabolic function becomes impaired, producing less energy and, consequently, the system is unable to sustain the neuronal energetic needs. We propose a multiscale model of electro-metabolic coupling in a reconstructed rat neocortex. This combines an electro-morphologically reconstructed electrophysiological model with a detailed NGV metabolic model. Our results demonstrate that the large-scale model effectively captures electro-metabolic processes at the circuit level, highlighting the importance of heterogeneity within the circuit, where energetic demands vary according to neuronal characteristics. Finally, in metabolic ageing, our model indicates that the middle cortical layers are particularly vulnerable to energy impairment.

Understanding the intracellular dynamics of brain cells entails performing three-dimensional molecular simulations incorporating ultrastructural models that can capture cellular membrane geometries at nanometer scales. While there is an abundance of neuronal morphologies available online, e.g. from NeuroMorpho.Org, converting those fairly abstract point-and-diameter representations into geometrically realistic and simulation-ready, i.e. watertight, manifolds is challenging. Many neuronal mesh reconstruction methods have been proposed; however, their resulting meshes are either biologically unplausible or non-watertight. We present an effective and unconditionally robust method capable of generating geometrically realistic and watertight surface manifolds of spiny cortical neurons from their morphological descriptions. The robustness of our method is assessed based on a mixed dataset of cortical neurons with a wide variety of morphological classes. The implementation is seamlessly extended and applied to synthetic astrocytic morphologies that are also plausibly biological in detail. Resulting meshes are ultimately used to create volumetric meshes with tetrahedral domains to perform scalable in silico reaction-diffusion simulations for revealing cellular structure-function relationships. Availability and implementation: Our method is implemented in NeuroMorphoVis, a neuroscience-specific open source Blender add-on, making it freely accessible for neuroscience researchers.

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.

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.

The mouse brain contains a rich diversity of inhibitory neuron types that have been characterized by their patterns of gene expression. However, it is still unclear how these cell types are distributed across the mouse brain. We developed a computational method to estimate the densities of different inhibitory neuron types across the mouse brain. Our method allows the unbiased integration of diverse and disparate datasets into one framework to predict inhibitory neuron densities for uncharted brain regions. We constrained our estimates based on previously computed brain-wide neuron densities, gene expression data from in situ hybridization image stacks together with a wide range of values reported in the literature. Using constrained optimization, we derived coherent estimates of cell densities for the different inhibitory neuron types. We estimate that 20.3% of all neurons in the mouse brain are inhibitory. Among all inhibitory neurons, 18% predominantly express parvalbumin (PV), 16% express somatostatin (SST), 3% express vasoactive intestinal peptide (VIP), and the remainder 63% belong to the residual GABAergic population. We find that our density estimations improve as more literature values are integrated. Our pipeline is extensible, allowing new cell types or data to be integrated as they become available. The data, algorithms, software, and results of our pipeline are publicly available and update the Blue Brain Cell Atlas. This work therefore leverages the research community to collectively converge on the numbers of each cell type in each brain region.

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.

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.

Recorded potentials in the extracellular space (ECS) of the brain is a standard measure of population activity in neural tissue. Computational models that simulate the relationship between the ECS potential and its underlying neurophysiological processes are commonly used in the interpretation of such measurements. Standard methods, such as volume-conductor theory and current-source density theory, assume that diffusion has a negligible effect on the ECS potential, at least in the range of frequencies picked up by most recording systems. This assumption remains to be verified. We here present a hybrid simulation framework that accounts for diffusive effects on the ECS potential. The framework uses (1) the NEURON simulator to compute the activity and ionic output currents from multicompartmental neuron models, and (2) the electrodiffusive Kirchhoff-Nernst-Planck framework to simulate the resulting dynamics of the potential and ion concentrations in the ECS, accounting for the effect of electrical migration as well as diffusion. Using this framework, we explore the effect that ECS diffusion has on the electrical potential surrounding a small population of 10 pyramidal neurons. The neural model was tuned so that simulations over ∼100 seconds of biological time led to shifts in ECS concentrations by a few millimolars, similar to what has been seen in experiments. By comparing simulations where ECS diffusion was absent with simulations where ECS diffusion was included, we made the following key findings: (i) ECS diffusion shifted the local potential by up to ∼0.2 mV. (ii) The power spectral density (PSD) of the diffusion-evoked potential shifts followed a 1/f2 power law. (iii) Diffusion effects dominated the PSD of the ECS potential for frequencies up to several hertz. In scenarios with large, but physiologically realistic ECS concentration gradients, diffusion was thus found to affect the ECS potential well within the frequency range picked up in experimental recordings.

The spatial arrangement of Ca2+ channels and vesicles remains unknown for most CNS synapses, despite of the crucial importance of this geometrical parameter for the Ca2+ control of transmitter release. At a large model synapse, the calyx of Held, transmitter release is controlled by several Ca2+ channels in a "domain overlap" mode, at least in young animals. To study the geometrical constraints of Ca2+ channel placement in domain overlap control of release, we used stochastic MCell modelling, at active zones for which the position of docked vesicles was derived from electron microscopy (EM). We found that random placement of Ca2+ channels was unable to produce high slope values between release and presynaptic Ca2+ entry, a hallmark of domain overlap, and yielded excessively large release probabilities. The simple assumption that Ca2+ channels can be located anywhere at active zones, except below a critical distance of ~ 30 nm away from docked vesicles ("exclusion zone"), rescued high slope values and low release probabilities. Alternatively, high slope values can also be obtained by placing all Ca2+ channels into a single supercluster, which however results in significantly higher heterogeneity of release probabilities. We also show experimentally that high slope values, and the sensitivity to the slow Ca2+ chelator EGTA-AM, are maintained with developmental maturation of the calyx synapse. Taken together, domain overlap control of release represents a highly organized active zone architecture in which Ca2+ channels must obey a certain distance to docked vesicles. Furthermore, domain overlap can be employed by near-mature, fast-releasing synapses.