Sparse recurrent excitatory connectivity in the microcircuit of the adult mouse and human cortex

Spine loss in depression impairs dendritic signal integration in human cortical microcircuit models

Major depressive disorder (depression) is associated with altered dendritic structure and function of cortical pyramidal neurons, due to decreased inhibition from somatostatin (SST) interneurons and loss of spines and associated synapses, as indicated in postmortem human studies. Dendrites mediate signal processing through synaptic integration and nonlinear properties including backpropagating action potentials and dendritic Na+ spikes that enhance the neuron's computational power. However, it is currently unclear how depression-related dendritic changes impact signal integration. Here, we integrated human neuronal data of active dendritic properties and spine loss in depression into detailed computational models of human cortical microcircuits. We show that spine loss dampens signal response, worsening signal detection impairment than due to reduced SST interneuron inhibition alone. Furthermore, altered intrinsic properties due to spine loss abolished nonlinear dendritic signal integration and impaired recurrent microcircuit activity. Our study mechanistically links cellular changes in depression to impaired dendritic processing in human cortical microcircuits.

Neural heterogeneity enhances reliable neural information processing: Local sensitivity and globally input-slaved transient dynamics

Cortical neuronal activity varies over time and across repeated trials, yet consistently represents stimulus features. The dynamical mechanism underlying this reliable representation and computation remains elusive. This study uncovers a mechanism for reliable neural information processing, leveraging a biologically plausible network model incorporating neural heterogeneity. First, we investigate neuronal timescale diversity, revealing that it disrupts intrinsic coherent spatiotemporal patterns, induces firing rate heterogeneity, enhances local responsive sensitivity, and aligns network activity closely with input. The system exhibits globally input-slaved transient dynamics, essential for reliable neural information processing. Other neural heterogeneities, such as nonuniform input connections, spike threshold heterogeneity, and network in-degree heterogeneity, play similar roles, highlighting the importance of neural heterogeneity in shaping consistent stimulus representation. This mechanism offers a potentially general framework for understanding neural heterogeneity in reliable computation and informs the design of reservoir computing models endowed with liquid wave reservoirs for neuromorphic computing.

Computational protocol for modeling and analyzing synaptic dynamics using SRPlasticity

Transient changes in synaptic strength, known as short-term plasticity (STP), play a fundamental role in neuronal communication. Here, we present a protocol for using SRPlasticity, a software package that implements a computational model of STP. SRPlasticity supports automatic characterization of electrophysiological data and simulation of synaptic responses. We describe steps for installing and utilizing SRPlasticity, preprocessing data, fitting models, and simulating responses. We then detail procedures for analyzing spike response plasticity (SRP) model parameters to infer functional groupings of STP. For complete details on the use and execution of this protocol, please refer to Rossbroich et al.1 and Beninger et al.2.

Subthreshold moment analysis of neuronal populations driven by synchronous synaptic inputs

Even when driven by the same stimulus, neuronal responses are well-known to exhibit a striking level of spiking variability. In-vivo electrophysiological recordings also reveal a surprisingly large degree of variability at the subthreshold level. In prior work, we considered biophysically relevant neuronal models to account for the observed magnitude of membrane voltage fluctuations. We found that accounting for these fluctuations requires weak but nonzero synchrony in the spiking activity, in amount that are consistent with experimentally measured spiking correlations. Here we investigate whether such synchrony can explain additional statistical features of the measured neural activity, including neuronal voltage covariability and voltage skewness. Addressing this question involves conducting a generalized moment analysis of conductance-based neurons in response to input drives modeled as correlated jump processes. Technically, we perform such an analysis using fixed-point techniques from queuing theory that are applicable in the stationary regime of activity. We found that weak but nonzero synchrony can consistently explain the experimentally reported voltage covariance and skewness. This confirms the role of synchrony as a primary driver of cortical variability and supports that physiological neural activity emerges as a population-level phenomenon, especially in the spontaneous regime.

Subthreshold variability of neuronal populations driven by synchronous synaptic inputs

Even when driven by the same stimulus, neuronal responses are well-known to exhibit a striking level of spiking variability. In-vivo electrophysiological recordings also reveal a surprisingly large degree of variability at the subthreshold level. In prior work, we considered biophysically relevant neuronal models to account for the observed magnitude of membrane voltage fluctuations. We found that accounting for these fluctuations requires weak but nonzero synchrony in the spiking activity, in amount that are consistent with experimentally measured spiking correlations. Here we investigate whether such synchrony can explain additional statistical features of the measured neural activity, including neuronal voltage covariability and voltage skewness. Addressing this question involves conducting a generalized moment analysis of conductance-based neurons in response to input drives modeled as correlated jump processes. Technically, we perform such an analysis using fixed-point techniques from queuing theory that are applicable in the stationary regime of activity. We found that weak but nonzero synchrony can consistently explain the experimentally reported voltage covariance and skewness. This confirms the role of synchrony as a primary driver of cortical variability and supports that physiological neural activity emerges as a population-level phenomenon, especially in the spontaneous regime.

Efficient coding in biophysically realistic excitatory-inhibitory spiking networks

The principle of efficient coding posits that sensory cortical networks are designed to encode maximal sensory information with minimal metabolic cost. Despite the major influence of efficient coding in neuroscience, it has remained unclear whether fundamental empirical properties of neural network activity can be explained solely based on this normative principle. Here, we derive the structural, coding, and biophysical properties of excitatory-inhibitory recurrent networks of spiking neurons that emerge directly from imposing that the network minimizes an instantaneous loss function and a time-averaged performance measure enacting efficient coding. We assumed that the network encodes a number of independent stimulus features varying with a time scale equal to the membrane time constant of excitatory and inhibitory neurons. The optimal network has biologically plausible biophysical features, including realistic integrate-and-fire spiking dynamics, spike-triggered adaptation, and a non-specific excitatory external input. The excitatory-inhibitory recurrent connectivity between neurons with similar stimulus tuning implements feature-specific competition, similar to that recently found in visual cortex. Networks with unstructured connectivity cannot reach comparable levels of coding efficiency. The optimal ratio of excitatory vs inhibitory neurons and the ratio of mean inhibitory-to-inhibitory vs excitatory-to-inhibitory connectivity are comparable to those of cortical sensory networks. The efficient network solution exhibits an instantaneous balance between excitation and inhibition. The network can perform efficient coding even when external stimuli vary over multiple time scales. Together, these results suggest that key properties of biological neural networks may be accounted for by efficient coding.

