Reconstruction and Simulation of Neocortical Microcircuitry

Cortical parvalbumin-expressing interneurons sample network oscillations in their synaptic activity

Synaptic activity is thought to be the primary input of the frequency bands conveyed in the electroencephalogram (EEG) and local field potentials (LFPs) recorded on the cortex. Here we ask whether synaptic activity observed in parvalbumin expressing (PV + ) neurons recorded in isolated cortical tissue bear these frequency bands. The muscarinic agonist carbachol (CCh) was used to increase cortical excitability. PV + neurons play a significant role in perisomatic inhibition and the synchronization of cortical ensembles to generate gamma (γ) oscillations during cholinergic modulation. γ-oscillations associate with cognitive activities co-existing with slower rhythms. While CCh induces depolarization and firing in pyramidal neurons, it only causes barrages of synaptic potentials without firing in most PV + neurons. We show that the frequency spectra of CCh-induced synaptic events recorded onto layer 5 PV + neurons display the various frequency bands generated by cortical networks: from δ to γ. Isolation of inhibitory events shows potency increases in the δ band and decreases in other bands. Isolated excitatory events exhibit a decrease in the β-band. Excitatory potentials appear to drive the circuitry while inhibitory ones appear to regulate events frequency. Muscarinic M1-class receptors are mainly responsible for the synaptic activity from which oscillatory bands emerge. These results demonstrate that PV + interneurons "sample" network activity through the ligand-gated synaptic events that receive from it. We conclude that random synaptic events recorded in single neurons contain the wide range of brain oscillations as revealed by frequency spectra and power density analyses.

Adding Space to Random Networks of Spiking Neurons: A Method Based on Scaling the Network Size

Many spiking neural network models are based on random graphs that do not include topological and structural properties featured in real brain networks. To turn these models into spatial networks that describe the topographic arrangement of connections is a challenging task because one has to deal with neurons at the spatial network boundary. Addition of space may generate spurious network behavior like oscillations introduced by periodic boundary conditions or unbalanced neuronal spiking due to lack or excess of connections. Here, we introduce a boundary solution method for networks with added spatial extension that prevents the occurrence of spurious spiking behavior. The method is based on a recently proposed technique for scaling the network size that preserves first- and second-order statistics.

On the validity of electric brain signal predictions based on population firing rates

Neural activity at the population level is commonly studied experimentally through measurements of electric brain signals like local field potentials (LFPs), or electroencephalography (EEG) signals. To allow for comparison between observed and simulated neural activity it is therefore important that simulations of neural activity can accurately predict these brain signals. Simulations of neural activity at the population level often rely on point-neuron network models or firing-rate models. While these simplified representations of neural activity are computationally efficient, they lack the explicit spatial information needed for calculating LFP/EEG signals. Different heuristic approaches have been suggested for overcoming this limitation, but the accuracy of these approaches has not fully been assessed. One such heuristic approach, the so-called kernel method, has previously been applied with promising results and has the additional advantage of being well-grounded in the biophysics underlying electric brain signal generation. It is based on calculating rate-to-LFP/EEG kernels for each synaptic pathway in a network model, after which LFP/EEG signals can be obtained directly from population firing rates. This amounts to a massive reduction in the computational effort of calculating brain signals because the brain signals are calculated for each population instead of for each neuron. Here, we investigate how and when the kernel method can be expected to work, and present a theoretical framework for predicting its accuracy. We show that the relative error of the brain signal predictions is a function of the single-cell kernel heterogeneity and the spike-train correlations. Finally, we demonstrate that the kernel method is most accurate for contributions which are also dominating the brain signals: spatially clustered and correlated synaptic input to large populations of pyramidal cells. We thereby further establish the kernel method as a promising approach for calculating electric brain signals from large-scale neural simulations.

An unsupervised map of excitatory neuron dendritic morphology in the mouse visual cortex

Neurons in the neocortex exhibit astonishing morphological diversity, which is critical for properly wiring neural circuits and giving neurons their functional properties. However, the organizational principles underlying this morphological diversity remain an open question. Here, we took a data-driven approach using graph-based machine learning methods to obtain a low-dimensional morphological "bar code" describing more than 30,000 excitatory neurons in mouse visual areas V1, AL, and RL that were reconstructed from the millimeter scale MICrONS serial-section electron microscopy volume. Contrary to previous classifications into discrete morphological types (m-types), our data-driven approach suggests that the morphological landscape of cortical excitatory neurons is better described as a continuum, with a few notable exceptions in layers 5 and 6. Dendritic morphologies in layers 2-3 exhibited a trend towards a decreasing width of the dendritic arbor and a smaller tuft with increasing cortical depth. Inter-area differences were most evident in layer 4, where V1 contained more atufted neurons than higher visual areas. Moreover, we discovered neurons in V1 on the border to layer 5, which avoided deeper layers with their dendrites. In summary, we suggest that excitatory neurons' morphological diversity is better understood by considering axes of variation than using distinct m-types.

Integrating multimodal data to understand cortical circuit architecture and function

In recent years there has been a tremendous growth in new technologies that allow large-scale investigation of different characteristics of the nervous system at an unprecedented level of detail. There is a growing trend to use combinations of these new techniques to determine direct links between different modalities. In this Perspective, we focus on the mouse visual cortex, as this is one of the model systems in which much progress has been made in the integration of multimodal data to advance understanding. We review several approaches that allow integration of data regarding various properties of cortical cell types, connectivity at the level of brain areas, cell types and individual cells, and functional neural activity in vivo. The increasingly crucial contributions of computation and theory in analyzing and systematically modeling data are also highlighted. Together with open sharing of data, tools and models, integrative approaches are essential tools in modern neuroscience for improving our understanding of the brain architecture, mechanisms and function.

Connectomics of predicted Sst transcriptomic types in mouse visual cortex

Neural circuit function is shaped both by the cell types that comprise the circuit and the connections between them1. Neural cell types have previously been defined by morphology2,3, electrophysiology4, transcriptomic expression5,6, connectivity7-9 or a combination of such modalities10-12. The Patch-seq technique enables the characterization of morphology, electrophysiology and transcriptomic properties from individual cells13-15. These properties were integrated to define 28 inhibitory, morpho-electric-transcriptomic (MET) cell types in mouse visual cortex16, which do not include synaptic connectivity. Conversely, large-scale electron microscopy (EM) enables morphological reconstruction and a near-complete description of a neuron's local synaptic connectivity, but does not include transcriptomic or electrophysiological information. Here, we leveraged morphological information from Patch-seq to predict the transcriptomically defined cell subclass and/or MET-type of inhibitory neurons within a large-scale EM dataset. We further analysed Martinotti cells-a somatostatin (Sst)-positive17 morphological cell type18,19-which were classified successfully into Sst MET-types with distinct axon myelination and synaptic output connectivity patterns. We demonstrate that morphological features can be used to link cell types across experimental modalities, enabling further comparison of connectivity to gene expression and electrophysiology. We observe unique connectivity rules for predicted Sst cell types.

A connectome manipulation framework for the systematic and reproducible study of structure-function relationships through simulations

