Reconstruction and Simulation of Neocortical Microcircuitry

The role of calcium and CaMKII in sleep

Sleep is an evolutionarily conserved phenotype shared by most of the animals on the planet. Prolonged wakefulness will result in increased sleep need or sleep pressure. However, its mechanisms remain elusive. Recent findings indicate that Ca²⁺ signaling, known to control diverse physiological functions, also regulates sleep. This review intends to summarize research advances in Ca²⁺ and Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) in sleep regulation. Significant changes in sleep phenotype have been observed through calcium-related channels, receptors, and pumps. Mathematical modeling for neuronal firing patterns during NREM sleep suggests that these molecules compose a Ca²⁺-dependent hyperpolarization mechanism. The intracellular Ca²⁺ may then trigger sleep induction and maintenance through the activation of CaMKII, one of the sleep-promoting kinases. CaMKII and its multisite phosphorylation status may provide a link between transient calcium dynamics typically observed in neurons and sleep-wake dynamics observed on the long-time scale.

A method to estimate the cellular composition of the mouse brain from heterogeneous datasets

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.

A functional neuron maturation device provides convenient application on microelectrode array for neural network measurement

BACKGROUND

Microelectrode array (MEA) systems are valuable for in vitro assessment of neurotoxicity and drug efficiency. However, several difficulties such as protracted functional maturation and high experimental costs hinder the use of MEA analysis requiring human induced pluripotent stem cells (hiPSCs). Neural network functional parameters are also needed for in vitro to in vivo extrapolation.

METHODS

In the present study, we produced a cost effective nanofiber culture platform, the SCAD device, for long-term culture of hiPSC-derived neurons and primary peripheral neurons. The notable advantage of SCAD device is convenient application on multiple MEA systems for neuron functional analysis.

RESULTS

We showed that the SCAD device could promote functional maturation of cultured hiPSC-derived neurons, and neurons responded appropriately to convulsant agents. Furthermore, we successfully analyzed parameters for in vitro to in vivo extrapolation, i.e., low-frequency components and synaptic propagation velocity of the signal, potentially reflecting neural network functions from neurons cultured on SCAD device. Finally, we measured the axonal conduction velocity of peripheral neurons.

CONCLUSIONS

Neurons cultured on SCAD devices might constitute a reliable in vitro platform to investigate neuron functions, drug efficacy and toxicity, and neuropathological mechanisms by MEA.

Closed-loop direct control of seizure focus in a rodent model of temporal lobe epilepsy via localized electric fields applied sequentially

Direct electrical stimulation of the seizure focus can achieve the early termination of epileptic oscillations. However, direct intervention of the hippocampus, the most prevalent seizure focus in temporal lobe epilepsy is thought to be not practicable due to its large size and elongated shape. Here, in a rat model, we report a sequential narrow-field stimulation method for terminating seizures, while focusing stimulus energy at the spatially extensive hippocampal structure. The effects and regional specificity of this method were demonstrated via electrophysiological and biological responses. Our proposed modality demonstrates spatiotemporal preciseness and selectiveness for modulating the pathological target region which may have potential for further investigation as a therapeutic approach.

Enhancement of calculation accuracy of the integrated photonic tensor flow processer by global optical power allocation

We present a global optical power allocation architecture, which can enhance the calculation accuracy of the integrated photonic tensor flow processor (PTFP). By adjusting the optical power splitting ratio according to the weight value and loss of each calculating unit, this architecture can efficiently use optical power so that the signal-to-noise ratio of the PTFP is enhanced. In the case of considering the on-chip optical delay line and spectral loss, the calculation accuracy measured in the experiment is enhanced by more than 1 bit compared with the fixed optical power allocation architecture.

Long-Term Synaptic Plasticity Tunes the Gain of Information Channels through the Cerebellum Granular Layer

A central hypothesis on brain functioning is that long-term potentiation (LTP) and depression (LTD) regulate the signals transfer function by modifying the efficacy of synaptic transmission. In the cerebellum, granule cells have been shown to control the gain of signals transmitted through the mossy fiber pathway by exploiting synaptic inhibition in the glomeruli. However, the way LTP and LTD control signal transformation at the single-cell level in the space, time and frequency domains remains unclear. Here, the impact of LTP and LTD on incoming activity patterns was analyzed by combining patch-clamp recordings in acute cerebellar slices and mathematical modeling. LTP reduced the delay, increased the gain and broadened the frequency bandwidth of mossy fiber burst transmission, while LTD caused opposite changes. These properties, by exploiting NMDA subthreshold integration, emerged from microscopic changes in spike generation in individual granule cells such that LTP anticipated the emission of spikes and increased their number and precision, while LTD sorted the opposite effects. Thus, akin with the expansion recoding process theoretically attributed to the cerebellum granular layer, LTP and LTD could implement selective filtering lines channeling information toward the molecular and Purkinje cell layers for further processing.

A faster way to model neuronal circuitry

Artificial neural networks could pave the way for efficiently simulating large-scale models of neuronal networks in the nervous system.

Quantitative spatial analysis reveals that the local axons of lamina I projection neurons and interneurons exhibit distributions that predict distinct roles in spinal sensory processing

Our knowledge about the detailed wiring of neuronal circuits in the spinal dorsal horn (DH), where initial sensory processing takes place, is still very sparse. While a substantial amount of data is available on the somatodendritic morphology of DH neurons, the laminar and segmental distribution patterns and consequential function of individual axons are much less characterized. In the present study, we fully reconstructed the axonal and dendritic processes of 10 projection neurons (PNs) and 15 interneurons (INs) in lamina I of the rat, to reveal quantitative differences in their distribution. We also performed whole-cell patch-clamp recordings to test the predicted function of certain axon collaterals. In line with our earlier qualitative description, we found that lamina I INs in the lateral aspect of the superficial DH send axon collaterals toward the medial part and occupy mostly laminae I-III, providing anatomical basis for a lateromedial flow of information within the DH. Local axon collaterals of PNs were more extensively distributed including dorsal commissural axon collaterals that might refer to those reported earlier linking the lateral aspect of the left and right DHs. PN collaterals dominated the dorsolateral funiculus and laminae IV-VI, suggesting propriospinal and ventral connections. Indeed, patch-clamp recordings confirmed the existence of a dorsoventral excitatory drive upon activation of neurokinin-1 receptors that, although being expressed in various lamina I neurons, are specifically enriched in PNs. In summary, lamina I PNs and INs have almost identical dendritic input fields, while their segmental axon collateral distribution patterns are distinct. INs, whose somata reside in lamina I, establish local connections, may show asymmetry, and contribute to bridging the medial and lateral halves of the DH. PNs, on the other hand, preferably relay their integrated dendritic input to deeper laminae of the spinal gray matter where it might be linked to other ascending pathways or the premotor network, resulting in a putative direct contribution to the nociceptive withdrawal reflex.

Quasicriticality explains variability of human neural dynamics across life span

