Reconstruction and Simulation of Neocortical Microcircuitry

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.

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.

Scaling up cell-counting efforts in neuroscience through semi-automated methods

Quantifying how the cellular composition of brain regions vary across development, aging, sex, and disease, is crucial in experimental neuroscience, and the accuracy of different counting methods is continuously debated. Due to the tedious nature of most counting procedures, studies are often restricted to one or a few brain regions. Recently, there have been considerable methodological advances in combining semi-automated feature extraction with brain atlases for cell quantification. Such methods hold great promise for scaling up cell-counting efforts. However, little focus has been paid to how these methods should be implemented and reported to support reproducibility. Here, we provide an overview of practices for conducting and reporting cell counting in mouse and rat brains, showing that critical details for interpretation are typically lacking. We go on to discuss how novel methods may increase efficiency and reproducibility of cell counting studies. Lastly, we provide practical recommendations for researchers planning cell counting.

Shared and divergent principles of synaptic transmission between cortical excitatory neurons in rodent and human brain

Information transfer between principal neurons in neocortex occurs through (glutamatergic) synaptic transmission. In this focussed review, we provide a detailed overview on the strength of synaptic neurotransmission between pairs of excitatory neurons in human and laboratory animals with a specific focus on data obtained using patch clamp electrophysiology. We reach two major conclusions: (1) the synaptic strength, measured as unitary excitatory postsynaptic potential (or uEPSP), is remarkably consistent across species, cortical regions, layers and/or cell-types (median 0.5 mV, interquartile range 0.4-1.0 mV) with most variability associated with the cell-type specific connection studied (min 0.1-max 1.4 mV), (2) synaptic function cannot be generalized across human and rodent, which we exemplify by discussing the differences in anatomical and functional properties of pyramidal-to-pyramidal connections within human and rodent cortical layers 2 and 3. With only a handful of studies available on synaptic transmission in human, it is obvious that much remains unknown to date. Uncovering the shared and divergent principles of synaptic transmission across species however, will almost certainly be a pivotal step toward understanding human cognitive ability and brain function in health and disease.

Identification of interacting neural populations: methods and statistical considerations

As improved recording technologies have created new opportunities for neurophysiological investigation, emphasis has shifted from individual neurons to multiple populations that form circuits, and it has become important to provide evidence of cross-population coordinated activity. We review various methods for doing so, placing them in six major categories while avoiding technical descriptions and instead focusing on high-level motivations and concerns. Our aim is to indicate what the methods can achieve and the circumstances under which they are likely to succeed. Toward this end, we include a discussion of four cross-cutting issues: the definition of neural populations, trial-to-trial variability and Poisson-like noise, time-varying dynamics, and causality.

A whole-task brain model of associative recognition that accounts for human behavior and neuroimaging data

Brain models typically focus either on low-level biological detail or on qualitative behavioral effects. In contrast, we present a biologically-plausible spiking-neuron model of associative learning and recognition that accounts for both human behavior and low-level brain activity across the whole task. Based on cognitive theories and insights from machine-learning analyses of M/EEG data, the model proceeds through five processing stages: stimulus encoding, familiarity judgement, associative retrieval, decision making, and motor response. The results matched human response times and source-localized MEG data in occipital, temporal, prefrontal, and precentral brain regions; as well as a classic fMRI effect in prefrontal cortex. This required two main conceptual advances: a basal-ganglia-thalamus action-selection system that relies on brief thalamic pulses to change the functional connectivity of the cortex, and a new unsupervised learning rule that causes very strong pattern separation in the hippocampus. The resulting model shows how low-level brain activity can result in goal-directed cognitive behavior in humans.

Dynamics of neocortical networks: connectivity beyond the canonical microcircuit

