Computational synthesis of cortical dendritic morphologies

Non-homogenous axonal bouton distribution in whole-brain single-cell neuronal networks

We examined the distribution of pre-synaptic contacts in axons of mouse neurons and constructed whole-brain single-cell neuronal networks using an extensive dataset of 1,891 fully reconstructed neurons. We found that bouton locations were not homogeneous throughout the axon and among brain regions. As our algorithm was able to generate whole-brain single-cell connectivity matrices from full morphology reconstruction datasets, we further found that non-homogeneous bouton locations have a significant impact on network wiring, including degree distribution, triad census, and community structure. By perturbing neuronal morphology, we further explored the link between anatomical details and network topology. In our in silico exploration, we found that dendritic and axonal tree span would have the greatest impact on network wiring, followed by synaptic contact deletion. Our results suggest that neuroanatomical details must be carefully addressed in studies of whole-brain networks at the single-cell level.

Mathematical generation of data-driven hippocampal CA1 pyramidal neurons and interneurons copies via A-GLIF models for large-scale networks covering the experimental variability range

Efficient and accurate large-scale networks are a fundamental tool in modeling brain areas, to advance our understanding of neuronal dynamics. However, their implementation faces two key issues: computational efficiency and heterogeneity. Computational efficiency is achieved using simplified neurons, whereas there are no practical solutions available to solve the problem of reproducing in a large-scale network the experimentally observed heterogeneity of the intrinsic properties of neurons. This is important, because the use of identical nodes in a network can generate artifacts which can hinder an adequate representation of the properties of a real network. To this aim, we introduce a mathematical procedure to generate an arbitrary large number of copies of simplified hippocampal CA1 pyramidal neurons and interneurons models, which exhibit the full range of firing dynamics observed in these cells - including adapting, non-adapting and bursting. For this purpose, we rely on a recently published adaptive generalized leaky integrate-and-fire (A-GLIF) modeling approach, leveraging on its ability to reproduce the rich set of electrophysiological behaviors of these types of neurons under a variety of different stimulation currents. The generation procedure is based on a perturbation of model's parameters related to the initial data, firing block, and internal dynamics, and suitably validated against experimental data to ensure that the firing dynamics of any given cell copy remains within the experimental range. A classification procedure confirmed that the firing behavior of most of the pyramidal/interneuron copies was consistent with the experimental data. This approach allows to obtain heterogeneous copies with mathematically controlled firing properties. A full set of heterogeneous neurons composing the CA1 region of a rat hippocampus (approximately 1.2 million neurons), are provided in a database freely available in the live paper section of the EBRAINS platform. By adapting the underlying A-GLIF framework, it will be possible to extend the numerical approach presented here to create, in a mathematically controlled manner, an arbitrarily large number of non-identical copies of cell populations with firing properties related to other brain areas.

Accelerating the continuous community sharing of digital neuromorphology data

The tree-like morphology of neurons and glia is a key cellular determinant of circuit connectivity and metabolic function in the nervous system of essentially all animals. To elucidate the contribution of specific cell types to both physiological and pathological brain states, it is important to access detailed neuroanatomy data for quantitative analysis and computational modeling. NeuroMorpho.Org is the largest online collection of freely available digital neural reconstructions and related metadata and is continuously updated with new uploads. Earlier in the project, we released multiple datasets together yearly, but this process caused an average delay of several months in making the data public. Moreover, in the past 5 years, >80% of invited authors agreed to share their data with the community via NeuroMorpho.Org, up from <20% in the first 5 years of the project. In the same period, the average number of reconstructions per publication increased 600%, creating the need for automatic processing to release more reconstructions in less time. The progressive automation of our pipeline enabled the transition to agile releases of individual datasets as soon as they are ready. The overall time from data identification to public sharing decreased by 63.7%; 78% of the datasets are now released in less than 3 months with an average workflow duration below 40 days. Furthermore, the mean processing time per reconstruction dropped from 3 hours to 2 minutes. With these continuous improvements, NeuroMorpho.Org strives to forge a positive culture of open data. Most importantly, the new, original research enabled through reuse of datasets across the world has a multiplicative effect on science discovery, benefiting both authors and users.

Efficient, continual, and generalized learning in the brain - neural mechanism of Mental Schema 2.0