Aging impacts the brain's structural and functional organization and over time leads to various disorders, such as Alzheimer's disease and cognitive impairment. The process also impacts sensory function, bringing about a general slowing in various perceptual and cognitive functions. Here, we analyze the Cambridge Centre for Ageing and Neuroscience (Cam-CAN) resting-state magnetoencephalography (MEG) dataset-the largest aging cohort available-in light of the quasicriticality framework, a novel organizing principle for brain functionality which relates information processing and scaling properties of brain activity to brain connectivity and stimulus. Examination of the data using this framework reveals interesting correlations with age and gender of test subjects. Using simulated data as verification, our results suggest a link between changes to brain connectivity due to aging and increased dynamical fluctuations of neuronal firing rates. Our findings suggest a platform to develop biomarkers of neurological health.

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.

Sims and Vulnerability: On the Ethics of Creating Emulated Minds

It might become possible to build artificial minds with the capacity for experience. This raises a plethora of ethical issues, explored, among others, in the context of whole brain emulations (WBE). In this paper, I will take up the problem of vulnerability - given, for various reasons, less attention in the literature - that the conscious emulations will likely exhibit. Specifically, I will examine the role that vulnerability plays in generating ethical issues that may arise when dealing with WBEs. I will argue that concerns about vulnerability are more matters of institutional design than individual ethics, both when it comes to creating humanlike brain emulations, and when animal-like emulations are concerned. Consequently, the article contains reflection on some institutional measures that can be taken to protect the sims' interests. It concludes that an institutional framework more likely to succeed in this task is competitive and poly-centric, rather than monopolistic and centralized.

Multiscale imaging informs translational mouse modeling of neurological disease

Multiscale neurophysiology reveals that simple motor actions are associated with changes in neuronal firing in virtually every brain region studied. Accordingly, the assessment of focal pathology such as stroke or progressive neurodegenerative diseases must also extend widely across brain areas. To derive mechanistic information through imaging, multiple resolution scales and multimodal factors must be included, such as the structure and function of specific neurons and glial cells and the dynamics of specific neurotransmitters. Emerging multiscale methods in preclinical animal studies that span micro- to macroscale examinations fill this gap, allowing a circuit-based understanding of pathophysiological mechanisms. Combined with high-performance computation and open-source data repositories, these emerging multiscale and large field-of-view techniques include live functional ultrasound, multi- and single-photon wide-scale light microscopy, video-based miniscopes, and tissue-penetrating fiber photometry, as well as variants of post-mortem expansion microscopy. We present these technologies and outline use cases and data pipelines to uncover new knowledge within animal models of stroke, Alzheimer's disease, and movement disorders.

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

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

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

BACKGROUND

Whole-cell patch-clamp electrophysiology is an essential technique for understanding how single neurons translate their diverse inputs into a functional output. The relative inaccessibility of live human cortical neurons for experimental manipulation has made it difficult to determine the unique features of how human cortical neurons differ from their counterparts in other species.

FINDINGS

We present a curated repository of whole-cell patch-clamp recordings from surgically resected human cortical tissue, encompassing 118 neurons from 35 individuals (age range, 21-59 years; 17 male, 18 female). Recorded human cortical neurons derive from layers 2 and 3 (L2&3), deep layer 3 (L3c), or layer 5 (L5) and are annotated with a rich set of subject and experimental metadata. For comparison, we also provide a limited set of comparable recordings from 21-day-old mice (11 cells from 5 mice). All electrophysiological recordings are provided in the Neurodata Without Borders (NWB) format and are available for further analysis via the Distributed Archives for Neurophysiology Data Integration online repository. The associated data conversion code is made publicly available and can help others in converting electrophysiology datasets to the open NWB standard for general reuse.

CONCLUSION

These data can be used for novel analyses of biophysical characteristics of human cortical neurons, including in cross-species or cross-lab comparisons or in building computational models of individual human neurons.

Model simulations unveil the structure-function-dynamics relationship of the cerebellar cortical microcircuit

The cerebellar network is renowned for its regular architecture that has inspired foundational computational theories. However, the relationship between circuit structure, function and dynamics remains elusive. To tackle the issue, we developed an advanced computational modeling framework that allows us to reconstruct and simulate the structure and function of the mouse cerebellar cortex using morphologically realistic multi-compartmental neuron models. The cerebellar connectome is generated through appropriate connection rules, unifying a collection of scattered experimental data into a coherent construct and providing a new model-based ground-truth about circuit organization. Naturalistic background and sensory-burst stimulation are used for functional validation against recordings in vivo, monitoring the impact of cellular mechanisms on signal propagation, inhibitory control, and long-term synaptic plasticity. Our simulations show how mossy fibers entrain the local neuronal microcircuit, boosting the formation of columns of activity travelling from the granular to the molecular layer providing a new resource for the investigation of local microcircuit computation and of the neural correlates of behavior.

Ultrafast simulation of large-scale neocortical microcircuitry with biophysically realistic neurons

Understanding the activity of the mammalian brain requires an integrative knowledge of circuits at distinct scales, ranging from ion channel gating to circuit connectomics. Computational models are regularly employed to understand how multiple parameters contribute synergistically to circuit behavior. However, traditional models of anatomically and biophysically realistic neurons are computationally demanding, especially when scaled to model local circuits. To overcome this limitation, we trained several artificial neural network (ANN) architectures to model the activity of realistic multicompartmental cortical neurons. We identified an ANN architecture that accurately predicted subthreshold activity and action potential firing. The ANN could correctly generalize to previously unobserved synaptic input, including in models containing nonlinear dendritic properties. When scaled, processing times were orders of magnitude faster compared with traditional approaches, allowing for rapid parameter-space mapping in a circuit model of Rett syndrome. Thus, we present a novel ANN approach allowing for rapid, detailed network experiments using inexpensive and commonly available computational resources.

Computational models of neurotransmission at cerebellar synapses unveil the impact on network computation

The neuroscientific field benefits from the conjoint evolution of experimental and computational techniques, allowing for the reconstruction and simulation of complex models of neurons and synapses. Chemical synapses are characterized by presynaptic vesicle cycling, neurotransmitter diffusion, and postsynaptic receptor activation, which eventually lead to postsynaptic currents and subsequent membrane potential changes. These mechanisms have been accurately modeled for different synapses and receptor types (AMPA, NMDA, and GABA) of the cerebellar cortical network, allowing simulation of their impact on computation. Of special relevance is short-term synaptic plasticity, which generates spatiotemporal filtering in local microcircuits and controls burst transmission and information flow through the network. Here, we present how data-driven computational models recapitulate the properties of neurotransmission at cerebellar synapses. The simulation of microcircuit models is starting to reveal how diverse synaptic mechanisms shape the spatiotemporal profiles of circuit activity and computation.

Impact of somatostatin interneurons on interactions between barrels in plasticity induced by whisker deprivation