Synaptic connectivity at the neuronal level is characterized by highly nonrandom features. Hypotheses about their role can be developed by correlating structural metrics to functional features. But, to prove causation, manipulations of connectivity would have to be studied. However, the fine-grained scale at which nonrandom trends are expressed makes this approach challenging to pursue experimentally. Simulations of neuronal networks provide an alternative route to study arbitrarily complex manipulations in morphologically and biophysically detailed models. Here, we present Connectome-Manipulator, a Python framework for rapid connectome manipulations of large-scale network models in Scalable Open Network Architecture TemplAte (SONATA) format. In addition to creating or manipulating the connectome of a model, it provides tools to fit parameters of stochastic connectivity models against existing connectomes. This enables rapid replacement of any existing connectome with equivalent connectomes at different levels of complexity, or transplantation of connectivity features from one connectome to another, for systematic study. We employed the framework in the detailed model of the rat somatosensory cortex in two exemplary use cases: transplanting interneuron connectivity trends from electron microscopy data and creating simplified connectomes of excitatory connectivity. We ran a series of network simulations and found diverse shifts in the activity of individual neuron populations causally linked to these manipulations.

Localization of the epileptogenic network from scalp EEG using a patient-specific whole-brain model

Computational modeling is a key tool for elucidating the neuronal mechanisms underlying epileptic activity. Despite considerable progress, existing models often lack realistic accuracy in representing electrophysiological epileptic activity. In this study, we used a comprehensive human brain model based on a neural mass model, which is tailored to the layered structure of the neocortex and incorporates patient-specific imaging data. This approach allowed the simulation of scalp EEGs in an epileptic patient suffering from type 2 focal cortical dysplasia (FCD). The simulation specifically addressed epileptic activity induced by FCD, faithfully reproducing intracranial interictal epileptiform discharges (IEDs) recorded with electrocorticography. For constructing the patient-specific scalp EEG, we carefully defined a clear delineation of the epileptogenic zone by numerical simulations to ensure fidelity to the topography, polarity, and diffusion characteristics of IEDs. This nuanced approach improves the accuracy of the simulated EEG signal, provides a more accurate representation of epileptic activity, and enhances our understanding of the mechanism behind the epileptogenic networks. The accuracy of the model was confirmed by a postoperative reevaluation with a secondary EEG simulation that was consistent with the lesion's removal. Ultimately, this personalized approach may prove instrumental in optimizing and tailoring epilepsy treatment strategies.

Genotype-function-phenotype correlations for SCN1A variants identified by clinical genetic testing

OBJECTIVE

Interpretation of clinical genetic testing, which identifies a potential genetic etiology in 25% of children with epilepsy, is limited by variants of uncertain significance. Understanding functional consequences of variants can help distinguish pathogenic from benign alleles. We combined automated patch clamp recording with neurophysiological simulations to discern genotype-function-phenotype correlations in a real-world cohort of children with SCN1A-associated epilepsy.

METHODS

Clinical data were extracted for children with SCN1A variants identified by clinical genetic testing. Functional properties of non-truncating NaV1.1 variant channels were determined using automated patch clamp recording. Functional data were incorporated into a parvalbumin-positive (PV+) interneuron computer model to predict variant effects on neuron firing and were compared with longitudinal clinical data describing epilepsy types, neurocognitive outcomes, and medication response.

RESULTS

Twelve SCN1A variants were identified (nine non-truncating). Six non-truncating variants exhibited no measurable sodium current in heterologous cells consistent with complete loss of function (LoF). Two variants caused either partial LoF (L479P) or a mixture of gain and loss of function (I1356M). The remaining non-truncating variant (T1250M) exhibited normal function. Functional data changed classification of pathogenicity for six variants. Complete LoF variants were universally associated with seizure onset before one year of age and febrile seizures, and were often associated with drug resistant epilepsy and below average cognitive outcomes. Simulations demonstrated abnormal firing in heterozygous model neurons containing dysfunctional variants.

INTERPRETATION

In SCN1A-associated epilepsy, functional analysis and neuron simulation studies resolved variants of uncertain significance and correlated with aspects of phenotype and medication response.

Statistical method accounts for microscopic electric field distortions around neurons when simulating activation thresholds

INTRODUCTION

Notwithstanding advances in computational models of neuromodulation, there are mismatches between simulated and experimental activation thresholds. Transcranial Magnetic Stimulation (TMS) of the primary motor cortex generates motor evoked potentials (MEPs). At the threshold of MEP generation, whole-head models predict macroscopic (at millimeter scale) electric fields (50-70 V/m) which are considerably below conventionally simulated cortical neuron thresholds (175-350 V/m).

METHODS

We hypothesize that this apparent contradiction is in part a consequence of electrical field warping by brain microstructure. Classical neuronal models ignore the physical presence of neighboring neurons and microstructure and assume that the macroscopic field directly acts on the neurons. In previous work, we performed advanced numerical calculations considering realistic microscopic compartments (e.g., cells, blood vessels), resulting in locally inhomogeneous (micrometer scale) electric field and altered neuronal activation thresholds. Here we combine detailed neural threshold simulations under homogeneous field assumptions with microscopic field calculations, leveraging a novel statistical approach.

RESULTS

We show that, provided brain-region specific microstructure metrics, a single statistically derived scaling factor between microscopic and macroscopic electric fields can be applied in predicting neuronal thresholds. For the cortical sample considered, the statistical method matches TMS experimental thresholds.

CONCLUSIONS

Our approach can be broadly applied to neuromodulation models, where fully coupled microstructure scale simulations may not be computationally tractable.

Organizational Principles of the Primate Cerebral Cortex at the Single-Cell Level

The primate cerebral cortex, the major organ for cognition, consists of an immense number of neurons. However, the organizational principles governing these neurons remain unclear. By accessing the single-cell spatial transcriptome of over 25 million neuron cells across the entire macaque cortex, it is discovered that the distribution of neurons within cortical layers is highly non-random. Strikingly, over three-quarters of these neurons are located in distinct neuronal clusters. Within these clusters, different cell types tend to collaborate rather than function independently. Typically, excitatory neuron clusters mainly consist of excitatory-excitatory combinations, while inhibitory clusters primarily contain excitatory-inhibitory combinations. Both cluster types have roughly equal numbers of neurons in each layer. Importantly, most excitatory and inhibitory neuron clusters form spatial partnerships, indicating a balanced local neuronal network and correlating with specific functional regions. These organizational principles are conserved across mouse cortical regions. These findings suggest that different brain regions of the cortex may exhibit similar mechanisms at the neuronal population level.

Subthreshold electric fields bidirectionally modulate neurotransmitter release through axon polarization

