Physiology and anatomy of synaptic connections between thick tufted pyramidal neurones in the developing rat neocortex

Directional intermodular coupling enriches functional complexity in biological neuronal networks

Hierarchically modular organization is a canonical network topology that is evolutionarily conserved in the nervous systems of animals. Within the network, neurons form directional connections defined by the growth of their axonal terminals. However, this topology is dissimilar to the network formed by dissociated neurons in culture because they form randomly connected networks on homogeneous substrates. In this study, we fabricated microfluidic devices to reconstitute hierarchically modular neuronal networks in culture (in vitro) and investigated how non-random structures, such as directional connectivity between modules, affect global network dynamics. Embedding directional connections in a pseudo-feedforward manner suppressed excessive synchrony in cultured neuronal networks and enhanced the integration-segregation balance. Modeling the behavior of biological neuronal networks using spiking neural networks (SNNs) further revealed that modularity and directionality cooperate to shape such network dynamics. Finally, we demonstrate that for a given network topology, the statistics of network dynamics, such as global network activation, correlation coefficient, and functional complexity, can be analytically predicted based on eigendecomposition of the transition matrix in the state-transition model. Hence, the integration of bioengineering and cell culture technologies enables us not only to reconstitute complex network circuitry in the nervous system but also to understand the structure-function relationships in biological neuronal networks by bridging theoretical modeling with in vitro experiments.

A dynamic bifurcation mechanism explains cortex-wide neural correlates of conscious access

Conscious access is suggested to involve "ignition," an all-or-none activation across cortical areas. To elucidate this phenomenon, we carry out computer simulations of a detection task using a mesoscale connectome-based model for the multiregional macaque cortex. The model uncovers a dynamic bifurcation mechanism that gives rise to ignition in a network of associative regions. A hierarchical N-methyl-D-aspartate (NMDA)/α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor gradient plays a critical role: fast AMPA receptors drive feedforward signal propagation, while slow NMDA receptors in feedback pathways shape and sustain the ignited network. Intriguingly, the model suggests higher NMDA-to-AMPA receptor ratios in sensory areas compared to association areas, a prediction supported by in vitro autoradiography data. Furthermore, the model accounts for diverse behavioral and physiological phenomena linked to consciousness. This work sheds light on how receptor gradients along the cortical hierarchy enable distributed cognitive functions and provides a biologically constrained computational framework for investigating the neurophysiological basis of conscious access.

Optogenetic estimation of synaptic connections in brain slices

BACKGROUND

Detection of synaptic connections is essential for understanding neural circuits. By using optogenetics, current injection, and glutamate uncaging to activate presynaptic cells and simultaneously recording the subsequent response of postsynaptic cells, the presence of synaptic connections can be confirmed. However, these methods present throughput challenges, such as the need for simultaneous multicellular patch-clamp recording and two-photon microscopy. These challenges lead to a trade-off between sacrificing resolution and experimental throughput.

NEW METHOD

We adopted the laser, typically used for local field ablation, and combined this with post hoc analysis. We successfully approximated the synaptic connection probabilities using only an epi-fluorescence microscope and single-cell recordings.

RESULTS

We sequentially stimulated the channelrhodopsin 2-expressing cells surrounding the recorded cell and approximated the synaptic connection probabilities. This probability value was comparable to that obtained from simultaneous multi-cell patch-clamp recordings, which included more than 600 pairs.

COMPARISON WITH EXISTING METHODS

Our setup allows us to estimate connection probabilities within 100 s, outperforming existing methods. We successfully estimated synaptic connection probabilities using only the optical path typically used by an epi-fluorescence microscope and single-cell recordings. It may also be suitable for dendritic ablation experiments.

CONCLUSIONS

The proposed method simplifies the estimation of connection probabilities, which is expected to advance the study of neural circuits in conditions such as autism and schizophrenia where connection probabilities vary. Furthermore, this approach is applicable not only to local circuits but also to long-range connections, thus increasing experimental throughput.

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.

Synaptic interactions between stellate cells and parvalbumin interneurons in layer 2 of the medial entorhinal cortex are organized at the scale of grid cell clusters

Interactions between excitatory and inhibitory neurons are critical to computations in cortical circuits but their organization is difficult to assess with standard electrophysiological approaches. Within the medial entorhinal cortex, representation of location by grid and other spatial cells involves circuits in layer 2 in which excitatory stellate cells interact with each other via inhibitory parvalbumin expressing interneurons. Whether this connectivity is structured to support local circuit computations is unclear. Here, we introduce strategies to address the functional organization of excitatory-inhibitory interactions using crossed Cre- and Flp-driver mouse lines to direct targeted presynaptic optogenetic activation and postsynaptic cell identification. We then use simultaneous patch-clamp recordings from postsynaptic neurons to assess their shared input from optically activated presynaptic populations. We find that extensive axonal projections support spatially organized connectivity between stellate cells and parvalbumin interneurons, such that direct connections are often, but not always, shared by nearby neurons, whereas multisynaptic interactions coordinate inputs to neurons with greater spatial separation. We suggest that direct excitatory-inhibitory synaptic interactions may operate at the scale of grid cell clusters, with local modules defined by excitatory-inhibitory connectivity, while indirect interactions may coordinate activity at the scale of grid cell modules.

A dendritic mechanism for balancing synaptic flexibility and stability

Biological and artificial neural networks learn by modifying synaptic weights, but it is unclear how these systems retain previous knowledge and also acquire new information. Here, we show that cortical pyramidal neurons can solve this plasticity-versus-stability dilemma by differentially regulating synaptic plasticity at distinct dendritic compartments. Oblique dendrites of adult mouse layer 5 cortical pyramidal neurons selectively receive monosynaptic thalamic input, integrate linearly, and lack burst-timing synaptic potentiation. In contrast, basal dendrites, which do not receive thalamic input, exhibit conventional NMDA receptor (NMDAR)-mediated supralinear integration and synaptic potentiation. Congruently, spiny synapses on oblique branches show decreased structural plasticity in vivo. The selective decline in NMDAR activity and expression at synapses on oblique dendrites is controlled by a critical period of visual experience. Our results demonstrate a biological mechanism for how single neurons can safeguard a set of inputs from ongoing plasticity by altering synaptic properties at distinct dendritic domains.

Clustered synapses develop in distinct dendritic domains in visual cortex before eye opening

