Reconstruction and Simulation of Neocortical Microcircuitry

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.

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.

Enabling Electric Field Modeling of Microscopically Realistic Brain

Across all domains of brain stimulation (neuromodulation), conventional analysis of neuron activation involves two discrete steps: i) prediction of macroscopic electric field, ignoring presence of cells and; ii) prediction of cell activation from tissue electric fields. The first step assumes that current flow is not distorted by the dense tortuous network of cell structures. The deficiencies of this assumption have long been recognized, but - except for trivial geometries - ignored, because it presented intractable computation hurdles. This study introduces a novel approach for analyzing electric fields within a microscopically realistic brain volume. Our pipeline overcomes the technical intractability that prevented such analysis while also showing significant implications for brain stimulation. Contrary to the standard finite element method (FEM), we suggest using a nested iterative boundary element method (BEM) coupled with the fast multipole method (FMM). This approach allows for solving problems with multiple length scales more efficiently. A target application is a subvolume of the L2/3 P36 mouse primary visual cortex containing approximately 400 detailed densely packed neuronal cells at a resolution of 100 nm, which is obtained from scanning electron microscopy data. Our immediate result is a reduction of the stimulation field strength necessary for neuron activation by a factor of 0.85-0.55 (by 15%-45%) as compared to macroscopic predictions. This is in line with modern experimental data stating that existing macroscopic theories substantially overestimate electric field levels necessary for brain stimulation.

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.

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.

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.

A Modular Implementation to Handle and Benchmark Drift Correction for High-Density Extracellular Recordings

High-density neural devices are now offering the possibility to record from neuronal populations in vivo at unprecedented scale. However, the mechanical drifts often observed in these recordings are currently a major issue for "spike sorting," an essential analysis step to identify the activity of single neurons from extracellular signals. Although several strategies have been proposed to compensate for such drifts, the lack of proper benchmarks makes it hard to assess the quality and effectiveness of motion correction. In this paper, we present a benchmark study to precisely and quantitatively evaluate the performance of several state-of-the-art motion correction algorithms introduced in the literature. Using simulated recordings with induced drifts, we dissect the origins of the errors performed while applying a motion correction algorithm as a preprocessing step in the spike sorting pipeline. We show how important it is to properly estimate the positions of the neurons from extracellular traces in order to correctly estimate the probe motion, compare several interpolation procedures, and highlight what are the current limits for motion correction approaches.

On the physiological and structural contributors to the overall balance of excitation and inhibition in local cortical networks

Overall balance of excitation and inhibition in cortical networks is central to their functionality and normal operation. Such orchestrated co-evolution of excitation and inhibition is established through convoluted local interactions between neurons, which are organized by specific network connectivity structures and are dynamically controlled by modulating synaptic activities. Therefore, identifying how such structural and physiological factors contribute to establishment of overall balance of excitation and inhibition is crucial in understanding the homeostatic plasticity mechanisms that regulate the balance. We use biologically plausible mathematical models to extensively study the effects of multiple key factors on overall balance of a network. We characterize a network's baseline balanced state by certain functional properties, and demonstrate how variations in physiological and structural parameters of the network deviate this balance and, in particular, result in transitions in spontaneous activity of the network to high-amplitude slow oscillatory regimes. We show that deviations from the reference balanced state can be continuously quantified by measuring the ratio of mean excitatory to mean inhibitory synaptic conductances in the network. Our results suggest that the commonly observed ratio of the number of inhibitory to the number of excitatory neurons in local cortical networks is almost optimal for their stability and excitability. Moreover, the values of inhibitory synaptic decay time constants and density of inhibitory-to-inhibitory network connectivity are critical to overall balance and stability of cortical networks. However, network stability in our results is sufficiently robust against modulations of synaptic quantal conductances, as required by their role in learning and memory. Our study based on extensive bifurcation analyses thus reveal the functional optimality and criticality of structural and physiological parameters in establishing the baseline operating state of local cortical networks.

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.