Augmenting flexibility: mutual inhibition between inhibitory neurons expands functional diversity

Recent advances in microcircuit analysis of nervous systems have revealed a plethora of mutual connections between inhibitory interneurons across many different species and brain regions. The abundance of these mutual connections has not been fully explained. Strikingly, we show that neural circuits with mutually inhibitory connections are able to rapidly and flexibly switch between distinct functions. That is, multiple functions coexist for a single set of synaptic weights. Here, we develop a theoretical framework to explain how inhibitory recurrent circuits give rise to this flexibility and show that mutual inhibition doubles the number of cusp bifurcations in small neural circuits. As a concrete example, we study a class of functional motifs we call coupled recurrent inhibitory and recurrent excitatory loops (CRIRELs). These CRIRELs have the advantage of being both multi-functional and controllable, performing a plethora of functions, including decisions, memory, toggle, and so forth. Finally, we demonstrate how mutual inhibition maximizes storage capacity for larger networks.

Bridging Neuroscience and AI: Environmental Enrichment as a model for forward knowledge transfer in continual learning

Continual learning (CL) refers to an agent's capability to learn from a continuous stream of data and transfer knowledge without forgetting old information. One crucial aspect of CL is forward transfer, i.e., improved and faster learning on a new task by leveraging information from prior knowledge. While this ability comes naturally to biological brains, it poses a significant challenge for artificial intelligence (AI). Here, we suggest that environmental enrichment (EE) can be used as a biological model for studying forward transfer, inspiring human-like AI development. EE refers to animal studies that enhance cognitive, social, motor, and sensory stimulation and is a model for what, in humans, is referred to as 'cognitive reserve'. Enriched animals show significant improvement in learning speed and performance on new tasks, typically exhibiting forward transfer. We explore anatomical, molecular, and neuronal changes post-EE and discuss how artificial neural networks (ANNs) can be used to predict neural computation changes after enriched experiences. Finally, we provide a synergistic way of combining neuroscience and AI research that paves the path toward developing AI capable of rapid and efficient new task learning.

Human hippocampal CA3 uses specific functional connectivity rules for efficient associative memory

Our brain has remarkable computational power, generating sophisticated behaviors, storing memories over an individual's lifetime, and producing higher cognitive functions. However, little of our neuroscience knowledge covers the human brain. Is this organ truly unique, or is it a scaled version of the extensively studied rodent brain? Combining multicellular patch-clamp recording with expansion-based superresolution microscopy and full-scale modeling, we determined the cellular and microcircuit properties of the human hippocampal CA3 region, a fundamental circuit for memory storage. In contrast to neocortical networks, human hippocampal CA3 displayed sparse connectivity, providing a circuit architecture that maximizes associational power. Human synapses showed unique reliability, high precision, and long integration times, exhibiting both species- and circuit-specific properties. Together with expanded neuronal numbers, these circuit characteristics greatly enhanced the memory storage capacity of CA3. Our results reveal distinct microcircuit properties of the human hippocampus and begin to unravel the inner workings of our most complex organ.

Efficient coding in biophysically realistic excitatory-inhibitory spiking networks

The principle of efficient coding posits that sensory cortical networks are designed to encode maximal sensory information with minimal metabolic cost. Despite the major influence of efficient coding in neuroscience, it has remained unclear whether fundamental empirical properties of neural network activity can be explained solely based on this normative principle. Here, we derive the structural, coding, and biophysical properties of excitatory-inhibitory recurrent networks of spiking neurons that emerge directly from imposing that the network minimizes an instantaneous loss function and a time-averaged performance measure enacting efficient coding. We assumed that the network encodes a number of independent stimulus features varying with a time scale equal to the membrane time constant of excitatory and inhibitory neurons. The optimal network has biologically-plausible biophysical features, including realistic integrate-and-fire spiking dynamics, spike-triggered adaptation, and a non-specific excitatory external input. The excitatory-inhibitory recurrent connectivity between neurons with similar stimulus tuning implements feature-specific competition, similar to that recently found in visual cortex. Networks with unstructured connectivity cannot reach comparable levels of coding efficiency. The optimal ratio of excitatory vs inhibitory neurons and the ratio of mean inhibitory-to-inhibitory vs excitatory-to-inhibitory connectivity are comparable to those of cortical sensory networks. The efficient network solution exhibits an instantaneous balance between excitation and inhibition. The network can perform efficient coding even when external stimuli vary over multiple time scales. Together, these results suggest that key properties of biological neural networks may be accounted for by efficient coding.

Harnessing the potential of human induced pluripotent stem cells, functional assays and machine learning for neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) affect 4.7% of the global population and are associated with delays in brain development and a spectrum of impairments that can lead to lifelong disability and even mortality. Identification of biomarkers for accurate diagnosis and medications for effective treatment are lacking, in part due to the historical use of preclinical model systems that do not translate well to the clinic for neurological disorders, such as rodents and heterologous cell lines. Human-induced pluripotent stem cells (hiPSCs) are a promising in vitro system for modeling NDDs, providing opportunities to understand mechanisms driving NDDs in human neurons. Functional assays, including patch clamping, multielectrode array, and imaging-based assays, are popular tools employed with hiPSC disease models for disease investigation. Recent progress in machine learning (ML) algorithms also presents unprecedented opportunities to advance the NDD research process. In this review, we compare two-dimensional and three-dimensional hiPSC formats for disease modeling, discuss the applications of functional assays, and offer insights on incorporating ML into hiPSC-based NDD research and drug screening.

Principles of visual cortex excitatory microcircuit organization

Synapse-specific connectivity and dynamics determine microcircuit function but are challenging to explore with classic paired recordings due to their low throughput. We therefore implemented optomapping, a ∼100-fold faster two-photon optogenetic method. In mouse primary visual cortex (V1), we optomapped 30,454 candidate inputs to reveal 1,790 excitatory inputs to pyramidal, basket, and Martinotti cells. Across these cell types, log-normal distribution of synaptic efficacies emerged as a principle. For pyramidal cells, optomapping reproduced the canonical circuit but unexpectedly uncovered that the excitation of basket cells concentrated to layer 5 and that of Martinotti cells dominated in layer 2/3. The excitation of basket cells was stronger and reached farther than the excitation of pyramidal cells, which may promote stability. Short-term plasticity surprisingly depended on cortical layer in addition to target cell. Finally, optomapping revealed an overrepresentation of shared inputs for interconnected layer-6 pyramidal cells. Thus, by resolving the throughput problem, optomapping uncovered hitherto unappreciated principles of V1 structure.