Subthreshold electric fields modulate brain activity and demonstrate potential in several therapeutic applications. Electric fields are known to generate heterogenous membrane polarization within neurons due to their complex morphologies. While the effects of somatic and dendritic polarization in postsynaptic neurons have been characterized, the functional consequences of axonal polarization on neurotransmitter release from the presynapse are unknown. Here, we combined noninvasive optogenetic indicators of voltage, calcium and neurotransmitter release to study the subcellular response within single neurons to subthreshold electric fields. We first captured the detailed spatiotemporal polarization profile produced by uniform electric fields within individual neurons. Small polarization of presynaptic boutons produces rapid and powerful modulation of neurotransmitter release, with the direction - facilitation or inhibition - depending on the direction of polarization. We determined that subthreshold electric fields drive this effect by rapidly altering the number of synaptic vesicles participating in neurotransmission, producing effects which resemble short-term plasticity akin to presynaptic homeostatic plasticity. These results provide key insights into the mechanisms of subthreshold electric fields at the cellular level.

Abstract Figure

Postsynaptic spiking determines anti-Hebbian LTD in visual cortex basket cells

Long-term plasticity at pyramidal cell to basket cell (PC → BC) synapses is important for the functioning of cortical microcircuits. It is well known that at neocortical PC → PC synapses, dendritic calcium (Ca2+) dynamics signal coincident pre-and postsynaptic spiking which in turn triggers long-term potentiation (LTP). However, the link between dendritic Ca2+ dynamics and long-term plasticity at PC → BC synapses of primary visual cortex (V1) is not as well known. Here, we explored if PC → BC synaptic plasticity in developing V1 is sensitive to postsynaptic spiking. Two-photon (2P) Ca2+ imaging revealed that action potentials (APs) in dendrites of V1 layer-5 (L5) BCs back-propagated decrementally but actively to the location of PC → BC putative synaptic contacts. Pairing excitatory inputs with postsynaptic APs elicited dendritic Ca2+ supralinearities for pre-before-postsynaptic but not post-before-presynaptic temporal ordering, suggesting that APs could impact synaptic plasticity. In agreement, extracellular stimulation as well as high-throughput 2P optogenetic mapping of plasticity both revealed that pre-before-postsynaptic but not post-before-presynaptic pairing resulted in anti-Hebbian long-term depression (LTD). Our results demonstrate that V1 BC dendritic Ca2+ nonlinearities and synaptic plasticity at PC → BC connections are both sensitive to somatic spiking.

Digital twin brain simulator for real-time consciousness monitoring and virtual intervention using primate electrocorticogram data

At the forefront of bridging computational brain modeling with personalized medicine, this study introduces a novel, real-time, electrocorticogram (ECoG) simulator, based on the digital twin brain concept. Utilizing advanced data assimilation techniques, specifically a Variational Bayesian Recurrent Neural Network model with hierarchical latent units, the simulator dynamically predicts ECoG signals reflecting real-time brain latent states. By assimilating broad ECoG signals from macaque monkeys across awake and anesthetized conditions, the model successfully updated its latent states in real-time, enhancing precision of ECoG signal simulations. Behind successful data assimilation, self-organization of latent states in the model was observed, reflecting brain states and individuality. This self-organization facilitated simulation of virtual drug administration and uncovered functional networks underlying changes in brain function during anesthesia. These results show that the proposed model can simulate brain signals in real-time with high accuracy and is also useful for revealing underlying information processing dynamics.

A detailed theory of thalamic and cortical microcircuits for predictive visual inference

Understanding cortical microcircuitry requires theoretical models that can tease apart their computational logic from biological details. Although Bayesian inference serves as an abstract framework of cortical computation, precisely mapping concrete instantiations of computational models to biology under real-world tasks is necessary to produce falsifiable neural models. On the basis of a recent generative model, recursive cortical networks, that demonstrated excellent performance on vision benchmarks, we derive a theoretical cortical microcircuit by placing the requirements of the computational model within biological constraints. The derived model suggests precise algorithmic roles for the columnar and laminar feed-forward, feedback, and lateral connections, the thalamic pathway, blobs and interblobs, and the innate lineage-specific interlaminar connectivity within cortical columns. The model also explains several visual phenomena, including the subjective contour effect and neon-color spreading effect, with circuit-level precision. Our model and methodology provides a path forward in understanding cortical and thalamic computations.

Frequency-dependent phase entrainment of cortical cell types during tACS: computational modeling evidence

Objective. Transcranial alternating current stimulation (tACS) enables non-invasive modulation of brain activity, holding promise for clinical and research applications. Yet, it remains unclear how the stimulation frequency differentially impacts various neuron types. Here, we aimed to quantify the frequency-dependent behavior of key neocortical cell types.Approach. We used both detailed (anatomical multicompartments) and simplified (three compartments) single-cell modeling approaches based on the Hodgkin-Huxley formalism to study neocortical excitatory and inhibitory cells under various tACS intensities and frequencies within the 5-50 Hz range at rest and during basal 10 Hz activity.Main results. L5 pyramidal cells (PCs) exhibited the highest polarizability at direct current, ranging from 0.21 to 0.25 mm and decaying exponentially with frequency. Inhibitory neurons displayed membrane resonance in the 5-15 Hz range with lower polarizability, although bipolar cells had higher polarizability. Layer 5 PC demonstrated the highest entrainment close to 10 Hz, which decayed with frequency. In contrast, inhibitory neurons entrainment increased with frequency, reaching levels akin to PC. Results from simplified models could replicate phase preferences, while amplitudes tended to follow opposite trends in PC.Significance. tACS-induced membrane polarization is frequency-dependent, revealing observable resonance behavior. Whilst optimal phase entrainment of sustained activity is achieved in PC when tACS frequency matches endogenous activity, inhibitory neurons tend to be entrained at higher frequencies. Consequently, our results highlight the potential for precise, cell-specific targeting for tACS.

Contributions of action potentials to scalp EEG: Theory and biophysical simulations

Differences in the apparent 1/f component of neural power spectra require correction depending on the underlying neural mechanisms, which remain incompletely understood. Past studies suggest that neuronal spiking produces broadband signals and shapes the spectral trend of invasive macroscopic recordings, but it is unclear to what extent action potentials (APs) influence scalp EEG. Here, we combined biophysical simulations with statistical modelling to examine the amplitude and spectral content of scalp potentials generated by the electric fields from spiking activity. In physiological parameter regimes, we found that APs contribute negligibly to the EEG spectral trend. Consistent with this, comparing our biophysical simulations with previously published data from pharmacologically paralyzed subjects suggested that the EEG spectral trend can be explained by a combination of synaptic timescales and electromyogram contamination. We also modelled rhythmic EEG generation, finding that APs can generate detectable narrowband power between approximately 60 and 1000 Hz, reaching frequencies much faster than would be possible from synaptic currents. Finally, we show that different spectral detrending strategies are required for AP generated oscillations compared to synaptically generated oscillations, suggesting that existing detrending methods for EEG spectra need to be modified for high frequency signals.

Heterogeneous and higher-order cortical connectivity undergirds efficient, robust, and reliable neural codes

We hypothesized that the heterogeneous architecture of biological neural networks provides a substrate to regulate the well-known tradeoff between robustness and efficiency, thereby allowing different subpopulations of the same network to optimize for different objectives. To distinguish between subpopulations, we developed a metric based on the mathematical theory of simplicial complexes that captures the complexity of their connectivity by contrasting its higher-order structure to a random control and confirmed its relevance in several openly available connectomes. Using a biologically detailed cortical model and an electron microscopic dataset, we showed that subpopulations with low simplicial complexity exhibit efficient activity. Conversely, subpopulations of high simplicial complexity play a supporting role in boosting the reliability of the network as a whole, softening the robustness-efficiency tradeoff. Crucially, we found that both types of subpopulations can and do coexist within a single connectome in biological neural networks, due to the heterogeneity of their connectivity.