Topological data analysis of the firings of a network of stochastic spiking neurons

Topological data analysis is becoming more and more popular in recent years. It has found various applications in many different fields, for its convenience in analyzing and understanding the structure and dynamic of complex systems. We used topological data analysis to analyze the firings of a network of stochastic spiking neurons, which can be in a sub-critical, critical, or super-critical state depending on the value of the control parameter. We calculated several topological features regarding Betti curves and then analyzed the behaviors of these features, using them as inputs for machine learning to discriminate the three states of the network.

Accelerating Neural ODEs Using Model Order Reduction

Embedding nonlinear dynamical systems into artificial neural networks is a powerful new formalism for machine learning. By parameterizing ordinary differential equations (ODEs) as neural network layers, these Neural ODEs are memory-efficient to train, process time series naturally, and incorporate knowledge of physical systems into deep learning (DL) models. However, the practical applications of Neural ODEs are limited due to long inference times because the outputs of the embedded ODE layers are computed numerically with differential equation solvers that can be computationally demanding. Here, we show that mathematical model order reduction (MOR) methods can be used for compressing and accelerating Neural ODEs by accurately simulating the continuous nonlinear dynamics in low-dimensional subspaces. We implement our novel compression method by developing Neural ODEs that integrate the necessary subspace-projection and interpolation operations as layers of the neural network. We validate our approach by comparing it to neuron pruning and singular value decomposition (SVD)-based weight truncation methods from the literature in image and time-series classification tasks. The methods are evaluated by acceleration versus accuracy when adjusting the level of compression. On this spectrum, we achieve a favorable balance over existing methods by using MOR when compressing a convolutional Neural ODE. In compressing a recurrent Neural ODE, SVD-based weight truncation yields good performance. Based on our results, our integration of MOR with Neural ODEs can facilitate efficient, dynamical system-driven DL in resource-constrained applications.

Evolution of Consciousness

The natural evolution of consciousness in different animal species mandates that conscious experiences are causally potent in order to confer any advantage in the struggle for survival. Any endeavor to construct a physical theory of consciousness based on emergence within the framework of classical physics, however, leads to causally impotent conscious experiences in direct contradiction to evolutionary theory since epiphenomenal consciousness cannot evolve through natural selection. Here, we review recent theoretical advances in describing sentience and free will as fundamental aspects of reality granted by quantum physical laws. Modern quantum information theory considers quantum states as a physical resource that endows quantum systems with the capacity to perform physical tasks that are classically impossible. Reductive identification of conscious experiences with the quantum information comprised in quantum brain states allows for causally potent consciousness that is capable of performing genuine choices for future courses of physical action. The consequent evolution of brain cortical networks contributes to increased computational power, memory capacity, and cognitive intelligence of the living organisms.

Data-driven multiscale model of macaque auditory thalamocortical circuits reproduces in vivo dynamics

We developed a detailed model of macaque auditory thalamocortical circuits, including primary auditory cortex (A1), medial geniculate body (MGB), and thalamic reticular nucleus, utilizing the NEURON simulator and NetPyNE tool. The A1 model simulates a cortical column with over 12,000 neurons and 25 million synapses, incorporating data on cell-type-specific neuron densities, morphology, and connectivity across six cortical layers. It is reciprocally connected to the MGB thalamus, which includes interneurons and core and matrix-layer-specific projections to A1. The model simulates multiscale measures, including physiological firing rates, local field potentials (LFPs), current source densities (CSDs), and electroencephalography (EEG) signals. Laminar CSD patterns, during spontaneous activity and in response to broadband noise stimulus trains, mirror experimental findings. Physiological oscillations emerge spontaneously across frequency bands comparable to those recorded in vivo. We elucidate population-specific contributions to observed oscillation events and relate them to firing and presynaptic input patterns. The model offers a quantitative theoretical framework to integrate and interpret experimental data and predict its underlying cellular and circuit mechanisms.