The activity of inhibitory interneurons has a profound role in shaping cortical plasticity. Somatostatin-expressing interneurons (SOM-INs) are involved in several aspects of experience-dependent cortical rewiring. We addressed the question of the barrel cortex SOM-IN engagement in plasticity formation induced by sensory deprivation in adult mice (2-3 months old). We used a spared vibrissa paradigm, resulting in a massive sensory map reorganization. Using chemogenetic manipulation, the activity of barrel cortex SOM-INs was blocked or activated by continuous clozapine N-oxide (CNO) administration during one-week-long deprivation. To visualize the deprivation-induced plasticity, [¹⁴C]-2-deoxyglucose mapping of cortical functional representation of the spared whisker was performed at the end of the deprivation. The plasticity was manifested as an extension of cortical activation in response to spared vibrissae stimulation. We found that SOM-IN inhibition in the cortical column of the spared whisker did not influence the areal extent of the cortex activated by the spared whisker. However, blocking the activity of SOM-INs in the deprived column, adjacent to the spared one, decreased the plasticity of the spared whisker representation. SOM-IN activation did not affect plasticity. These data show that SOM-IN activity is part of cortical circuitry that affects interbarrel interactions underlying deprivation-induced plasticity in adult mice.

In vivo extracellular recordings of thalamic and cortical visual responses reveal V1 connectivity rules

The brain's connectome provides the scaffold for canonical neural computations. However, a comparison of connectivity studies in the mouse primary visual cortex (V1) reveals that the average number and strength of connections between specific neuron types can vary. Can variability in V1 connectivity measurements coexist with canonical neural computations? We developed a theory-driven approach to deduce V1 network connectivity from visual responses in mouse V1 and visual thalamus (dLGN). Our method revealed that the same recorded visual responses were captured by multiple connectivity configurations. Remarkably, the magnitude and selectivity of connectivity weights followed a specific order across most of the inferred connectivity configurations. We argue that this order stems from the specific shapes of the recorded contrast response functions and contrast invariance of orientation tuning. Remarkably, despite variability across connectivity studies, connectivity weights computed from individual published connectivity reports followed the order we identified with our method, suggesting that the relations between the weights, rather than their magnitudes, represent a connectivity motif supporting canonical V1 computations.

The present and future of neural interfaces

The 2020's decade will likely witness an unprecedented development and deployment of neurotechnologies for human rehabilitation, personalized use, and cognitive or other enhancement. New materials and algorithms are already enabling active brain monitoring and are allowing the development of biohybrid and neuromorphic systems that can adapt to the brain. Novel brain-computer interfaces (BCIs) have been proposed to tackle a variety of enhancement and therapeutic challenges, from improving decision-making to modulating mood disorders. While these BCIs have generally been developed in an open-loop modality to optimize their internal neural decoders, this decade will increasingly witness their validation in closed-loop systems that are able to continuously adapt to the user's mental states. Therefore, a proactive ethical approach is needed to ensure that these new technological developments go hand in hand with the development of a sound ethical framework. In this perspective article, we summarize recent developments in neural interfaces, ranging from neurohybrid synapses to closed-loop BCIs, and thereby identify the most promising macro-trends in BCI research, such as simulating vs. interfacing the brain, brain recording vs. brain stimulation, and hardware vs. software technology. Particular attention is devoted to central nervous system interfaces, especially those with application in healthcare and human enhancement. Finally, we critically assess the possible futures of neural interfacing and analyze the short- and long-term implications of such neurotechnologies.

Characteristic columnar connectivity caters to cortical computation: Replication, simulation, and evaluation of a microcircuit model

The neocortex, and with it the mammalian brain, achieves a level of computational efficiency like no other existing computational engine. A deeper understanding of its building blocks (cortical microcircuits), and their underlying computational principles is thus of paramount interest. To this end, we need reproducible computational models that can be analyzed, modified, extended and quantitatively compared. In this study, we further that aim by providing a replication of a seminal cortical column model. This model consists of noisy Hodgkin-Huxley neurons connected by dynamic synapses, whose connectivity scheme is based on empirical findings from intracellular recordings. Our analysis confirms the key original finding that the specific, data-based connectivity structure enhances the computational performance compared to a variety of alternatively structured control circuits. For this comparison, we use tasks based on spike patterns and rates that require the systems not only to have simple classification capabilities, but also to retain information over time and to be able to compute nonlinear functions. Going beyond the scope of the original study, we demonstrate that this finding is independent of the complexity of the neuron model, which further strengthens the argument that it is the connectivity which is crucial. Finally, a detailed analysis of the memory capabilities of the circuits reveals a stereotypical memory profile common across all circuit variants. Notably, the circuit with laminar structure does not retain stimulus any longer than any other circuit type. We therefore conclude that the model's computational advantage lies in a sharper representation of the stimuli.

The EBRAINS Hodgkin-Huxley Neuron Builder: An online resource for building data-driven neuron models

In the last decades, brain modeling has been established as a fundamental tool for understanding neural mechanisms and information processing in individual cells and circuits at different scales of observation. Building data-driven brain models requires the availability of experimental data and analysis tools as well as neural simulation environments and, often, large scale computing facilities. All these components are rarely found in a comprehensive framework and usually require ad hoc programming. To address this, we developed the EBRAINS Hodgkin-Huxley Neuron Builder (HHNB), a web resource for building single cell neural models via the extraction of activity features from electrophysiological traces, the optimization of the model parameters via a genetic algorithm executed on high performance computing facilities and the simulation of the optimized model in an interactive framework. Thanks to its inherent characteristics, the HHNB facilitates the data-driven model building workflow and its reproducibility, hence fostering a collaborative approach to brain modeling.

Objective Supervised Machine Learning-Based Classification and Inference of Biological Neuronal Networks

The classification of biological neuron types and networks poses challenges to the full understanding of the human brain’s organisation and functioning. In this paper, we develop a novel objective classification model of biological neuronal morphology and electrical types and their networks, based on the attributes of neuronal communication using supervised machine learning solutions. This presents advantages compared to the existing approaches in neuroinformatics since the data related to mutual information or delay between neurons obtained from spike trains are more abundant than conventional morphological data. We constructed two open-access computational platforms of various neuronal circuits from the Blue Brain Project realistic models, named Neurpy and Neurgen. Then, we investigated how we could perform network tomography with cortical neuronal circuits for the morphological, topological and electrical classification of neurons. We extracted the simulated data of 10,000 network topology combinations with five layers, 25 morphological type (m-type) cells, and 14 electrical type (e-type) cells. We applied the data to several different classifiers (including Support Vector Machine (SVM), Decision Trees, Random Forest, and Artificial Neural Networks). We achieved accuracies of up to 70%, and the inference of biological network structures using network tomography reached up to 65% of accuracy. Objective classification of biological networks can be achieved with cascaded machine learning methods using neuron communication data. SVM methods seem to perform better amongst used techniques. Our research not only contributes to existing classification efforts but sets the road-map for future usage of brain–machine interfaces towards an in vivo objective classification of neurons as a sensing mechanism of the brain’s structure.