Therapeutic dose prediction of α5-GABA receptor modulation from simulated EEG of depression severity

Treatment for major depressive disorder (depression) often has partial efficacy and a large portion of patients are treatment resistant. Recent studies implicate reduced somatostatin (SST) interneuron inhibition in depression, and new pharmacology boosting this inhibition via positive allosteric modulators of α5-GABAA receptors (α5-PAM) offers a promising effective treatment. However, testing the effect of α5-PAM on human brain activity is limited, meriting the use of detailed simulations. We utilized our previous detailed computational models of human depression microcircuits with reduced SST interneuron inhibition and α5-PAM effects, to simulate EEG of individual microcircuits across depression severity and α5-PAM doses. We developed machine learning models that predicted optimal dose from EEG with high accuracy and recovered microcircuit activity and EEG. This study provides dose prediction models for α5-PAM administration based on EEG biomarkers of depression severity. Given limitations in doing the above in the living human brain, the results and tools we developed will facilitate translation of α5-PAM treatment to clinical use.

Neural dynamics and seizure correlations: Insights from neural mass models in a Tetanus Toxin rat model of epilepsy

This study focuses on the use of a neural mass model to investigate potential relationships between functional connectivity and seizure frequency in epilepsy. We fitted a three-layer neural mass model of a cortical column to intracranial EEG (iEEG) data from a Tetanus Toxin rat model of epilepsy, which also included responses to periodic electrical stimulation. Our results show that some of the connectivity weights between different neural populations correlate significantly with the number of seizures each day, offering valuable insights into the dynamics of neural circuits during epileptogenesis. We also simulated single-pulse electrical stimulation of the neuronal populations to observe their responses after the connectivity weights were optimized to fit background (non-seizure) EEG data. The recovery time, defined as the time from stimulation until the membrane potential returns to baseline, was measured as a representation of the critical slowing down phenomenon observed in nonlinear systems operating near a bifurcation boundary. The results revealed that recovery times in the responses of the computational model fitted to the EEG data were longer during 5 min periods preceding seizures compared to 1 hr before seizures in four out of six rats. Analysis of the iEEG recorded in response to electrical stimulation revealed results similar to the computational model in four out of six rats. This study supports the potential use of this computational model as a model-based biomarker for seizure prediction when direct electrical stimulation to the brain is not feasible.

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

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

Structural influences on synaptic plasticity: The role of presynaptic connectivity in the emergence of E/I co-tuning

Cortical neurons are versatile and efficient coding units that develop strong preferences for specific stimulus characteristics. The sharpness of tuning and coding efficiency is hypothesized to be controlled by delicately balanced excitation and inhibition. These observations suggest a need for detailed co-tuning of excitatory and inhibitory populations. Theoretical studies have demonstrated that a combination of plasticity rules can lead to the emergence of excitation/inhibition (E/I) co-tuning in neurons driven by independent, low-noise signals. However, cortical signals are typically noisy and originate from highly recurrent networks, generating correlations in the inputs. This raises questions about the ability of plasticity mechanisms to self-organize co-tuned connectivity in neurons receiving noisy, correlated inputs. Here, we study the emergence of input selectivity and weight co-tuning in a neuron receiving input from a recurrent network via plastic feedforward connections. We demonstrate that while strong noise levels destroy the emergence of co-tuning in the readout neuron, introducing specific structures in the non-plastic pre-synaptic connectivity can re-establish it by generating a favourable correlation structure in the population activity. We further show that structured recurrent connectivity can impact the statistics in fully plastic recurrent networks, driving the formation of co-tuning in neurons that do not receive direct input from other areas. Our findings indicate that the network dynamics created by simple, biologically plausible structural connectivity patterns can enhance the ability of synaptic plasticity to learn input-output relationships in higher brain areas.

Coordinated changes in a cortical circuit sculpt effects of novelty on neural dynamics

Recent studies have found dramatic cell-type-specific responses to stimulus novelty, highlighting the importance of analyzing the cortical circuitry at this granularity to understand brain function. Although initial work characterized activity by cell type, the alterations in cortical circuitry due to interacting novelty effects remain unclear. We investigated circuit mechanisms underlying the observed neural dynamics in response to novel stimuli using a large-scale public dataset of electrophysiological recordings in behaving mice and a population network model. The model was constrained by multi-patch synaptic physiology and electron microscopy data. We found generally weaker connections under novel stimuli, with shifts in the balance between somatostatin (SST) and vasoactive intestinal polypeptide (VIP) populations and increased excitatory influences on parvalbumin (PV) and SST populations. These findings systematically characterize how cortical circuits adapt to stimulus novelty.

Modulation of large rhythmic depolarizations in human large basket cells by norepinephrine and acetylcholine

Rhythmic brain activity is critical to many brain functions and is sensitive to neuromodulation, but so far very few studies have investigated this activity on the cellular level in vitro in human brain tissue samples. This study reveals and characterizes a novel rhythmic network activity in the human neocortex. Using intracellular patch-clamp recordings of human cortical neurons, we identify large rhythmic depolarizations (LRDs) driven by glutamate release but not by GABA. These LRDs are intricate events made up of multiple depolarizing phases, occurring at ~0.3 Hz, have large amplitudes and long decay times. Unlike human tissue, rat neocortex layers 2/3 exhibit no such activity under identical conditions. LRDs are mainly observed in a subset of L2/3 interneurons that receive substantial excitatory inputs and are likely large basket cells based on their morphology. LRDs are highly sensitive to norepinephrine (NE) and acetylcholine (ACh), two neuromodulators that affect network dynamics. NE increases LRD frequency through β-adrenergic receptor activity while ACh decreases it via M4 muscarinic receptor activation. Multi-electrode array recordings show that NE enhances and synchronizes oscillatory network activity, whereas ACh causes desynchronization. Thus, NE and ACh distinctly modulate LRDs, exerting specific control over human neocortical activity.

Reduced Synaptic Heterogeneity in a Tetanus Toxin Model of Epilepsy: Insights from Computational Modeling