The neocortical network consists of two types of excitatory neurons and a variety of GABAergic inhibitory interneurons, which are organized in distinct microcircuits providing feedforward, feedback, lateral inhibition, and disinhibition. This network is activated by layer- and cell-type specific inputs from first and higher order thalamic nuclei, other subcortical regions, and by cortico-cortical projections. Parallel and serial information processing occurs simultaneously in different intracortical subnetworks and is influenced by neuromodulatory inputs arising from the basal forebrain (cholinergic), raphe nuclei (serotonergic), locus coeruleus (noradrenergic), and ventral tegmentum (dopaminergic). Neocortical neurons differ in their intrinsic firing pattern, in their local and global synaptic connectivity, and in the dynamics of their synaptic interactions. During repetitive stimulation, synaptic connections between distinct neuronal cell types show short-term facilitation or depression, thereby activating or inactivating intracortical microcircuits. Specific networks are capable to generate local and global activity patterns (e.g., synchronized oscillations), which contribute to higher cognitive function and behavior. This review article aims to give a brief overview on our current understanding of the structure and function of the neocortical network.

Inhibitory stabilized network behaviour in a balanced neural mass model of a cortical column

Strong inhibitory recurrent connections can reduce the tendency for a neural network to become unstable. This is known as inhibitory stabilization; networks that are unstable in the absence of strong inhibitory feedback because of their unstable excitatory recurrent connections are known as Inhibition Stabilized Networks (ISNs). One of the characteristics of ISNs is their "paradoxical response", where perturbing the inhibitory neurons with additional excitatory input results in a decrease in their activity after a temporal delay instead of increasing their activity. Here, we develop a model of populations of neurons across different layers of cortex. Within each layer, there is one population of inhibitory neurons and one population of excitatory neurons. The connectivity weights across different populations in the model are derived from a synaptic physiology database provided by the Allen Institute. The model shows a gradient of excitation-inhibition balance across different layers in the cortex, where superficial layers are more inhibitory dominated compared to deeper layers. To investigate the presence of ISNs across different layers, we measured the membrane potentials of neural populations in the model after perturbing inhibitory populations. The results show that layer 2/3 in the model does not operate in the ISN regime but layers 4 and 5 do operate in the ISN regime. These results accord with neurophysiological findings that explored the presence of ISNs across different layers in the cortex. The results show that there may be a systematic macroscopic gradient of inhibitory stabilization across different layers in the cortex that depends on the level of excitation-inhibition balance, and that the strength of the paradoxical response increases as the model moves closer to bifurcation points.

Architectural organization of ∼1,500-neuron modular minicolumnar disinhibitory circuits in healthy and Alzheimer's cortices

Acquisition of neuronal circuit architectures, central to understanding brain function and dysfunction, remains prohibitively challenging. Here I report the development of a simultaneous and sequential octuple-sexdecuple whole-cell patch-clamp recording system that enables architectural reconstruction of complex cortical circuits. The method unveils the canonical layer 1 single bouquet cell (SBC)-led disinhibitory neuronal circuits across the mouse somatosensory, motor, prefrontal, and medial entorhinal cortices. The ∼1,500-neuron modular circuits feature the translaminar, unidirectional, minicolumnar, and independent disinhibition and optimize cortical complexity, subtlety, plasticity, variation, and redundancy. Moreover, architectural reconstruction uncovers age-dependent deficits at SBC-disinhibited synapses in the senescence-accelerated mouse prone 8, an animal model of Alzheimer's disease. The deficits exhibit the characteristic Alzheimer's-like cortical spread and correlation with cognitive impairments. These findings decrypt operations of the elementary processing units in healthy and Alzheimer's mouse cortices and validate the efficacy of octuple-sexdecuple patch-clamp recordings for architectural reconstruction of complex neuronal circuits.

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. Simulating polarization by uniform local E-field approximated the spatial distribution of tDCS polarization but produced large errors in some regions.

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.

3D synaptic organization of layer III of the human anterior cingulate and temporopolar cortex