Synaptic inputs to cortical neurons are highly structured in adult sensory systems, such that neighboring synapses along dendrites are activated by similar stimuli. This organization of synaptic inputs, called synaptic clustering, is required for high-fidelity signal processing, and clustered synapses can already be observed before eye opening. However, how clustered inputs emerge during development is unknown. Here, we employed concurrent in vivo whole-cell patch-clamp and dendritic calcium imaging to map spontaneous synaptic inputs to dendrites of layer 2/3 neurons in the mouse primary visual cortex during the second postnatal week until eye opening. We found that the number of functional synapses and the frequency of transmission events increase several fold during this developmental period. At the beginning of the second postnatal week, synapses assemble specifically in confined dendritic segments, whereas other segments are devoid of synapses. By the end of the second postnatal week, just before eye opening, dendrites are almost entirely covered by domains of co-active synapses. Finally, co-activity with their neighbor synapses correlates with synaptic stabilization and potentiation. Thus, clustered synapses form in distinct functional domains presumably to equip dendrites with computational modules for high-capacity sensory processing when the eyes open.

Differences in the consolidation by spontaneous and evoked ripples in the presence of active dendrites

Ripples are a typical form of neural activity in hippocampal neural networks associated with the replay of episodic memories during sleep as well as sleep-related plasticity and memory consolidation. The emergence of ripples has been observed both dependent as well as independent of input from other brain areas and often coincides with dendritic spikes. Yet, it is unclear how input-evoked and spontaneous ripples as well as dendritic excitability affect plasticity and consolidation. Here, we use mathematical modeling to compare these cases. We find that consolidation as well as the emergence of spontaneous ripples depends on a reliable propagation of activity in feed-forward structures which constitute memory representations. This propagation is facilitated by excitable dendrites, which entail that a few strong synapses are sufficient to trigger neuronal firing. In this situation, stimulation-evoked ripples lead to the potentiation of weak synapses within the feed-forward structure and, thus, to a consolidation of a more general sequence memory. However, spontaneous ripples that occur without stimulation, only consolidate a sparse backbone of the existing strong feed-forward structure. Based on this, we test a recently hypothesized scenario in which the excitability of dendrites is transiently enhanced after learning, and show that such a transient increase can strengthen, restructure and consolidate even weak hippocampal memories, which would be forgotten otherwise. Hence, a transient increase in dendritic excitability would indeed provide a mechanism for stabilizing memories.

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.

Synchronization of delayed coupled neurons with multiple synaptic connections

Synchronization is a key feature of the brain dynamics and is necessary for information transmission across brain regions and in higher brain functions like cognition, learning and memory. Experimental findings demonstrated that in cortical microcircuits there are multiple synapses between pairs of connected neurons. Synchronization of neurons in the presence of multiple synaptic connections may be relevant for optimal learning and memory, however, its effect on the dynamics of the neurons is not adequately studied. Here, we address the question that how changes in the strength of the synaptic connections and transmission delays between neurons impact synchronization in a two-neuron system with multiple synapses. To this end, we analytically and computationally investigated synchronization dynamics by considering both phase oscillator model and conductance-based Hodgkin-Huxley (HH) model. Our results show that symmetry/asymmetry of feedforward and feedback connections crucially determines stability of the phase locking of the system based on the strength of connections and delays. In both models, the two-neuron system with multiple synapses achieves in-phase synchrony in the presence of small and large delays, whereas an anti-phase synchronization state is favored for median delays. Our findings can expand the understanding of the functional role of multisynaptic contacts in neuronal synchronization and may shed light on the dynamical consequences of pathological multisynaptic connectivity in a number of brain disorders.

Data-driven multiscale computational models of cortical and subcortical regions

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

Enhanced Release Probability without Changes in Synaptic Delay during Analogue-Digital Facilitation

Neuronal timing with millisecond precision is critical for many brain functions such as sensory perception, learning and memory formation. At the level of the chemical synapse, the synaptic delay is determined by the presynaptic release probability (Pr) and the waveform of the presynaptic action potential (AP). For instance, paired-pulse facilitation or presynaptic long-term potentiation are associated with reductions in the synaptic delay, whereas paired-pulse depression or presynaptic long-term depression are associated with an increased synaptic delay. Parallelly, the AP broadening that results from the inactivation of voltage gated potassium (Kv) channels responsible for the repolarization phase of the AP delays the synaptic response, and the inactivation of sodium (Nav) channels by voltage reduces the synaptic latency. However, whether synaptic delay is modulated during depolarization-induced analogue-digital facilitation (d-ADF), a form of context-dependent synaptic facilitation induced by prolonged depolarization of the presynaptic neuron and mediated by the voltage-inactivation of presynaptic Kv1 channels, remains unclear. We show here that despite Pr being elevated during d-ADF at pyramidal L5-L5 cell synapses, the synaptic delay is surprisingly unchanged. This finding suggests that both Pr- and AP-dependent changes in synaptic delay compensate for each other during d-ADF. We conclude that, in contrast to other short- or long-term modulations of presynaptic release, synaptic timing is not affected during d-ADF because of the opposite interaction of Pr- and AP-dependent modulations of synaptic delay.

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

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

Dorsal hippocampus to nucleus accumbens projections drive reinforcement via activation of accumbal dynorphin neurons

The hippocampus is pivotal in integrating emotional processing, learning, memory, and reward-related behaviors. The dorsal hippocampus (dHPC) is particularly crucial for episodic, spatial, and associative memory, and has been shown to be necessary for context- and cue-associated reward behaviors. The nucleus accumbens (NAc), a central structure in the mesolimbic reward pathway, integrates the salience of aversive and rewarding stimuli. Despite extensive research on dHPC→NAc direct projections, their sufficiency in driving reinforcement and reward-related behavior remains to be determined. Our study establishes that activating excitatory neurons in the dHPC is sufficient to induce reinforcing behaviors through its direct projections to the dorso-medial subregion of the NAc shell (dmNAcSh). Notably, dynorphin-containing neurons specifically contribute to dHPC-driven reinforcing behavior, even though both dmNAcSh dynorphin- and enkephalin-containing neurons are activated with dHPC stimulation. Our findings unveil a pathway governing reinforcement, advancing our understanding of the hippocampal circuity's role in reward-seeking behaviors.

Understanding the genetic mechanisms and cognitive impairments in Down syndrome: towards a holistic approach

The most common genetic cause of intellectual disability is Down syndrome (DS), trisomy 21. It commonly results from three copies of human chromosome 21 (HC21). There are no mutations or deletions involved in DS. Instead, the phenotype is caused by altered transcription of the genes on HC21. These transcriptional variations are responsible for a myriad of symptoms affecting every organ system. A very debilitating aspect of DS is intellectual disability (ID). Although tremendous advances have been made to try and understand the underlying mechanisms of ID, there is a lack of a unified, holistic view to defining the cause and managing the cognitive impairments. In this literature review, we discuss the mechanisms of neuronal over-inhibition, abnormal morphology, and other genetic factors in contributing to the development of ID in DS patients and to gain a holistic understanding of ID in DS patients. We also highlight potential therapeutic approaches to improve the quality of life of DS patients.

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.

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.