A neural mass model was used to assess connectivity strength across diverse populations by fitting the model to background EEG data obtained from a Tetanus Toxin rat model of epilepsy. Our findings reveal a notable decline in the variability of estimated parameters when using EEG data recorded from rats in the Tetanus Toxin group compared with the control group. A detailed comparison of standard deviations in estimated parameters between day 1 and day 20 recordings, coinciding with a heightened number of seizures, underscores the impact of Tetanus Toxin on diminishing synaptic strength variability across recordings. This study supports electrophysiological studies suggesting that epileptogenesis induces a reduction in biophysical heterogeneity, potentially leading to an increase in network synchrony associated with epilepsy. Furthermore, our computational model establishes a foundation for future explorations of the implications of this diminished variability.

Variation and convergence in the morpho-functional properties of the mammalian neocortex

Man's natural inclination to classify and hierarchize the living world has prompted neurophysiologists to explore possible differences in brain organisation between mammals, with the aim of understanding the diversity of their behavioural repertoires. But what really distinguishes the human brain from that of a platypus, an opossum or a rodent? In this review, we compare the structural and electrical properties of neocortical neurons in the main mammalian radiations and examine their impact on the functioning of the networks they form. We discuss variations in overall brain size, number of neurons, length of their dendritic trees and density of spines, acknowledging their increase in humans as in most large-brained species. Our comparative analysis also highlights a remarkable consistency, particularly pronounced in marsupial and placental mammals, in the cell typology, intrinsic and synaptic electrical properties of pyramidal neuron subtypes, and in their organisation into functional circuits. These shared cellular and network characteristics contribute to the emergence of strikingly similar large-scale physiological and pathological brain dynamics across a wide range of species. These findings support the existence of a core set of neural principles and processes conserved throughout mammalian evolution, from which a number of species-specific adaptations appear, likely allowing distinct functional needs to be met in a variety of environmental contexts.

Directed and acyclic synaptic connectivity in the human layer 2-3 cortical microcircuit

The computational capabilities of neuronal networks are fundamentally constrained by their specific connectivity. Previous studies of cortical connectivity have mostly been carried out in rodents; whether the principles established therein also apply to the evolutionarily expanded human cortex is unclear. We studied network properties within the human temporal cortex using samples obtained from brain surgery. We analyzed multineuron patch-clamp recordings in layer 2-3 pyramidal neurons and identified substantial differences compared with rodents. Reciprocity showed random distribution, synaptic strength was independent from connection probability, and connectivity of the supragranular temporal cortex followed a directed and mostly acyclic graph topology. Application of these principles in neuronal models increased dimensionality of network dynamics, suggesting a critical role for cortical computation.

Multimodal Nature of the Single-cell Primate Brain Atlas: Morphology, Transcriptome, Electrophysiology, and Connectivity

Primates exhibit complex brain structures that augment cognitive function. The neocortex fulfills high-cognitive functions through billions of connected neurons. These neurons have distinct transcriptomic, morphological, and electrophysiological properties, and their connectivity principles vary. These features endow the primate brain atlas with a multimodal nature. The recent integration of next-generation sequencing with modified patch-clamp techniques is revolutionizing the way to census the primate neocortex, enabling a multimodal neuronal atlas to be established in great detail: (1) single-cell/single-nucleus RNA-seq technology establishes high-throughput transcriptomic references, covering all major transcriptomic cell types; (2) patch-seq links the morphological and electrophysiological features to the transcriptomic reference; (3) multicell patch-clamp delineates the principles of local connectivity. Here, we review the applications of these technologies in the primate neocortex and discuss the current advances and tentative gaps for a comprehensive understanding of the primate neocortex.

Emergence of co-tuning in inhibitory neurons as a network phenomenon mediated by randomness, correlations, and homeostatic plasticity

Cortical excitatory neurons show clear tuning to stimulus features, but the tuning properties of inhibitory interneurons are ambiguous. While inhibitory neurons have been considered to be largely untuned, some studies show that some parvalbumin-expressing (PV) neurons do show feature selectivity and participate in co-tuned subnetworks with pyramidal neurons. In this study, we first use mean-field theory to demonstrate that a combination of homeostatic plasticity governing the synaptic dynamics of the connections from PV to excitatory neurons, heterogeneity in the excitatory postsynaptic potentials that impinge on PV neurons, and shared correlated input from layer 4 results in the functional and structural self-organization of PV subnetworks. Second, we show that structural and functional feature tuning of PV neurons emerges more clearly at the network level, i.e., that population-level measures identify functional and structural co-tuning of PV neurons that are not evident in pairwise individual-level measures. Finally, we show that such co-tuning can enhance network stability at the cost of reduced feature selectivity.

Functional subtypes of synaptic dynamics in mouse and human

Synapses preferentially respond to particular temporal patterns of activity with a large degree of heterogeneity that is informally or tacitly separated into classes. Yet, the precise number and properties of such classes are unclear. Do they exist on a continuum and, if so, when is it appropriate to divide that continuum into functional regions? In a large dataset of glutamatergic cortical connections, we perform model-based characterization to infer the number and characteristics of functionally distinct subtypes of synaptic dynamics. In rodent data, we find five clusters that partially converge with transgenic-associated subtypes. Strikingly, the application of the same clustering method in human data infers a highly similar number of clusters, supportive of stable clustering. This nuanced dictionary of functional subtypes shapes the heterogeneity of cortical synaptic dynamics and provides a lens into the basic motifs of information transmission in the brain.

In-silico testing of new pharmacology for restoring inhibition and human cortical function in depression

Reduced inhibition by somatostatin-expressing interneurons is associated with depression. Administration of positive allosteric modulators of α5 subunit-containing GABAA receptor (α5-PAM) that selectively target this lost inhibition exhibit antidepressant and pro-cognitive effects in rodent models of chronic stress. However, the functional effects of α5-PAM on the human brain in vivo are unknown, and currently cannot be assessed experimentally. We modeled the effects of α5-PAM on tonic inhibition as measured in human neurons, and tested in silico α5-PAM effects on detailed models of human cortical microcircuits in health and depression. We found that α5-PAM effectively recovered impaired cortical processing as quantified by stimulus detection metrics, and also recovered the power spectral density profile of the microcircuit EEG signals. We performed an α5-PAM dose-response and identified simulated EEG biomarker candidates. Our results serve to de-risk and facilitate α5-PAM translation and provide biomarkers in non-invasive brain signals for monitoring target engagement and drug efficacy.

Exact Analysis of the Subthreshold Variability for Conductance-Based Neuronal Models with Synchronous Synaptic Inputs