Computational models reveal how chloride dynamics determine the optimal distribution of inhibitory synapses to minimise dendritic excitability

Many neurons in the mammalian central nervous system have complex dendritic arborisations and active dendritic conductances that enable these cells to perform sophisticated computations. How dendritically targeted inhibition affects local dendritic excitability is not fully understood. Here we use computational models of branched dendrites to investigate where GABAergic synapses should be placed to minimise dendritic excitability over time. To do so, we formulate a metric we term the “Inhibitory Level” (IL), which quantifies the effectiveness of synaptic inhibition for reducing the depolarising effect of nearby excitatory input. GABAergic synaptic inhibition is dependent on the reversal potential for GABA_A receptors (EGABA), which is primarily set by the transmembrane chloride ion (Cl^-) concentration gradient. We, therefore, investigated how variable EGABA and dynamic chloride affects dendritic inhibition. We found that the inhibitory effectiveness of dendritic GABAergic synapses combines at an encircled branch junction. The extent of this inhibitory accumulation is dependent on the number of branches and location of synapses but is independent of EGABA. This inhibitory accumulation occurs even for very distally placed inhibitory synapses when they are hyperpolarising–but not when they are shunting. When accounting for Cl^- fluxes and dynamics in Cl^- concentration, we observed that Cl^- loading is detrimental to inhibitory effectiveness. This enabled us to determine the most inhibitory distribution of GABAergic synapses which is close to–but not at–a shared branch junction. This distribution balances a trade-off between a stronger combined inhibitory influence when synapses closely encircle a branch junction with the deleterious effects of increased Cl^- by loading that occurs when inhibitory synapses are co-located.

Dendritic branches allow for a rich repertoire of computational capabilities for neurons within the brain. Inhibitory synaptic inputs, which utilise the neurotransmitter GABA, refine and enhance dendritic computations. They are traditionally viewed with regards to their inhibitory effect on action potential generation at the neuronal cell body. Here, we studied the local effects of inhibitory synapses on excitability in dendrites. We also considered the dynamic nature of inhibition that deteriorates the longer it is active due to intracellular chloride ion loading. The central goal of our investigation was to find the best locations for multiple inhibitory synapses to maximise their combined inhibitory effectiveness on nearby excitation in the dendritic tree. We found that the optimal distribution is when inhibitory synapses closely encircle a branch junction, without being co-located at the junction itself. This maximises how their inhibitory influence combines whilst minimising the deleterious effects of chloride loading.

Large-scale biophysically detailed model of somatosensory thalamocortical circuits in NetPyNE

The primary somatosensory cortex (S1) of mammals is critically important in the perception of touch and related sensorimotor behaviors. In 2015, the Blue Brain Project (BBP) developed a groundbreaking rat S1 microcircuit simulation with over 31,000 neurons with 207 morpho-electrical neuron types, and 37 million synapses, incorporating anatomical and physiological information from a wide range of experimental studies. We have implemented this highly detailed and complex S1 model in NetPyNE, using the data available in the Neocortical Microcircuit Collaboration Portal. NetPyNE provides a Python high-level interface to NEURON and allows defining complicated multiscale models using an intuitive declarative standardized language. It also facilitates running parallel simulations, automates the optimization and exploration of parameters using supercomputers, and provides a wide range of built-in analysis functions. This will make the S1 model more accessible and simpler to scale, modify and extend in order to explore research questions or interconnect to other existing models. Despite some implementation differences, the NetPyNE model preserved the original cell morphologies, electrophysiological responses and spatial distribution for all 207 cell types; and the connectivity properties of all 1941 pathways, including synaptic dynamics and short-term plasticity (STP). The NetPyNE S1 simulations produced reasonable physiological firing rates and activity patterns across all populations. When STP was included, the network generated a 1 Hz oscillation comparable to the original model in vitro-like state. By then reducing the extracellular calcium concentration, the model reproduced the original S1 in vivo-like states with asynchronous activity. These results validate the original study using a new modeling tool. Simulated local field potentials (LFPs) exhibited realistic oscillatory patterns and features, including distance- and frequency-dependent attenuation. The model was extended by adding thalamic circuits, including 6 distinct thalamic populations with intrathalamic, thalamocortical (TC) and corticothalamic connectivity derived from experimental data. The thalamic model reproduced single known cell and circuit-level dynamics, including burst and tonic firing modes and oscillatory patterns, providing a more realistic input to cortex and enabling study of TC interactions. Overall, our work provides a widely accessible, data-driven and biophysically-detailed model of the somatosensory TC circuits that can be employed as a community tool for researchers to study neural dynamics, function and disease.

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

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

Constructing functional models from biophysically-detailed neurons

Improving biological plausibility and functional capacity are two important goals for brain models that connect low-level neural details to high-level behavioral phenomena. We develop a method called “oracle-supervised Neural Engineering Framework” (osNEF) to train biologically-detailed spiking neural networks that realize a variety of cognitively-relevant dynamical systems. Specifically, we train networks to perform computations that are commonly found in cognitive systems (communication, multiplication, harmonic oscillation, and gated working memory) using four distinct neuron models (leaky-integrate-and-fire neurons, Izhikevich neurons, 4-dimensional nonlinear point neurons, and 4-compartment, 6-ion-channel layer-V pyramidal cell reconstructions) connected with various synaptic models (current-based synapses, conductance-based synapses, and voltage-gated synapses). We show that osNEF networks exhibit the target dynamics by accounting for nonlinearities present within the neuron models: performance is comparable across all four systems and all four neuron models, with variance proportional to task and neuron model complexity. We also apply osNEF to build a model of working memory that performs a delayed response task using a combination of pyramidal cells and inhibitory interneurons connected with NMDA and GABA synapses. The baseline performance and forgetting rate of the model are consistent with animal data from delayed match-to-sample tasks (DMTST): we observe a baseline performance of 95% and exponential forgetting with time constant τ = 8.5s, while a recent meta-analysis of DMTST performance across species observed baseline performances of 58 − 99% and exponential forgetting with time constants of τ = 2.4 − 71s. These results demonstrate that osNEF can train functional brain models using biologically-detailed components and open new avenues for investigating the relationship between biophysical mechanisms and functional capabilities.

Computational models of biologically realistic neural networks help scientists understand and recreate a wide variety of brain processes, responsible for everything from fish locomotion to human cognition. To be useful, these models must both recreate features of the brain, such as the electrical, chemical, and geometric properties of neurons, and perform useful functional operations, such as storing and retrieving information from a short term memory. Here, we develop a new method for training networks built from biologically detailed components. We simulate networks that contain a variety of complex neurons and synapses, then show that our method successfully trains them to perform a variety of cognitive operations. Most notably, we train a working memory model that contains detailed reconstructions of cortical neurons, and demonstrate that it performs a memory task with performance that is comparable to simple animals. Researchers can use our method to train detailed brain models and investigate how biological features (or deficits thereof) relate to cognition, which may provide insights into the biological basis of mental disorders such as Parkinson’s disease.

