Reconstruction and Simulation of Neocortical Microcircuitry

Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons

Understanding the brain requires studying its multiscale interactions from molecules to networks. The increasing availability of large-scale datasets detailing brain circuit composition, connectivity, and activity is transforming neuroscience. However, integrating and interpreting this data remains challenging. Concurrently, advances in supercomputing and sophisticated modeling tools now enable the development of highly detailed, large-scale biophysical circuit models. These mechanistic multiscale models offer a method to systematically integrate experimental data, facilitating investigations into brain structure, function, and disease. This review, based on a Society for Neuroscience 2024 MiniSymposium, aims to disseminate recent advances in large-scale mechanistic modeling to the broader community. It highlights (1) examples of current models for various brain regions developed through experimental data integration; (2) their predictive capabilities regarding cellular and circuit mechanisms underlying experimental recordings (e.g., membrane voltage, spikes, local-field potential, electroencephalography/magnetoencephalography) and brain function; and (3) their use in simulating biomarkers for brain diseases like epilepsy, depression, schizophrenia, and Parkinson's, aiding in understanding their biophysical underpinnings and developing novel treatments. The review showcases state-of-the-art models covering hippocampus, somatosensory, visual, motor, auditory cortical, and thalamic circuits across species. These models predict neural activity at multiple scales and provide insights into the biophysical mechanisms underlying sensation, motor behavior, brain signals, neural coding, disease, pharmacological interventions, and neural stimulation. Collaboration with experimental neuroscientists and clinicians is essential for the development and validation of these models, particularly as datasets grow. Hence, this review aims to foster interest in detailed brain circuit models, leading to cross-disciplinary collaborations that accelerate brain research.

Theoretical framework for learning through structural plasticity

A growing body of research indicates that structural plasticity mechanisms are crucial for learning and memory consolidation. Starting from a simple phenomenological model, we exploit a mean-field approach to develop a theoretical framework of learning through this kind of plasticity, capable of taking into account several features of the connectivity and pattern of activity of biological neural networks, including probability distributions of neuron firing rates, selectivity of the responses of single neurons to multiple stimuli, probabilistic connection rules, and noisy stimuli. More importantly, it describes the effects of stabilization, pruning, and reorganization of synaptic connections. This framework is used to compute the values of some relevant quantities used to characterize the learning and memory capabilities of the neuronal network in training and testing procedures as the number of training patterns and other model parameters vary. The results are then compared with those obtained through simulations with firing-rate-based neuronal network models.

Areal specializations in the morpho-electric and transcriptomic properties of primate layer 5 extratelencephalic projection neurons

Large-scale analysis of single-cell gene expression has revealed transcriptomically defined cell subclasses present throughout the primate neocortex with gene expression profiles that differ depending upon neocortical region. Here, we test whether the interareal differences in gene expression translate to regional specializations in the physiology and morphology of infragranular glutamatergic neurons by performing Patch-seq experiments in brain slices from the temporal cortex (TCx) and motor cortex (MCx) of the macaque. We confirm that transcriptomically defined extratelencephalically projecting neurons of layer 5 (L5 ET neurons) include retrogradely labeled corticospinal neurons in the MCx and find multiple physiological properties and ion channel genes that distinguish L5 ET from non-ET neurons in both areas. Additionally, while infragranular ET and non-ET neurons retain distinct neuronal properties across multiple regions, there are regional morpho-electric and gene expression specializations in the L5 ET subclass, providing mechanistic insights into the specialized functional architecture of the primate neocortex.

Biostasis: A Roadmap for Research in Preservation and Potential Revival of Humans

Human biostasis, the preservation of a human when all other contemporary options for extension of quality life are exhausted, offers the speculative potential for survival via continuation of life in the future. While provably reversible preservation, also known as suspended animation, is not yet possible for humans, the primary justification for contemporary biostasis is the preservation of the brain, which is broadly considered the seat of memories, personality, and identity. By preserving the information contained within the brain's structures, it may be possible to resuscitate a healthy whole individual using advanced future technologies. There are numerous challenges in biostasis, including inadequacies in current preservation techniques, methods to evaluate the quality of preservation, and potential future revival technologies. In this report, we describe a roadmap that attempts to delineate research directions that could improve the field of biostasis, focusing on optimizing preservation protocols and establishing metrics for querying preservation quality, as well as pre- and post-cardiac arrest factors, stabilization strategies, and methods for long-term preservation. We acknowledge the highly theoretical nature of future revival technologies and the importance of achieving high-fidelity brain preservation to maximize the potential of future repair technologies. We plan to update the research roadmap biennially. Our goal is to encourage multidisciplinary communication and collaboration in this field.

Information Thermodynamics: From Physics to Neuroscience

This paper provides a perspective on applying the concepts of information thermodynamics, developed recently in non-equilibrium statistical physics, to problems in theoretical neuroscience. Historically, information and energy in neuroscience have been treated separately, in contrast to physics approaches, where the relationship of entropy production with heat is a central idea. It is argued here that also in neural systems, information and energy can be considered within the same theoretical framework. Starting from basic ideas of thermodynamics and information theory on a classic Brownian particle, it is shown how noisy neural networks can infer its probabilistic motion. The decoding of the particle motion by neurons is performed with some accuracy, and it has some energy cost, and both can be determined using information thermodynamics. In a similar fashion, we also discuss how neural networks in the brain can learn the particle velocity and maintain that information in the weights of plastic synapses from a physical point of view. Generally, it is shown how the framework of stochastic and information thermodynamics can be used practically to study neural inference, learning, and information storing.

DendroTweaks: An interactive approach for unraveling dendritic dynamics

Neurons rely on the interplay between dendritic morphology and ion channels to transform synaptic inputs into a sequence of somatic spikes. Detailed biophysical models with active dendrites have been instrumental in exploring this interaction. However, such models can be challenging to understand and validate due to the large number of parameters involved. In this work, we introduce DendroTweaks - a toolbox designed to illuminate how morpho-electric properties map to dendritic events and how these dendritic events shape neuronal output. DendroTweaks features a web-based graphical interface, where users can explore single-cell neuronal models and adjust their morphological and biophysical parameters with real-time visual feedback. In particular, DendroTweaks is tailored to interactive fine-tuning of subcellular properties, such as kinetics and distributions of ion channels, as well as the dynamics and allocation of synaptic inputs. It offers an automated approach for standardization and refinement of voltage-gated ion channel models to make them more comprehensible and reusable. The toolbox allows users to run various experimental protocols and record data from multiple dendritic and somatic locations, thereby enhancing model validation. Finally, it aims to deepen our understanding of which dendritic properties are essential for neuronal input-output transformation. Using this knowledge, one can simplify models through a built-in morphology reduction algorithm and export them for further use in faster, more interpretable networks. With DendroTweaks, users can gain better control and understanding of their models, advancing research on dendritic input-output transformations and their role in network computations.

A layered microcircuit model of somatosensory cortex with three interneuron types and cell-type-specific short-term plasticity

Three major types of GABAergic interneurons, parvalbumin-, somatostatin-, and vasoactive intestinal peptide-expressing (PV, SOM, VIP) cells, play critical but distinct roles in the cortical microcircuitry. Their specific electrophysiology and connectivity shape their inhibitory functions. To study the network dynamics and signal processing specific to these cell types in the cerebral cortex, we developed a multi-layer model incorporating biologically realistic interneuron parameters from rodent somatosensory cortex. The model is fitted to in vivo data on cell-type-specific population firing rates. With a protocol of cell-type-specific stimulation, network responses when activating different neuron types are examined. The model reproduces the experimentally observed inhibitory effects of PV and SOM cells and disinhibitory effect of VIP cells on excitatory cells. We further create a version of the model incorporating cell-type-specific short-term synaptic plasticity (STP). While the ongoing activity with and without STP is similar, STP modulates the responses of Exc, SOM, and VIP cells to cell-type-specific stimulation, presumably by changing the dominant inhibitory pathways. With slight adjustments, the model also reproduces sensory responses of specific interneuron types recorded in vivo. Our model provides predictions on network dynamics involving cell-type-specific short-term plasticity and can serve to explore the computational roles of inhibitory interneurons in sensory functions.