The spiking activity of neocortical neurons exhibits a striking level of variability, even when these networks are driven by identical stimuli. The approximately Poisson firing of neurons has led to the hypothesis that these neural networks operate in the asynchronous state. In the asynchronous state, neurons fire independently from one another, so that the probability that a neuron experience synchronous synaptic inputs is exceedingly low. While the models of asynchronous neurons lead to observed spiking variability, it is not clear whether the asynchronous state can also account for the level of subthreshold membrane potential variability. We propose a new analytical framework to rigorously quantify the subthreshold variability of a single conductance-based neuron in response to synaptic inputs with prescribed degrees of synchrony. Technically, we leverage the theory of exchangeability to model input synchrony via jump-process-based synaptic drives; we then perform a moment analysis of the stationary response of a neuronal model with all-or-none conductances that neglects postspiking reset. As a result, we produce exact, interpretable closed forms for the first two stationary moments of the membrane voltage, with explicit dependence on the input synaptic numbers, strengths, and synchrony. For biophysically relevant parameters, we find that the asynchronous regime yields realistic subthreshold variability (voltage variance ≃4-9 mV2) only when driven by a restricted number of large synapses, compatible with strong thalamic drive. By contrast, we find that achieving realistic subthreshold variability with dense cortico-cortical inputs requires including weak but nonzero input synchrony, consistent with measured pairwise spiking correlations. We also show that, without synchrony, the neural variability averages out to zero for all scaling limits with vanishing synaptic weights, independent of any balanced state hypothesis. This result challenges the theoretical basis for mean-field theories of the asynchronous state.

Exact analysis of the subthreshold variability for conductance-based neuronal models with synchronous synaptic inputs

The spiking activity of neocortical neurons exhibits a striking level of variability, even when these networks are driven by identical stimuli. The approximately Poisson firing of neurons has led to the hypothesis that these neural networks operate in the asynchronous state. In the asynchronous state neurons fire independently from one another, so that the probability that a neuron experience synchronous synaptic inputs is exceedingly low. While the models of asynchronous neurons lead to observed spiking variability, it is not clear whether the asynchronous state can also account for the level of subthreshold membrane potential variability. We propose a new analytical framework to rigorously quantify the subthreshold variability of a single conductance-based neuron in response to synaptic inputs with prescribed degrees of synchrony. Technically we leverage the theory of exchangeability to model input synchrony via jump-process-based synaptic drives; we then perform a moment analysis of the stationary response of a neuronal model with all-or-none conductances that neglects post-spiking reset. As a result, we produce exact, interpretable closed forms for the first two stationary moments of the membrane voltage, with explicit dependence on the input synaptic numbers, strengths, and synchrony. For biophysically relevant parameters, we find that the asynchronous regime only yields realistic subthreshold variability (voltage variance ≃ 4 - 9 m V 2 ) when driven by a restricted number of large synapses, compatible with strong thalamic drive. By contrast, we find that achieving realistic subthreshold variability with dense cortico-cortical inputs requires including weak but nonzero input synchrony, consistent with measured pairwise spiking correlations. We also show that without synchrony, the neural variability averages out to zero for all scaling limits with vanishing synaptic weights, independent of any balanced state hypothesis. This result challenges the theoretical basis for mean-field theories of the asynchronous state.

How deep is the brain? The shallow brain hypothesis

Deep learning and predictive coding architectures commonly assume that inference in neural networks is hierarchical. However, largely neglected in deep learning and predictive coding architectures is the neurobiological evidence that all hierarchical cortical areas, higher or lower, project to and receive signals directly from subcortical areas. Given these neuroanatomical facts, today's dominance of cortico-centric, hierarchical architectures in deep learning and predictive coding networks is highly questionable; such architectures are likely to be missing essential computational principles the brain uses. In this Perspective, we present the shallow brain hypothesis: hierarchical cortical processing is integrated with a massively parallel process to which subcortical areas substantially contribute. This shallow architecture exploits the computational capacity of cortical microcircuits and thalamo-cortical loops that are not included in typical hierarchical deep learning and predictive coding networks. We argue that the shallow brain architecture provides several critical benefits over deep hierarchical structures and a more complete depiction of how mammalian brains achieve fast and flexible computational capabilities.

Coordinated changes in a cortical circuit sculpt effects of novelty on neural dynamics

Recent studies have found dramatic cell-type specific responses to stimulus novelty, highlighting the importance of analyzing the cortical circuitry at the cell-type specific level of granularity to understand brain function. Although initial work classified and characterized activity for each cell type, the specific alterations in cortical circuitry-particularly when multiple novelty effects interact-remain unclear. To address this gap, we employed a large-scale public dataset of electrophysiological recordings in the visual cortex of awake, behaving mice using Neuropixels probes and designed population network models to investigate the observed changes in neural dynamics in response to a combination of distinct forms of novelty. The model parameters were rigorously constrained by publicly available structural datasets, including multi-patch synaptic physiology and electron microscopy data. Our systematic optimization approach identified tens of thousands of model parameter sets that replicate the observed neural activity. Analysis of these solutions revealed generally weaker connections under novel stimuli, as well as a shift in the balance e between SST and VIP populations. Along with this, PV and SST populations experienced overall more excitatory influences compared to excitatory and VIP populations. Our results also highlight the role of VIP neurons in multiple aspects of visual stimulus processing and altering gain and saturation dynamics under novel conditions. In sum, our findings provide a systematic characterization of how the cortical circuit adapts to stimulus novelty by combining multiple rich public datasets.

Structural and functional specializations of human fast-spiking neurons support fast cortical signaling

Fast-spiking interneurons (FSINs) provide fast inhibition that synchronizes neuronal activity and is critical for cognitive function. Fast synchronization frequencies are evolutionary conserved in the expanded human neocortex despite larger neuron-to-neuron distances that challenge fast input-output transfer functions of FSINs. Here, we test in human neurons from neurosurgery tissue, which mechanistic specializations of human FSINs explain their fast-signaling properties in human cortex. With morphological reconstructions, multipatch recordings, and biophysical modeling, we find that despite threefold longer dendritic path, human FSINs maintain fast inhibition between connected pyramidal neurons through several mechanisms: stronger synapse strength of excitatory inputs, larger dendrite diameter with reduced complexity, faster AP initiation, and faster and larger inhibitory output, while Na+ current activation/inactivation properties are similar. These adaptations underlie short input-output delays in fast inhibition of human pyramidal neurons through FSINs, explaining how cortical synchronization frequencies are conserved despite expanded and sparse network topology of human cortex.