Cortical origin of theta error signals

A multi-scale approach elucidated the origin of the error-related-negativity (ERN), with its associated theta-rhythm, and the post-error-positivity (Pe) in macaque supplementary eye field (SEF). Using biophysical modeling, synaptic inputs to a subpopulation of layer-3 (L3) and layer-5 (L5) pyramidal cells (PCs) were optimized to reproduce error-related spiking modulation and inter-spike intervals. The intrinsic dynamics of dendrites in L5 but not L3 error PCs generate theta rhythmicity with random phases. Saccades synchronized the phases of the theta-rhythm, which was magnified on errors. Contributions from error PCs to the laminar current source density (CSD) observed in SEF were negligible and could not explain the observed association between error-related spiking modulation in L3 PCs and scalp-EEG. CSD from recorded laminar field potentials in SEF was comprised of multipolar components, with monopoles indicating strong electro-diffusion, dendritic/axonal electrotonic current leakage outside SEF, or violations of the model assumptions. Our results also demonstrate the involvement of secondary cortical regions, in addition to SEF, particularly for the later Pe component. The dipolar component from the observed CSD paralleled the ERN dynamics, while the quadrupolar component paralleled the Pe. These results provide the most advanced explanation to date of the cellular mechanisms generating the ERN.

Connecting single-cell transcriptomes to projectomes in mouse visual cortex

The mammalian brain is composed of diverse neuron types that play different functional roles. Recent single-cell RNA sequencing approaches have led to a whole brain taxonomy of transcriptomically-defined cell types, yet cell type definitions that include multiple cellular properties can offer additional insights into a neuron's role in brain circuits. While the Patch-seq method can investigate how transcriptomic properties relate to the local morphological and electrophysiological properties of cell types, linking transcriptomic identities to long-range projections is a major unresolved challenge. To address this, we collected coordinated Patch-seq and whole brain morphology data sets of excitatory neurons in mouse visual cortex. From the Patch-seq data, we defined 16 integrated morpho-electric-transcriptomic (MET)-types; in parallel, we reconstructed the complete morphologies of 300 neurons. We unified the two data sets with a multi-step classifier, to integrate cell type assignments and interrogate cross-modality relationships. We find that transcriptomic variations within and across MET-types correspond with morphological and electrophysiological phenotypes. In addition, this variation, along with the anatomical location of the cell, can be used to predict the projection targets of individual neurons. We also shed new light on infragranular cell types and circuits, including cell-type-specific, interhemispheric projections. With this approach, we establish a comprehensive, integrated taxonomy of excitatory neuron types in mouse visual cortex and create a system for integrated, high-dimensional cell type classification that can be extended to the whole brain and potentially across species.

Disruption of the autism gene and chromatin regulator KDM5A alters hippocampal cell identity

Chromatin regulation plays a pivotal role in establishing and maintaining cellular identity and is one of the top pathways disrupted in autism spectrum disorder (ASD). The hippocampus, composed of distinct cell types, is often affected in patients with ASD. However, the specific hippocampal cell types and their transcriptional programs that are dysregulated in ASD are unknown. Using single-nucleus RNA sequencing, we show that the ASD gene, lysine demethylase 5A (KDM5A), regulates the development of specific subtypes of excitatory and inhibitory neurons. We found that KDM5A is essential for establishing hippocampal cell identity by controlling a differentiation switch early in development. Our findings define a role for the chromatin regulator KDM5A in establishing hippocampal cell identity and contribute to the emerging convergent mechanisms across ASD.

Controlling morpho-electrophysiological variability of neurons with detailed biophysical models

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

Axons of cortical basket cells originating from dendrites develop higher local complexity than axons emerging from basket cell somata