Cell-specific effects of temporal interference stimulation on cortical function

Temporal interference (TI) stimulation is a popular non-invasive neurostimulation technique that utilizes the following salient neural behavior: pure sinusoid (generated in off-target brain regions) appears to cause no stimulation, whereas modulated sinusoid (generated in target brain regions) does. To understand its effects and mechanisms, we examine responses of different cell types, excitatory pyramidal (Pyr) and inhibitory parvalbumin-expressing (PV) neurons, to pure and modulated sinusoids, in intact network as well as in isolation. In intact network, we present data showing that PV neurons are much less likely than Pyr neurons to exhibit TI stimulation. Remarkably, in isolation, our data shows that almost all Pyr neurons stop exhibiting TI stimulation. We conclude that TI stimulation is largely a network phenomenon. Indeed, PV neurons actively inhibit Pyr neurons in the off-target regions due to pure sinusoids (in off-target regions) generating much higher PV firing rates than modulated sinusoids in the target regions. Additionally, we use computational studies to support and extend our experimental observations.

NeoCoMM: A neocortical neuroinspired computational model for the reconstruction and simulation of epileptiform events

BACKGROUND

Understanding the pathophysiological dynamics that underline Interictal Epileptiform Events (IEEs) such as epileptic spikes, spike-and-waves or High-Frequency Oscillations (HFOs) is of major importance in the context of neocortical refractory epilepsy, as it paves the way for the development of novel therapies. Typically, these events are detected in Local Field Potential (LFP) recordings obtained through depth electrodes during pre-surgical investigations. Although essential, the underlying pathophysiological mechanisms for the generation of these epileptic neuromarkers remain unclear. The aim of this paper is to propose a novel neurophysiologically relevant reconstruction of the neocortical microcircuitry in the context of epilepsy. This reconstruction intends to facilitate the analysis of a comprehensive set of parameters encompassing physiological, morphological, and biophysical aspects that directly impact the generation and recording of different IEEs.

METHOD

a novel microscale computational model of an epileptic neocortical column was introduced. This model incorporates the intricate multilayered structure of the cortex and allows for the simulation of realistic interictal epileptic signals. The proposed model was validated through comparisons with real IEEs recorded using intracranial stereo-electroencephalography (SEEG) signals from both humans and animals. Using the model, the user can recreate epileptiform patterns observed in different species (human, rodent, and mouse) and study the intracellular activity associated with these patterns.

RESULTS

Our model allowed us to unravel the relationship between glutamatergic and GABAergic synaptic transmission of the epileptic neural network and the type of generated IEE. Moreover, sensitivity analyses allowed for the exploration of the pathophysiological parameters responsible for the transitions between these events. Finally, the presented modeling framework also provides an Electrode Tissue Model (ETI) that adds realism to the simulated signals and offers the possibility of studying their sensitivity to the electrode characteristics.

CONCLUSION

The model (NeoCoMM) presented in this work can be of great use in different applications since it offers an in silico framework for sensitivity analysis and hypothesis testing. It can also be used as a starting point for more complex studies.

Rosetta stone for the population dynamics of spiking neuron networks

Populations of spiking neuron models have densities of their microscopic variables (e.g., single-cell membrane potentials) whose evolution fully capture the collective dynamics of biological networks, even outside equilibrium. Despite its general applicability, the Fokker-Planck equation governing such evolution is mainly studied within the borders of the linear response theory, although alternative spectral expansion approaches offer some advantages in the study of the out-of-equilibrium dynamics. This is mainly due to the difficulty in computing the state-dependent coefficients of the expanded system of differential equations. Here, we address this issue by deriving analytic expressions for such coefficients by pairing perturbative solutions of the Fokker-Planck approach with their counterparts from the spectral expansion. A tight relationship emerges between several of these coefficients and the Laplace transform of the interspike interval density (i.e., the distribution of first-passage times). "Coefficients" like the current-to-rate gain function, the eigenvalues of the Fokker-Planck operator and its eigenfunctions at the boundaries are derived without resorting to integral expressions. For the leaky integrate-and-fire neurons, the coupling terms between stationary and nonstationary modes are also worked out paving the way to accurately characterize the critical points and the relaxation timescales in networks of interacting populations.

Distribution of calbindin-positive neurons across areas and layers of the marmoset cerebral cortex

The diversity of the mammalian cerebral cortex demands technical approaches to map the spatial distribution of neurons with different biochemical identities. This issue is magnified in the case of the primate cortex, characterized by a large number of areas with distinctive cytoarchitectures. To date, no full map of the distribution of cells expressing a specific protein has been reported for the cortex of any primate. Here we have charted the 3-dimensional distribution of neurons expressing the calcium-binding protein calbindin (CB+ neurons) across the entire marmoset cortex, using a combination of immunohistochemistry, automated cell identification, computerized reconstruction, and cytoarchitecture-aware registration. CB+ neurons formed a heterogeneous population, which together corresponded to 10-20% of the cortical neurons. They occurred in higher proportions in areas corresponding to low hierarchical levels of processing, such as sensory cortices. Although CB+ neurons were concentrated in the supragranular and granular layers, there were clear global trends in their laminar distribution. For example, their relative density in infragranular layers increased with hierarchical level along sensorimotor processing streams, and their density in layer 4 was lower in areas involved in sensorimotor integration, action planning and motor control. These results reveal new quantitative aspects of the cytoarchitectural organization of the primate cortex, and demonstrate an approach to mapping the full distribution of neurochemically distinct cells throughout the brain which is readily applicable to most other mammalian species.

Neuronal firing rate diversity lowers the dimension of population covariability

Populations of neurons produce activity with two central features. First, neuronal responses are very diverse - specific stimuli or behaviors prompt some neurons to emit many action potentials, while other neurons remain relatively silent. Second, the trial-to-trial fluctuations of neuronal response occupy a low dimensional space, owing to significant correlations between the activity of neurons. These two features define the quality of neuronal representation. We link these two aspects of population response using a recurrent circuit model and derive the following relation: the more diverse the firing rates of neurons in a population, the lower the effective dimension of population trial-to-trial covariability. This surprising prediction is tested and validated using simultaneously recorded neuronal populations from numerous brain areas in mice, non-human primates, and in the motor cortex of human participants. Using our relation we present a theory where a more diverse neuronal code leads to better fine discrimination performance from population activity. In line with this theory, we show that neuronal populations across the brain exhibit both more diverse mean responses and lower-dimensional fluctuations when the brain is in more heightened states of information processing. In sum, we present a key organizational principle of neuronal population response that is widely observed across the nervous system and acts to synergistically improve population representation.

Cortical dynamics of perception as trains of coherent gamma oscillations, with the pulvinar as central coordinator