Compressive sensing of functional connectivity maps from patterned optogenetic stimulation of neuronal ensembles

Mapping functional connectivity between neurons is an essential step toward probing the neural computations mediating behavior. Accurately determining synaptic connectivity maps in populations of neurons is challenging in terms of yield, accuracy, and experimental time. Here, we developed a compressive sensing approach to reconstruct synaptic connectivity maps based on random two-photon cell-targeted optogenetic stimulation and membrane voltage readout of many putative postsynaptic neurons. Using a biophysical network model of interconnected populations of excitatory and inhibitory neurons, we characterized mapping recall and precision as a function of network observability, sparsity, number of neurons stimulated, off-target stimulation, synaptic reliability, propagation latency, and network topology. We found that mapping can be achieved with far fewer measurements than the standard pairwise sequential approach, with network sparsity and synaptic reliability serving as primary determinants of the performance. Our results suggest a rapid and efficient method to reconstruct functional connectivity of sparsely connected neuronal networks.

The influence of inter-regional delays in generating large-scale brain networks of phase synchronization

Large-scale networks of phase synchronization are considered to regulate the communication between brain regions fundamental to cognitive function, but the mapping to their structural substrates, i.e., the structure-function relationship, remains poorly understood. Biophysical Network Models (BNMs) have demonstrated the influences of local oscillatory activity and inter-regional anatomical connections in generating alpha-band (8-12 Hz) networks of phase synchronization observed with Electroencephalography (EEG) and Magnetoencephalography (MEG). Yet, the influence of inter-regional conduction delays remains unknown. In this study, we compared a BNM with standard "distance-dependent delays", which assumes constant conduction velocity, to BNMs with delays specified by two alternative methods accounting for spatially varying conduction velocities, "isochronous delays" and "mixed delays". We followed the Approximate Bayesian Computation (ABC) workflow, i) specifying neurophysiologically informed prior distributions of BNM parameters, ii) verifying the suitability of the prior distributions with Prior Predictive Checks, iii) fitting each of the three BNMs to alpha-band MEG resting-state data (N = 75) with Bayesian optimization for Likelihood-Free Inference (BOLFI), and iv) choosing between the fitted BNMs with ABC model comparison on a separate MEG dataset (N = 30). Prior Predictive Checks revealed the range of dynamics generated by each of the BNMs to encompass those seen in the MEG data, suggesting the suitability of the prior distributions. Fitting the models to MEG data yielded reliable posterior distributions of the parameters of each of the BNMs. Finally, model comparison revealed the BNM with "distance-dependent delays", as the most probable to describe the generation of alpha-band networks of phase synchronization seen in MEG. These findings suggest that distance-dependent delays might contribute to the neocortical architecture of human alpha-band networks of phase synchronization. Hence, our study illuminates the role of inter-regional delays in generating the large-scale networks of phase synchronization that might subserve the communication between regions vital to cognition.

Modeling functional cell types in spike train data

A major goal of computational neuroscience is to build accurate models of the activity of neurons that can be used to interpret their function in circuits. Here, we explore using functional cell types to refine single-cell models by grouping them into functionally relevant classes. Formally, we define a hierarchical generative model for cell types, single-cell parameters, and neural responses, and then derive an expectation-maximization algorithm with variational inference that maximizes the likelihood of the neural recordings. We apply this "simultaneous" method to estimate cell types and fit single-cell models from simulated data, and find that it accurately recovers the ground truth parameters. We then apply our approach to in vitro neural recordings from neurons in mouse primary visual cortex, and find that it yields improved prediction of single-cell activity. We demonstrate that the discovered cell-type clusters are well separated and generalizable, and thus amenable to interpretation. We then compare discovered cluster memberships with locational, morphological, and transcriptomic data. Our findings reveal the potential to improve models of neural responses by explicitly allowing for shared functional properties across neurons.

Artificial Neural Network Model with Astrocyte-Driven Short-Term Memory

In this study, we introduce an innovative hybrid artificial neural network model incorporating astrocyte-driven short-term memory. The model combines a convolutional neural network with dynamic models of short-term synaptic plasticity and astrocytic modulation of synaptic transmission. The model's performance was evaluated using simulated data from visual change detection experiments conducted on mice. Comparisons were made between the proposed model, a recurrent neural network simulating short-term memory based on sustained neural activity, and a feedforward neural network with short-term synaptic depression (STPNet) trained to achieve the same performance level as the mice. The results revealed that incorporating astrocytic modulation of synaptic transmission enhanced the model's performance.

Exact analysis of the subthreshold variability for conductance-based neuronal models with synchronous synaptic inputs

The spiking activity of neocortical neurons exhibits a striking level of variability, even when these networks are driven by identical stimuli. The approximately Poisson firing of neurons has led to the hypothesis that these neural networks operate in the asynchronous state. In the asynchronous state neurons fire independently from one another, so that the probability that a neuron experience synchronous synaptic inputs is exceedingly low. While the models of asynchronous neurons lead to observed spiking variability, it is not clear whether the asynchronous state can also account for the level of subthreshold membrane potential variability. We propose a new analytical framework to rigorously quantify the subthreshold variability of a single conductance-based neuron in response to synaptic inputs with prescribed degrees of synchrony. Technically we leverage the theory of exchangeability to model input synchrony via jump-process-based synaptic drives; we then perform a moment analysis of the stationary response of a neuronal model with all-or-none conductances that neglects post-spiking reset. As a result, we produce exact, interpretable closed forms for the first two stationary moments of the membrane voltage, with explicit dependence on the input synaptic numbers, strengths, and synchrony. For biophysically relevant parameters, we find that the asynchronous regime only yields realistic subthreshold variability (voltage variance ≅ 4-9mV 2 ) when driven by a restricted number of large synapses, compatible with strong thalamic drive. By contrast, we find that achieving realistic subthreshold variability with dense cortico-cortical inputs requires including weak but nonzero input synchrony, consistent with measured pairwise spiking correlations. We also show that without synchrony, the neural variability averages out to zero for all scaling limits with vanishing synaptic weights, independent of any balanced state hypothesis. This result challenges the theoretical basis for mean-field theories of the asynchronous state.

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.

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.

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.

Species-specific wiring of cortical circuits for small-world networks in the primary visual cortex