Self-organized and self-sustained ensemble activity patterns in simulation of mouse primary motor cortex

The idea of self-organized signal processing in the cerebral cortex has become a focus of research since Beggs and Plentz1 reported avalanches in local field potential recordings from organotypic cultures and acute slices of rat somatosensory cortex. How the cortex intrinsically organizes signals remains unknown. A current hypothesis was proposed by the condensed matter physicists Bak, Tang, and Wiesenfeld2 when they conjectured that if neuronal avalanche activity followed inverse power law distributions, then brain activity may be set around phase transitions within self-organized signals. We asked if we would observe self-organized signals in an isolated slice of our data driven detailed simulation of the mouse primary motor cortex? If we did, would we observe avalanches with power-law distributions in size and duration and what would they look like? Our results demonstrate that a brief unstructured stimulus (100ms, 57μA current) to a small subset of neurons (about 181 of more than 10,000) in a simulated mouse primary motor cortex slice results in self-organized and self-sustained avalanches with power-law size and duration distributions and values similar to those reported from in vivo and in vitro experiments. We observed 4 cross-layer and cross-neuron population patterns, 3 of which displayed a dominant rhythmic component. Avalanches were each composed of one or more of the 4 population patterns.

Computational Generation of Long-range Axonal Morphologies

Long-range axons are fundamental to brain connectivity and functional organization, enabling communication between different brain regions. Recent advances in experimental techniques have yielded a substantial number of whole-brain axonal reconstructions. While previous computational generative models of neurons have predominantly focused on dendrites, generating realistic axonal morphologies is more challenging due to their distinct targeting. In this study, we present a novel algorithm for axon synthesis that combines algebraic topology with the Steiner tree algorithm, an extension of the minimum spanning tree, to generate both the local and long-range compartments of axons. We demonstrate that our computationally generated axons closely replicate experimental data in terms of their morphological properties. This approach enables the generation of biologically accurate long-range axons that span large distances and connect multiple brain regions, advancing the digital reconstruction of the brain. Ultimately, our approach opens up new possibilities for large-scale in-silico simulations, advancing research into brain function and disorders.

The NeuroML ecosystem for standardized multi-scale modeling in neuroscience

Data-driven models of neurons and circuits are important for understanding how the properties of membrane conductances, synapses, dendrites, and the anatomical connectivity between neurons generate the complex dynamical behaviors of brain circuits in health and disease. However, the inherent complexity of these biological processes makes the construction and reuse of biologically detailed models challenging. A wide range of tools have been developed to aid their construction and simulation, but differences in design and internal representation act as technical barriers to those who wish to use data-driven models in their research workflows. NeuroML, a model description language for computational neuroscience, was developed to address this fragmentation in modeling tools. Since its inception, NeuroML has evolved into a mature community standard that encompasses a wide range of model types and approaches in computational neuroscience. It has enabled the development of a large ecosystem of interoperable open-source software tools for the creation, visualization, validation, and simulation of data-driven models. Here, we describe how the NeuroML ecosystem can be incorporated into research workflows to simplify the construction, testing, and analysis of standardized models of neural systems, and supports the FAIR (Findability, Accessibility, Interoperability, and Reusability) principles, thus promoting open, transparent and reproducible science.

Combined statistical-biophysical modeling links ion channel genes to physiology of cortical neuron types

Neural cell types have classically been characterized by their anatomy and electrophysiology. More recently, single-cell transcriptomics has enabled an increasingly fine genetically defined taxonomy of cortical cell types, but the link between the gene expression of individual cell types and their physiological and anatomical properties remains poorly understood. Here, we develop a hybrid modeling approach to bridge this gap. Our approach combines statistical and mechanistic models to predict cells' electrophysiological activity from their gene expression pattern. To this end, we fit biophysical Hodgkin-Huxley-based models for a wide variety of cortical cell types using simulation-based inference, while overcoming the challenge posed by the mismatch between the mathematical model and the data. Using multimodal Patch-seq data, we link the estimated model parameters to gene expression using an interpretable sparse linear regression model. Our approach recovers specific ion channel gene expressions as predictive of biophysical model parameters including ion channel densities, directly implicating their mechanistic role in determining neural firing.

Enabling electric field model of microscopically realistic brain

BACKGROUND

Modeling brain stimulation at the microscopic scale may reveal new paradigms for various stimulation modalities.

OBJECTIVE

We present the largest map to date of extracellular electric field distributions within a layer L2/L3 mouse primary visual cortex brain sample. This was enabled by the automated analysis of serial section electron microscopy images with improved handling of image defects, covering a volume of 250 × 140 × 90 μm³.

METHODS

The map was obtained by applying a uniform brain stimulation electric field at three different polarizations and accurately computing microscopic field perturbations using the boundary element fast multipole method. We used the map to identify the effect of microscopic field perturbations on the activation thresholds of individual neurons. Previous relevant studies modeled a macroscopically homogeneous cortical volume.

RESULT

Our result shows that the microscopic field perturbations - an 'electric field spatial noise' with a mean value of zero - only modestly influence the macroscopically predicted stimulation field strengths necessary for neuronal activation. The thresholds do not change by more than 10 % on average.

CONCLUSION

Under the stated limitations and assumptions of our method, this result essentially justifies the conventional theory of "invisible" neurons embedded in a macroscopic brain model for transcranial magnetic and transcranial electrical stimulation. However, our result is solely sample-specific and is only relevant to this relatively small sample with 396 neurons. It largely neglects the effect of the microcapillary network. Furthermore, we only considered the uniform impressed field and a single-pulse stimulation time course.

Shapes of direct cortical responses vs. short-range axono-cortical evoked potentials: The effects of direct electrical stimulation applied to the human brain

OBJECTIVE

Direct cortical responses (DCR) and axono-cortical evoked potentials (ACEP) are generated by electrically stimulating the cortex either directly or indirectly through white matter pathways, potentially leading to different electrogenic processes. For ACEP, the slow conduction velocity of axons (median ≈ 4 m.s-1) is anticipated to induce a delay. For DCR, direct electrical stimulation (DES) of the cortex is expected to elicit additional cortical activity involving smaller and slower non-myelinated axons. We tried to validate these hypotheses.

METHODS

DES was administered either directly on the cortex or to white matter fascicles within the resection cavity, while recording DCR or ACEP at the cortical level in nine patients.

RESULTS

Short but significant delays (≈ 2 ms) were measurable for ACEP immediately following the initial component (≈ 7 ms). Subsequent activities (≈ 40 ms) exhibited notable differences between DCR and ACEP, suggesting the presence of additional cortical activities for DCR.

CONCLUSION

Distinctions between ACEPs and DCRs can be made based on a delay at the onset of early components and the dissimilarity in the shape of the later components (>40 ms after the DES artifact).

SIGNIFICANCE