The human anterior cingulate and temporopolar cortices have been proposed as highly connected nodes involved in high-order cognitive functions, but their synaptic organization is still basically unknown due to the difficulties involved in studying the human brain. Using Focused Ion Beam/Scanning Electron Microscopy (FIB/SEM) to study the synaptic organization of the human brain obtained with a short post-mortem delay allows excellent results to be obtained. We have used this technology to analyze layer III of the anterior cingulate cortex (Brodmann area 24) and the temporopolar cortex, including the temporal pole (Brodmann area 38 ventral and dorsal) and anterior middle temporal gyrus (Brodmann area 21). Our results, based on 6695 synaptic junctions fully reconstructed in 3D, revealed that Brodmann areas 24, 21 and ventral area 38 showed similar synaptic density and synaptic size, whereas dorsal area 38 displayed the highest synaptic density and the smallest synaptic size. However, the proportion of the different types of synapses (excitatory and inhibitory), the postsynaptic targets, and the shapes of excitatory and inhibitory synapses were similar, regardless of the region examined. These observations indicate that certain aspects of the synaptic organization are rather homogeneous, whereas others show specific variations across cortical regions.

Single-cell spatial transcriptome reveals cell-type organization in the macaque cortex

Elucidating the cellular organization of the cerebral cortex is critical for understanding brain structure and function. Using large-scale single-nucleus RNA sequencing and spatial transcriptomic analysis of 143 macaque cortical regions, we obtained a comprehensive atlas of 264 transcriptome-defined cortical cell types and mapped their spatial distribution across the entire cortex. We characterized the cortical layer and region preferences of glutamatergic, GABAergic, and non-neuronal cell types, as well as regional differences in cell-type composition and neighborhood complexity. Notably, we discovered a relationship between the regional distribution of various cell types and the region's hierarchical level in the visual and somatosensory systems. Cross-species comparison of transcriptomic data from human, macaque, and mouse cortices further revealed primate-specific cell types that are enriched in layer 4, with their marker genes expressed in a region-dependent manner. Our data provide a cellular and molecular basis for understanding the evolution, development, aging, and pathogenesis of the primate brain.

Quantitative analysis of the optogenetic excitability of CA1 neurons

Introduction

Optogenetics has emerged as a promising technique for modulating neuronal activity and holds potential for the treatment of neurological disorders such as temporal lobe epilepsy (TLE). However, clinical translation still faces many challenges. This in-silico study aims to enhance the understanding of optogenetic excitability in CA1 cells and to identify strategies for improving stimulation protocols.

Methods

Employing state-of-the-art computational models coupled with Monte Carlo simulated light propagation, the optogenetic excitability of four CA1 cells, two pyramidal and two interneurons, expressing ChR2(H134R) is investigated.

Results and discussion

The results demonstrate that confining the opsin to specific neuronal membrane compartments significantly improves excitability. An improvement is also achieved by focusing the light beam on the most excitable cell region. Moreover, the perpendicular orientation of the optical fiber relative to the somato-dendritic axis yields superior results. Inter-cell variability is observed, highlighting the importance of considering neuron degeneracy when designing optogenetic tools. Opsin confinement to the basal dendrites of the pyramidal cells renders the neuron the most excitable. A global sensitivity analysis identified opsin location and expression level as having the greatest impact on simulation outcomes. The error reduction of simulation outcome due to coupling of neuron modeling with light propagation is shown. The results promote spatial confinement and increased opsin expression levels as important improvement strategies. On the other hand, uncertainties in these parameters limit precise determination of the irradiance thresholds. This study provides valuable insights on optogenetic excitability of CA1 cells useful for the development of improved optogenetic stimulation protocols for, for instance, TLE treatment.

Brain criticality predicts individual levels of inter-areal synchronization in human electrophysiological data

Neuronal oscillations and their synchronization between brain areas are fundamental for healthy brain function. Yet, synchronization levels exhibit large inter-individual variability that is associated with behavioral variability. We test whether individual synchronization levels are predicted by individual brain states along an extended regime of critical-like dynamics - the Griffiths phase (GP). We use computational modelling to assess how synchronization is dependent on brain criticality indexed by long-range temporal correlations (LRTCs). We analyze LRTCs and synchronization of oscillations from resting-state magnetoencephalography and stereo-electroencephalography data. Synchronization and LRTCs are both positively linearly and quadratically correlated among healthy subjects, while in epileptogenic areas they are negatively linearly correlated. These results show that variability in synchronization levels is explained by the individual position along the GP with healthy brain areas operating in its subcritical and epileptogenic areas in its supercritical side. We suggest that the GP is fundamental for brain function allowing individual variability while retaining functional advantages of criticality.