Connectivity concepts in neuronal network modeling

Sustainable research on computational models of neuronal networks requires published models to be understandable, reproducible, and extendable. Missing details or ambiguities about mathematical concepts and assumptions, algorithmic implementations, or parameterizations hinder progress. Such flaws are unfortunately frequent and one reason is a lack of readily applicable standards and tools for model description. Our work aims to advance complete and concise descriptions of network connectivity but also to guide the implementation of connection routines in simulation software and neuromorphic hardware systems. We first review models made available by the computational neuroscience community in the repositories ModelDB and Open Source Brain, and investigate the corresponding connectivity structures and their descriptions in both manuscript and code. The review comprises the connectivity of networks with diverse levels of neuroanatomical detail and exposes how connectivity is abstracted in existing description languages and simulator interfaces. We find that a substantial proportion of the published descriptions of connectivity is ambiguous. Based on this review, we derive a set of connectivity concepts for deterministically and probabilistically connected networks and also address networks embedded in metric space. Beside these mathematical and textual guidelines, we propose a unified graphical notation for network diagrams to facilitate an intuitive understanding of network properties. Examples of representative network models demonstrate the practical use of the ideas. We hope that the proposed standardizations will contribute to unambiguous descriptions and reproducible implementations of neuronal network connectivity in computational neuroscience.

Signal processing and computational modeling for interpretation of SEEG-recorded interictal epileptiform discharges in epileptogenic and non-epileptogenic zones

OBJECTIVE

In partial epilepsies, interictal epileptiform discharges (IEDs) are paroxysmal events observed in epileptogenic and non-epileptogenic zones. IEDs' generation and recurrence are subject to different hypotheses: they appear through glutamatergic and GABAergic processes; they may trigger seizures or prevent seizure propagation. This paper focuses on a specific class of IEDs, spike-waves (SWs), characterized by a short-duration spike followed by a longer duration wave, both of the same polarity. Signal analysis and neurophysiological mathematical models are used to interpret puzzling IED generation.

APPROACH

Interictal activity was recorded by intracranial stereo-electroencephalography (SEEG) electrodes in five different patients. SEEG experts identified the epileptic and non-epileptic zones in which IEDs were detected. After quantifying spatial and temporal features of the detected IEDs, the most significant features for classifying epileptic and non-epileptic zones were determined. A neurophysiologically-plausible mathematical model was then introduced to simulate the IEDs and understand the underlying differences observed in epileptic and non-epileptic zone IEDs.

MAIN RESULTS

Two classes of SWs were identified according to subtle differences in morphology and timing of the spike and wave component. Results showed that type-1 SWs were generated in epileptogenic regions also involved at seizure onset, while type-2 SWs were produced in the propagation or non-involved areas. The modeling study indicated that synaptic kinetics, cortical organization, and network interactions determined the morphology of the simulated SEEG signals. Modeling results suggested that the IED morphologies were linked to the degree of preserved inhibition.

SIGNIFICANCE

This work contributes to the understanding of different mechanisms generating IEDs in epileptic networks. The combination of signal analysis and computational models provides an efficient framework for exploring IEDs in partial epilepsies and classifying epileptogenic and non-epileptogenic zones.

Repeated methamphetamine administration produces cognitive deficits through augmentation of GABAergic synaptic transmission in the prefrontal cortex

Methamphetamine (METH) abuse is associated with the emergence of cognitive deficits and hypofrontality, a pathophysiological marker of many neuropsychiatric disorders that is produced by altered balance of local excitatory and inhibitory synaptic transmission. However, there is a dearth of information regarding the cellular and synaptic mechanisms underlying METH-induced cognitive deficits and associated hypofrontal states. Using PV-Cre transgenic rats that went through a METH sensitization regime or saline (SAL) followed by 7-10 days of home cage abstinence combined with cognitive tests, chemogenetic experiments, and whole-cell patch recordings on the prelimbic prefrontal cortex (PFC), we investigated the cellular and synaptic mechanisms underlying METH-induce hypofrontality. We report here that repeated METH administration in rats produces deficits in working memory and increases in inhibitory synaptic transmission onto pyramidal neurons in the PFC. The increased PFC inhibition is detected by an increase in spontaneous and evoked inhibitory postsynaptic synaptic currents (IPSCs), an increase in GABAergic presynaptic function, and a shift in the excitatory-inhibitory balance onto PFC deep-layer pyramidal neurons. We find that pharmacological blockade of D1 dopamine receptor function reduces the METH-induced augmentation of IPSCs, suggesting a critical role for D1 dopamine signaling in METH-induced hypofrontality. In addition, repeated METH administration increases the intrinsic excitability of parvalbumin-positive fast spiking interneurons (PV + FSIs), a key local interneuron population in PFC that contributes to the control of inhibitory tone. Using a cell type-specific chemogenetic approach, we show that increasing PV + FSIs activity in the PFC is necessary and sufficient to cause deficits in temporal order memory similar to those induced by METH. Conversely, reducing PV + FSIs activity in the PFC of METH-exposed rats rescues METH-induced temporal order memory deficits. Together, our findings reveal that repeated METH exposure increases PFC inhibitory tone through a D1 dopamine signaling-dependent potentiation of inhibitory synaptic transmission, and that reduction of PV + FSIs activity can rescue METH-induced cognitive deficits, suggesting a potential therapeutic approach to treating cognitive symptoms in patients suffering from METH use disorder.

Buice M. MouseNet: A biologically constrained convolutional neural network model for the mouse visual cortex

Convolutional neural networks trained on object recognition derive inspiration from the neural architecture of the visual system in mammals, and have been used as models of the feedforward computation performed in the primate ventral stream. In contrast to the deep hierarchical organization of primates, the visual system of the mouse has a shallower arrangement. Since mice and primates are both capable of visually guided behavior, this raises questions about the role of architecture in neural computation. In this work, we introduce a novel framework for building a biologically constrained convolutional neural network model of the mouse visual cortex. The architecture and structural parameters of the network are derived from experimental measurements, specifically the 100-micrometer resolution interareal connectome, the estimates of numbers of neurons in each area and cortical layer, and the statistics of connections between cortical layers. This network is constructed to support detailed task-optimized models of mouse visual cortex, with neural populations that can be compared to specific corresponding populations in the mouse brain. Using a well-studied image classification task as our working example, we demonstrate the computational capability of this mouse-sized network. Given its relatively small size, MouseNet achieves roughly 2/3rds the performance level on ImageNet as VGG16. In combination with the large scale Allen Brain Observatory Visual Coding dataset, we use representational similarity analysis to quantify the extent to which MouseNet recapitulates the neural representation in mouse visual cortex. Importantly, we provide evidence that optimizing for task performance does not improve similarity to the corresponding biological system beyond a certain point. We demonstrate that the distributions of some physiological quantities are closer to the observed distributions in the mouse brain after task training. We encourage the use of the MouseNet architecture by making the code freely available.