Synaptic weights that correlate with presynaptic selectivity increase decoding performance

The activity of neurons in the visual cortex is often characterized by tuning curves, which are thought to be shaped by Hebbian plasticity during development and sensory experience. This leads to the prediction that neural circuits should be organized such that neurons with similar functional preference are connected with stronger weights. In support of this idea, previous experimental and theoretical work have provided evidence for a model of the visual cortex characterized by such functional subnetworks. A recent experimental study, however, have found that the postsynaptic preferred stimulus was defined by the total number of spines activated by a given stimulus and independent of their individual strength. While this result might seem to contradict previous literature, there are many factors that define how a given synaptic input influences postsynaptic selectivity. Here, we designed a computational model in which postsynaptic functional preference is defined by the number of inputs activated by a given stimulus. Using a plasticity rule where synaptic weights tend to correlate with presynaptic selectivity, and is independent of functional-similarity between pre- and postsynaptic activity, we find that this model can be used to decode presented stimuli in a manner that is comparable to maximum likelihood inference.

Spike timing-dependent plasticity alters electrosensory neuron synaptic strength in vitro but does not consistently predict changes in sensory tuning in vivo

How do sensory systems optimize detection of behaviorally relevant stimuli when the sensory environment is constantly changing? We addressed the role of spike timing-dependent plasticity (STDP) in driving changes in synaptic strength in a sensory pathway and whether those changes in synaptic strength could alter sensory tuning. It is challenging to precisely control temporal patterns of synaptic activity in vivo and replicate those patterns in vitro in behaviorally relevant ways. This makes it difficult to make connections between STDP-induced changes in synaptic physiology and plasticity in sensory systems. Using the mormyrid species Brevimyrus niger and Brienomyrus brachyistius, which produce electric organ discharges for electrolocation and communication, we can precisely control the timing of synaptic input in vivo and replicate these same temporal patterns of synaptic input in vitro. In central electrosensory neurons in the electric communication pathway, using whole cell intracellular recordings in vitro, we paired presynaptic input with postsynaptic spiking at different delays. Using whole cell intracellular recordings in awake, behaving fish, we paired sensory stimulation with postsynaptic spiking using the same delays. We found that Hebbian STDP predictably alters sensory tuning in vitro and is mediated by NMDA receptors. However, the change in synaptic responses induced by sensory stimulation in vivo did not adhere to the direction predicted by the STDP observed in vitro. Further analysis suggests that this difference is influenced by polysynaptic activity, including inhibitory interneurons. Our findings suggest that STDP rules operating at identified synapses may not drive predictable changes in sensory responses at the circuit level.NEW & NOTEWORTHY We replicated behaviorally relevant temporal patterns of synaptic activity in vitro and used the same patterns during sensory stimulation in vivo. There was a Hebbian spike timing-dependent plasticity (STDP) pattern in vitro, but sensory responses in vivo did not shift according to STDP predictions. Analysis suggests that this disparity is influenced by differences in polysynaptic activity, including inhibitory interneurons. These results suggest that STDP rules at synapses in vitro do not necessarily apply to circuits in vivo.

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.

GluN2B-NMDAR subunit contribution on synaptic plasticity: A phenomenological model for CA3-CA1 synapses

Synaptic plasticity is believed to be a key mechanism underlying learning and memory. We developed a phenomenological N-methyl-D-aspartate (NMDA) receptor-based voltage-dependent synaptic plasticity model for synaptic modifications at hippocampal CA3-CA1 synapses on a hippocampal CA1 pyramidal neuron. The model incorporates the GluN2A-NMDA and GluN2B-NMDA receptor subunit-based functions and accounts for the synaptic strength dependence on the postsynaptic NMDA receptor composition and functioning without explicitly modeling the NMDA receptor-mediated intracellular calcium, a local trigger of synaptic plasticity. We embedded the model into a two-compartmental model of a hippocampal CA1 pyramidal cell and validated it against experimental data of spike-timing-dependent synaptic plasticity (STDP), high and low-frequency stimulation. The developed model predicts altered learning rules in synapses formed on the apical dendrites of the detailed compartmental model of CA1 pyramidal neuron in the presence of the GluN2B-NMDA receptor hypofunction and can be used in hippocampal networks to model learning in health and disease.

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.

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.

Uncovering hidden network architecture from spiking activities using an exact statistical input-output relation of neurons

Identifying network architecture from observed neural activities is crucial in neuroscience studies. A key requirement is knowledge of the statistical input-output relation of single neurons in vivo. By utilizing an exact analytical solution of the spike-timing for leaky integrate-and-fire neurons under noisy inputs balanced near the threshold, we construct a framework that links synaptic type, strength, and spiking nonlinearity with the statistics of neuronal population activity. The framework explains structured pairwise and higher-order interactions of neurons receiving common inputs under different architectures. We compared the theoretical predictions with the activity of monkey and mouse V1 neurons and found that excitatory inputs given to pairs explained the observed sparse activity characterized by strong negative triple-wise interactions, thereby ruling out the alternative explanation by shared inhibition. Moreover, we showed that the strong interactions are a signature of excitatory rather than inhibitory inputs whenever the spontaneous rate is low. We present a guide map of neural interactions that help researchers to specify the hidden neuronal motifs underlying observed interactions found in empirical data.

Cellular composition and circuit organization of the locus coeruleus of adult mice

The locus coeruleus (LC) houses the vast majority of noradrenergic neurons in the brain and regulates many fundamental functions, including fight and flight response, attention control, and sleep/wake cycles. While efferent projections of the LC have been extensively investigated, little is known about its local circuit organization. Here, we performed large-scale multipatch recordings of noradrenergic neurons in adult mouse LC to profile their morpho-electric properties while simultaneously examining their interactions. LC noradrenergic neurons are diverse and could be classified into two major morpho-electric types. While fast excitatory synaptic transmission among LC noradrenergic neurons was not observed in our preparation, these mature LC neurons connected via gap junction at a rate similar to their early developmental stage and comparable to other brain regions. Most electrical connections form between dendrites and are restricted to narrowly spaced pairs or small clusters of neurons of the same type. In addition, more than two electrically coupled cell pairs were often identified across a cohort of neurons from individual multicell recording sets that followed a chain-like organizational pattern. The assembly of LC noradrenergic neurons thus follows a spatial and cell-type-specific wiring principle that may be imposed by a unique chain-like rule.

Reconstruction of sparse recurrent connectivity and inputs from the nonlinear dynamics of neuronal networks