The comparison of different types of evoked potentials allows to better understand the effects of DES.

The Neurobehavioral State hypothesis

Since the early attempts to understand the brain made by Greek philosophers more than 2000 years ago, one of the main questions in neuroscience has been how the brain perceives all the stimuli in the environment and uses this information to implement a response. Recent hypotheses of the neural code rely on the existence of an ideal observer, whether on specific areas of the cerebral cortex or distributed network composed of cortical and subcortical elements. The Neurobehavioral State hypothesis stipulates that neurons are in a quasi-stable state due to the dynamic interaction of their molecular components. This increases their computational capabilities and electrophysiological behavior further than a binary active/inactive state. Together, neuronal populations across the brain learn to identify and associate internal and external stimuli with actions and emotions. Furthermore, such associations can be stored through the regulation of neuronal components as new quasi-stable states. Using this framework, behavior arises as the result of the dynamic interaction between internal and external stimuli together with previously established quasi-stable states that delineate the behavioral response. Finally, the Neurobehavioral State hypothesis is firmly grounded on present evidence of the complex dynamics within the brain, from the molecular to the network level, and avoids the need for a central observer by proposing the brain configures itself through experience-driven associations.

Hemispheric Asymmetry of Intracortical Myelin Orientation in the Mouse Auditory Cortex

Communication sound processing in mouse AC is lateralized. Both left and right AC are highly specialised and differ in auditory stimulus representation, functional connectivity and field topography. Previous studies have highlighted intracortical functional circuits that explain hemispheric stimulus preference. However, the underlying microstructure remains poorly understood. In this study, we examine structural lateralization of AC on the basis of immunohistochemically stained and tissue-cleared adult mouse brains (n = 11). We found hemispheric asymmetries of intracortical myelination, most prominently in layer 2/3, which featured more intercolumnar connections in the right AC. Furthermore, we found a larger structural asymmetry in the right AC. We also investigated sex differences. In male mice, myelination direction in the right AC is tilted to the anterior side. This pattern is inverted in female mice. However, the spatial distribution of neuronal cell bodies in the left and right AC along the laminar axis of the cortex was remarkably symmetric in all samples. These results suggest that basic developmentally defined structures such as cortical columns remain untouched by lateral specialisation, but more plastic myelinated axons show diverse hemispheric asymmetries. These asymmetries may contribute to specialisation on lateralized tasks such as vocal communication processing or specialisation on spectral or temporal complexity of stimuli.

Modeling of Blood Flow Dynamics in Rat Somatosensory Cortex

Background: The cerebral microvasculature forms a dense network of interconnected blood vessels where flow is modulated partly by astrocytes. Increased neuronal activity stimulates astrocytes to release vasoactive substances at the endfeet, altering the diameters of connected vessels. Methods: Our study simulated the coupling between blood flow variations and vessel diameter changes driven by astrocytic activity in the rat somatosensory cortex. We developed a framework with three key components: coupling between the vasculature and synthesized astrocytic morphologies, a fluid dynamics model to compute flow in each vascular segment, and a stochastic process replicating the effect of astrocytic endfeet on vessel radii. Results: The model was validated against experimental flow values from the literature across cortical depths. We found that local vasodilation from astrocyte activity increased blood flow, especially in capillaries, exhibiting a layer-specific response in deeper cortical layers. Additionally, the highest blood flow variability occurred in capillaries, emphasizing their role in cerebral perfusion regulation. We discovered that astrocytic activity impacted blood flow dynamics in a localized, clustered manner, with most vascular segments influenced by two to three neighboring endfeet. Conclusions: These insights enhance our understanding of neurovascular coupling and guide future research on blood flow-related diseases.

What makes human cortical pyramidal neurons functionally complex

Humans exhibit unique cognitive abilities within the animal kingdom, but the neural mechanisms driving these advanced capabilities remain poorly understood. Human cortical neurons differ from those of other species, such as rodents, in both their morphological and physiological characteristics. Could the distinct properties of human cortical neurons help explain the superior cognitive capabilities of humans? Understanding this relationship requires a metric to quantify how neuronal properties contribute to the functional complexity of single neurons, yet no such standardized measure currently exists. Here, we propose the Functional Complexity Index (FCI), a generalized, deep learning-based framework to assess the input-output complexity of neurons. By comparing the FCI of cortical pyramidal neurons from different layers in rats and humans, we identified key morpho-electrical factors that underlie functional complexity. Human cortical pyramidal neurons were found to be significantly more functionally complex than their rat counterparts, primarily due to differences in dendritic membrane area and branching pattern, as well as density and nonlinearity of NMDA-mediated synaptic receptors. These findings reveal the structural-biophysical basis for the enhanced functional properties of human neurons.

Of mice and men: Dendritic architecture differentiates human from mice neuronal networks

The organizational principles that distinguish the human brain from other species have been a long-standing enigma in neuroscience. Focusing on the uniquely evolved human cortical layers 2 and 3, we computationally reconstruct the cortical architecture for mice and humans. We show that human pyramidal cells form highly complex networks, demonstrated by the increased number and simplex dimension compared to mice. This is surprising because human pyramidal cells are much sparser in the cortex. We show that the number and size of neurons fail to account for this increased network complexity, suggesting that another morphological property is a key determinant of network connectivity. Topological comparison of dendritic structure reveals much higher perisomatic (basal and oblique) branching density in human pyramidal cells. Using topological tools we quantitatively show that this neuronal structural property directly impacts network complexity, including the formation of a rich subnetwork structure. We conclude that greater dendritic complexity, a defining attribute of human L2 and 3 neurons, may provide the human cortex with enhanced computational capacity and cognitive flexibility.

The impact of Parkinson's disease on striatal network connectivity and corticostriatal drive: An in silico study

Striatum, the input stage of the basal ganglia, is important for sensory-motor integration, initiation and selection of behavior, as well as reward learning. Striatum receives glutamatergic inputs from mainly cortex and thalamus. In rodents, the striatal projection neurons (SPNs), giving rise to the direct and the indirect pathway (dSPNs and iSPNs, respectively), account for 95% of the neurons, and the remaining 5% are GABAergic and cholinergic interneurons. Interneuron axon terminals as well as local dSPN and iSPN axon collaterals form an intricate striatal network. Following chronic dopamine depletion as in Parkinson's disease (PD), both morphological and electrophysiological striatal neuronal features have been shown to be altered in rodent models. Our goal with this in silico study is twofold: (a) to predict and quantify how the intrastriatal network connectivity structure becomes altered as a consequence of the morphological changes reported at the single-neuron level and (b) to investigate how the effective glutamatergic drive to the SPNs would need to be altered to account for the activity level seen in SPNs during PD. In summary, we predict that the richness of the connectivity motifs in the striatal network is significantly decreased during PD while, at the same time, a substantial enhancement of the effective glutamatergic drive to striatum is present.

Predicting modular functions and neural coding of behavior from a synaptic wiring diagram