Translational neuronal ensembles: Neuronal microcircuits in psychology, physiology, pharmacology and pathology

Multi-recording techniques show evidence that neurons coordinate their firing forming ensembles and that brain networks are made by connections between ensembles. While “canonical” microcircuits are composed of interconnected principal neurons and interneurons, it is not clear how they participate in recorded neuronal ensembles: “groups of neurons that show spatiotemporal co-activation”. Understanding synapses and their plasticity has become complex, making hard to consider all details to fill the gap between cellular-synaptic and circuit levels. Therefore, two assumptions became necessary: First, whatever the nature of the synapses these may be simplified by “functional connections”. Second, whatever the mechanisms to achieve synaptic potentiation or depression, the resultant synaptic weights are relatively stable. Both assumptions have experimental basis cited in this review, and tools to analyze neuronal populations are being developed based on them. Microcircuitry processing followed with multi-recording techniques show temporal sequences of neuronal ensembles resembling computational routines. These sequences can be aligned with the steps of behavioral tasks and behavior can be modified upon their manipulation, supporting the hypothesis that they are memory traces. In vitro, recordings show that these temporal sequences can be contained in isolated tissue of histological scale. Sequences found in control conditions differ from those recorded in pathological tissue obtained from animal disease models and those recorded after the actions of clinically useful drugs to treat disease states, setting the basis for new bioassays to test drugs with potential clinical use. These findings make the neuronal ensembles theoretical framework a dynamic neuroscience paradigm.

From descriptive connectome to mechanistic connectome: Generative modeling in functional magnetic resonance imaging analysis

As a newly emerging field, connectomics has greatly advanced our understanding of the wiring diagram and organizational features of the human brain. Generative modeling-based connectome analysis, in particular, plays a vital role in deciphering the neural mechanisms of cognitive functions in health and dysfunction in diseases. Here we review the foundation and development of major generative modeling approaches for functional magnetic resonance imaging (fMRI) and survey their applications to cognitive or clinical neuroscience problems. We argue that conventional structural and functional connectivity (FC) analysis alone is not sufficient to reveal the complex circuit interactions underlying observed neuroimaging data and should be supplemented with generative modeling-based effective connectivity and simulation, a fruitful practice that we term "mechanistic connectome." The transformation from descriptive connectome to mechanistic connectome will open up promising avenues to gain mechanistic insights into the delicate operating principles of the human brain and their potential impairments in diseases, which facilitates the development of effective personalized treatments to curb neurological and psychiatric disorders.

A realistic morpho-anatomical connection strategy for modelling full-scale point-neuron microcircuits

The modeling of extended microcircuits is emerging as an effective tool to simulate the neurophysiological correlates of brain activity and to investigate brain dysfunctions. However, for specific networks, a realistic modeling approach based on the combination of available physiological, morphological and anatomical data is still an open issue. One of the main problems in the generation of realistic networks lies in the strategy adopted to build network connectivity. Here we propose a method to implement a neuronal network at single cell resolution by using the geometrical probability volumes associated with pre- and postsynaptic neurites. This allows us to build a network with plausible connectivity properties without the explicit use of computationally intensive touch detection algorithms using full 3D neuron reconstructions. The method has been benchmarked for the mouse hippocampus CA1 area, and the results show that this approach is able to generate full-scale brain networks at single cell resolution that are in good agreement with experimental findings. This geometric reconstruction of axonal and dendritic occupancy, by effectively reflecting morphological and anatomical constraints, could be integrated into structured simulators generating entire circuits of different brain areas facilitating the simulation of different brain regions with realistic models.

Brain signal predictions from multi-scale networks using a linearized framework

Simulations of neural activity at different levels of detail are ubiquitous in modern neurosciences, aiding the interpretation of experimental data and underlying neural mechanisms at the level of cells and circuits. Extracellular measurements of brain signals reflecting transmembrane currents throughout the neural tissue remain commonplace. The lower frequencies (≲ 300Hz) of measured signals generally stem from synaptic activity driven by recurrent interactions among neural populations and computational models should also incorporate accurate predictions of such signals. Due to limited computational resources, large-scale neuronal network models (≳ 10⁶ neurons or so) often require reducing the level of biophysical detail and account mainly for times of action potentials (‘spikes’) or spike rates. Corresponding extracellular signal predictions have thus poorly accounted for their biophysical origin.

Here we propose a computational framework for predicting spatiotemporal filter kernels for such extracellular signals stemming from synaptic activity, accounting for the biophysics of neurons, populations, and recurrent connections. Signals are obtained by convolving population spike rates by appropriate kernels for each connection pathway and summing the contributions. Our main results are that kernels derived via linearized synapse and membrane dynamics, distributions of cells, conduction delay, and volume conductor model allow for accurately capturing the spatiotemporal dynamics of ground truth extracellular signals from conductance-based multicompartment neuron networks. One particular observation is that changes in the effective membrane time constants caused by persistent synapse activation must be accounted for.

The work also constitutes a major advance in computational efficiency of accurate, biophysics-based signal predictions from large-scale spike and rate-based neuron network models drastically reducing signal prediction times compared to biophysically detailed network models. This work also provides insight into how experimentally recorded low-frequency extracellular signals of neuronal activity may be approximately linearly dependent on spiking activity. A new software tool LFPykernels serves as a reference implementation of the framework.

Understanding the brain’s function and activity in healthy and pathological states across spatial scales and times spanning entire lives is one of humanity’s great undertakings. In experimental and clinical work probing the brain’s activity, a variety of electric and magnetic measurement techniques are routinely applied. However interpreting the extracellularly measured signals remains arduous due to multiple factors, mainly the large number of neurons contributing to the signals and complex interactions occurring in recurrently connected neuronal circuits. To understand how neurons give rise to such signals, mechanistic modeling combined with forward models derived using volume conductor theory has proven to be successful, but this approach currently does not scale to the systems level (encompassing millions of neurons or more) where simplified or abstract neuron representations typically are used. Motivated by experimental findings implying approximately linear relationships between times of neuronal action potentials and extracellular population signals, we provide a biophysics-based method for computing causal filters relating spikes and extracellular signals that can be applied with spike times or rates of large-scale neuronal network models for predictions of population signals without relying on ad hoc approximations.

The quest for multiscale brain modeling

Addressing the multiscale organization of the brain, which is fundamental to the dynamic repertoire of the organ, remains challenging. In principle, it should be possible to model neurons and synapses in detail and then connect them into large neuronal assemblies to explain the relationship between microscopic phenomena, large-scale brain functions, and behavior. It is more difficult to infer neuronal functions from ensemble measurements such as those currently obtained with brain activity recordings. In this article we consider theories and strategies for combining bottom-up models, generated from principles of neuronal biophysics, with top-down models based on ensemble representations of network activity and on functional principles. These integrative approaches are hoped to provide effective multiscale simulations in virtual brains and neurorobots, and pave the way to future applications in medicine and information technologies.