Long-range horizontal connections (LRCs) are conspicuous anatomical structures in the primary visual cortex (V1) of mammals, yet their detailed functions in relation to visual processing are not fully understood. Here, we show that LRCs are key components to organize a "small-world network" optimized for each size of the visual cortex, enabling the cost-efficient integration of visual information. Using computational simulations of a biologically inspired model neural network, we found that sparse LRCs added to networks, combined with dense local connections, compose a small-world network and significantly enhance image classification performance. We confirmed that the performance of the network appeared to be strongly correlated with the small-world coefficient of the model network under various conditions. Our theoretical model demonstrates that the amount of LRCs to build a small-world network depends on each size of cortex and that LRCs are beneficial only when the size of the network exceeds a certain threshold. Our model simulation of various sizes of cortices validates this prediction and provides an explanation of the species-specific existence of LRCs in animal data. Our results provide insight into a biological strategy of the brain to balance functional performance and resource cost.

Target cell-specific synaptic dynamics of excitatory to inhibitory neuron connections in supragranular layers of human neocortex

Rodent studies have demonstrated that synaptic dynamics from excitatory to inhibitory neuron types are often dependent on the target cell type. However, these target cell-specific properties have not been well investigated in human cortex, where there are major technical challenges in reliably obtaining healthy tissue, conducting multiple patch-clamp recordings on inhibitory cell types, and identifying those cell types. Here, we take advantage of newly developed methods for human neurosurgical tissue analysis with multiple patch-clamp recordings, post-hoc fluorescent in situ hybridization (FISH), machine learning-based cell type classification and prospective GABAergic AAV-based labeling to investigate synaptic properties between pyramidal neurons and PVALB- vs. SST-positive interneurons. We find that there are robust molecular differences in synapse-associated genes between these neuron types, and that individual presynaptic pyramidal neurons evoke postsynaptic responses with heterogeneous synaptic dynamics in different postsynaptic cell types. Using molecular identification with FISH and classifiers based on transcriptomically identified PVALB neurons analyzed by Patch-seq, we find that PVALB neurons typically show depressing synaptic characteristics, whereas other interneuron types including SST-positive neurons show facilitating characteristics. Together, these data support the existence of target cell-specific synaptic properties in human cortex that are similar to rodent, thereby indicating evolutionary conservation of local circuit connectivity motifs from excitatory to inhibitory neurons and their synaptic dynamics.

Targeting operational regimes of interest in recurrent neural networks

Neural computations emerge from local recurrent neural circuits or computational units such as cortical columns that comprise hundreds to a few thousand neurons. Continuous progress in connectomics, electrophysiology, and calcium imaging require tractable spiking network models that can consistently incorporate new information about the network structure and reproduce the recorded neural activity features. However, for spiking networks, it is challenging to predict which connectivity configurations and neural properties can generate fundamental operational states and specific experimentally reported nonlinear cortical computations. Theoretical descriptions for the computational state of cortical spiking circuits are diverse, including the balanced state where excitatory and inhibitory inputs balance almost perfectly or the inhibition stabilized state (ISN) where the excitatory part of the circuit is unstable. It remains an open question whether these states can co-exist with experimentally reported nonlinear computations and whether they can be recovered in biologically realistic implementations of spiking networks. Here, we show how to identify spiking network connectivity patterns underlying diverse nonlinear computations such as XOR, bistability, inhibitory stabilization, supersaturation, and persistent activity. We establish a mapping between the stabilized supralinear network (SSN) and spiking activity which allows us to pinpoint the location in parameter space where these activity regimes occur. Notably, we find that biologically-sized spiking networks can have irregular asynchronous activity that does not require strong excitation-inhibition balance or large feedforward input and we show that the dynamic firing rate trajectories in spiking networks can be precisely targeted without error-driven training algorithms.

Transcriptomic cell type structures in vivo neuronal activity across multiple timescales

Cell type is hypothesized to be a key determinant of a neuron's role within a circuit. Here, we examine whether a neuron's transcriptomic type influences the timing of its activity. We develop a deep-learning architecture that learns features of interevent intervals across timescales (ms to >30 min). We show that transcriptomic cell-class information is embedded in the timing of single neuron activity in the intact brain of behaving animals (calcium imaging and extracellular electrophysiology) as well as in a bio-realistic model of the visual cortex. Further, a subset of excitatory cell types are distinguishable but can be classified with higher accuracy when considering cortical layer and projection class. Finally, we show that computational fingerprints of cell types may be universalizable across structured stimuli and naturalistic movies. Our results indicate that transcriptomic class and type may be imprinted in the timing of single neuron activity across diverse stimuli.

How synaptic strength, short-term plasticity, and input synchrony contribute to neuronal spike output

Neurons integrate from thousands of synapses whose strengths span an order of magnitude. Intriguingly, in mouse neocortex, the few 'strong' synapses are formed between similarly tuned cells, suggesting they determine spiking output. This raises the question of how other computational primitives, including 'background' activity from the many 'weak' synapses, short-term plasticity, and temporal factors contribute to spiking. We used paired recordings and extracellular stimulation experiments to map excitatory postsynaptic potential (EPSP) amplitudes and paired-pulse ratios of synaptic connections formed between pyramidal neurons in layer 2/3 (L2/3) of barrel cortex. While net short-term plasticity was weak, strong synaptic connections were exclusively depressing. Importantly, we found no evidence for clustering of synaptic properties on individual neurons. Instead, EPSPs and paired-pulse ratios of connections converging onto the same cells spanned the full range observed across L2/3, which critically constrains theoretical models of cortical filtering. To investigate how different computational primitives of synaptic information processing interact to shape spiking, we developed a computational model of a pyramidal neuron in the excitatory L2/3 circuitry, which was constrained by our experiments and published in vivo data. We found that strong synapses were substantially depressed during ongoing activation and their ability to evoke correlated spiking primarily depended on their high temporal synchrony and high firing rates observed in vivo. However, despite this depression, their larger EPSP amplitudes strongly amplified information transfer and responsiveness. Thus, our results contribute to a nuanced framework of how cortical neurons exploit synergies between temporal coding, synaptic properties, and noise to transform synaptic inputs into spikes.

Layer-dependent stability of intracortical recordings and neuronal cell loss