A long-standing goal in neuroscience is to understand how a circuit's form influences its function. Here, we reconstruct and analyze a synaptic wiring diagram of the larval zebrafish brainstem to predict key functional properties and validate them through comparison with physiological data. We identify modules of strongly connected neurons that turn out to be specialized for different behavioral functions, the control of eye and body movements. The eye movement module is further organized into two three-block cycles that support the positive feedback long hypothesized to underlie low-dimensional attractor dynamics in oculomotor control. We construct a neural network model based directly on the reconstructed wiring diagram that makes predictions for the cellular-resolution coding of eye position and neural dynamics. These predictions are verified statistically with calcium imaging-based neural activity recordings. This work demonstrates how connectome-based brain modeling can reveal previously unknown anatomical structure in a neural circuit and provide insights linking network form to function.

Beyond-local neural information processing in neuronal networks

While there is much knowledge about local neuronal circuitry, considerably less is known about how neuronal input is integrated and combined across neuronal networks to encode higher order brain functions. One challenge lies in the large number of complex neural interactions. Neural networks use oscillating activity for information exchange between distributed nodes. To better understand building principles underlying the observation of synchronized oscillatory activity in a large-scale network, we developed a reductionistic neuronal network model. Fundamental building principles are laterally and temporally interconnected virtual nodes (microcircuits), wherein each node was modeled as a local oscillator. By this building principle, the neuronal network model can integrate information in time and space. The simulation gives rise to a wave interference pattern that spreads over all simulated columns in form of a travelling wave. The model design stabilizes states of efficient information processing across all participating neuronal equivalents. Model-specific oscillatory patterns, generated by complex input stimuli, were similar to electrophysiological high-frequency signals that we could confirm in the primate visual cortex during a visual perception task. Important oscillatory model pre-runners, limitations and strength of our reductionistic model are discussed. Our simple scalable model shows unique integration properties and successfully reproduces a variety of biological phenomena such as harmonics, coherence patterns, frequency-speed relationships, and oscillatory activities. We suggest that our scalable model simulates aspects of a basic building principle underlying oscillatory, large-scale integration of information in small and large brains.

Evaluation and comparison of methods for neuronal parameter optimization using the Neuroptimus software framework

Finding optimal parameters for detailed neuronal models is a ubiquitous challenge in neuroscientific research. In recent years, manual model tuning has been gradually replaced by automated parameter search using a variety of different tools and methods. However, using most of these software tools and choosing the most appropriate algorithm for a given optimization task require substantial technical expertise, which prevents the majority of researchers from using these methods effectively. To address these issues, we developed a generic platform (called Neuroptimus) that allows users to set up neural parameter optimization tasks via a graphical interface, and to solve these tasks using a wide selection of state-of-the-art parameter search methods implemented by five different Python packages. Neuroptimus also offers several features to support more advanced usage, including the ability to run most algorithms in parallel, which allows it to take advantage of high-performance computing architectures. We used the common interface provided by Neuroptimus to conduct a detailed comparison of more than twenty different algorithms (and implementations) on six distinct benchmarks that represent typical scenarios in neuronal parameter search. We quantified the performance of the algorithms in terms of the best solutions found and in terms of convergence speed. We identified several algorithms, including covariance matrix adaptation evolution strategy and particle swarm optimization, that consistently, without any fine-tuning, found good solutions in all of our use cases. By contrast, some other algorithms including all local search methods provided good solutions only for the simplest use cases, and failed completely on more complex problems. We also demonstrate the versatility of Neuroptimus by applying it to an additional use case that involves tuning the parameters of a subcellular model of biochemical pathways. Finally, we created an online database that allows uploading, querying and analyzing the results of optimization runs performed by Neuroptimus, which enables all researchers to update and extend the current benchmarking study. The tools and analysis we provide should aid members of the neuroscience community to apply parameter search methods more effectively in their research.

An integrative data-driven model simulating C. elegans brain, body and environment interactions

The behavior of an organism is influenced by the complex interplay between its brain, body and environment. Existing data-driven models focus on either the brain or the body-environment. Here we present BAAIWorm, an integrative data-driven model of Caenorhabditis elegans, which consists of two submodels: the brain model and the body-environment model. The brain model was built by multicompartment models with realistic morphology, connectome and neural population dynamics based on experimental data. Simultaneously, the body-environment model used a lifelike body and a three-dimensional physical environment. Through the closed-loop interaction between the two submodels, BAAIWorm reproduced the realistic zigzag movement toward attractors observed in C. elegans. Leveraging this model, we investigated the impact of neural system structure on both neural activities and behaviors. Consequently, BAAIWorm can enhance our understanding of how the brain controls the body to interact with its surrounding environment.

Circuits and mechanisms for TMS-induced corticospinal waves: Connecting sensitivity analysis to the network graph

Transcranial magnetic stimulation (TMS) is a non-invasive, FDA-cleared treatment for neuropsychiatric disorders with broad potential for new applications, but the neural circuits that are engaged during TMS are still poorly understood. Recordings of neural activity from the corticospinal tract provide a direct readout of the response of motor cortex to TMS, and therefore a new opportunity to model neural circuit dynamics. The study goal was to use epidural recordings from the cervical spine of human subjects to develop a computational model of a motor cortical macrocolumn through which the mechanisms underlying the response to TMS, including direct and indirect waves, could be investigated. An in-depth sensitivity analysis was conducted to identify important pathways, and machine learning was used to identify common circuit features among these pathways. Sensitivity analysis identified neuron types that preferentially contributed to single corticospinal waves. Single wave preference could be predicted using the average connection probability of all possible paths between the activated neuron type and L5 pyramidal tract neurons (PTNs). For these activations, the total conduction delay of the shortest path to L5 PTNs determined the latency of the corticospinal wave. Finally, there were multiple neuron type activations that could preferentially modulate a particular corticospinal wave. The results support the hypothesis that different pathways of circuit activation contribute to different corticospinal waves with participation of both excitatory and inhibitory neurons. Moreover, activation of both afferents to the motor cortex as well as specific neuron types within the motor cortex initiated different I-waves, and the results were interpreted to propose the cortical origins of afferents that may give rise to certain I-waves. The methodology provides a workflow for performing computationally tractable sensitivity analyses on complex models and relating the results to the network structure to both identify and understand mechanisms underlying the response to acute stimulation.

Sparse Spike Feature Learning to Recognize Traceable Interictal Epileptiform Spikes

Interictal epileptiform spikes (spikes) and epileptogenic focus are strongly correlated. However, partial spikes are insensitive to epileptogenic focus, which restricts epilepsy neurosurgery. Therefore, identifying spike subtypes that are strongly associated with epileptogenic focus (traceable spikes) could facilitate their use as reliable signal sources for accurately tracing epileptogenic focus. However, the sparse firing phenomenon in the transmission of intracranial neuronal discharges leads to differences within spikes that cannot be observed visually. Therefore, neuro-electro-physiologists are unable to identify traceable spikes that could accurately locate epileptogenic focus. Herein, we propose a novel sparse spike feature learning method to recognize traceable spikes and extract discrimination information related to epileptogenic focus. First, a multilevel eigensystem feature representation was determined based on a multilevel feature representation module to express the intrinsic properties of a spike. Second, the sparse feature learning module expressed the sparse spike multi-domain context feature representation to extract sparse spike feature representations. Among them, a sparse spike encoding strategy was implemented to effectively simulate the sparse firing phenomenon for the accurate encoding of the activity of intracranial neurosources. The sensitivity of the proposed method was 97.1%, demonstrating its effectiveness and significant efficiency relative to other state-of-the-art methods.