Neural Algorithms and Circuits for Motor Planning

The brain plans and executes volitional movements. The underlying patterns of neural population activity have been explored in the context of movements of the eyes, limbs, tongue, and head in nonhuman primates and rodents. How do networks of neurons produce the slow neural dynamics that prepare specific movements and the fast dynamics that ultimately initiate these movements? Recent work exploits rapid and calibrated perturbations of neural activity to test specific dynamical systems models that are capable of producing the observed neural activity. These joint experimental and computational studies show that cortical dynamics during motor planning reflect fixed points of neural activity (attractors). Subcortical control signals reshape and move attractors over multiple timescales, causing commitment to specific actions and rapid transitions to movement execution. Experiments in rodents are beginning to reveal how these algorithms are implemented at the level of brain-wide neural circuits.

Combining hypothesis- and data-driven neuroscience modeling in FAIR workflows

Modeling in neuroscience occurs at the intersection of different points of view and approaches. Typically, hypothesis-driven modeling brings a question into focus so that a model is constructed to investigate a specific hypothesis about how the system works or why certain phenomena are observed. Data-driven modeling, on the other hand, follows a more unbiased approach, with model construction informed by the computationally intensive use of data. At the same time, researchers employ models at different biological scales and at different levels of abstraction. Combining these models while validating them against experimental data increases understanding of the multiscale brain. However, a lack of interoperability, transparency, and reusability of both models and the workflows used to construct them creates barriers for the integration of models representing different biological scales and built using different modeling philosophies. We argue that the same imperatives that drive resources and policy for data - such as the FAIR (Findable, Accessible, Interoperable, Reusable) principles - also support the integration of different modeling approaches. The FAIR principles require that data be shared in formats that are Findable, Accessible, Interoperable, and Reusable. Applying these principles to models and modeling workflows, as well as the data used to constrain and validate them, would allow researchers to find, reuse, question, validate, and extend published models, regardless of whether they are implemented phenomenologically or mechanistically, as a few equations or as a multiscale, hierarchical system. To illustrate these ideas, we use a classical synaptic plasticity model, the Bienenstock-Cooper-Munro rule, as an example due to its long history, different levels of abstraction, and implementation at many scales.

Beyond LIF Neurons on Neuromorphic Hardware

Neuromorphic systems aim to provide accelerated low-power simulation of Spiking Neural Networks (SNNs), typically featuring simple and efficient neuron models such as the Leaky Integrate-and-Fire (LIF) model. Biologically plausible neuron models developed by neuroscientists are largely ignored in neuromorphic computing due to their increased computational costs. This work bridges this gap through implementation and evaluation of a single compartment Hodgkin-Huxley (HH) neuron and a multi-compartment neuron incorporating dendritic computation on the SpiNNaker, and SpiNNaker2 prototype neuromorphic systems. Numerical accuracy of the model implementations is benchmarked against reference models in the NEURON simulation environment, with excellent agreement achieved by both the fixed- and floating-point SpiNNaker implementations. The computational cost is evaluated in terms of timing measurements profiling neural state updates. While the additional model complexity understandably increases computation times relative to LIF models, it was found a wallclock time increase of only 8× was observed for the HH neuron (11× for the mutlicompartment model), demonstrating the potential of hardware accelerators in the next-generation neuromorphic system to optimize implementation of complex neuron models. The benefits of models directly corresponding to biophysiological data are demonstrated: HH neurons are able to express a range of output behaviors not captured by LIF neurons; and the dendritic compartment provides the first implementation of a spiking multi-compartment neuron model with XOR-solving capabilities on neuromorphic hardware. The work paves the way for inclusion of more biologically representative neuron models in neuromorphic systems, and showcases the benefits of hardware accelerators included in the next-generation SpiNNaker2 architecture.

Hyperbolic disc embedding of functional human brain connectomes using resting-state fMRI

The brain presents a real complex network of modular, small-world, and hierarchical nature, which are features of non-Euclidean geometry. Using resting-state functional magnetic resonance imaging, we constructed a scale-free binary graph for each subject, using internodal time series correlation of regions of interest as a proximity measure. The resulting network could be embedded onto manifolds of various curvatures and dimensions. While maintaining the fidelity of embedding (low distortion, high mean average precision), functional brain networks were found to be best represented in the hyperbolic disc. Using the 𝕊¹/ℍ² model, we reduced the dimension of the network into two-dimensional hyperbolic space and were able to efficiently visualize the internodal connections of the brain, preserving proximity as distances and angles on the hyperbolic discs. Each individual disc revealed relevance with its anatomic counterpart and absence of center-spaced node. Using the hyperbolic distance on the 𝕊¹/ℍ² model, we could detect the anomaly of network in autism spectrum disorder subjects. This procedure of embedding grants us a reliable new framework for studying functional brain networks and the possibility of detecting anomalies of the network in the hyperbolic disc on an individual scale.

Normalizing the brain connectome for communication through synchronization

Networks in neuroscience determine how brain function unfolds, and their perturbations lead to psychiatric disorders and brain disease. Brain networks are characterized by their connectomes, which comprise the totality of all connections, and are commonly described by graph theory. This approach is deeply rooted in a particle view of information processing, based on the quantification of informational bits such as firing rates. Oscillations and brain rhythms demand, however, a wave perspective of information processing based on synchronization. We extend traditional graph theory to a dual, particle-wave, perspective, integrate time delays due to finite transmission speeds, and derive a normalization of the connectome. When applied to the database of the Human Connectome Project, it explains the emergence of frequency-specific network cores including the visual and default mode networks. These findings are robust across human subjects (N = 100) and are a fundamental network property within the wave picture. The normalized connectome comprises the particle view in the limit of infinite transmission speeds and opens the applicability of graph theory to a wide range of novel network phenomena, including physiological and pathological brain rhythms. These two perspectives are orthogonal, but not incommensurable, when understood within the novel, here-proposed, generalized framework of structural connectivity.

The Neuron Phenotype Ontology: A FAIR Approach to Proposing and Classifying Neuronal Types