Reconstructing the recurrent structural connectivity of neuronal networks is a challenge crucial to address in characterizing neuronal computations. While directly measuring the detailed connectivity structure is generally prohibitive for large networks, we develop a novel framework for reverse-engineering large-scale recurrent network connectivity matrices from neuronal dynamics by utilizing the widespread sparsity of neuronal connections. We derive a linear input-output mapping that underlies the irregular dynamics of a model network composed of both excitatory and inhibitory integrate-and-fire neurons with pulse coupling, thereby relating network inputs to evoked neuronal activity. Using this embedded mapping and experimentally feasible measurements of the firing rate as well as voltage dynamics in response to a relatively small ensemble of random input stimuli, we efficiently reconstruct the recurrent network connectivity via compressive sensing techniques. Through analogous analysis, we then recover high dimensional natural stimuli from evoked neuronal network dynamics over a short time horizon. This work provides a generalizable methodology for rapidly recovering sparse neuronal network data and underlines the natural role of sparsity in facilitating the efficient encoding of network data in neuronal dynamics.

Bayesian reconstruction of memories stored in neural networks from their connectivity

The advent of comprehensive synaptic wiring diagrams of large neural circuits has created the field of connectomics and given rise to a number of open research questions. One such question is whether it is possible to reconstruct the information stored in a recurrent network of neurons, given its synaptic connectivity matrix. Here, we address this question by determining when solving such an inference problem is theoretically possible in specific attractor network models and by providing a practical algorithm to do so. The algorithm builds on ideas from statistical physics to perform approximate Bayesian inference and is amenable to exact analysis. We study its performance on three different models, compare the algorithm to standard algorithms such as PCA, and explore the limitations of reconstructing stored patterns from synaptic connectivity.

New insights into binocular rivalry from the reconstruction of evolving percepts using model network dynamics

When the two eyes are presented with highly distinct stimuli, the resulting visual percept generally switches every few seconds between the two monocular images in an irregular fashion, giving rise to a phenomenon known as binocular rivalry. While a host of theoretical studies have explored potential mechanisms for binocular rivalry in the context of evoked model dynamics in response to simple stimuli, here we investigate binocular rivalry directly through complex stimulus reconstructions based on the activity of a two-layer neuronal network model with competing downstream pools driven by disparate monocular stimuli composed of image pixels. To estimate the dynamic percept, we derive a linear input-output mapping rooted in the non-linear network dynamics and iteratively apply compressive sensing techniques for signal recovery. Utilizing a dominance metric, we are able to identify when percept alternations occur and use data collected during each dominance period to generate a sequence of percept reconstructions. We show that despite the approximate nature of the input-output mapping and the significant reduction in neurons downstream relative to stimulus pixels, the dominant monocular image is well-encoded in the network dynamics and improvements are garnered when realistic spatial receptive field structure is incorporated into the feedforward connectivity. Our model demonstrates gamma-distributed dominance durations and well obeys Levelt's four laws for how dominance durations change with stimulus strength, agreeing with key recurring experimental observations often used to benchmark rivalry models. In light of evidence that individuals with autism exhibit relatively slow percept switching in binocular rivalry, we corroborate the ubiquitous hypothesis that autism manifests from reduced inhibition in the brain by systematically probing our model alternation rate across choices of inhibition strength. We exhibit sufficient conditions for producing binocular rivalry in the context of natural scene stimuli, opening a clearer window into the dynamic brain computations that vary with the generated percept and a potential path toward further understanding neurological disorders.

Functional Applications of Stable Tau Oligomers in Cell Biology and Electrophysiology Studies

Aggregated tau protein plays a key role in the pathogenesis of neurodegenerative tauopathies including Alzheimer's disease. Soluble, low-molecular-weight tau oligomers are formed early in disease processes and are thought to have toxic functions that disrupt neuronal function. The dynamic and transient nature of tau oligomers complicates in vitro functional studies to explore the mechanistic links between oligomer formation and neurodegeneration. We have previously described a method of producing stable and structurally characterized oligomers that maintain their oligomeric conformation and prevent further aggregation. This method allows for the flexibility of stabilizing tau oligomers by specifically labelling cysteine residues with fluorescent or colorless maleimide conjugates. Here, we describe the functional applications of these preformed stable tau oligomers in cell biology and electrophysiological studies. These investigations allow real-time insights into the cellular uptake of exogenous tau oligomers and their functional effects in the recipient cells.

Processing of cell assemblies in the lateral entorhinal cortex

There is evidence that olfactory cortex responds to its afferent input with the generation of cell assemblies: collections of principal neurons that fire together over a time scale of tens of ms. If such assemblies form an odor representation, then a fundamental question is how each assembly then induces neuronal activity in downstream structures. We have addressed this question in a detailed model of superficial layers of lateral entorhinal cortex, a recipient of input from olfactory cortex and olfactory bulb. Our results predict that the response of the fan cell subpopulation can be approximated by a relatively simple Boolean process, somewhat along the lines of the McCulloch/Pitts scheme; this is the case because of the sparsity of recurrent excitation amongst fan cells. However, because of recurrent excitatory connections between layer 2 and layer 3 pyramidal cells, synaptic and probably also gap junctional, the response of pyramidal cell subnetworks cannot be so approximated. Because of the highly structured anatomy of entorhinal output projections, our model suggests that downstream targets of entorhinal cortex (dentate gyrus, hippocampal CA3, CA1, piriform cortex, olfactory bulb) receive differentially processed information.

Repetitively burst-spiking neurons in reeler mice show conserved but also highly variable morphological features of layer Vb-fated “thick-tufted” pyramidal cells

Reelin is a large extracellular glycoprotein that is secreted by Cajal-Retzius cells during embryonic development to regulate neuronal migration and cell proliferation but it also seems to regulate ion channel distribution and synaptic vesicle release properties of excitatory neurons well into adulthood. Mouse mutants with a compromised reelin signaling cascade show a highly disorganized neocortex but the basic connectional features of the displaced excitatory principal cells seem to be relatively intact. Very little is known, however, about the intrinsic electrophysiological and morphological properties of individual cells in the reeler cortex. Repetitive burst-spiking (RB) is a unique property of large, thick-tufted pyramidal cells of wild-type layer Vb exclusively, which project to several subcortical targets. In addition, they are known to possess sparse but far-reaching intracortical recurrent collaterals. Here, we compared the electrophysiological properties and morphological features of neurons in the reeler primary somatosensory cortex with those of wild-type controls. Whereas in wild-type mice, RB pyramidal cells were only detected in layer Vb, and the vast majority of reeler RB pyramidal cells were found in the superficial third of the cortical depth. There were no obvious differences in the intrinsic electrophysiological properties and basic morphological features (such as soma size or the number of dendrites) were also well preserved. However, the spatial orientation of the entire dendritic tree was highly variable in the reeler neocortex, whereas it was completely stereotyped in wild-type mice. It seems that basic quantitative features of layer Vb-fated RB pyramidal cells are well conserved in the highly disorganized mutant neocortex, whereas qualitative morphological features vary, possibly to properly orient toward the appropriate input pathways, which are known to show an atypical oblique path through the reeler cortex. The oblique dendritic orientation thus presumably reflects a re-orientation of dendritic input domains toward spatially highly disorganized afferent projections.