How to incorporate biological insights into network models and why it matters

Due to the staggering complexity of the brain and its neural circuitry, neuroscientists rely on the analysis of mathematical models to elucidate its function. From Hodgkin and Huxley's detailed description of the action potential in 1952 to today, new theories and increasing computational power have opened up novel avenues to study how neural circuits implement the computations that underlie behaviour. Computational neuroscientists have developed many models of neural circuits that differ in complexity, biological realism or emergent network properties. With recent advances in experimental techniques for detailed anatomical reconstructions or large-scale activity recordings, rich biological data have become more available. The challenge when building network models is to reflect experimental results, either through a high level of detail or by finding an appropriate level of abstraction. Meanwhile, machine learning has facilitated the development of artificial neural networks, which are trained to perform specific tasks. While they have proven successful at achieving task-oriented behaviour, they are often abstract constructs that differ in many features from the physiology of brain circuits. Thus, it is unclear whether the mechanisms underlying computation in biological circuits can be investigated by analysing artificial networks that accomplish the same function but differ in their mechanisms. Here, we argue that building biologically realistic network models is crucial to establishing causal relationships between neurons, synapses, circuits and behaviour. More specifically, we advocate for network models that consider the connectivity structure and the recorded activity dynamics while evaluating task performance.

Simulations of cortical networks using spatially extended conductance-based neuronal models

The Hodgkin-Huxley model of action potential generation and propagation, published in the Journal of Physiology in 1952, initiated the field of biophysically detailed computational modelling in neuroscience, which has expanded to encompass a variety of species and components of the nervous system. Here we review the developments in this area with a focus on efforts in the community towards modelling the mammalian neocortex using spatially extended conductance-based neuronal models. The Hodgkin-Huxley formalism and related foundational contributions, such as Rall's cable theory, remain widely used in these efforts to the current day. We argue that at present the field is undergoing a qualitative change due to new very rich datasets describing the composition, connectivity and functional activity of cortical circuits, which are being integrated systematically into large-scale network models. This trend, combined with the accelerating development of convenient software tools supporting such complex modelling projects, is giving rise to highly detailed models of the cortex that are extensively constrained by the data, enabling computational investigation of a multitude of questions about cortical structure and function.

From single-neuron dynamics to higher-order circuit motifs in control and pathological brain networks

The convergence of advanced single-cell in vivo functional imaging techniques, computational modelling tools and graph-based network analytics has heralded new opportunities to study single-cell dynamics across large-scale networks, providing novel insights into principles of brain communication and pointing towards potential new strategies for treating neurological disorders. A major recent finding has been the identification of unusually richly connected hub cells that have capacity to synchronize networks and may also be critical in network dysfunction. While hub neurons are traditionally defined by measures that consider solely the number and strength of connections, novel higher-order graph analytics now enables the mining of massive networks for repeating subgraph patterns called motifs. As an illustration of the power offered by higher-order analysis of neuronal networks, we highlight how recent methodological advances uncovered a new functional cell type, the superhub, that is predicted to play a major role in regulating network dynamics. Finally, we discuss open questions that will be critical for assessing the importance of higher-order cellular-scale network analytics in understanding brain function in health and disease.

Differential expression of GABAA receptor subunits δ and α6 mediates tonic inhibition in parvalbumin and somatostatin interneurons in the mouse hippocampus