Neuronal differentiation is regulated by neuronal activity. Here, we analyzed dendritic and axonal growth of Basket cells (BCs) and non-Basket cells (non-BCs) using sparse transfection of channelrhodopsin-YFP and repetitive optogenetic stimulation in slice cultures of rat visual cortex. Neocortical interneurons often display axon-carrying dendrites (AcDs). We found that the AcDs of BCs and non-BCs were, on average, the most complex dendrites. Further, the AcD configuration had an influence on BC axonal development. Axons originating from an AcD formed denser arborizations with more terminal endings within the dendritic field of the parent cell. Intriguingly, this occurred already in unstimulated BCs, and complexity was not increased further by optogenetic stimulation. However, optogenetic stimulation exerted a growth-promoting effect on axons emerging from BC somata. The axons of non-BCs neither responded to the AcD configuration nor to the optogenetic stimulation. The results suggest that the formation of locally dense BC plexuses is regulated by spontaneous activity. Moreover, in the AcD configuration, the AcD and the axon it carries mutually support each other's growth.

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

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

Laminar neural dynamics of auditory evoked responses: Computational modeling of local field potentials in auditory cortex of non-human primates

Evoked neural responses to sensory stimuli have been extensively investigated in humans and animal models both to enhance our understanding of brain function and to aid in clinical diagnosis of neurological and neuropsychiatric conditions. Recording and imaging techniques such as electroencephalography (EEG), magnetoencephalography (MEG), local field potentials (LFPs), and calcium imaging provide complementary information about different aspects of brain activity at different spatial and temporal scales. Modeling and simulations provide a way to integrate these different types of information to clarify underlying neural mechanisms. In this study, we aimed to shed light on the neural dynamics underlying auditory evoked responses by fitting a rate-based model to LFPs recorded via multi-contact electrodes which simultaneously sampled neural activity across cortical laminae. Recordings included neural population responses to best-frequency (BF) and non-BF tones at four representative sites in primary auditory cortex (A1) of awake monkeys. The model considered major neural populations of excitatory, parvalbumin-expressing (PV), and somatostatin-expressing (SOM) neurons across layers 2/3, 4, and 5/6. Unknown parameters, including the connection strength between the populations, were fitted to the data. Our results revealed similar population dynamics, fitted model parameters, predicted equivalent current dipoles (ECD), tuning curves, and lateral inhibition profiles across recording sites and animals, in spite of quite different extracellular current distributions. We found that PV firing rates were higher in BF than in non-BF responses, mainly due to different strengths of tonotopic thalamic input, whereas SOM firing rates were higher in non-BF than in BF responses due to lateral inhibition. In conclusion, we demonstrate the feasibility of the model-fitting approach in identifying the contributions of cell-type specific population activity to stimulus-evoked LFPs across cortical laminae, providing a foundation for further investigations into the dynamics of neural circuits underlying cortical sensory processing.

Multi-scale model of axonal and dendritic polarization by transcranial direct current stimulation in realistic head geometry

BACKGROUND

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation modality that can alter cortical excitability. However, it remains unclear how the subcellular elements of different neuron types are polarized by specific electric field (E-field) distributions.

OBJECTIVE

To quantify neuronal polarization generated by tDCS in a multi-scale computational model.

METHODS

We embedded layer-specific, morphologically-realistic cortical neuron models in a finite element model of the E-field in a human head and simulated steady-state polarization generated by conventional primary-motor-cortex-supraorbital (M1-SO) and 4 × 1 high-definition (HD) tDCS. We quantified somatic, axonal, and dendritic polarization of excitatory pyramidal cells in layers 2/3, 5, and 6, as well as inhibitory interneurons in layers 1 and 4 of the hand knob.

RESULTS