Synchronization of spikes carried by the visual streams is strategic for the proper binding of cortical assemblies, hence for the perception of visual objects as coherent units. Perception of a complex visual scene involves multiple trains of gamma oscillations, coexisting at each stage in visual and associative cortex. Here, we analyze how this synchrony is managed, so that the perception of each visual object can emerge despite this complex interweaving of cortical activations. After a brief review of structural and temporal facts, we analyze the interactions which make the oscillations coherent for the visual elements related to the same object. We continue with the propagation of these gamma oscillations within the sensory chain. The dominant role of the pulvinar and associated reticular thalamic nucleus as cortical coordinator is the common thread running through this step-by-step description. Synchronization mechanisms are analyzed in the context of visual perception, although the present considerations are not limited to this sense. A simple experiment is described, with the aim of assessing the validity of the elements developed here. A first set of results is provided, together with a proposed method to go further in this investigation.

Genetic mechanisms for impaired synaptic plasticity in schizophrenia revealed by computational modeling

Schizophrenia phenotypes are suggestive of impaired cortical plasticity in the disease, but the mechanisms of these deficits are unknown. Genomic association studies have implicated a large number of genes that regulate neuromodulation and plasticity, indicating that the plasticity deficits have a genetic origin. Here, we used biochemically detailed computational modeling of postsynaptic plasticity to investigate how schizophrenia-associated genes regulate long-term potentiation (LTP) and depression (LTD). We combined our model with data from postmortem RNA expression studies (CommonMind gene-expression datasets) to assess the consequences of altered expression of plasticity-regulating genes for the amplitude of LTP and LTD. Our results show that the expression alterations observed post mortem, especially those in the anterior cingulate cortex, lead to impaired protein kinase A (PKA)-pathway-mediated LTP in synapses containing GluR1 receptors. We validated these findings using a genotyped electroencephalogram (EEG) dataset where polygenic risk scores for synaptic and ion channel-encoding genes as well as modulation of visual evoked potentials were determined for 286 healthy controls. Our results provide a possible genetic mechanism for plasticity impairments in schizophrenia, which can lead to improved understanding and, ultimately, treatment of the disorder.

Constructing neural networks with pre-specified dynamics

A main goal in neuroscience is to understand the computations carried out by neural populations that give animals their cognitive skills. Neural network models allow to formulate explicit hypotheses regarding the algorithms instantiated in the dynamics of a neural population, its firing statistics, and the underlying connectivity. Neural networks can be defined by a small set of parameters, carefully chosen to procure specific capabilities, or by a large set of free parameters, fitted with optimization algorithms that minimize a given loss function. In this work we alternatively propose a method to make a detailed adjustment of the network dynamics and firing statistic to better answer questions that link dynamics, structure, and function. Our algorithm-termed generalised Firing-to-Parameter (gFTP)-provides a way to construct binary recurrent neural networks whose dynamics strictly follows a user pre-specified transition graph that details the transitions between population firing states triggered by stimulus presentations. Our main contribution is a procedure that detects when a transition graph is not realisable in terms of a neural network, and makes the necessary modifications in order to obtain a new transition graph that is realisable and preserves all the information encoded in the transitions of the original graph. With a realisable transition graph, gFTP assigns values to the network firing states associated with each node in the graph, and finds the synaptic weight matrices by solving a set of linear separation problems. We test gFTP performance by constructing networks with random dynamics, continuous attractor-like dynamics that encode position in 2-dimensional space, and discrete attractor dynamics. We then show how gFTP can be employed as a tool to explore the link between structure, function, and the algorithms instantiated in the network dynamics.

Heterozygous expression of a Kcnt1 gain-of-function variant has differential effects on SST- and PV-expressing cortical GABAergic neurons

More than twenty recurrent missense gain-of-function (GOF) mutations have been identified in the sodium-activated potassium (KNa) channel gene KCNT1 in patients with severe developmental and epileptic encephalopathies (DEEs), most of which are resistant to current therapies. Defining the neuron types most vulnerable to KCNT1 GOF will advance our understanding of disease mechanisms and provide refined targets for precision therapy efforts. Here, we assessed the effects of heterozygous expression of a Kcnt1 GOF variant (Y777H) on KNa currents and neuronal physiology among cortical glutamatergic and GABAergic neurons in mice, including those expressing vasoactive intestinal polypeptide (VIP), somatostatin (SST), and parvalbumin (PV), to identify and model the pathogenic mechanisms of autosomal dominant KCNT1 GOF variants in DEEs. Although the Kcnt1-Y777H variant had no effects on glutamatergic or VIP neuron function, it increased subthreshold KNa currents in both SST and PV neurons but with opposite effects on neuronal output; SST neurons became hypoexcitable with a higher rheobase current and lower action potential (AP) firing frequency, whereas PV neurons became hyperexcitable with a lower rheobase current and higher AP firing frequency. Further neurophysiological and computational modeling experiments showed that the differential effects of the Y777H variant on SST and PV neurons are not likely due to inherent differences in these neuron types, but to an increased persistent sodium current in PV, but not SST, neurons. The Y777H variant also increased excitatory input onto, and chemical and electrical synaptic connectivity between, SST neurons. Together, these data suggest differential pathogenic mechanisms, both direct and compensatory, contribute to disease phenotypes, and provide a salient example of how a pathogenic ion channel variant can cause opposite functional effects in closely related neuron subtypes due to interactions with other ionic conductances.

A comprehensive data-driven model of cat primary visual cortex

Knowledge integration based on the relationship between structure and function of the neural substrate is one of the main targets of neuroinformatics and data-driven computational modeling. However, the multiplicity of data sources, the diversity of benchmarks, the mixing of observables of different natures, and the necessity of a long-term, systematic approach make such a task challenging. Here we present a first snapshot of a long-term integrative modeling program designed to address this issue in the domain of the visual system: a comprehensive spiking model of cat primary visual cortex. The presented model satisfies an extensive range of anatomical, statistical and functional constraints under a wide range of visual input statistics. In the presence of physiological levels of tonic stochastic bombardment by spontaneous thalamic activity, the modeled cortical reverberations self-generate a sparse asynchronous ongoing activity that quantitatively matches a range of experimentally measured statistics. When integrating feed-forward drive elicited by a high diversity of visual contexts, the simulated network produces a realistic, quantitatively accurate interplay between visually evoked excitatory and inhibitory conductances; contrast-invariant orientation-tuning width; center surround interactions; and stimulus-dependent changes in the precision of the neural code. This integrative model offers insights into how the studied properties interact, contributing to a better understanding of visual cortical dynamics. It provides a basis for future development towards a comprehensive model of low-level perception.

Microinstrumentation for Brain Organoids

Brain organoids are three-dimensional aggregates of self-organized differentiated stem cells that mimic the structure and function of human brain regions. Organoids bridge the gaps between conventional drug screening models such as planar mammalian cell culture, animal studies, and clinical trials. They can revolutionize the fields of developmental biology, neuroscience, toxicology, and computer engineering. Conventional microinstrumentation for conventional cellular engineering, such as planar microfluidic chips; microelectrode arrays (MEAs); and optical, magnetic, and acoustic techniques, has limitations when applied to three-dimensional (3D) organoids, primarily due to their limits with inherently two-dimensional geometry and interfacing. Hence, there is an urgent need to develop new instrumentation compatible with live cell culture techniques and with scalable 3D formats relevant to organoids. This review discusses conventional planar approaches and emerging 3D microinstrumentation necessary for advanced organoid-machine interfaces. Specifically, this article surveys recently developed microinstrumentation, including 3D printed and curved microfluidics, 3D and fast-scan optical techniques, buckling and self-folding MEAs, 3D interfaces for electrochemical measurements, and 3D spatially controllable magnetic and acoustic technologies relevant to two-way information transfer with brain organoids. This article highlights key challenges that must be addressed for robust organoid culture and reliable 3D spatiotemporal information transfer.