Enabling Electric Field Model of Microscopically Realistic Brain

Modeling brain stimulation at the microscopic scale may reveal new paradigms for various stimulation modalities. We present the largest map to date of extracellular electric field distributions within a layer L2/L3 mouse primary visual cortex brain sample. This was enabled by the automated analysis of serial section electron microscopy images with improved handling of image defects, covering a volume of 250 × 140 × 90 μm 3 . The map was obtained by applying a uniform brain stimulation electric field at three different polarizations and accurately computing microscopic field perturbations using the boundary element fast multipole method. We used the map to identify the effect of microscopic field perturbations on the activation thresholds of individual neurons. Previous relevant studies modeled a macroscopically homogeneous cortical volume. Our result shows that the microscopic field perturbations - an 'electric field spatial noise' with a mean value of zero - only modestly influence the macroscopically predicted stimulation field strengths necessary for neuronal activation. The thresholds do not change by more than 10% on average. Under the stated limitations and assumptions of our method, this result justifies the conventional theory of "invisible" neurons embedded in a macroscopic brain model for transcranial magnetic and transcranial electrical stimulation. However, our result is solely sample-specific and largely neglects the effect of the microcapillary network. Furthermore, we only considered the uniform impressed field and a single- pulse stimulation time course.

Significance statement

This study is arguably the first attempt to model brain stimulation at the microscopic scale, enabled by automated analysis of modern scanning electron microscopy images of the brain. It concentrates on modeling microscopic perturbations of the extracellular electric field caused by the physical cell structure and is applicable to any type of brain stimulation.

Data availability statement

Post-processed cell CAD models (383, stl format), microcapillary CAD models (34, stl format), post-processed neuron morphologies (267, swc format), extracellular electric field and potential distributions at different polarizations (267x3, MATLAB format), *.ses projects files for biophysical modeling with Neuron software (267x2, Neuron format), and computed neuron activating thresholds at different conditions (267x8, Excel tables, without the sample polarization correction from Section 2.8) are made available online through BossDB , a volumetric open-source database for 3D and 4D neuroscience data.

Future projections for mammalian whole-brain simulations based on technological trends in related fields

Large-scale brain simulation allows us to understand the interaction of vast numbers of neurons having nonlinear dynamics to help understand the information processing mechanisms in the brain. The scale of brain simulations continues to rise as computer performance improves exponentially. However, a simulation of the human whole brain has not yet been achieved as of 2024 due to insufficient computational performance and brain measurement data. This paper examines technological trends in supercomputers, cell type classification, connectomics, and large-scale activity measurements relevant to whole-brain simulation. Based on these trends, we attempt to predict the feasible timeframe for mammalian whole-brain simulation. Our estimates suggest that mouse whole-brain simulation at the cellular level could be realized around 2034, marmoset around 2044, and human likely later than 2044.

Real-time multicompartment Hodgkin-Huxley neuron emulation on SoC FPGA

Advanced computational models and simulations to unravel the complexities of brain function have known a growing interest in recent years in the field of neurosciences, driven by significant technological progress in computing platforms. Multicompartment models, which capture the detailed morphological and functional properties of neural circuits, represent a significant advancement in this area providing more biological coherence than single compartment modeling. These models serve as a cornerstone for exploring the neural basis of sensory processing, learning paradigms, adaptive behaviors, and neurological disorders. Yet, the high complexity of these models presents a challenge for their real-time implementation, which is essential for exploring alternative therapies for neurological disorders such as electroceutics that rely on biohybrid interaction. Here, we present an accessible, user-friendly, and real-time emulator for multicompartment Hodgkin-Huxley neurons on SoC FPGA. Our system enables real-time emulation of multicompartment neurons while emphasizing cost-efficiency, flexibility, and ease of use. We showcase an implementation utilizing a technology that remains underrepresented in the current literature for this specific application. We anticipate that our system will contribute to the enhancement of computation platforms by presenting an alternative architecture for multicompartment computation. Additionally, it constitutes a step toward developing neuromorphic-based neuroprostheses for bioelectrical therapeutics through an embedded real-time platform running at a similar timescale to biological networks.

Rectified activity-dependent population plasticity implicates cortical adaptation for memory and cognitive functions

Cortical network undergoes rewiring everyday due to learning and memory events. To investigate the trends of population adaptation in neocortex overtime, we record cellular activity of large-scale cortical populations in response to neutral environments and conditioned contexts and identify a general intrinsic cortical adaptation mechanism, naming rectified activity-dependent population plasticity (RAPP). Comparing each adjacent day, the previously activated neurons reduce activity, but remain with residual potentiation, and increase population variability in proportion to their activity during previous recall trials. RAPP predicts both the decay of context-induced activity patterns and the emergence of sparse memory traces. Simulation analysis reveal that the local inhibitory connections might account for the residual potentiation in RAPP. Intriguingly, introducing the RAPP phenomenon in the artificial neural network show promising improvement in small sample size pattern recognition tasks. Thus, RAPP represents a phenomenon of cortical adaptation, contributing to the emergence of long-lasting memory and high cognitive functions.

Community-based reconstruction and simulation of a full-scale model of the rat hippocampus CA1 region

The CA1 region of the hippocampus is one of the most studied regions of the rodent brain, thought to play an important role in cognitive functions such as memory and spatial navigation. Despite a wealth of experimental data on its structure and function, it has been challenging to integrate information obtained from diverse experimental approaches. To address this challenge, we present a community-based, full-scale in silico model of the rat CA1 that integrates a broad range of experimental data, from synapse to network, including the reconstruction of its principal afferents, the Schaffer collaterals, and a model of the effects that acetylcholine has on the system. We tested and validated each model component and the final network model, and made input data, assumptions, and strategies explicit and transparent. The unique flexibility of the model allows scientists to potentially address a range of scientific questions. In this article, we describe the methods used to set up simulations to reproduce in vitro and in vivo experiments. Among several applications in the article, we focus on theta rhythm, a prominent hippocampal oscillation associated with various behavioral correlates and use our computer model to reproduce experimental findings. Finally, we make data, code, and model available through the hippocampushub.eu portal, which also provides an extensive set of analyses of the model and a user-friendly interface to facilitate adoption and usage. This community-based model represents a valuable tool for integrating diverse experimental data and provides a foundation for further research into the complex workings of the hippocampal CA1 region.

Structural and molecular characterization of astrocyte and vasculature connectivity in the mouse hippocampus and cortex