The challenge of defining and cataloging the building blocks of the brain requires a standardized approach to naming neurons and organizing knowledge about their properties. The US Brain Initiative Cell Census Network, Human Cell Atlas, Blue Brain Project, and others are generating vast amounts of data and characterizing large numbers of neurons throughout the nervous system. The neuroscientific literature contains many neuron names (e.g. parvalbumin-positive interneuron or layer 5 pyramidal cell) that are commonly used and generally accepted. However, it is often unclear how such common usage types relate to many evidence-based types that are proposed based on the results of new techniques. Further, comparing different types across labs remains a significant challenge. Here, we propose an interoperable knowledge representation, the Neuron Phenotype Ontology (NPO), that provides a standardized and automatable approach for naming cell types and normalizing their constituent phenotypes using identifiers from community ontologies as a common language. The NPO provides a framework for systematically organizing knowledge about cellular properties and enables interoperability with existing neuron naming schemes. We evaluate the NPO by populating a knowledge base with three independent cortical neuron classifications derived from published data sets that describe neurons according to molecular, morphological, electrophysiological, and synaptic properties. Competency queries to this knowledge base demonstrate that the NPO knowledge model enables interoperability between the three test cases and neuron names commonly used in the literature.

State-dependent activity dynamics of hypothalamic stress effector neurons

The stress response necessitates an immediate boost in vital physiological functions from their homeostatic operation to an elevated emergency response. However, the neural mechanisms underlying this state-dependent change remain largely unknown. Using a combination of in vivo and ex vivo electrophysiology with computational modeling, we report that corticotropin releasing hormone (CRH) neurons in the paraventricular nucleus of the hypothalamus (PVN), the effector neurons of hormonal stress response, rapidly transition between distinct activity states through recurrent inhibition. Specifically, in vivo optrode recording shows that under non-stress conditions, CRH_PVN neurons often fire with rhythmic brief bursts (RB), which, somewhat counterintuitively, constrains firing rate due to long (~2 s) interburst intervals. Stressful stimuli rapidly switch RB to continuous single spiking (SS), permitting a large increase in firing rate. A spiking network model shows that recurrent inhibition can control this activity-state switch, and more broadly the gain of spiking responses to excitatory inputs. In biological CRH_PVN neurons ex vivo, the injection of whole-cell currents derived from our computational model recreates the in vivo-like switch between RB and SS, providing direct evidence that physiologically relevant network inputs enable state-dependent computation in single neurons. Together, we present a novel mechanism for state-dependent activity dynamics in CRH_PVN neurons.

Modernizing the NEURON Simulator for Sustainability, Portability, and Performance

The need for reproducible, credible, multiscale biological modeling has led to the development of standardized simulation platforms, such as the widely-used NEURON environment for computational neuroscience. Developing and maintaining NEURON over several decades has required attention to the competing needs of backwards compatibility, evolving computer architectures, the addition of new scales and physical processes, accessibility to new users, and efficiency and flexibility for specialists. In order to meet these challenges, we have now substantially modernized NEURON, providing continuous integration, an improved build system and release workflow, and better documentation. With the help of a new source-to-source compiler of the NMODL domain-specific language we have enhanced NEURON's ability to run efficiently, via the CoreNEURON simulation engine, on a variety of hardware platforms, including GPUs. Through the implementation of an optimized in-memory transfer mechanism this performance optimized backend is made easily accessible to users, providing training and model-development paths from laptop to workstation to supercomputer and cloud platform. Similarly, we have been able to accelerate NEURON's reaction-diffusion simulation performance through the use of just-in-time compilation. We show that these efforts have led to a growing developer base, a simpler and more robust software distribution, a wider range of supported computer architectures, a better integration of NEURON with other scientific workflows, and substantially improved performance for the simulation of biophysical and biochemical models.

Modeling cortical synaptic effects of anesthesia and their cholinergic reversal

General anesthetics work through a variety of molecular mechanisms while resulting in the common end point of sedation and loss of consciousness. Generally, the administration of common anesthetics induces reduction in synaptic excitation while promoting synaptic inhibition. Exogenous modulation of the anesthetics' synaptic effects can help determine the neuronal pathways involved in anesthesia. For example, both animal and human studies have shown that exogenously induced increases in acetylcholine in the brain can elicit wakeful-like behavior despite the continued presence of the anesthetic. However, the underlying mechanisms of anesthesia reversal at the cellular level have not been investigated. Here we apply a computational model of a network of excitatory and inhibitory neurons to simulate the network-wide effects of anesthesia, due to changes in synaptic inhibition and excitation, and their reversal by cholinergic activation through muscarinic receptors. We use a differential evolution algorithm to fit model parameters to match measures of spiking activity, neuronal connectivity, and network dynamics recorded in the visual cortex of rodents during anesthesia with desflurane in vivo. We find that facilitating muscarinic receptor effects of acetylcholine on top of anesthetic-induced synaptic changes predicts the reversal of anesthetic suppression of neurons' spiking activity, functional connectivity, as well as pairwise and population interactions. Thus, our model predicts a specific neuronal mechanism for the cholinergic reversal of anesthesia consistent with experimental behavioral observations.

Biophysical Kv3 channel alterations dampen excitability of cortical PV interneurons and contribute to network hyperexcitability in early Alzheimer's

In Alzheimer’s disease (AD), a multitude of genetic risk factors and early biomarkers are known. Nevertheless, the causal factors responsible for initiating cognitive decline in AD remain controversial. Toxic plaques and tangles correlate with progressive neuropathology, yet disruptions in circuit activity emerge before their deposition in AD models and patients. Parvalbumin (PV) interneurons are potential candidates for dysregulating cortical excitability as they display altered action potential (AP) firing before neighboring excitatory neurons in prodromal AD. Here, we report a novel mechanism responsible for PV hypoexcitability in young adult familial AD mice. We found that biophysical modulation of K_v3 channels, but not changes in their mRNA or protein expression, were responsible for dampened excitability in young 5xFAD mice. These K⁺ conductances could efficiently regulate near-threshold AP firing, resulting in gamma-frequency-specific network hyperexcitability. Thus, biophysical ion channel alterations alone may reshape cortical network activity prior to changes in their expression levels. Our findings demonstrate an opportunity to design a novel class of targeted therapies to ameliorate cortical circuit hyperexcitability in early AD.

Scaling and Benchmarking an Evolutionary Algorithm for Constructing Biophysical Neuronal Models

Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.

Morphological pseudotime ordering and fate mapping reveal diversification of cerebellar inhibitory interneurons

Understanding how diverse neurons are assembled into circuits requires a framework for describing cell types and their developmental trajectories. Here we combine genetic fate-mapping, pseudotemporal profiling of morphogenesis, and dual morphology and RNA labeling to resolve the diversification of mouse cerebellar inhibitory interneurons. Molecular layer interneurons (MLIs) derive from a common progenitor population but comprise diverse dendritic-, somatic-, and axon initial segment-targeting interneurons. Using quantitative morphology from 79 mature MLIs, we identify two discrete morphological types and presence of extensive within-class heterogeneity. Pseudotime trajectory inference using 732 developmental morphologies indicate the emergence of distinct MLI types during migration, before reaching their final positions. By comparing MLI identities from morphological and transcriptomic signatures, we demonstrate the dissociation between these modalities and that subtype divergence can be resolved from axonal morphogenesis prior to marker gene expression. Our study illustrates the utility of applying single-cell methods to quantify morphology for defining neuronal diversification.