Volume electron microscopy analysis of synapses in primary regions of the human cerebral cortex

Functional and structural studies investigating macroscopic connectivity in the human cerebral cortex suggest that high-order associative regions exhibit greater connectivity compared to primary ones. However, the synaptic organization of these brain regions remains unexplored. In the present work, we conducted volume electron microscopy to investigate the synaptic organization of the human brain obtained at autopsy. Specifically, we examined layer III of Brodmann areas 17, 3b, and 4, as representative areas of primary visual, somatosensorial, and motor cortex. Additionally, we conducted comparative analyses with our previous datasets of layer III from temporopolar and anterior cingulate associative cortical regions (Brodmann areas 24, 38, and 21). 9,690 synaptic junctions were 3D reconstructed, showing that certain synaptic characteristics are specific to particular regions. The number of synapses per volume, the proportion of the postsynaptic targets, and the synaptic size may distinguish one region from another, regardless of whether they are associative or primary cortex. By contrast, other synaptic characteristics were common to all analyzed regions, such as the proportion of excitatory and inhibitory synapses, their shapes, their spatial distribution, and a higher proportion of synapses located on dendritic spines. The present results provide further insights into the synaptic organization of the human cerebral cortex.

High-throughput analysis of dendrite and axonal arbors reveals transcriptomic correlates of neuroanatomy

Neuronal anatomy is central to the organization and function of brain cell types. However, anatomical variability within apparently homogeneous populations of cells can obscure such insights. Here, we report large-scale automation of neuronal morphology reconstruction and analysis on a dataset of 813 inhibitory neurons characterized using the Patch-seq method, which enables measurement of multiple properties from individual neurons, including local morphology and transcriptional signature. We demonstrate that these automated reconstructions can be used in the same manner as manual reconstructions to understand the relationship between some, but not all, cellular properties used to define cell types. We uncover gene expression correlates of laminar innervation on multiple transcriptomically defined neuronal subclasses and types. In particular, our results reveal correlates of the variability in Layer 1 (L1) axonal innervation in a transcriptomically defined subpopulation of Martinotti cells in the adult mouse neocortex.

Synthesis of geometrically realistic and watertight neuronal ultrastructure manifolds for in silico modeling

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

Morphological classification of neurons based on Sugeno fuzzy integration and multi-classifier fusion

In order to extract more important morphological features of neuron images and achieve accurate classification of the neuron type, a method is proposed that uses Sugeno fuzzy integral integration of three optimized deep learning models, namely AlexNet, VGG11_bn, and ResNet-50. Firstly, using the pre-trained model of AlexNet and the output layer is fine-tuned to improve the model's performance. Secondly, in the VGG11_bn network, Global Average Pooling (GAP) is adopted to replace the traditional fully connected layer to reduce the number of parameters. Additionally, the generalization ability of the model is improved by transfer learning. Thirdly, the SE(squeeze and excitation) module is added to the ResNet-50 variant ResNeXt-50 to adjust the channel weight and capture the key information of the input data. The GELU activation function is used to better fit the data distribution. Finally, Sugeno fuzzy integral is used to fuse the output of each model to get the final classification result. The experimental results showed that on the Img_raw, Img_resample and Img_XYalign dataset, the accuracy of 4-category classification reached 98.04%, 91.75% and 93.13%, respectively, and the accuracy of 12-category classification reached 97.82%, 85.68% and 87.60%, respectively. The proposed method has good classification performance in the morphological classification of neurons.

Phasor-based analysis of a neuromorphic architecture for microwave sensing

This article presents a design procedure for implementing artificial neural networks (ANNs) using conventional microwave components at the hardware level with potential applications in radar and remote sensing. The main objective is to develop structured hardware design methods for implementing artificial neurons, utilizing microwave devices to create neuromorphic devices compatible with high-frequency electromagnetic waves. The research aims to address the challenge of encoding and modulating information in electromagnetic waves into a format suitable for the neuromorphic device by using frequency-modulated information instead of intensity-modulated information. It also proposes a method for integrating principal component analysis as a dimensionality reduction technique with the implementation of ANNs on a single hardware. As a dummy task, the process outlined here is used to implement an artificial neural network at the hardware level, with a specific emphasis on creating hardware that is capable of performing matrix multiplications in the form of dot products while also being able to extract the resulting data in an interpretable manner. The proposed implementation involves the use of directional couplers to implement weights and sample the resulting signal at specific intervals to obtain the multiplication result.

Parallel Synapses with Transmission Nonlinearities Enhance Neuronal Classification Capacity

Cortical neurons often establish multiple synaptic contacts with the same postsynaptic neuron. To avoid functional redundancy of these parallel synapses, it is crucial that each synapse exhibits distinct computational properties. Here we model the current to the soma contributed by each synapse as a sigmoidal transmission function of its presynaptic input, with learnable parameters such as amplitude, slope, and threshold. We evaluate the classification capacity of a neuron equipped with such nonlinear parallel synapses, and show that with a small number of parallel synapses per axon, it substantially exceeds that of the Perceptron. Furthermore, the number of correctly classified data points can increase superlinearly as the number of presynaptic axons grows. When training with an unrestricted number of parallel synapses, our model neuron can effectively implement an arbitrary aggregate transmission function for each axon, constrained only by monotonicity. Nevertheless, successful learning in the model neuron often requires only a small number of parallel synapses. We also apply these parallel synapses in a feedforward neural network trained to classify MNIST images, and show that they can increase the test accuracy. This demonstrates that multiple nonlinear synapses per input axon can substantially enhance a neuron's computational power.

Modeling the impact of neuromorphological alterations in Down syndrome on fast neural oscillations

Cognitive disorders, including Down syndrome (DS), present significant morphological alterations in neuron architectural complexity. However, the relationship between neuromorphological alterations and impaired brain function is not fully understood. To address this gap, we propose a novel computational model that accounts for the observed cell deformations in DS. The model consists of a cross-sectional layer of the mouse motor cortex, composed of 3000 neurons. The network connectivity is obtained by accounting explicitly for two single-neuron morphological parameters: the mean dendritic tree radius and the spine density in excitatory pyramidal cells. We obtained these values by fitting reconstructed neuron data corresponding to three mouse models: wild-type (WT), transgenic (TgDyrk1A), and trisomic (Ts65Dn). Our findings reveal a dynamic interplay between pyramidal and fast-spiking interneurons leading to the emergence of gamma activity (∼40 Hz). In the DS models this gamma activity is diminished, corroborating experimental observations and validating our computational methodology. We further explore the impact of disrupted excitation-inhibition balance by mimicking the reduction recurrent inhibition present in DS. In this case, gamma power exhibits variable responses as a function of the external input to the network. Finally, we perform a numerical exploration of the morphological parameter space, unveiling the direct influence of each structural parameter on gamma frequency and power. Our research demonstrates a clear link between changes in morphology and the disruption of gamma oscillations in DS. This work underscores the potential of computational modeling to elucidate the relationship between neuron architecture and brain function, and ultimately improve our understanding of cognitive disorders.

Multitask learning of a biophysically-detailed neuron model