The relation of astrocytic endfeet to the vasculature plays a key functional role in the neuro-glia-vasculature unit. We characterize the spatial organization of astrocytes and the structural aspects that facilitate their involvement in molecular exchanges. Using double transgenic mice, we performed co-immunostaining, confocal microscopy, and three-dimensional digital segmentation to investigate the biophysical and molecular organization of astrocytes and their intricate endfoot network at the micrometer level in the isocortex and hippocampus. The results showed that hippocampal astrocytes had smaller territories, reduced endfoot dimensions, and fewer contacts with blood vessels compared with those in the isocortex. Additionally, we found that both connexins 43 and 30 have a higher density in the endfoot and the former is overexpressed relative to the latter. However, due to the limitations of the method, further studies are needed to determine the exact localization on the endfoot. The quantitative information obtained in this study will be useful for modeling the interactions of astrocytes with the vasculature.

Gather your neurons and model together: Community times ahead

Bottom-up, data-driven, large-scale models provide a mechanistic understanding of neuronal functions. A new study in PLOS Biology builds a biologically realistic model of the rodent CA1 region that aims to become an accessible tool for the whole hippocampal community.

Beyond Barrels: Diverse Thalamocortical Projection Motifs in the Mouse Ventral Posterior Complex

Thalamocortical pathways from the rodent ventral posterior (VP) thalamic complex to the somatosensory cerebral cortex areas are a key model in modern neuroscience. However, beyond the intensively studied projection from medial VP (VPM) to the primary somatosensory area (S1), the wiring of these pathways remains poorly characterized. We combined micropopulation tract-tracing and single-cell transfection experiments to map the pathways arising from different portions of the VP complex in male mice. We found that pathways originating from different VP regions show differences in area/lamina arborization pattern and axonal varicosity size. Neurons from the rostral VPM subnucleus innervate trigeminal S1 in point-to-point fashion. In contrast, a caudal VPM subnucleus innervates heavily and topographically second somatosensory area (S2), but not S1. Neurons in a third, intermediate VPM subnucleus innervate through branched axons both S1 and S2, with markedly different laminar patterns in each area. A small anterodorsal subnucleus selectively innervates dysgranular S1. The parvicellular VPM subnucleus selectively targets the insular cortex and adjacent portions of S1 and S2. Neurons in the rostral part of the lateral VP nucleus (VPL) innervate spinal S1, while caudal VPL neurons simultaneously target S1 and S2. Rostral and caudal VP nuclei show complementary patterns of calcium-binding protein expression. In addition to the cortex, neurons in caudal VP subnuclei target the sensorimotor striatum. Our finding of a massive projection from VP to S2 separate from the VP projections to S1 adds critical anatomical evidence to the notion that different somatosensory submodalities are processed in parallel in S1 and S2.

Layer-specific dynamics of local field potentials in monkey V1 during electrical stimulation

The mammalian neocortex, organized into six cellular layers or laminae, forms a cortical network within layers. Layer specific computations are crucial for sensory processing of visual stimuli within primary visual cortex. Laminar recordings of local field potentials (LFPs) are a powerful tool to study neural activity within cortical layers. Electric brain stimulation is widely used in basic neuroscience and in a large range of clinical applications. However, the layer-specific effects of electric stimulation on LFPs remain unclear. To address this gap, we conducted laminar LFP recordings of the primary visual cortex in monkeys while presenting a flash visual stimulus. Simultaneously, we applied a low frequency sinusoidal current to the occipital lobe with offset frequency to the flash stimulus repetition rate. We analyzed the modulation of visual-evoked potentials with respect to the applied phase of the electric stimulation. Our results reveal that only the deeper layers, but not the superficial layers, show phase-dependent changes in LFP components with respect to the applied current. Employing a cortical column model, we demonstrate that these in vivo observations can be explained by phase-dependent changes in the driving force within neurons of deeper layers. Our findings offer crucial insight into the selective modulation of cortical layers through electric stimulation, thus advancing approaches for more targeted neuromodulation.

Synchrony dynamics underlie irregular neocortical spiking

Cortical neurons are characterized by their variable spiking patterns. We challenge prevalent theories for the origin of spiking variability. We examine the specific hypothesis that cortical synchrony drives spiking variability in vivo . Using dynamic clamp, we demonstrate that intrinsic neuronal properties do not contribute substantially to spiking variability, but rather spiking variability emerges from weakly synchronous network drive. With large-scale electrophysiology we quantify the degree of synchrony and its time scale in cortical networks in vivo . We demonstrate that physiological levels of synchrony are sufficient to generate irregular responses found in vivo . Further, this synchrony shifts over timescales ranging from 25 to 200 ms, depending on the presence of external sensory input. Such shifts occur when the network moves from spontaneous to driven modes, leading naturally to a decline in response variability as observed across cortical areas. Finally, while individual neurons exhibit reliable responses to physiological drive, different neurons respond in a distinct fashion according to their intrinsic properties, contributing to stable synchrony across the neural network.

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

More than 20 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 (Kcnt1Y777H) 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 Kcnt1Y777H 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 Kcnt1Y777H 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 Kcnt1Y777H 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.

Multi-scale spiking network model of human cerebral cortex

Although the structure of cortical networks provides the necessary substrate for their neuronal activity, the structure alone does not suffice to understand the activity. Leveraging the increasing availability of human data, we developed a multi-scale, spiking network model of human cortex to investigate the relationship between structure and dynamics. In this model, each area in one hemisphere of the Desikan-Killiany parcellation is represented by a $1\,\mathrm{mm^{2}}$ column with a layered structure. The model aggregates data across multiple modalities, including electron microscopy, electrophysiology, morphological reconstructions, and diffusion tensor imaging, into a coherent framework. It predicts activity on all scales from the single-neuron spiking activity to the area-level functional connectivity. We compared the model activity with human electrophysiological data and human resting-state functional magnetic resonance imaging (fMRI) data. This comparison reveals that the model can reproduce aspects of both spiking statistics and fMRI correlations if the inter-areal connections are sufficiently strong. Furthermore, we study the propagation of a single-spike perturbation and macroscopic fluctuations through the network. The open-source model serves as an integrative platform for further refinements and future in silico studies of human cortical structure, dynamics, and function.

Reconciliation of weak pairwise spike-train correlations and highly coherent local field potentials across space

Multi-electrode arrays covering several square millimeters of neural tissue provide simultaneous access to population signals such as extracellular potentials and spiking activity of one hundred or more individual neurons. The interpretation of the recorded data calls for multiscale computational models with corresponding spatial dimensions and signal predictions. Multi-layer spiking neuron network models of local cortical circuits covering about $1\,{\text{mm}^{2}}$ have been developed, integrating experimentally obtained neuron-type-specific connectivity data and reproducing features of observed in-vivo spiking statistics. Local field potentials can be computed from the simulated spiking activity. We here extend a local network and local field potential model to an area of $4\times 4\,{\text{mm}^{2}}$, preserving the neuron density and introducing distance-dependent connection probabilities and conduction delays. We find that the upscaling procedure preserves the overall spiking statistics of the original model and reproduces asynchronous irregular spiking across populations and weak pairwise spike-train correlations in agreement with experimental recordings from sensory cortex. Also compatible with experimental observations, the correlation of local field potential signals is strong and decays over a distance of several hundred micrometers. Enhanced spatial coherence in the low-gamma band around $50\,\text{Hz}$ may explain the recent report of an apparent band-pass filter effect in the spatial reach of the local field potential.