3D Printed Skull Cap and Benchtop Fabricated Microwire-Based Microelectrode Array for Custom Rat Brain Recordings

Microwire microelectrode arrays (MEAs) have been a popular low-cost tool for chronic electrophysiological recordings and are an inexpensive means to record the electrical dynamics crucial to brain function. However, both the fabrication and implantation procedures for multi-MEAs on a single rodent are time-consuming and the accuracy and quality are highly manual skill-dependent. To address the fabrication and implantation challenges for microwire MEAs, (1) a computer-aided designed and 3D printed skull cap for the pre-determined implantation locations of each MEA and (2) a benchtop fabrication approach for low-cost custom microwire MEAs were developed. A proof-of-concept design of a 32-channel 4-MEA (8-wire each) recording system was prototyped and tested through Sprague Dawley rat recordings. The skull cap design, based on the CT-scan of a single rat conforms well with multiple Sprague Dawley rats of various sizes, ages, and weight with a minimal bregma alignment error (A/P axis standard error of the mean = 0.25 mm, M/L axis standard error of the mean = 0.07 mm, n = 6). The prototyped 32-channel system was able to record the spiking activities over five months. The developed benchtop fabrication method and the 3D printed skull cap implantation platform would enable neuroscience groups to conduct in-house design, fabrication, and implantation of customizable microwire MEAs at a lower cost than the current commercial options and experience a shorter lead time for the design modifications and iterations.

Synaptic alterations in visual cortex reshape contrast-dependent gamma oscillations and inhibition-excitation ratio in a genetic mouse model of migraine

Background

Migraine affects a significant fraction of the world population, yet its etiology is not completely understood. In vitro results highlighted thalamocortical and intra-cortical glutamatergic synaptic gain-of-function associated with a monogenic form of migraine (familial-hemiplegic-migraine-type-1: FHM1). However, how these alterations reverberate on cortical activity remains unclear. As altered responsivity to visual stimuli and abnormal processing of visual sensory information are common hallmarks of migraine, herein we investigated the effects of FHM1-driven synaptic alterations in the visual cortex of awake mice.

Methods

We recorded extracellular field potentials from the primary visual cortex (V1) of head-fixed awake FHM1 knock-in (n = 12) and wild type (n = 12) mice in response to square-wave gratings with different visual contrasts. Additionally, we reproduced in silico the obtained experimental results with a novel spiking neurons network model of mouse V1, by implementing in the model both the synaptic alterations characterizing the FHM1 genetic mouse model adopted.

Results

FHM1 mice displayed similar amplitude but slower temporal evolution of visual evoked potentials. Visual contrast stimuli induced a lower increase of multi-unit activity in FHM1 mice, while the amount of information content about contrast level remained, however, similar to WT.

Spectral analysis of the local field potentials revealed an increase in the β/low γ range of WT mice following the abrupt reversal of contrast gratings. Such frequency range transitioned to the high γ range in FHM1 mice. Despite this change in the encoding channel, these oscillations preserved the amount of information conveyed about visual contrast. The computational model showed how these network effects may arise from a combination of changes in thalamocortical and intra-cortical synaptic transmission, with the former inducing a lower cortical activity and the latter inducing the higher frequencies ɣ oscillations.

Conclusions

Contrast-driven ɣ modulation in V1 activity occurs at a much higher frequency in FHM1. This is likely to play a role in the altered processing of visual information. Computational studies suggest that this shift is specifically due to enhanced cortical excitatory transmission. Our network model can help to shed light on the relationship between cellular and network levels of migraine neural alterations.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s10194-022-01495-9.

Dynamics of phase oscillator networks with synaptic weight and structural plasticity

We study the dynamics of Kuramoto oscillator networks with two distinct adaptation processes, one varying the coupling strengths and the other altering the network structure. Such systems model certain networks of oscillatory neurons where the neuronal dynamics, synaptic weights, and network structure interact with and shape each other. We model synaptic weight adaptation with spike-timing-dependent plasticity (STDP) that runs on a longer time scale than neuronal spiking. Structural changes that include addition and elimination of contacts occur at yet a longer time scale than the weight adaptations. First, we study the steady-state dynamics of Kuramoto networks that are bistable and can settle in synchronized or desynchronized states. To compare the impact of adding structural plasticity, we contrast the network with only STDP to one with a combination of STDP and structural plasticity. We show that the inclusion of structural plasticity optimizes the synchronized state of a network by allowing for synchronization with fewer links than a network with STDP alone. With non-identical units in the network, the addition of structural plasticity leads to the emergence of correlations between the oscillators' natural frequencies and node degrees. In the desynchronized regime, the structural plasticity decreases the number of contacts, leading to a sparse network. In this way, adding structural plasticity strengthens both synchronized and desynchronized states of a network. Second, we use desynchronizing coordinated reset stimulation and synchronizing periodic stimulation to induce desynchronized and synchronized states, respectively. Our findings indicate that a network with a combination of STDP and structural plasticity may require stronger and longer stimulation to switch between the states than a network with STDP only.

Pre- and postsynaptically expressed spike-timing-dependent plasticity contribute differentially to neuronal learning

A plethora of experimental studies have shown that long-term synaptic plasticity can be expressed pre- or postsynaptically depending on a range of factors such as developmental stage, synapse type, and activity patterns. The functional consequences of this diversity are not clear, although it is understood that whereas postsynaptic expression of plasticity predominantly affects synaptic response amplitude, presynaptic expression alters both synaptic response amplitude and short-term dynamics. In most models of neuronal learning, long-term synaptic plasticity is implemented as changes in connective weights. The consideration of long-term plasticity as a fixed change in amplitude corresponds more closely to post- than to presynaptic expression, which means theoretical outcomes based on this choice of implementation may have a postsynaptic bias. To explore the functional implications of the diversity of expression of long-term synaptic plasticity, we adapted a model of long-term plasticity, more specifically spike-timing-dependent plasticity (STDP), such that it was expressed either independently pre- or postsynaptically, or in a mixture of both ways. We compared pair-based standard STDP models and a biologically tuned triplet STDP model, and investigated the outcomes in a minimal setting, using two different learning schemes: in the first, inputs were triggered at different latencies, and in the second a subset of inputs were temporally correlated. We found that presynaptic changes adjusted the speed of learning, while postsynaptic expression was more efficient at regulating spike timing and frequency. When combining both expression loci, postsynaptic changes amplified the response range, while presynaptic plasticity allowed control over postsynaptic firing rates, potentially providing a form of activity homeostasis. Our findings highlight how the seemingly innocuous choice of implementing synaptic plasticity by single weight modification may unwittingly introduce a postsynaptic bias in modelling outcomes. We conclude that pre- and postsynaptically expressed plasticity are not interchangeable, but enable complimentary functions.