Inhibitory γ-aminobutyric acid (GABA)-ergic interneurons mediate inhibition in neuronal circuitry and support normal brain function. Consequently, dysregulation of inhibition is implicated in various brain disorders. Parvalbumin (PV) and somatostatin (SST) interneurons, the two major types of GABAergic inhibitory interneurons in the hippocampus, exhibit distinct morpho-physiological properties and coordinate information processing and memory formation. However, the molecular mechanisms underlying the specialized properties of PV and SST interneurons remain unclear. This study aimed to compare the transcriptomic differences between these two classes of interneurons in the hippocampus using the ribosome tagging approach. The results revealed distinct expressions of genes such as voltage-gated ion channels and GABAA receptor subunits between PV and SST interneurons. Gabrd and Gabra6 were identified as contributors to the contrasting tonic GABAergic inhibition observed in PV and SST interneurons. Moreover, some of the differentially expressed genes were associated with schizophrenia and epilepsy. In conclusion, our results provide molecular insights into the distinct roles of PV and SST interneurons in health and disease.

Efficient simulation of neural development using shared memory parallelization

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

Code available at

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

Application of quantum machine learning using quantum kernel algorithms on multiclass neuron M-type classification

The functional characterization of different neuronal types has been a longstanding and crucial challenge. With the advent of physical quantum computers, it has become possible to apply quantum machine learning algorithms to translate theoretical research into practical solutions. Previous studies have shown the advantages of quantum algorithms on artificially generated datasets, and initial experiments with small binary classification problems have yielded comparable outcomes to classical algorithms. However, it is essential to investigate the potential quantum advantage using real-world data. To the best of our knowledge, this study is the first to propose the utilization of quantum systems to classify neuron morphologies, thereby enhancing our understanding of the performance of automatic multiclass neuron classification using quantum kernel methods. We examined the influence of feature engineering on classification accuracy and found that quantum kernel methods achieved similar performance to classical methods, with certain advantages observed in various configurations.

Rapid estimation of cortical neuron activation thresholds by transcranial magnetic stimulation using convolutional neural networks

BACKGROUND

Transcranial magnetic stimulation (TMS) can modulate neural activity by evoking action potentials in cortical neurons. TMS neural activation can be predicted by coupling subject-specific head models of the TMS-induced electric field (E-field) to populations of biophysically realistic neuron models; however, the significant computational cost associated with these models limits their utility and eventual translation to clinically relevant applications.

OBJECTIVE

To develop computationally efficient estimators of the activation thresholds of multi-compartmental cortical neuron models in response to TMS-induced E-field distributions.

METHODS

Multi-scale models combining anatomically accurate finite element method (FEM) simulations of the TMS E-field with layer-specific representations of cortical neurons were used to generate a large dataset of activation thresholds. 3D convolutional neural networks (CNNs) were trained on these data to predict thresholds of model neurons given their local E-field distribution. The CNN estimator was compared to an approach using the uniform E-field approximation to estimate thresholds in the non-uniform TMS-induced E-field.

RESULTS

The 3D CNNs estimated thresholds with mean absolute percent error (MAPE) on the test dataset below 2.5% and strong correlation between the CNN predicted and actual thresholds for all cell types (R2 > 0.96). The CNNs estimated thresholds with a 2-4 orders of magnitude reduction in the computational cost of the multi-compartmental neuron models. The CNNs were also trained to predict the median threshold of populations of neurons, speeding up computation further.

CONCLUSION

3D CNNs can estimate rapidly and accurately the TMS activation thresholds of biophysically realistic neuron models using sparse samples of the local E-field, enabling simulating responses of large neuron populations or parameter space exploration on a personal computer.

Correlation Analysis of Molecularly-Defined Cortical Interneuron Populations with Morpho-Electric Properties in Layer V of Mouse Neocortex