The human brain operates at multiple levels, from molecules to circuits, and understanding these complex processes requires integrated research efforts. Simulating biophysically-detailed neuron models is a computationally expensive but effective method for studying local neural circuits. Recent innovations have shown that artificial neural networks (ANNs) can accurately predict the behavior of these detailed models in terms of spikes, electrical potentials, and optical readouts. While these methods have the potential to accelerate large network simulations by several orders of magnitude compared to conventional differential equation based modelling, they currently only predict voltage outputs for the soma or a select few neuron compartments. Our novel approach, based on enhanced state-of-the-art architectures for multitask learning (MTL), allows for the simultaneous prediction of membrane potentials in each compartment of a neuron model, at a speed of up to two orders of magnitude faster than classical simulation methods. By predicting all membrane potentials together, our approach not only allows for comparison of model output with a wider range of experimental recordings (patch-electrode, voltage-sensitive dye imaging), it also provides the first stepping stone towards predicting local field potentials (LFPs), electroencephalogram (EEG) signals, and magnetoencephalography (MEG) signals from ANN-based simulations. While LFP and EEG are an important downstream application, the main focus of this paper lies in predicting dendritic voltages within each compartment to capture the entire electrophysiology of a biophysically-detailed neuron model. It further presents a challenging benchmark for MTL architectures due to the large amount of data involved, the presence of correlations between neighbouring compartments, and the non-Gaussian distribution of membrane potentials.

Tracing nerve fibers with volume electron microscopy to quantitatively analyze brain connectivity

The highly complex structure of the brain requires an approach that can unravel its connectivity. Using volume electron microscopy and a dedicated software we can trace and measure all nerve fibers present within different samples of brain tissue. With this software tool, individual dendrites and axons are traced, obtaining a simplified "skeleton" of each fiber, which is linked to its corresponding synaptic contacts. The result is an intricate meshwork of axons and dendrites interconnected by a cloud of synaptic junctions. To test this methodology, we apply it to the stratum radiatum of the hippocampus and layers 1 and 3 of the somatosensory cortex of the mouse. We find that nerve fibers are densely packed in the neuropil, reaching up to 9 kilometers per cubic mm. We obtain the number of synapses, the number and lengths of dendrites and axons, the linear densities of synapses established by dendrites and axons, and their location on dendritic spines and shafts. The quantitative data obtained through this method enable us to identify subtle traits and differences in the synaptic organization of the samples, which might have been overlooked in a qualitative analysis.

A lightweight data-driven spiking neuronal network model of Drosophila olfactory nervous system with dedicated hardware support

Data-driven spiking neuronal network (SNN) models enable in-silico analysis of the nervous system at the cellular and synaptic level. Therefore, they are a key tool for elucidating the information processing principles of the brain. While extensive research has focused on developing data-driven SNN models for mammalian brains, their complexity poses challenges in achieving precision. Network topology often relies on statistical inference, and the functions of specific brain regions and supporting neuronal activities remain unclear. Additionally, these models demand huge computing facilities and their simulation speed is considerably slower than real-time. Here, we propose a lightweight data-driven SNN model that strikes a balance between simplicity and reproducibility. The model is built using a qualitative modeling approach that can reproduce key dynamics of neuronal activity. We target the Drosophila olfactory nervous system, extracting its network topology from connectome data. The model was successfully implemented on a small entry-level field-programmable gate array and simulated the activity of a network in real-time. In addition, the model reproduced olfactory associative learning, the primary function of the olfactory system, and characteristic spiking activities of different neuron types. In sum, this paper propose a method for building data-driven SNN models from biological data. Our approach reproduces the function and neuronal activities of the nervous system and is lightweight, acceleratable with dedicated hardware, making it scalable to large-scale networks. Therefore, our approach is expected to play an important role in elucidating the brain's information processing at the cellular and synaptic level through an analysis-by-construction approach. In addition, it may be applicable to edge artificial intelligence systems in the future.

Brain organoid as a model to study the role of mitochondria in neurodevelopmental disorders: achievements and weaknesses

Mitochondrial diseases are a group of severe pathologies that cause complex neurodegenerative disorders for which, in most cases, no therapy or treatment is available. These organelles are critical regulators of both neurogenesis and homeostasis of the neurological system. Consequently, mitochondrial damage or dysfunction can occur as a cause or consequence of neurodevelopmental or neurodegenerative diseases. As genetic knowledge of neurodevelopmental disorders advances, associations have been identified between genes that encode mitochondrial proteins and neurological symptoms, such as neuropathy, encephalomyopathy, ataxia, seizures, and developmental delays, among others. Understanding how mitochondrial dysfunction can alter these processes is essential in researching rare diseases. Three-dimensional (3D) cell cultures, which self-assemble to form specialized structures composed of different cell types, represent an accessible manner to model organogenesis and neurodevelopmental disorders. In particular, brain organoids are revolutionizing the study of mitochondrial-based neurological diseases since they are organ-specific and model-generated from a patient's cell, thereby overcoming some of the limitations of traditional animal and cell models. In this review, we have collected which neurological structures and functions recapitulate in the different types of reported brain organoids, focusing on those generated as models of mitochondrial diseases. In addition to advancements in the generation of brain organoids, techniques, and approaches for studying neuronal structures and physiology, drug screening and drug repositioning studies performed in brain organoids with mitochondrial damage and neurodevelopmental disorders have also been reviewed. This scope review will summarize the evidence on limitations in studying the function and dynamics of mitochondria in brain organoids.

The Long Journey from Animal Electricity to the Discovery of Ion Channels and the Modelling of the Human Brain

This retrospective begins with Galvani's experiments on frogs at the end of the 18th century and his discovery of 'animal electricity'. It goes on to illustrate the numerous contributions to the field of physical chemistry in the second half of the 19th century (Nernst's equilibrium potential, based on the work of Wilhelm Ostwald, Max Planck's ion electrodiffusion, Einstein's studies of Brownian motion) which led Bernstein to propose his membrane theory in the early 1900s as an explanation of Galvani's findings and cell excitability. These processes were fully elucidated by Hodgkin and Huxley in 1952 who detailed the ionic basis of resting and action potentials, but without addressing the question of where these ions passed. The emerging question of the existence of ion channels, widely debated over the next two decades, was finally accepted and, a decade later, many of them began to be cloned. This led to the possibility of modelling the activity of individual neurons in the brain and then that of simple circuits. Taking advantage of the remarkable advances in computer science in the new millennium, together with a much deeper understanding of brain architecture, more ambitious scientific goals were dreamed of to understand the brain and how it works. The retrospective concludes by reviewing the main efforts in this direction, namely the construction of a digital brain, an in silico copy of the brain that would run on supercomputers and behave just like a real brain.

A Multimodal Fitting Approach to Construct Single-Neuron Models With Patch Clamp and High-Density Microelectrode Arrays

In computational neuroscience, multicompartment models are among the most biophysically realistic representations of single neurons. Constructing such models usually involves the use of the patch-clamp technique to record somatic voltage signals under different experimental conditions. The experimental data are then used to fit the many parameters of the model. While patching of the soma is currently the gold-standard approach to build multicompartment models, several studies have also evidenced a richness of dynamics in dendritic and axonal sections. Recording from the soma alone makes it hard to observe and correctly parameterize the activity of nonsomatic compartments. In order to provide a richer set of data as input to multicompartment models, we here investigate the combination of somatic patch-clamp recordings with recordings of high-density microelectrode arrays (HD-MEAs). HD-MEAs enable the observation of extracellular potentials and neural activity of neuronal compartments at subcellular resolution. In this work, we introduce a novel framework to combine patch-clamp and HD-MEA data to construct multicompartment models. We first validate our method on a ground-truth model with known parameters and show that the use of features extracted from extracellular signals, in addition to intracellular ones, yields models enabling better fits than using intracellular features alone. We also demonstrate our procedure using experimental data by constructing cell models from in vitro cell cultures. The proposed multimodal fitting procedure has the potential to augment the modeling efforts of the computational neuroscience community and provide the field with neuronal models that are more realistic and can be better validated.

An integrative, multiscale view on neural theories of consciousness