Differences between functional properties of pre- or postsynaptically expressed long-term plasticity have not yet been explored in much detail. In this paper, we used minimalist models of STDP with different expression loci, in search of fundamental functional consequences. Biologically, presynaptic expression acts mostly on neurotransmitter release, thereby altering short-term synaptic dynamics, whereas postsynaptic expression affects mainly synaptic gain. We compared models where plasticity was expressed only presynaptically or postsynaptically, or in both ways. We found that postsynaptic plasticity had a bigger impact over response times, while both pre- and postsynaptic plasticity were similarly capable of detecting correlated inputs. A model with biologically tuned expression of plasticity achieved the same outcome over a range of frequencies. Also, postsynaptic spiking frequency was not directly affected by presynaptic plasticity of short-term plasticity alone, however in combination with a postsynaptic component, it helped restrain positive feedback, contributing to activity homeostasis. In conclusion, expression locus may determine affinity for distinct coding schemes while also contributing to keep activity within bounds. Our findings highlight the importance of carefully implementing expression of plasticity in biological modelling, since the locus of expression may affect functional outcomes in simulations.

A calcium-based plasticity model for predicting long-term potentiation and depression in the neocortex

Pyramidal cells (PCs) form the backbone of the layered structure of the neocortex, and plasticity of their synapses is thought to underlie learning in the brain. However, such long-term synaptic changes have been experimentally characterized between only a few types of PCs, posing a significant barrier for studying neocortical learning mechanisms. Here we introduce a model of synaptic plasticity based on data-constrained postsynaptic calcium dynamics, and show in a neocortical microcircuit model that a single parameter set is sufficient to unify the available experimental findings on long-term potentiation (LTP) and long-term depression (LTD) of PC connections. In particular, we find that the diverse plasticity outcomes across the different PC types can be explained by cell-type-specific synaptic physiology, cell morphology and innervation patterns, without requiring type-specific plasticity. Generalizing the model to in vivo extracellular calcium concentrations, we predict qualitatively different plasticity dynamics from those observed in vitro. This work provides a first comprehensive null model for LTP/LTD between neocortical PC types in vivo, and an open framework for further developing models of cortical synaptic plasticity.

Intrinsic excitability mechanisms of neuronal ensemble formation

Neuronal ensembles are coactive groups of cortical neurons, found in spontaneous and evoked activity, that can mediate perception and behavior. To understand the mechanisms that lead to the formation of ensembles, we co-activated layer 2/3 pyramidal neurons in brain slices from mouse visual cortex, in animals of both sexes, replicating in vitro an optogenetic protocol to generate ensembles in vivo. Using whole-cell and perforated patch-clamp pair recordings we found that, after optogenetic or electrical stimulation, coactivated neurons increased their correlated activity, a hallmark of ensemble formation. Coactivated neurons showed small biphasic changes in presynaptic plasticity, with an initial depression followed by a potentiation after a recovery period. Optogenetic and electrical stimulation also induced significant increases in frequency and amplitude of spontaneous EPSPs, even after single-cell stimulation. In addition, we observed unexpected strong and persistent increases in neuronal excitability after stimulation, with increases in membrane resistance and reductions in spike threshold. A pharmacological agent that blocks changes in membrane resistance reverted this effect. These significant increases in excitability can explain the observed biphasic synaptic plasticity. We conclude that cell-intrinsic changes in excitability are involved in the formation of neuronal ensembles. We propose an ‘iceberg’ model, by which increased neuronal excitability makes subthreshold connections suprathreshold, enhancing the effect of already existing synapses, and generating a new neuronal ensemble.

In the brain, groups of neurons that are activated together – also known as neuronal ensembles – are the basic units that underpin perception and behavior. Yet, exactly how these coactive circuits are established remains under investigation.

In 1949, Canadian psychologist Donald Hebb proposed that, when brains learn something new, the neurons which are activated together connect to form ensembles, and their connections become stronger each time this specific piece of knowledge is recalled. This idea that ‘neurons that fire together, wire together’ can explain how memories are acquired and recalled, by strengthening their wiring.

However, recent studies have questioned whether strengthening connections is the only mechanism by which neural ensembles can be created. Changes in the excitability of neurons (how easily they are to fire and become activated) may also play a role. In other words, ensembles could emerge because certain neurons become more excitable and fire more readily.

To solve this conundrum, Alejandre-García et al. examined both hypotheses in the same system. Neurons in slices of the mouse visual cortex were stimulated electrically or optically, via a technique that controls neural activity with light. The activity of individual neurons and their connections was then measured with electrodes.

Spontaneous activity among connected neurons increased after stimulation, indicative of the formation of neuronal ensembles. Connected neurons also showed small changes in the strength of their connections, which first decreased and then rebounded after an initial recovery period.

Intriguingly, cells also showed unexpected strong and persistent increases in neuronal excitability after stimulation, such that neurons fired more readily to the same stimulus. In other words, neurons maintained a cellular memory of having been stimulated. The authors conclude that ensembles form because connected neurons become more excitable, which in turn, may strengthen connections of the circuit at a later stage.

These results provide fresh insights about the neural circuits underpinning learning and memory. In time, the findings could also help to understand disorders such as Alzheimer’s disease and schizophrenia, which are characterised by memory impairments and disordered thinking.

Storage capacity of networks with discrete synapses and sparsely encoded memories

Attractor neural networks are one of the leading theoretical frameworks for the formation and retrieval of memories in networks of biological neurons. In this framework, a pattern imposed by external inputs to the network is said to be learned when this pattern becomes a fixed point attractor of the network dynamics. The storage capacity is the maximum number of patterns that can be learned by the network. In this paper, we study the storage capacity of fully connected and sparsely connected networks with a binarized Hebbian rule, for arbitrary coding levels. Our results show that a network with discrete synapses has a similar storage capacity as the model with continuous synapses, and that this capacity tends asymptotically towards the optimal capacity, in the space of all possible binary connectivity matrices, in the sparse coding limit. We also derive finite coding level corrections for the asymptotic solution in the sparse coding limit. The result indicates the capacity of networks with Hebbian learning rules converges to the optimal capacity extremely slowly when the coding level becomes small. Our results also show that in networks with sparse binary connectivity matrices, the information capacity per synapse is larger than in the fully connected case, and thus such networks store information more efficiently.

Computational Concepts for Reconstructing and Simulating Brain Tissue