Axonal and dendritic terminals were polarized more than the soma in all neurons, with peak axonal and dendritic polarization of 0.92 mV and 0.21 mV, respectively, compared to peak somatic polarization of 0.07 mV for 1.8 mA M1-SO stimulation. Both montages generated regions of depolarization and hyperpolarization beneath the M1 anode; M1-SO produced slightly stronger, more diffuse polarization peaking in the central sulcus, while 4 × 1 HD produced higher peak polarization in the gyral crown. The E-field component normal to the cortical surface correlated strongly with pyramidal neuron somatic polarization (R2>0.9), but exhibited weaker correlations with peak pyramidal axonal and dendritic polarization (R2:0.5-0.9) and peak polarization in all subcellular regions of interneurons (R2:0.3-0.6). Simulating polarization by uniform local E-field extracted at the soma approximated the spatial distribution of tDCS polarization but produced large errors in some regions (median absolute percent error: 7.9 %).

CONCLUSIONS

Polarization of pre- and postsynaptic compartments of excitatory and inhibitory cortical neurons may play a significant role in tDCS neuromodulation. These effects cannot be predicted from the E-field distribution alone but rather require calculation of the neuronal response.

The energy challenges of artificial superintelligence

We argue here that contemporary semiconductor computing technology poses a significant if not insurmountable barrier to the emergence of any artificial general intelligence system, let alone one anticipated by many to be "superintelligent". This limit on artificial superintelligence (ASI) emerges from the energy requirements of a system that would be more intelligent but orders of magnitude less efficient in energy use than human brains. An ASI would have to supersede not only a single brain but a large population given the effects of collective behavior on the advancement of societies, further multiplying the energy requirement. A hypothetical ASI would likely consume orders of magnitude more energy than what is available in highly-industrialized nations. We estimate the energy use of ASI with an equation we term the "Erasi equation", for the Energy Requirement for Artificial SuperIntelligence. Additional efficiency consequences will emerge from the current unfocussed and scattered developmental trajectory of AI research. Taken together, these arguments suggest that the emergence of an ASI is highly unlikely in the foreseeable future based on current computer architectures, primarily due to energy constraints, with biomimicry or other new technologies being possible solutions.

Morphoelectric and transcriptomic divergence of the layer 1 interneuron repertoire in human versus mouse neocortex

Neocortical layer 1 (L1) is a site of convergence between pyramidal-neuron dendrites and feedback axons where local inhibitory signaling can profoundly shape cortical processing. Evolutionary expansion of human neocortex is marked by distinctive pyramidal neurons with extensive L1 branching, but whether L1 interneurons are similarly diverse is underexplored. Using Patch-seq recordings from human neurosurgical tissue, we identified four transcriptomic subclasses with mouse L1 homologs, along with distinct subtypes and types unmatched in mouse L1. Subclass and subtype comparisons showed stronger transcriptomic differences in human L1 and were correlated with strong morphoelectric variability along dimensions distinct from mouse L1 variability. Accompanied by greater layer thickness and other cytoarchitecture changes, these findings suggest that L1 has diverged in evolution, reflecting the demands of regulating the expanded human neocortical circuit.

Signature morphoelectric properties of diverse GABAergic interneurons in the human neocortex

Human cortex transcriptomic studies have revealed a hierarchical organization of γ-aminobutyric acid-producing (GABAergic) neurons from subclasses to a high diversity of more granular types. Rapid GABAergic neuron viral genetic labeling plus Patch-seq (patch-clamp electrophysiology plus single-cell RNA sequencing) sampling in human brain slices was used to reliably target and analyze GABAergic neuron subclasses and individual transcriptomic types. This characterization elucidated transitions between PVALB and SST subclasses, revealed morphological heterogeneity within an abundant transcriptomic type, identified multiple spatially distinct types of the primate-specialized double bouquet cells (DBCs), and shed light on cellular differences between homologous mouse and human neocortical GABAergic neuron types. These results highlight the importance of multimodal phenotypic characterization for refinement of emerging transcriptomic cell type taxonomies and for understanding conserved and specialized cellular properties of human brain cell types.