How is conscious experience related to material brain processes? A variety of theories aiming to answer this age-old question have emerged from the recent surge in consciousness research, and some are now hotly debated. Although most researchers have so far focused on the development and validation of their preferred theory in relative isolation, this article, written by a group of scientists representing different theories, takes an alternative approach. Noting that various theories often try to explain different aspects or mechanistic levels of consciousness, we argue that the theories do not necessarily contradict each other. Instead, several of them may converge on fundamental neuronal mechanisms and be partly compatible and complementary, so that multiple theories can simultaneously contribute to our understanding. Here, we consider unifying, integration-oriented approaches that have so far been largely neglected, seeking to combine valuable elements from various theories.

Multiscale neuro-inspired models for interpretation of EEG signals in patients with epilepsy

OBJECTIVE

The aim is to gain insight into the pathophysiological mechanisms underlying interictal epileptiform discharges observed in electroencephalographic (EEG) and stereo-EEG (SEEG, depth electrodes) recordings performed during pre-surgical evaluation of patients with drug-resistant epilepsy.

METHODS

We developed novel neuro-inspired computational models of the human cerebral cortex at three different levels of description: i) microscale (detailed neuron models), ii) mesoscale (neuronal mass models) and iii) macroscale (whole brain models). Although conceptually different, micro- and mesoscale models share some similar features, such as the typology of neurons (pyramidal cells and three types of interneurons), their spatial arrangement in cortical layers, and their synaptic connectivity (excitatory and inhibitory). The whole brain model consists of a large-scale network of interconnected neuronal masses, with connectivity based on the human connectome.

RESULTS

For these three levels of description, the fine-tuning of free parameters and the quantitative comparison with real data allowed us to reproduce interictal epileptiform discharges with a high degree of fidelity and to formulate hypotheses about the cell- and network-related mechanisms underlying the generation of fast ripples and SEEG-recorded epileptic spikes and spike-waves.

CONCLUSIONS

The proposed models provide valuable insights into the pathophysiological mechanisms underlying the generation of epileptic events. The knowledge gained from these models effectively complements the clinical analysis of SEEG data collected during the evaluation of patients with epilepsy.

SIGNIFICANCE

These models are likely to play a key role in the mechanistic interpretation of epileptiform activity.

Top-down modulation in canonical cortical circuits with short-term plasticity

Cortical dynamics and computations are strongly influenced by diverse GABAergic interneurons, including those expressing parvalbumin (PV), somatostatin (SST), and vasoactive intestinal peptide (VIP). Together with excitatory (E) neurons, they form a canonical microcircuit and exhibit counterintuitive nonlinear phenomena. One instance of such phenomena is response reversal, whereby SST neurons show opposite responses to top-down modulation via VIP depending on the presence of bottom-up sensory input, indicating that the network may function in different regimes under different stimulation conditions. Combining analytical and computational approaches, we demonstrate that model networks with multiple interneuron subtypes and experimentally identified short-term plasticity mechanisms can implement response reversal. Surprisingly, despite not directly affecting SST and VIP activity, PV-to-E short-term depression has a decisive impact on SST response reversal. We show how response reversal relates to inhibition stabilization and the paradoxical effect in the presence of several short-term plasticity mechanisms demonstrating that response reversal coincides with a change in the indispensability of SST for network stabilization. In summary, our work suggests a role of short-term plasticity mechanisms in generating nonlinear phenomena in networks with multiple interneuron subtypes and makes several experimentally testable predictions.

Physical and functional convergence of the autism risk genes Scn2a and Ank2 in neocortical pyramidal cell dendrites

Dysfunction in sodium channels and their ankyrin scaffolding partners have both been implicated in neurodevelopmental disorders, including autism spectrum disorder (ASD). In particular, the genes SCN2A, which encodes the sodium channel NaV1.2, and ANK2, which encodes ankyrin-B, have strong ASD association. Recent studies indicate that ASD-associated haploinsufficiency in Scn2a impairs dendritic excitability and synaptic function in neocortical pyramidal cells, but how NaV1.2 is anchored within dendritic regions is unknown. Here, we show that ankyrin-B is essential for scaffolding NaV1.2 to the dendritic membrane of mouse neocortical neurons and that haploinsufficiency of Ank2 phenocopies intrinsic dendritic excitability and synaptic deficits observed in Scn2a+/- conditions. These results establish a direct, convergent link between two major ASD risk genes and reinforce an emerging framework suggesting that neocortical pyramidal cell dendritic dysfunction can contribute to neurodevelopmental disorder pathophysiology.

The emergence of multiscale connectomics-based approaches in stroke recovery

Stroke is a leading cause of adult disability. Understanding stroke damage and recovery requires deciphering changes in complex brain networks across different spatiotemporal scales. While recent developments in brain readout technologies and progress in complex network modeling have revolutionized current understanding of the effects of stroke on brain networks at a macroscale, reorganization of smaller scale brain networks remains incompletely understood. In this review, we use a conceptual framework of graph theory to define brain networks from nano- to macroscales. Highlighting stroke-related brain connectivity studies at multiple scales, we argue that multiscale connectomics-based approaches may provide new routes to better evaluate brain structural and functional remapping after stroke and during recovery.

Data-driven multiscale computational models of cortical and subcortical regions

Data-driven computational models of neurons, synapses, microcircuits, and mesocircuits have become essential tools in modern brain research. The goal of these multiscale models is to integrate and synthesize information from different levels of brain organization, from cellular properties, dendritic excitability, and synaptic dynamics to microcircuits, mesocircuits, and ultimately behavior. This article surveys recent advances in the genesis of data-driven computational models of mammalian neural networks in cortical and subcortical areas. I discuss the challenges and opportunities in developing data-driven multiscale models, including the need for interdisciplinary collaborations, the importance of model validation and comparison, and the potential impact on basic and translational neuroscience research. Finally, I highlight future directions and emerging technologies that will enable more comprehensive and predictive data-driven models of brain function and dysfunction.

Modeling Instantaneous Firing Rate of Deep Brain Stimulation Target Neuronal Ensembles in the Basal Ganglia and Thalamus

OBJECTIVE

Deep brain stimulation (DBS) is an effective treatment for movement disorders, including Parkinson disease and essential tremor. However, the underlying mechanisms of DBS remain elusive. Despite the capability of existing models in interpreting experimental data qualitatively, there are very few unified computational models that quantitatively capture the dynamics of the neuronal activity of varying stimulated nuclei-including subthalamic nucleus (STN), substantia nigra pars reticulata (SNr), and ventral intermediate nucleus (Vim)-across different DBS frequencies.

MATERIALS AND METHODS

Both synthetic and experimental data were used in the model fitting; the synthetic data were generated by an established spiking neuron model that was reported in our previous work, and the experimental data were provided using single-unit microelectrode recordings (MERs) during DBS (microelectrode stimulation). Based on these data, we developed a novel mathematical model to represent the firing rate of neurons receiving DBS, including neurons in STN, SNr, and Vim-across different DBS frequencies. In our model, the DBS pulses were filtered through a synapse model and a nonlinear transfer function to formulate the firing rate variability. For each DBS-targeted nucleus, we fitted a single set of optimal model parameters consistent across varying DBS frequencies.

RESULTS

Our model accurately reproduced the firing rates observed and calculated from both synthetic and experimental data. The optimal model parameters were consistent across different DBS frequencies.

CONCLUSIONS

The result of our model fitting was in agreement with experimental single-unit MER data during DBS. Reproducing neuronal firing rates of different nuclei of the basal ganglia and thalamus during DBS can be helpful to further understand the mechanisms of DBS and to potentially optimize stimulation parameters based on their actual effects on neuronal activity.