It has previously been shown that it is possible to derive a new class of biophysically detailed brain tissue models when one computationally analyzes and exploits the interdependencies or the multi-modal and multi-scale organization of the brain. These reconstructions, sometimes referred to as digital twins, enable a spectrum of scientific investigations. Building such models has become possible because of increase in quantitative data but also advances in computational capabilities, algorithmic and methodological innovations. This chapter presents the computational science concepts that provide the foundation to the data-driven approach to reconstructing and simulating brain tissue as developed by the EPFL Blue Brain Project, which was originally applied to neocortical microcircuitry and extended to other brain regions. Accordingly, the chapter covers aspects such as a knowledge graph-based data organization and the importance of the concept of a dataset release. We illustrate algorithmic advances in finding suitable parameters for electrical models of neurons or how spatial constraints can be exploited for predicting synaptic connections. Furthermore, we explain how in silico experimentation with such models necessitates specific addressing schemes or requires strategies for an efficient simulation. The entire data-driven approach relies on the systematic validation of the model. We conclude by discussing complementary strategies that not only enable judging the fidelity of the model but also form the basis for its systematic refinements.

The impact of neuron morphology on cortical network architecture

The neurons in the cerebral cortex are not randomly interconnected. This specificity in wiring can result from synapse formation mechanisms that connect neurons, depending on their electrical activity and genetically defined identity. Here, we report that the morphological properties of the neurons provide an additional prominent source by which wiring specificity emerges in cortical networks. This morphologically determined wiring specificity reflects similarities between the neurons’ axo-dendritic projections patterns, the packing density, and the cellular diversity of the neuropil. The higher these three factors are, the more recurrent is the topology of the network. Conversely, the lower these factors are, the more feedforward is the network’s topology. These principles predict the empirically observed occurrences of clusters of synapses, cell type-specific connectivity patterns, and nonrandom network motifs. Thus, we demonstrate that wiring specificity emerges in the cerebral cortex at subcellular, cellular, and network scales from the specific morphological properties of its neuronal constituents.

•

Neuronal network architectures reflect the morphologies of their constituents

•

Morphology predicts nonrandom connectivity from subcellular to network scales

•

Morphology predicts connectivity patterns consistent with those observed empirically

•

Neuron morphology is a major source for wiring specificity in the cerebral cortex

Revealing directed effective connectivity of cortical neuronal networks from measurements

In the study of biological networks, one of the major challenges is to understand the relationships between network structure and dynamics. In this paper, we model in vitro cortical neuronal cultures as stochastic dynamical systems and apply a method that reconstructs directed networks from dynamics [Ching and Tam, Phys. Rev. E 95, 010301(R) (2017)2470-004510.1103/PhysRevE.95.010301] to reveal directed effective connectivity, namely, the directed links and synaptic weights, of the neuronal cultures from voltage measurements recorded by a multielectrode array. The effective connectivity so obtained reproduces several features of cortical regions in rats and monkeys and has similar network properties as the synaptic network of the nematode Caenorhabditis elegans, whose entire nervous system has been mapped out. The distribution of the incoming degree is bimodal and the distributions of the average incoming and outgoing synaptic strength are non-Gaussian with long tails. The effective connectivity captures different information from the commonly studied functional connectivity, estimated using statistical correlation between spiking activities. The average synaptic strengths of excitatory incoming and outgoing links are found to increase with the spiking activity in the estimated effective connectivity but not in the functional connectivity estimated using the same sets of voltage measurements. These results thus demonstrate that the reconstructed effective connectivity can capture the general properties of synaptic connections and better reveal relationships between network structure and dynamics.

Flexible Neural Network Realized by the Probabilistic SiOx Memristive Synaptic Array for Energy-Efficient Image Learning

The human brain's neural networks are sparsely connected via tunable and probabilistic synapses, which may be essential for performing energy-efficient cognitive and intellectual functions. In this sense, the implementation of a flexible neural network with probabilistic synapses is a first step toward realizing the ultimate energy-efficient computing framework. Here, inspired by the efficient threshold-tunable and probabilistic rod-to-rod bipolar synapses in the human visual system, a 16 × 16 crossbar array comprising the vertical form of gate-tunable probabilistic SiO_x memristive synaptic barristor utilizing the Si/graphene heterojunction is designed and fabricated. Controllable stochastic switching dynamics in this array are achieved via various input voltage pulse schemes. In particular, the threshold tunability via electrostatic gating enables the efficient in situ alteration of the probabilistic switching activation (P_Act ) from 0 to 1.0, and can even modulate the degree of the P_Act change. A drop-connected algorithm based on the P_Act is constructed and used to successfully classify the shapes of several fashion items. The suggested approach can decrease the learning energy by up to ≈2,116 times relative to that of the conventional all-to-all connected network while exhibiting a high recognition accuracy of ≈93 %.

Reconstruction of neocortex: Organelles, compartments, cells, circuits, and activity

We assembled a semi-automated reconstruction of L2/3 mouse primary visual cortex from ∼250 × 140 × 90 μm3 of electron microscopic images, including pyramidal and non-pyramidal neurons, astrocytes, microglia, oligodendrocytes and precursors, pericytes, vasculature, nuclei, mitochondria, and synapses. Visual responses of a subset of pyramidal cells are included. The data are publicly available, along with tools for programmatic and three-dimensional interactive access. Brief vignettes illustrate the breadth of potential applications relating structure to function in cortical circuits and neuronal cell biology. Mitochondria and synapse organization are characterized as a function of path length from the soma. Pyramidal connectivity motif frequencies are predicted accurately using a configuration model of random graphs. Pyramidal cells receiving more connections from nearby cells exhibit stronger and more reliable visual responses. Sample code shows data access and analysis.

Local connectivity and synaptic dynamics in mouse and human neocortex

We present a unique, extensive, and open synaptic physiology analysis platform and dataset. Through its application, we reveal principles that relate cell type to synaptic properties and intralaminar circuit organization in the mouse and human cortex. The dynamics of excitatory synapses align with the postsynaptic cell subclass, whereas inhibitory synapse dynamics partly align with presynaptic cell subclass but with considerable overlap. Synaptic properties are heterogeneous in most subclass-to-subclass connections. The two main axes of heterogeneity are strength and variability. Cell subclasses divide along the variability axis, whereas the strength axis accounts for substantial heterogeneity within the subclass. In the human cortex, excitatory-to-excitatory synaptic dynamics are distinct from those in the mouse cortex and vary with depth across layers 2 and 3.

Metastable dynamics of neural circuits and networks