Whole human-brain mapping of single cortical neurons for profiling morphological diversity and stereotypy

Quantifying neuron morphology and distribution at the whole-brain scale is essential to understand the structure and diversity of cell types. It is exceedingly challenging to reuse recent technologies of single-cell labeling and whole-brain imaging to study human brains. We propose adaptive cell tomography (ACTomography), a low-cost, high-throughput, and high-efficacy tomography approach, based on adaptive targeting of individual cells. We established a platform to inject dyes into cortical neurons in surgical tissues of 18 patients with brain tumors or other conditions and one donated fresh postmortem brain. We collected three-dimensional images of 1746 cortical neurons, of which 852 neurons were reconstructed to quantify local dendritic morphology, and mapped to standard atlases. In our data, human neurons are more diverse across brain regions than by subject age or gender. The strong stereotypy within cohorts of brain regions allows generating a statistical tensor field of neuron morphology to characterize anatomical modularity of a human brain.

How far neuroscience is from understanding brains

The cellular biology of brains is relatively well-understood, but neuroscientists have not yet generated a theory explaining how brains work. Explanations of how neurons collectively operate to produce what brains can do are tentative and incomplete. Without prior assumptions about the brain mechanisms, I attempt here to identify major obstacles to progress in neuroscientific understanding of brains and central nervous systems. Most of the obstacles to our understanding are conceptual. Neuroscience lacks concepts and models rooted in experimental results explaining how neurons interact at all scales. The cerebral cortex is thought to control awake activities, which contrasts with recent experimental results. There is ambiguity distinguishing task-related brain activities from spontaneous activities and organized intrinsic activities. Brains are regarded as driven by external and internal stimuli in contrast to their considerable autonomy. Experimental results are explained by sensory inputs, behavior, and psychological concepts. Time and space are regarded as mutually independent variables for spiking, post-synaptic events, and other measured variables, in contrast to experimental results. Dynamical systems theory and models describing evolution of variables with time as the independent variable are insufficient to account for central nervous system activities. Spatial dynamics may be a practical solution. The general hypothesis that measurements of changes in fundamental brain variables, action potentials, transmitter releases, post-synaptic transmembrane currents, etc., propagating in central nervous systems reveal how they work, carries no additional assumptions. Combinations of current techniques could reveal many aspects of spatial dynamics of spiking, post-synaptic processing, and plasticity in insects and rodents to start with. But problems defining baseline and reference conditions hinder interpretations of the results. Furthermore, the facts that pooling and averaging of data destroy their underlying dynamics imply that single-trial designs and statistics are necessary.

Distinctive biophysical features of human cell-types: insights from studies of neurosurgically resected brain tissue

Electrophysiological characterization of live human tissue from epilepsy patients has been performed for many decades. Although initially these studies sought to understand the biophysical and synaptic changes associated with human epilepsy, recently, it has become the mainstay for exploring the distinctive biophysical and synaptic features of human cell-types. Both epochs of these human cellular electrophysiological explorations have faced criticism. Early studies revealed that cortical pyramidal neurons obtained from individuals with epilepsy appeared to function "normally" in comparison to neurons from non-epilepsy controls or neurons from other species and thus there was little to gain from the study of human neurons from epilepsy patients. On the other hand, contemporary studies are often questioned for the "normalcy" of the recorded neurons since they are derived from epilepsy patients. In this review, we discuss our current understanding of the distinct biophysical features of human cortical neurons and glia obtained from tissue removed from patients with epilepsy and tumors. We then explore the concept of within cell-type diversity and its loss (i.e., "neural homogenization"). We introduce neural homogenization to help reconcile the epileptogenicity of seemingly "normal" human cortical cells and circuits. We propose that there should be continued efforts to study cortical tissue from epilepsy patients in the quest to understand what makes human cell-types "human".

Reproducibility of biophysical in silico neuron states and spikes from event-based partial histories