Inference of network connectivity from temporally binned spike trains

BACKGROUND

Processing neural activity to reconstruct network connectivity is a central focus of neuroscience, yet the spatiotemporal requisites of biological nervous systems are challenging for current neuronal sensing modalities. Consequently, methods that leverage limited data to successfully infer synaptic connections, predict activity at single unit resolution, and decipher their effect on whole systems, can uncover critical information about neural processing. Despite the emergence of powerful methods for inferring connectivity, network reconstruction based on temporally subsampled data remains insufficiently unexplored.

NEW METHOD

We infer synaptic weights by processing firing rates within variable time bins for a heterogeneous feed-forward network of excitatory, inhibitory, and unconnected units. We assess classification and optimize model parameters for postsynaptic spike train reconstruction. We test our method on a physiological network of leaky integrate-and-fire neurons displaying bursting patterns and assess prediction of postsynaptic activity from microelectrode array data.

RESULTS

Results reveal parameters for improved prediction and performance and suggest that lower resolution data and limited access to neurons can be preferred.

COMPARISON WITH EXISTING METHOD(S)

Recent computational methods demonstrate highly improved reconstruction of connectivity from networks of parallel spike trains by considering spike lag, time-varying firing rates, and other underlying dynamics. However, these methods insufficiently explore temporal subsampling representative of novel data types.

CONCLUSIONS

We provide a framework for reverse engineering neural networks from data with limited temporal quality, describing optimal parameters for each bin size, which can be further improved using non-linear methods and applied to more complicated readouts and connectivity distributions in multiple brain circuits.

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.

Biophysical modulation and robustness of itinerant complexity in neuronal networks

Transient synchronization of bursting activity in neuronal networks, which occurs in patterns of metastable itinerant phase relationships between neurons, is a notable feature of network dynamics observed in vivo. However, the mechanisms that contribute to this dynamical complexity in neuronal circuits are not well understood. Local circuits in cortical regions consist of populations of neurons with diverse intrinsic oscillatory features. In this study, we numerically show that the phenomenon of transient synchronization, also referred to as metastability, can emerge in an inhibitory neuronal population when the neurons' intrinsic fast-spiking dynamics are appropriately modulated by slower inputs from an excitatory neuronal population. Using a compact model of a mesoscopic-scale network consisting of excitatory pyramidal and inhibitory fast-spiking neurons, our work demonstrates a relationship between the frequency of pyramidal population oscillations and the features of emergent metastability in the inhibitory population. In addition, we introduce a method to characterize collective transitions in metastable networks. Finally, we discuss potential applications of this study in mechanistically understanding cortical network dynamics.

Modeling circuit mechanisms of opposing cortical responses to visual flow perturbations

In an ever-changing visual world, animals' survival depends on their ability to perceive and respond to rapidly changing motion cues. The primary visual cortex (V1) is at the forefront of this sensory processing, orchestrating neural responses to perturbations in visual flow. However, the underlying neural mechanisms that lead to distinct cortical responses to such perturbations remain enigmatic. In this study, our objective was to uncover the neural dynamics that govern V1 neurons' responses to visual flow perturbations using a biologically realistic computational model. By subjecting the model to sudden changes in visual input, we observed opposing cortical responses in excitatory layer 2/3 (L2/3) neurons, namely, depolarizing and hyperpolarizing responses. We found that this segregation was primarily driven by the competition between external visual input and recurrent inhibition, particularly within L2/3 and L4. This division was not observed in excitatory L5/6 neurons, suggesting a more prominent role for inhibitory mechanisms in the visual processing of the upper cortical layers. Our findings share similarities with recent experimental studies focusing on the opposing influence of top-down and bottom-up inputs in the mouse primary visual cortex during visual flow perturbations.

Neural mechanisms of the temporal response of cortical neurons to intracortical microstimulation

BACKGROUND

Intracortical microstimulation (ICMS) is used to map neuronal circuitry in the brain and restore lost sensory function, including vision, hearing, and somatosensation. The temporal response of cortical neurons to single pulse ICMS is remarkably stereotyped and comprises short latency excitation followed by prolonged inhibition and, in some cases, rebound excitation. However, the neural origin of the different response components to ICMS are poorly understood, and the interactions between the three response components during trains of ICMS pulses remains unclear.

OBJECTIVE

We used computational modeling to determine the mechanisms contributing to the temporal response to ICMS in model cortical neurons.

METHODS

We implemented a biophysically based computational model of a cortical column comprising neurons with realistic morphology and synapses and quantified the temporal response of cortical neurons to different ICMS protocols. We characterized the temporal responses to single pulse ICMS across stimulation intensities and inhibitory (GABA-B/GABA-A) synaptic strengths. To probe interactions between response components, we quantified the response to paired pulse ICMS at different inter-pulse intervals and the response to short trains at different stimulation frequencies. Finally, we evaluated the performance of biomimetic ICMS trains in evoking sustained neural responses.

RESULTS

Single pulse ICMS evoked short latency excitation followed by a period of inhibition, but model neurons did not exhibit post-inhibitory excitation. The strength of short latency excitation increased and the duration of inhibition increased with increased stimulation amplitude. Prolonged inhibition resulted from both after-hyperpolarization currents and GABA-B synaptic transmission. During the paired pulse protocol, the strength of short latency excitation evoked by a test pulse decreased marginally compared to those evoked by a single pulse for interpulse intervals (IPI) < 100 m s. Further, the duration of inhibition evoked by the test pulse was prolonged compared to single pulse for IPIs <50 m s and was not predicted by linear superposition of individual inhibitory responses. For IPIs>50 m s, the duration of inhibition evoked by the test pulse was comparable to those evoked by a single pulse. Short ICMS trains evoked repetitive excitatory responses against a background of inhibition. However, the strength of the repetitive excitatory response declined during ICMS at higher frequencies. Further, the duration of inhibition at the cessation of ICMS at higher frequencies was prolonged compared to the duration following a single pulse. Biomimetic pulse trains evoked comparable neural response between the onset and offset phases despite the presence of stimulation induced inhibition.

CONCLUSIONS

The cortical column model replicated the short latency excitation and long-lasting inhibitory components of the stereotyped neural response documented in experimental studies of ICMS. Both cellular and synaptic mechanisms influenced the response components generated by ICMS. The non-linear interactions between response components resulted in dynamic ICMS-evoked neural activity and may play an important role in mediating the ICMS-induced precepts.

Cortical cell assemblies and their underlying connectivity: An in silico study

Recent developments in experimental techniques have enabled simultaneous recordings from thousands of neurons, enabling the study of functional cell assemblies. However, determining the patterns of synaptic connectivity giving rise to these assemblies remains challenging. To address this, we developed a complementary, simulation-based approach, using a detailed, large-scale cortical network model. Using a combination of established methods we detected functional cell assemblies from the stimulus-evoked spiking activity of 186,665 neurons. We studied how the structure of synaptic connectivity underlies assembly composition, quantifying the effects of thalamic innervation, recurrent connectivity, and the spatial arrangement of synapses on dendrites. We determined that these features reduce up to 30%, 22%, and 10% of the uncertainty of a neuron belonging to an assembly. The detected assemblies were activated in a stimulus-specific sequence and were grouped based on their position in the sequence. We found that the different groups were affected to different degrees by the structural features we considered. Additionally, connectivity was more predictive of assembly membership if its direction aligned with the temporal order of assembly activation, if it originated from strongly interconnected populations, and if synapses clustered on dendritic branches. In summary, reversing Hebb's postulate, we showed how cells that are wired together, fire together, quantifying how connectivity patterns interact to shape the emergence of assemblies. This includes a qualitative aspect of connectivity: not just the amount, but also the local structure matters; from the subcellular level in the form of dendritic clustering to the presence of specific network motifs.