Intracortical recordings can be used to voluntarily control external devices via brain-machine interfaces (BMI). Multiple factors, including the foreign body response (FBR), limit the stability of these neural signals over time. Current clinically approved devices consist of multi-electrode arrays with a single electrode site at the tip of each shank, confining the recording interface to a single layer of the cortex. Advancements in manufacturing technology have led to the development of high-density electrodes that can record from multiple layers. However, the long-term stability of neural recordings and the extent of neuronal cell loss around the electrode across different cortical depths have yet to be explored. To answer these questions, we recorded neural signals from rats chronically implanted with a silicon-substrate microelectrode array spanning the layers of the cortex. Our results show the long-term stability of intracortical recordings varies across cortical depth, with electrode sites around L4-L5 having the highest stability. Using machine learning guided segmentation, our novel histological technique, DeepHisto, revealed that the extent of neuronal cell loss varies across cortical layers, with L2/3 and L4 electrodes having the largest area of neuronal cell loss. These findings suggest that interfacing depth plays a major role in the FBR and long-term performance of intracortical neuroprostheses.

Determinants of functional synaptic connectivity among amygdala-projecting prefrontal cortical neurons in male mice

The medial prefrontal cortex (mPFC) mediates a variety of complex cognitive functions via its vast and diverse connections with cortical and subcortical structures. Understanding the patterns of synaptic connectivity that comprise the mPFC local network is crucial for deciphering how this circuit processes information and relays it to downstream structures. To elucidate the synaptic organization of the mPFC, we developed a high-throughput optogenetic method for mapping large-scale functional synaptic connectivity in acute brain slices. We show that in male mice, mPFC neurons that project to the basolateral amygdala (BLA) display unique spatial patterns of local-circuit synaptic connectivity, which distinguish them from the general mPFC cell population. When considering synaptic connections between pairs of mPFC neurons, the intrinsic properties of the postsynaptic cell and the anatomical positions of both cells jointly account for ~7.5% of the variation in the probability of connection. Moreover, anatomical distance and laminar position explain most of this fraction in variation. Our findings reveal the factors determining connectivity in the mPFC and delineate the architecture of synaptic connections in the BLA-projecting subnetwork.

Strong and reliable synaptic communication between pyramidal neurons in adult human cerebral cortex

Synaptic transmission constitutes the primary mode of communication between neurons. It is extensively studied in rodent but not human neocortex. We characterized synaptic transmission between pyramidal neurons in layers 2 and 3 using neurosurgically resected human middle temporal gyrus (MTG, Brodmann area 21), which is part of the distributed language circuitry. We find that local connectivity is comparable with mouse layer 2/3 connections in the anatomical homologue (temporal association area), but synaptic connections in human are 3-fold stronger and more reliable (0% vs 25% failure rates, respectively). We developed a theoretical approach to quantify properties of spinous synapses showing that synaptic conductance and voltage change in human dendritic spines are 3-4-folds larger compared with mouse, leading to significant NMDA receptor activation in human unitary connections. This model prediction was validated experimentally by showing that NMDA receptor activation increases the amplitude and prolongs decay of unitary excitatory postsynaptic potentials in human but not in mouse connections. Since NMDA-dependent recurrent excitation facilitates persistent activity (supporting working memory), our data uncovers cortical microcircuit properties in human that may contribute to language processing in MTG.

Modeling functional cell types in spike train data

A major goal of computational neuroscience is to build accurate models of the activity of neurons that can be used to interpret their function in circuits. Here, we explore using functional cell types to refine single-cell models by grouping them into functionally relevant classes. Formally, we define a hierarchical generative model for cell types, single-cell parameters, and neural responses, and then derive an expectation-maximization algorithm with variational inference that maximizes the likelihood of the neural recordings. We apply this "simultaneous" method to estimate cell types and fit single-cell models from simulated data, and find that it accurately recovers the ground truth parameters. We then apply our approach to in vitro neural recordings from neurons in mouse primary visual cortex, and find that it yields improved prediction of single-cell activity. We demonstrate that the discovered cell-type clusters are well separated and generalizable, and thus amenable to interpretation. We then compare discovered cluster memberships with locational, morphological, and transcriptomic data. Our findings reveal the potential to improve models of neural responses by explicitly allowing for shared functional properties across neurons.

Developmental neuronal origin regulates neocortical map formation

Sensory neurons in the neocortex exhibit distinct functional selectivity to constitute the neural map. While neocortical map of the visual cortex in higher mammals is clustered, it displays a striking "salt-and-pepper" pattern in rodents. However, little is known about the origin and basis of the interspersed neocortical map. Here we report that the intricate excitatory neuronal kinship-dependent synaptic connectivity influences precise functional map organization in the mouse primary visual cortex. While sister neurons originating from the same neurogenic radial glial progenitors (RGPs) preferentially develop synapses, cousin neurons derived from amplifying RGPs selectively antagonize horizontal synapse formation. Accordantly, cousin neurons in similar layers exhibit clear functional selectivity differences, contributing to a salt-and-pepper architecture. Removal of clustered protocadherins (cPCDHs), the largest subgroup of the diverse cadherin superfamily, eliminates functional selectivity differences between cousin neurons and alters neocortical map organization. These results suggest that developmental neuronal origin regulates neocortical map formation via cPCDHs.

Layer-specific distribution and expression pattern of AMPA- and NMDA-type glutamate receptors in the barrel field of the adult rat somatosensory cortex: a quantitative electron microscopic analysis

AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) and NMDA (N-methyl-d-aspartate) glutamate receptors are driving forces for synaptic transmission and plasticity at neocortical synapses. However, their distribution pattern in the adult rat neocortex is largely unknown and was quantified using freeze fracture replication combined with postimmunogold-labeling. Both receptors were co-localized at layer (L)4 and L5 postsynaptic densities (PSDs). At L4 dendritic shaft and spine PSDs, the number of gold grains detecting AMPA was similar, whereas at L5 shaft PSDs AMPA-receptors outnumbered those on spine PSDs. Their number was significantly higher at L5 vs. L4 PSDs. At L4 and L5 dendritic shaft PSDs, the number of gold grains detecting GluN1 was ~2-fold higher than at spine PSDs. The number of gold grains detecting the GluN1-subunit was higher for both shaft and spine PSDs in L5 vs. L4. Both receptors showed a large variability in L4 and L5. A high correlation between the number of gold grains and PSD size for both receptors and targets was observed. Both receptors were distributed over the entire PSD but showed a layer- and target-specific distribution pattern. The layer- and target-specific distribution of AMPA and GluN1 glutamate receptors partially contribute to the observed functional differences in synaptic transmission and plasticity in the neocortex.