Biophysically detailed simulations of neuronal activity often rely on solving large systems of differential equations; in some models, these systems have tens of thousands of states per cell. Numerically solving these equations is computationally intensive and requires making assumptions about the initial cell states. Additional realism from incorporating more biological detail is achieved at the cost of increasingly more states, more computational resources, and more modeling assumptions. We show that for both a point and morphologically-detailed cell model, the presence and timing of future action potentials is probabilistically well-characterized by the relative timings of a moderate number of recent events alone. Knowledge of initial conditions or full synaptic input history is not required. While model time constants, etc. impact the specifics, we demonstrate that for both individual spikes and sustained cellular activity, the uncertainty in spike response decreases as the number of known input events increases, to the point of approximate determinism. Further, we show cellular model states are reconstructable from ongoing synaptic events, despite unknown initial conditions. We propose that a strictly event-based modeling framework is capable of representing the complexity of cellular dynamics of the differential-equations models with significantly less per-cell state variables, thus offering a pathway toward utilizing modern data-driven modeling to scale up to larger network models while preserving individual cellular biophysics.

Simulations predict differing phase responses to excitation vs. inhibition in theta-resonant pyramidal neurons

Rhythmic activity is ubiquitous in neural systems, with theta-resonant pyramidal neurons integrating rhythmic inputs in many cortical structures. Impedance analysis has been widely used to examine frequency-dependent responses of neuronal membranes to rhythmic inputs, but it assumes that the neuronal membrane is a linear system, requiring the use of small signals to stay in a near-linear regime. However, postsynaptic potentials are often large and trigger nonlinear mechanisms (voltage-gated ion channels). The goals of this work were to 1) develop an analysis method to evaluate membrane responses in this nonlinear domain and 2) explore phase relationships between rhythmic stimuli and subthreshold and spiking membrane potential (Vmemb) responses in models of theta-resonant pyramidal neurons. Responses in these output regimes were asymmetrical, with different phase shifts during hyperpolarizing and depolarizing half-cycles. Suprathreshold theta-rhythmic stimuli produced nonstationary Vmemb responses. Sinusoidal inputs produced "phase retreat": action potentials occurred progressively later in cycles of the input stimulus, resulting from adaptation. Sinusoidal current with increasing amplitude over cycles produced "phase advance": action potentials occurred progressively earlier. Phase retreat, phase advance, and subthreshold phase shifts were modulated by multiple ion channel conductances. Our results suggest differential responses of cortical neurons depending on the frequency of oscillatory input, which will play a role in neuronal responses to shifts in network state. We hypothesize that intrinsic cellular properties complement network properties and contribute to in vivo phase-shift phenomena such as phase precession, seen in place and grid cells, and phase roll, also observed in hippocampal CA1 neurons.NEW & NOTEWORTHY We augmented electrical impedance analysis to characterize phase shifts between large-amplitude current stimuli and nonlinear, asymmetric membrane potential responses. We predict different frequency-dependent phase shifts in response excitation vs. inhibition, as well as shifts in spike timing over multiple input cycles, in theta-resonant pyramidal neurons. We hypothesize that these effects contribute to navigation-related phenomena such as phase precession and phase roll. Our neuron-level hypothesis complements, rather than falsifies, prior network-level explanations of these phenomena.

MicroRNAs are necessary for the emergence of Purkinje cell identity

Brain computations are dictated by the unique morphology and connectivity of neuronal subtypes, features established by closely timed developmental events. MicroRNAs (miRNAs) are critical for brain development, but current technologies lack the spatiotemporal resolution to determine how miRNAs instruct the steps leading to subtype identity. Here, we developed new tools to tackle this major gap. Fast and reversible miRNA loss-of-function revealed that miRNAs are necessary for cerebellar Purkinje cell (PC) differentiation, which previously appeared miRNA-independent, and resolved distinct miRNA critical windows in PC dendritogenesis and climbing fiber synaptogenesis, key determinants of PC identity. To identify underlying mechanisms, we generated a mouse model, which enables precise mapping of miRNAs and their targets in rare cell types. With PC-specific maps, we found that the PC-enriched miR-206 drives exuberant dendritogenesis and modulates synaptogenesis. Our results showcase vastly improved approaches for dissecting miRNA function and reveal that many critical miRNA mechanisms remain largely unexplored.