Cortical neurons emit seemingly erratic trains of action potentials or "spikes," and neural network dynamics emerge from the coordinated spiking activity within neural circuits. These rich dynamics manifest themselves in a variety of patterns, which emerge spontaneously or in response to incoming activity produced by sensory inputs. In this Review, we focus on neural dynamics that is best understood as a sequence of repeated activations of a number of discrete hidden states. These transiently occupied states are termed "metastable" and have been linked to important sensory and cognitive functions. In the rodent gustatory cortex, for instance, metastable dynamics have been associated with stimulus coding, with states of expectation, and with decision making. In frontal, parietal, and motor areas of macaques, metastable activity has been related to behavioral performance, choice behavior, task difficulty, and attention. In this article, we review the experimental evidence for neural metastable dynamics together with theoretical approaches to the study of metastable activity in neural circuits. These approaches include (i) a theoretical framework based on non-equilibrium statistical physics for network dynamics; (ii) statistical approaches to extract information about metastable states from a variety of neural signals; and (iii) recent neural network approaches, informed by experimental results, to model the emergence of metastable dynamics. By discussing these topics, we aim to provide a cohesive view of how transitions between different states of activity may provide the neural underpinnings for essential functions such as perception, memory, expectation, or decision making, and more generally, how the study of metastable neural activity may advance our understanding of neural circuit function in health and disease.

Frequency- and spike-timing-dependent mitochondrial Ca2+ signaling regulates the metabolic rate and synaptic efficacy in cortical

Mitochondrial activity is crucial for the plasticity of central synapses, but how the firing pattern of pre- and postsynaptic neurons affects the mitochondria remains elusive. We recorded changes in the fluorescence of cytosolic and mitochondrial Ca²⁺ indicators in cell bodies, axons, and dendrites of cortical pyramidal neurons in mouse brain slices while evoking pre- and postsynaptic spikes. Postsynaptic spike firing elicited fast mitochondrial Ca²⁺ responses that were about threefold larger in the somas and apical dendrites than in basal dendrites and axons. The amplitude of these responses and metabolic activity were extremely sensitive to the firing frequency. Furthermore, while an EPSP alone caused no detectable Ca²⁺ elevation in the dendritic mitochondria, the coincidence of EPSP with a backpropagating spike produced prominent, highly localized mitochondrial Ca²⁺ hotspots. Our results indicate that mitochondria decode the spike firing frequency and the Hebbian temporal coincidences into the Ca²⁺ signals, which are further translated into the metabolic output and most probably lead to long-term changes in synaptic efficacy.

Input rate encoding and gain control in dendrites of neocortical pyramidal neurons

Elucidating how neurons encode network activity is essential to understanding how the brain processes information. Neocortical pyramidal cells receive excitatory input onto spines distributed along dendritic branches. Local dendritic branch nonlinearities can boost the response to spatially clustered and synchronous input, but how this translates into the integration of patterns of ongoing activity remains unclear. To examine dendritic integration under naturalistic stimulus regimes, we use two-photon glutamate uncaging to repeatedly activate multiple dendritic spines at random intervals. In the proximal dendrites of two populations of layer 5 pyramidal neurons in the mouse motor cortex, spatially restricted synchrony is not a prerequisite for dendritic boosting. Branches encode afferent inputs with distinct rate sensitivities depending upon cell and branch type. Thus, inputs distributed along a dendritic branch can recruit supralinear boosting and the window of this nonlinearity may provide a mechanism by which dendrites can preferentially amplify slow-frequency network oscillations.

Synaptic connectivity to L2/3 of primary visual cortex measured by two-photon optogenetic stimulation

Understanding cortical microcircuits requires thorough measurement of physiological properties of synaptic connections formed within and between diverse subclasses of neurons. Towards this goal, we combined spatially precise optogenetic stimulation with multicellular recording to deeply characterize intralaminar and translaminar monosynaptic connections to supragranular (L2/3) neurons in the mouse visual cortex. The reliability and specificity of multiphoton optogenetic stimulation were measured across multiple Cre lines, and measurements of connectivity were verified by comparison to paired recordings and targeted patching of optically identified presynaptic cells. With a focus on translaminar pathways, excitatory and inhibitory synaptic connections from genetically defined presynaptic populations were characterized by their relative abundance, spatial profiles, strength, and short-term dynamics. Consistent with the canonical cortical microcircuit, layer 4 excitatory neurons and interneurons within L2/3 represented the most common sources of input to L2/3 pyramidal cells. More surprisingly, we also observed strong excitatory connections from layer 5 intratelencephalic neurons and potent translaminar inhibition from multiple interneuron subclasses. The hybrid approach revealed convergence to and divergence from excitatory and inhibitory neurons within and across cortical layers. Divergent excitatory connections often spanned hundreds of microns of horizontal space. In contrast, divergent inhibitory connections were more frequently measured from postsynaptic targets near each other.

Metastable attractors explain the variable timing of stable behavioral action sequences

The timing of self-initiated actions shows large variability even when they are executed in stable, well-learned sequences. Could this mix of reliability and stochasticity arise within the same neural circuit? We trained rats to perform a stereotyped sequence of self-initiated actions and recorded neural ensemble activity in secondary motor cortex (M2), which is known to reflect trial-by-trial action-timing fluctuations. Using hidden Markov models, we established a dictionary between activity patterns and actions. We then showed that metastable attractors, representing activity patterns with a reliable sequential structure and large transition timing variability, could be produced by reciprocally coupling a high-dimensional recurrent network and a low-dimensional feedforward one. Transitions between attractors relied on correlated variability in this mesoscale feedback loop, predicting a specific structure of low-dimensional correlations that were empirically verified in M2 recordings. Our results suggest a novel mesoscale network motif based on correlated variability supporting naturalistic animal behavior.

Heterogeneity of glutamatergic synapses: cellular mechanisms and network consequences

Chemical synapses are commonly known as a structurally and functionally highly diverse class of cell-cell contacts specialized to mediate communication between neurons. They represent the smallest "computational" unit of the brain and are typically divided into excitatory and inhibitory as well as modulatory categories. These categories are subdivided into diverse types, each representing a different structure-function repertoire that in turn are thought to endow neuronal networks with distinct computational properties. The diversity of structure and function found among a given category of synapses is referred to as heterogeneity. The main building blocks for this heterogeneity are synaptic vesicles, the active zone, the synaptic cleft, the postsynaptic density, and glial processes associated with the synapse. Each of these five structural modules entails a distinct repertoire of functions, and their combination specifies the range of functional heterogeneity at mammalian excitatory synapses, which are the focus of this review. We describe synapse heterogeneity that is manifested on different levels of complexity ranging from the cellular morphology of the pre- and postsynaptic cells toward the expression of different protein isoforms at individual release sites. We attempt to define the range of structural building blocks that are used to vary the basic functional repertoire of excitatory synaptic contacts and discuss sources and general mechanisms of synapse heterogeneity. Finally, we explore the possible impact of synapse heterogeneity on neuronal network function.