There has been tremendous progress in artificial neural networks (ANNs) over the past decade; however, the gap between ANNs and the biological brain as a learning device remains large. With the goal of closing this gap, this paper reviews learning mechanisms in the brain by focusing on three important issues in ANN research: efficiency, continuity, and generalization. We first discuss the method by which the brain utilizes a variety of self-organizing mechanisms to maximize learning efficiency, with a focus on the role of spontaneous activity of the brain in shaping synaptic connections to facilitate spatiotemporal learning and numerical processing. Then, we examined the neuronal mechanisms that enable lifelong continual learning, with a focus on memory replay during sleep and its implementation in brain-inspired ANNs. Finally, we explored the method by which the brain generalizes learned knowledge in new situations, particularly from the mathematical generalization perspective of topology. Besides a systematic comparison in learning mechanisms between the brain and ANNs, we propose "Mental Schema 2.0," a new computational property underlying the brain's unique learning ability that can be implemented in ANNs.

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.

Grey and white matter metrics demonstrate distinct and complementary prediction of differences in cognitive performance in children: Findings from ABCD (N= 11 876)

Individual differences in cognitive performance in childhood are a key predictor of significant life outcomes such as educational attainment and mental health. Differences in cognitive ability are governed in part by variations in brain structure. However, studies commonly focus on either grey or white matter metrics in humans, leaving open the key question as to whether grey or white matter microstructure play distinct or complementary roles supporting cognitive performance. To compare the role of grey and white matter in supporting cognitive performance, we used regularized structural equation models to predict cognitive performance with grey and white matter measures. Specifically, we compared how grey matter (volume, cortical thickness and surface area) and white matter measures (volume, fractional anisotropy and mean diffusivity) predicted individual differences in cognitive performance. The models were tested in 11,876 children (ABCD Study, 5680 female; 6196 male) at 10 years old. We found that grey and white matter metrics bring partly non-overlapping information to predict cognitive performance. The models with only grey or white matter explained respectively 15.4% and 12.4% of the variance in cognitive performance, while the combined model explained 19.0%. Zooming in we additionally found that different metrics within grey and white matter had different predictive power, and that the tracts/regions that were most predictive of cognitive performance differed across metric. These results show that studies focusing on a single metric in either grey or white matter to study the link between brain structure and cognitive performance are missing a key part of the equation.

Efficient simulation of neural development using shared memory parallelization

The Neural Development Simulator, NeuroDevSim, is a Python module that simulates the most important aspects of brain development: morphological growth, migration, and pruning. It uses an agent-based modeling approach inherited from the NeuroMaC software. Each cycle has agents called fronts execute model-specific code. In the case of a growing dendritic or axonal front, this will be a choice between extension, branching, or growth termination. Somatic fronts can migrate to new positions and any front can be retracted to prune parts of neurons. Collision detection prevents new or migrating fronts from overlapping with existing ones. NeuroDevSim is a multi-core program that uses an innovative shared memory approach to achieve parallel processing without messaging. We demonstrate linear strong parallel scaling up to 96 cores for large models and have run these successfully on 128 cores. Most of the shared memory parallelism is achieved without memory locking. Instead, cores have only write privileges to private sections of arrays, while being able to read the entire shared array. Memory conflicts are avoided by a coding rule that allows only active fronts to use methods that need writing access. The exception is collision detection, which is needed to avoid the growth of physically overlapping structures. For collision detection, a memory-locking mechanism was necessary to control access to grid points that register the location of nearby fronts. A custom approach using a serialized lock broker was able to manage both read and write locking. NeuroDevSim allows easy modeling of most aspects of neural development for models simulating a few complex or thousands of simple neurons or a mixture of both.

Code available at

https://github.com/CNS-OIST/NeuroDevSim.

Models of Purkinje cell dendritic tree selection during early cerebellar development

We investigate the relationship between primary dendrite selection of Purkinje cells and migration of their presynaptic partner granule cells during early cerebellar development. During postnatal development, each Purkinje cell grows more than three dendritic trees, from which a primary tree is selected for development, whereas the others completely retract. Experimental studies suggest that this selection process is coordinated by physical and synaptic interactions with granule cells, which undergo a massive migration at the same time. However, technical limitations hinder continuous experimental observation of multiple cell populations. To explore possible mechanisms underlying this selection process, we constructed a computational model using a new computational framework, NeuroDevSim. The study presents the first computational model that simultaneously simulates Purkinje cell growth and the dynamics of granule cell migrations during the first two postnatal weeks, allowing exploration of the role of physical and synaptic interactions upon dendritic selection. The model suggests that interaction with parallel fibers is important to establish the distinct planar morphology of Purkinje cell dendrites. Specific rules to select which dendritic trees to keep or retract result in larger winner trees with more synaptic contacts than using random selection. A rule based on afferent synaptic activity was less effective than rules based on dendritic size or numbers of synapses.

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.