Synaptic and dendritic architecture of different types of hippocampal somatostatin interneurons

GABAergic inhibitory neurons fundamentally shape the activity and plasticity of cortical circuits. A major subset of these neurons contains somatostatin (SOM); these cells play crucial roles in neuroplasticity, learning, and memory in many brain areas including the hippocampus, and are implicated in several neuropsychiatric diseases and neurodegenerative disorders. Two main types of SOM-containing cells in area CA1 of the hippocampus are oriens-lacunosum-moleculare (OLM) cells and hippocampo-septal (HS) cells. These cell types show many similarities in their soma-dendritic architecture, but they have different axonal targets, display different activity patterns in vivo, and are thought to have distinct network functions. However, a complete understanding of the functional roles of these interneurons requires a precise description of their intrinsic computational properties and their synaptic interactions. In the current study we generated, analyzed, and make available several key data sets that enable a quantitative comparison of various anatomical and physiological properties of OLM and HS cells in mouse. The data set includes detailed scanning electron microscopy (SEM)-based 3D reconstructions of OLM and HS cells along with their excitatory and inhibitory synaptic inputs. Combining this core data set with other anatomical data, patch-clamp electrophysiology, and compartmental modeling, we examined the precise morphological structure, inputs, outputs, and basic physiological properties of these cells. Our results highlight key differences between OLM and HS cells, particularly regarding the density and distribution of their synaptic inputs and mitochondria. For example, we estimated that an OLM cell receives about 8,400, whereas an HS cell about 15,600 synaptic inputs, about 16% of which are GABAergic. Our data and models provide insight into the possible basis of the different functionality of OLM and HS cell types and supply essential information for more detailed functional models of these neurons and the hippocampal network.

Finite-element modeling of neuromodulation via controlled delivery of potassium ions using conductive polymer-coated microelectrodes

Objective. The controlled delivery of potassium is an interesting neuromodulation modality, being potassium ions involved in shaping neuron excitability, synaptic transmission, network synchronization, and playing a key role in pathological conditions like epilepsy and spreading depression. Despite many successful examples of pre-clinical devices able to influence the extracellular potassium concentration, computational frameworks capturing the corresponding impact on neuronal activity are still missing.Approach. We present a finite-element model describing a PEDOT:PSS-coated microelectrode (herein, simplyionic actuator) able to release potassium and thus modulate the activity of a cortical neuron in anin-vitro-like setting. The dynamics of ions in the ionic actuator, the neural membrane, and the cellular fluids are solved self-consistently.Main results. We showcase the capability of the model to describe on a physical basis the modulation of the intrinsic excitability of the cell and of the synaptic transmission following the electro-ionic stimulation produced by the actuator. We consider three case studies for the ionic actuator with different levels of selectivity to potassium: ideal selectivity, no selectivity, and selectivity achieved by embedding ionophores in the polymer.Significance. This work is the first step toward a comprehensive computational framework aimed to investigate novel neuromodulation devices targeting specific ionic species, as well as to optimize their design and performance, in terms of the induced modulation of neural activity.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Frequency-dependent membrane polarization across neocortical cell types and subcellular elements by transcranial alternating current stimulation

Objective.Transcranial alternating current stimulation (tACS) is a non-invasive brain stimulation technique that directly interacts with ongoing brain oscillations in a frequency-dependent manner. However, it remains largely unclear how the cellular effects of tACS vary between cell types and subcellular elements.Approach.In this study, we use a set of morphologically realistic models of neocortical neurons to simulate the cellular response to uniform oscillating electric fields (EFs). We systematically characterize the membrane polarization in the soma, axons, and dendrites with varying field directions, intensities, and frequencies.Main results.Pyramidal cells are more sensitive to axial EF that is roughly parallel to the cortical column, while interneurons are sensitive to axial EF and transverse EF that is tangent to the cortical surface. Membrane polarization in each subcellular element increases linearly with EF intensity, and its slope, i.e. polarization length, highly depends on the stimulation frequency. At each frequency, pyramidal cells are more polarized than interneurons. Axons usually experience the highest polarization, followed by the dendrites and soma. Moreover, a visible frequency resonance presents in the apical dendrites of pyramidal cells, while the other subcellular elements primarily exhibit low-pass filtering properties. In contrast, each subcellular element of interneurons exhibits complex frequency-dependent polarization. Polarization phase in each subcellular element of cortical neurons lags that of field and exhibits high-pass filtering properties. These results demonstrate that the membrane polarization is not only frequency-dependent, but also cell type- and subcellular element-specific. Through relating effective length and ion mechanism with polarization, we emphasize the crucial role of cell morphology and biophysics in determining the frequency-dependent membrane polarization.Significance.Our findings highlight the diverse polarization patterns across cell types as well as subcellular elements, which provide some insights into the tACS cellular effects and should be considered when understanding the neural spiking activity by tACS.

Discovering optimal features for neuron-type identification from extracellular recordings

Advancements in multichannel recordings of single-unit activity (SUA) in vivo present an opportunity to discover novel features of spatially-varying extracellularly-recorded action potentials (EAPs) that are useful for identifying neuron-types. Traditional approaches to classifying neuron-types often rely on computing EAP waveform features based on conventions of single-channel recordings and thus inherit their limitations. However, spatiotemporal EAP waveforms are the product of signals from underlying current sources being mixed within the extracellular space. We introduce a machine learning approach to demix the underlying sources of spatiotemporal EAP waveforms. Using biophysically realistic computational models, we simulate EAP waveforms and characterize them by the relative prevalence of these sources, which we use as features for identifying the neuron-types corresponding to recorded single units. These EAP sources have distinct spatial and multi-resolution temporal patterns that are robust to various sampling biases. EAP sources also are shared across many neuron-types, are predictive of gross morphological features, and expose underlying morphological domains. We then organize known neuron-types into a hierarchy of latent morpho-electrophysiological types based on differences in the source prevalences, which provides a multi-level classification scheme. We validate the robustness, accuracy, and interpretations of our demixing approach by analyzing simulated EAPs from morphologically detailed models with classification and clustering methods. This simulation-based approach provides a machine learning strategy for neuron-type identification.

Connectivity Matrix Seriation via Relaxation

Volume electron microscopy together with computer-based image analysis are yielding neural circuit diagrams of ever larger regions of the brain. These datasets are usually represented in a cell-to-cell connectivity matrix and contain important information about prevalent circuit motifs allowing to directly test various theories on the computation in that brain structure. Of particular interest are the detection of cell assemblies and the quantification of feedback, which can profoundly change circuit properties. While the ordering of cells along the rows and columns doesn't change the connectivity, it can make special connectivity patterns recognizable. For example, ordering the cells along the flow of information, feedback and feedforward connections are segregated above and below the main matrix diagonal, respectively. Different algorithms are used to renumber matrices such as to minimize a given cost function, but either their performance becomes unsatisfying at a given size of the circuit or the CPU time needed to compute them scales in an unfavorable way with increasing number of neurons. Based on previous ideas, I describe an algorithm which is effective in matrix reordering with respect to both its performance as well as to its scaling in computing time. Rather than trying to reorder the matrix in discrete steps, the algorithm transiently relaxes the integer program by assigning a real-valued parameter to each cell describing its location on a continuous axis ('smooth-index') and finds the parameter set that minimizes the cost. I find that the smooth-index algorithm outperforms all algorithms I compared it to, including those based on topological sorting.