Highlights

Fast miRNA loss-of-function with T6B impairs postnatal Purkinje cell developmentReversible T6B reveals critical miRNA windows for dendritogenesis and synaptogenesisConditional Spy3-Ago2 mouse line enables miRNA-target network mapping in rare cellsPurkinje cell-enriched miR-206 regulates its unique dendritic and synaptic morphology.

KCNQ2/3 Gain-of-Function Variants and Cell Excitability: Differential Effects in CA1 versus L2/3 Pyramidal Neurons

Gain-of-function (GOF) pathogenic variants in the potassium channels KCNQ2 and KCNQ3 lead to hyperexcitability disorders such as epilepsy and autism spectrum disorders. However, the underlying cellular mechanisms of how these variants impair forebrain function are unclear. Here, we show that the R201C variant in KCNQ2 has opposite effects on the excitability of two types of mouse pyramidal neurons of either sex, causing hyperexcitability in layer 2/3 (L2/3) pyramidal neurons and hypoexcitability in CA1 pyramidal neurons. Similarly, the homologous R231C variant in KCNQ3 leads to hyperexcitability in L2/3 pyramidal neurons and hypoexcitability in CA1 pyramidal neurons. However, the effects of KCNQ3 gain-of-function on excitability are specific to superficial CA1 pyramidal neurons. These findings reveal a new level of complexity in the function of KCNQ2 and KCNQ3 channels in the forebrain and provide a framework for understanding the effects of gain-of-function variants and potassium channels in the brain.SIGNIFICANCE STATEMENT KCNQ2/3 gain-of-function (GOF) variants lead to severe forms of neurodevelopmental disorders, but the mechanisms by which these channels affect neuronal activity are poorly understood. In this study, using a series of transgenic mice we demonstrate that the same KCNQ2/3 GOF variants can lead to either hyperexcitability or hypoexcitability in different types of pyramidal neurons [CA1 vs layer (L)2/3]. Additionally, we show that expression of the recurrent KCNQ2 GOF variant R201C in forebrain pyramidal neurons could lead to seizures and SUDEP. Our data suggest that the effects of KCNQ2/3 GOF variants depend on specific cell types and brain regions, possibly accounting for the diverse range of phenotypes observed in individuals with KCNQ2/3 GOF variants.

A GPU-based computational framework that bridges neuron simulation and artificial intelligence

Biophysically detailed multi-compartment models are powerful tools to explore computational principles of the brain and also serve as a theoretical framework to generate algorithms for artificial intelligence (AI) systems. However, the expensive computational cost severely limits the applications in both the neuroscience and AI fields. The major bottleneck during simulating detailed compartment models is the ability of a simulator to solve large systems of linear equations. Here, we present a novel Dendritic Hierarchical Scheduling (DHS) method to markedly accelerate such a process. We theoretically prove that the DHS implementation is computationally optimal and accurate. This GPU-based method performs with 2-3 orders of magnitude higher speed than that of the classic serial Hines method in the conventional CPU platform. We build a DeepDendrite framework, which integrates the DHS method and the GPU computing engine of the NEURON simulator and demonstrate applications of DeepDendrite in neuroscience tasks. We investigate how spatial patterns of spine inputs affect neuronal excitability in a detailed human pyramidal neuron model with 25,000 spines. Furthermore, we provide a brief discussion on the potential of DeepDendrite for AI, specifically highlighting its ability to enable the efficient training of biophysically detailed models in typical image classification tasks.