Cortical interneurons can be categorized into distinct populations based on multiple modalities, including molecular signatures and morpho-electrical (M/E) properties. Recently, many transcriptomic signatures based on single-cell RNA-seq have been identified in cortical interneurons. However, whether different interneuron populations defined by transcriptomic signature expressions correspond to distinct M/E subtypes is still unknown. Here, we applied the Patch-PCR approach to simultaneously obtain the M/E properties and messenger RNA (mRNA) expression of >600 interneurons in layer V of the mouse somatosensory cortex (S1). Subsequently, we identified 11 M/E subtypes, 9 neurochemical cell populations (NCs), and 20 transcriptomic cell populations (TCs) in this cortical lamina. Further analysis revealed that cells in many NCs and TCs comprised several M/E types and were difficult to clearly distinguish morpho-electrically. A similar analysis of layer V interneurons of mouse primary visual cortex (V1) and motor cortex (M1) gave results largely comparable to S1. Comparison between S1, V1, and M1 suggested that, compared to V1, S1 interneurons were morpho-electrically more similar to M1. Our study reveals the presence of substantial M/E variations in cortical interneuron populations defined by molecular expression.

Distinctive properties of biological neural networks and recent advances in bottom-up approaches toward a better biologically plausible neural network

Although it may appear infeasible and impractical, building artificial intelligence (AI) using a bottom-up approach based on the understanding of neuroscience is straightforward. The lack of a generalized governing principle for biological neural networks (BNNs) forces us to address this problem by converting piecemeal information on the diverse features of neurons, synapses, and neural circuits into AI. In this review, we described recent attempts to build a biologically plausible neural network by following neuroscientifically similar strategies of neural network optimization or by implanting the outcome of the optimization, such as the properties of single computational units and the characteristics of the network architecture. In addition, we proposed a formalism of the relationship between the set of objectives that neural networks attempt to achieve, and neural network classes categorized by how closely their architectural features resemble those of BNN. This formalism is expected to define the potential roles of top-down and bottom-up approaches for building a biologically plausible neural network and offer a map helping the navigation of the gap between neuroscience and AI engineering.

Multiscale model of primary motor cortex circuits predicts in vivo cell-type-specific, behavioral state-dependent dynamics

Understanding cortical function requires studying multiple scales: molecular, cellular, circuit, and behavioral. We develop a multiscale, biophysically detailed model of mouse primary motor cortex (M1) with over 10,000 neurons and 30 million synapses. Neuron types, densities, spatial distributions, morphologies, biophysics, connectivity, and dendritic synapse locations are constrained by experimental data. The model includes long-range inputs from seven thalamic and cortical regions and noradrenergic inputs. Connectivity depends on cell class and cortical depth at sublaminar resolution. The model accurately predicts in vivo layer- and cell-type-specific responses (firing rates and LFP) associated with behavioral states (quiet wakefulness and movement) and experimental manipulations (noradrenaline receptor blockade and thalamus inactivation). We generate mechanistic hypotheses underlying the observed activity and analyzed low-dimensional population latent dynamics. This quantitative theoretical framework can be used to integrate and interpret M1 experimental data and sheds light on the cell-type-specific multiscale dynamics associated with several experimental conditions and behaviors.

Optical Diffractive Convolutional Neural Networks Implemented in an All-Optical Way

Optical neural networks can effectively address hardware constraints and parallel computing efficiency issues inherent in electronic neural networks. However, the inability to implement convolutional neural networks at the all-optical level remains a hurdle. In this work, we propose an optical diffractive convolutional neural network (ODCNN) that is capable of performing image processing tasks in computer vision at the speed of light. We explore the application of the 4f system and the diffractive deep neural network (D2NN) in neural networks. ODCNN is then simulated by combining the 4f system as an optical convolutional layer and the diffractive networks. We also examine the potential impact of nonlinear optical materials on this network. Numerical simulation results show that the addition of convolutional layers and nonlinear functions improves the classification accuracy of the network. We believe that the proposed ODCNN model can be the basic architecture for building optical convolutional networks.

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

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 modelling of post-synaptic plasticity to investigate how schizophrenia-associated genes regulate long-term potentiation (LTP) and depression (LTD). We combined our model with data from post-mortem mRNA 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 anterior cingulate cortex, lead to impaired PKA-pathway-mediated LTP in synapses containing GluR1 receptors. We validated these findings using a genotyped EEG dataset where polygenic risk scores for synaptic and ion channel-encoding genes as well as modulation of visual evoked potentials (VEP) 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.