A quantitative description of short-term plasticity at excitatory synapses in layer 2/3 of rat primary visual cortex

A robust balancing mechanism for spiking neural networks

Dynamical balance of excitation and inhibition is usually invoked to explain the irregular low firing activity observed in the cortex. We propose a robust nonlinear balancing mechanism for a random network of spiking neurons, which works also in the absence of strong external currents. Biologically, the mechanism exploits the plasticity of excitatory-excitatory synapses induced by short-term depression. Mathematically, the nonlinear response of the synaptic activity is the key ingredient responsible for the emergence of a stable balanced regime. Our claim is supported by a simple self-consistent analysis accompanied by extensive simulations performed for increasing network sizes. The observed regime is essentially fluctuation driven and characterized by highly irregular spiking dynamics of all neurons.

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.

Interpretation of presynaptic phenotypes of synaptic plasticity in terms of a two-step priming process

Studies on synaptic proteins involved in neurotransmitter release often aim at distinguishing between their roles in vesicle priming (the docking of synaptic vesicles to the plasma membrane and the assembly of a release machinery) as opposed to the process of vesicle fusion. This has traditionally been done by estimating two parameters, the size of the pool of fusion-competent vesicles (the readily releasable pool, RRP) and the probability that such vesicles are released by an action potential, with the aim of determining how these parameters are affected by molecular perturbations. Here, it is argued that the assumption of a homogeneous RRP may be too simplistic and may blur the distinction between vesicle priming and fusion. Rather, considering priming as a dynamic and reversible multistep process allows alternative interpretations of mutagenesis-induced changes in synaptic transmission and suggests mechanisms for variability in synaptic strength and short-term plasticity among synapses, as well as for interactions between short- and long-term plasticity. In many cases, assigned roles of proteins or causes for observed phenotypes are shifted from fusion- to priming-related when considering multistep priming. Activity-dependent enhancement of priming is an essential element in this alternative view and its variation among synapse types can explain why some synapses show depression and others show facilitation at low to intermediate stimulation frequencies. Multistep priming also suggests a mechanism for frequency invariance of steady-state release, which can be observed in some synapses involved in sensory processing.

Silences, spikes and bursts: Three-part knot of the neural code

When a neuron breaks silence, it can emit action potentials in a number of patterns. Some responses are so sudden and intense that electrophysiologists felt the need to single them out, labelling action potentials emitted at a particularly high frequency with a metonym - bursts. Is there more to bursts than a figure of speech? After all, sudden bouts of high-frequency firing are expected to occur whenever inputs surge. The burst coding hypothesis advances that the neural code has three syllables: silences, spikes and bursts. We review evidence supporting this ternary code in terms of devoted mechanisms for burst generation, synaptic transmission and synaptic plasticity. We also review the learning and attention theories for which such a triad is beneficial.

High-Endurance Long-Term Potentiation in Neuromorphic Organic Electrochemical Transistors by PEDOT:PSS Electrochemical Polymerization on the Gate Electrode

The brain exhibits extraordinary information processing capabilities thanks to neural networks that can operate in parallel with minimal energy consumption. Memory and learning require the creation of new neural networks through the long-term modification of the structure of the synapses, a phenomenon called long-term plasticity. Here, we use an organic electrochemical transistor to simulate long-term potentiation and depotentiation processes. Similarly to what happens in a synapse, the polymerization of the 3,4-ethylenedioxythiophene (EDOT) on the gate electrode modifies the structure of the device and boosts the ability of the gate potential to modify the conductivity of the channel. Operando AFM measurements were carried out to demonstrate the correlation between neuromorphic behavior and modification of the gate electrode. Long-term enhancement depends on both the number of pulses used and the gate potential, which generates long-term potentiation when a threshold of +0.7 V is overcome. Long-term depotentiation occurs by applying a +3.0 V potential and exploits the overoxidation of the deposited PEDOT:PSS. The induced states are stable for at least 2 months. The developed device shows very interesting characteristics in the field of neuromorphic electronics.

Input-specific synaptic depression shapes temporal integration in mouse visual cortex

Efficient sensory processing requires the nervous system to adjust to ongoing features of the environment. In primary visual cortex (V1), neuronal activity strongly depends on recent stimulus history. Existing models can explain effects of prolonged stimulus presentation but remain insufficient for explaining effects observed after shorter durations commonly encountered under natural conditions. We investigated the mechanisms driving adaptation in response to brief (100 ms) stimuli in L2/3 V1 neurons by performing in vivo whole-cell recordings to measure membrane potential and synaptic inputs. We find that rapid adaptation is generated by stimulus-specific suppression of excitatory and inhibitory synaptic inputs. Targeted optogenetic experiments reveal that these synaptic effects are due to input-specific short-term depression of transmission between layers 4 and 2/3. Thus, brief stimulus presentation engages a distinct adaptation mechanism from that previously reported in response to prolonged stimuli, enabling flexible control of sensory encoding across a wide range of timescales.

Short-term neuronal and synaptic plasticity act in synergy for deviance detection in spiking networks

Sensory areas of cortex respond more strongly to infrequent stimuli when these violate previously established regularities, a phenomenon known as deviance detection (DD). Previous modeling work has mainly attempted to explain DD on the basis of synaptic plasticity. However, a large fraction of cortical neurons also exhibit firing rate adaptation, an underexplored potential mechanism. Here, we investigate DD in a spiking neuronal network model with two types of short-term plasticity, fast synaptic short-term depression (STD) and slower threshold adaptation (TA). We probe the model with an oddball stimulation paradigm and assess DD by evaluating the network responses. We find that TA is sufficient to elicit DD. It achieves this by habituating neurons near the stimulation site that respond earliest to the frequently presented standard stimulus (local fatigue), which diminishes the response and promotes the recovery (global fatigue) of the wider network. Further, we find a synergy effect between STD and TA, where they interact with each other to achieve greater DD than the sum of their individual effects. We show that this synergy is caused by the local fatigue added by STD, which inhibits the global response to the frequently presented stimulus, allowing greater recovery of TA-mediated global fatigue and making the network more responsive to the deviant stimulus. Finally, we show that the magnitude of DD strongly depends on the timescale of stimulation. We conclude that highly predictable information can be encoded in strong local fatigue, which allows greater global recovery and subsequent heightened sensitivity for DD.

Distinct competitive impacts of palatability of taste stimuli on sampling dynamics during a preference test

Food or taste preference tests are analogous to naturalistic decisions in which the animal selects which stimuli to sample and for how long to sample them. The data acquired in such tests, the relative amounts of the alternative stimuli that are sampled and consumed, indicate the preference for each. While such preferences are typically recorded as a single quantity, an analysis of the ongoing sampling dynamics producing the preference can reveal otherwise hidden aspects of the decision-making process that depend on its underlying neural circuit mechanisms. Here, we perform a dynamic analysis of two factors that give rise to preferences in a two-alternative task, namely the distribution of durations of sampling bouts of each stimulus and the likelihood of returning to the same stimulus or switching to the alternative-that is, the transition probability-following each bout. The results of our analysis support a specific computational model of decision making whereby an exponential distribution of bout durations has a mean that is positively correlated with the palatability of that stimulus but also negatively correlated with the palatability of the alternative. This impact of the alternative stimulus on the distribution of bout durations decays over a timescale of tens of seconds, even though the memory of the alternative stimulus lasts far longer-long enough to impact the transition probabilities upon ending bouts. Together, our findings support a state transition model for bout durations and suggest a separate memory mechanism for stimulus selection. (PsycInfo Database Record (c) 2023 APA, all rights reserved).

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.

Neural adaptation and fractional dynamics as a window to underlying neural excitability

The relationship between macroscale electrophysiological recordings and the dynamics of underlying neural activity remains unclear. We have previously shown that low frequency EEG activity (<1 Hz) is decreased at the seizure onset zone (SOZ), while higher frequency activity (1-50 Hz) is increased. These changes result in power spectral densities (PSDs) with flattened slopes near the SOZ, which are assumed to be areas of increased excitability. We wanted to understand possible mechanisms underlying PSD changes in brain regions of increased excitability. We hypothesized that these observations are consistent with changes in adaptation within the neural circuit. We developed a theoretical framework and tested the effect of adaptation mechanisms, such as spike frequency adaptation and synaptic depression, on excitability and PSDs using filter-based neural mass models and conductance-based models. We compared the contribution of single timescale adaptation and multiple timescale adaptation. We found that adaptation with multiple timescales alters the PSDs. Multiple timescales of adaptation can approximate fractional dynamics, a form of calculus related to power laws, history dependence, and non-integer order derivatives. Coupled with input changes, these dynamics changed circuit responses in unexpected ways. Increased input without synaptic depression increases broadband power. However, increased input with synaptic depression may decrease power. The effects of adaptation were most pronounced for low frequency activity (< 1Hz). Increased input combined with a loss of adaptation yielded reduced low frequency activity and increased higher frequency activity, consistent with clinical EEG observations from SOZs. Spike frequency adaptation and synaptic depression, two forms of multiple timescale adaptation, affect low frequency EEG and the slope of PSDs. These neural mechanisms may underlie changes in EEG activity near the SOZ and relate to neural hyperexcitability. Neural adaptation may be evident in macroscale electrophysiological recordings and provide a window to understanding neural circuit excitability.

Input-specific synaptic depression shapes temporal integration in mouse visual cortex

Efficient sensory processing requires the nervous system to adjust to ongoing features of the environment. In primary visual cortex (V1), neuronal activity strongly depends on recent stimulus history. Existing models can explain effects of prolonged stimulus presentation, but remain insufficient for explaining effects observed after shorter durations commonly encountered under natural conditions. We investigated the mechanisms driving adaptation in response to brief (100 ms) stimuli in L2/3 V1 neurons by performing in vivo whole-cell recordings to measure membrane potential and synaptic inputs. We find that rapid adaptation is generated by stimulus-specific suppression of excitatory and inhibitory synaptic inputs. Targeted optogenetic experiments reveal that these synaptic effects are due to input-specific short-term depression of transmission between layers 4 and 2/3. Thus, distinct mechanisms are engaged following brief and prolonged stimulus presentation and together enable flexible control of sensory encoding across a wide range of time scales.

Developmental disruption and restoration of brain synaptome architecture in the murine Pax6 neurodevelopmental disease model

Neurodevelopmental disorders of genetic origin delay the acquisition of normal abilities and cause disabling phenotypes. Nevertheless, spontaneous attenuation and even complete amelioration of symptoms in early childhood and adolescence can occur in many disorders, suggesting that brain circuits possess an intrinsic capacity to overcome the deficits arising from some germline mutations. We examined the molecular composition of almost a trillion excitatory synapses on a brain-wide scale between birth and adulthood in mice carrying a mutation in the homeobox transcription factor Pax6, a neurodevelopmental disorder model. Pax6 haploinsufficiency had no impact on total synapse number at any age. By contrast, the molecular composition of excitatory synapses, the postnatal expansion of synapse diversity and the acquisition of normal synaptome architecture were delayed in all brain regions, interfering with networks and electrophysiological simulations of cognitive functions. Specific excitatory synapse types and subtypes were affected in two key developmental age-windows. These phenotypes were reversed within 2-3 weeks of onset, restoring synapse diversity and synaptome architecture to the normal developmental trajectory. Synapse subtypes with rapid protein turnover mediated the synaptome remodeling. This brain-wide capacity for remodeling of synapse molecular composition to recover and maintain the developmental trajectory of synaptome architecture may help confer resilience to neurodevelopmental genetic disorders.

Dynamic branching in a neural network model for probabilistic prediction of sequences

An important function of the brain is to predict which stimulus is likely to occur based on the perceived cues. The present research studied the branching behavior of a computational network model of populations of excitatory and inhibitory neurons, both analytically and through simulations. Results show how synaptic efficacy, retroactive inhibition and short-term synaptic depression determine the dynamics of selection between different branches predicting sequences of stimuli of different probabilities. Further results show that changes in the probability of the different predictions depend on variations of neuronal gain. Such variations allow the network to optimize the probability of its predictions to changing probabilities of the sequences without changing synaptic efficacy.

Classical-Contextual Interactions in V1 May Rely on Dendritic Computations

A signature feature of the neocortex is the dense network of horizontal connections (HCs) through which pyramidal neurons (PNs) exchange "contextual" information. In primary visual cortex (V1), HCs are thought to facilitate boundary detection, a crucial operation for object recognition, but how HCs modulate PN responses to boundary cues within their classical receptive fields (CRF) remains unknown. We began by "asking" natural images, through a structured data collection and ground truth labeling process, what function a V1 cell should use to compute boundary probability from aligned edge cues within and outside its CRF. The "answer" was an asymmetric 2-D sigmoidal function, whose nonlinear form provides the first normative account for the "multiplicative" center-flanker interactions previously reported in V1 neurons (Kapadia et al., 1995, 2000; Polat et al., 1998). Using a detailed compartmental model, we then show that this boundary-detecting classical-contextual interaction function can be computed by NMDAR-dependent spatial synaptic interactions within PN dendrites - the site where classical and contextual inputs first converge in the cortex. In additional simulations, we show that local interneuron circuitry activated by HCs can powerfully leverage the nonlinear spatial computing capabilities of PN dendrites, providing the cortex with a highly flexible substrate for integration of classical and contextual information.

Dynamics of Neuromuscular Transmission Reproduced by Calcium-Dependent and Reversible Serial Transitions in the Vesicle Fusion Complex

Neuromuscular transmission, from spontaneous release to facilitation and depression, was accurately reproduced by a mechanistic kinetic model of sequential maturation transitions in the molecular fusion complex. The model incorporates three predictions. First, calcium-dependent forward transitions take vesicles from docked to preprimed to primed states, followed by fusion. Second, prepriming and priming are reversible. Third, fusion and recycling are unidirectional. The model was fed with experimental data from previous studies, whereas the backward (β) and recycling (ρ) rate constant values were fitted. Classical experiments were successfully reproduced with four transition states in the model when every forward (α) rate constant had the same value, and both backward rate constants were 50–100 times larger. Such disproportion originated an abruptly decreasing gradient of resting vesicles from docked to primed states. By contrast, a three-state version of the model failed to reproduce the dynamics of transmission by using the same set of parameters. Simulations predict the following: (1) Spontaneous release reflects primed to fusion spontaneous transitions. (2) Calcium elevations synchronize the series of forward transitions that lead to fusion. (3) Facilitation reflects a transient increase of priming following the calcium-dependent maturation transitions. (4) The calcium sensors that produce facilitation are those that evoke the transitions form docked to primed states. (5) Backward transitions and recycling restore the resting state. (6) Depression reflects backward transitions and slow recycling after intense release. Altogether, our results predict that fusion is produced by one calcium sensor, whereas the modulation of the number of vesicles that fuse depends on the calcium sensors that promote the early transition states. Such finely tuned kinetics offers a mechanism for collective non-linear transitional adaptations of a homogeneous vesicle pool to the ever-changing pattern of electrical activity in the neuromuscular junction.

Paired Feed-Forward Excitation With Delayed Inhibition Allows High Frequency Computations Across Brain Regions

The transmission of high frequency temporal information across brain regions is critical to perception, but the mechanisms underlying such transmission remain unclear. Long-range projection patterns across brain areas are often comprised of paired feed-forward excitation followed closely by delayed inhibition, including the thalamic triad synapse, thalamic projections to cortex, and projections within the hippocampus. Previous studies have shown that these joint projections produce a shortened period of depolarization, sharpening the timing window over which the postsynaptic neuron can fire. Here we show that these projections can facilitate the transmission of high frequency computations even at frequencies that are highly filtered by neuronal membranes. This temporal facilitation occurred over a range of synaptic parameter values, including variations in synaptic strength, synaptic time constants, short-term synaptic depression, and the delay between excitation and inhibition. Further, these projections can coordinate computations across multiple network levels, even amid ongoing local activity. We suggest that paired feed-forward excitation and inhibition provide a hybrid signal-carrying both a value and a clock-like trigger-to allow circuits to be responsive to input whenever it arrives.

Novel stimuli evoke excess activity in the mouse primary visual cortex

Rapid detection and processing of stimulus novelty are key elements of adaptive behavior. Predictive coding theories postulate that novel stimuli should be encoded differently from familiar stimuli. Here, we show that the majority of neurons in layer 2/3 of the mouse primary visual cortex exhibit a significant excess response to novel visual stimuli. The distinction between novel and familiar images developed rapidly, requiring only a few repeated presentations. We show that this phenomenon can be described by a model of cascading adaptation. This ubiquitous mechanism makes it likely that similar computations could be carried out in many brain areas.

To explore how neural circuits represent novel versus familiar inputs, we presented mice with repeated sets of images with novel images sparsely substituted. Using two-photon calcium imaging to record from layer 2/3 neurons in the mouse primary visual cortex, we found that novel images evoked excess activity in the majority of neurons. This novelty response rapidly emerged, arising with a time constant of 2.6 ± 0.9 s. When a new image set was repeatedly presented, a majority of neurons had similarly elevated activity for the first few presentations, which decayed to steady state with a time constant of 1.4 ± 0.4 s. When we increased the number of images in the set, the novelty response’s amplitude decreased, defining a capacity to store ∼15 familiar images under our conditions. These results could be explained quantitatively using an adaptive subunit model in which presynaptic neurons have individual tuning and gain control. This result shows that local neural circuits can create different representations for novel versus familiar inputs using generic, widely available mechanisms.

Improved alpha-beta power reduction via combined electrical and ultrasonic stimulation in a parkinsonian cortex-basal ganglia-thalamus computational model

Objective. To investigate computationally the interaction of combined electrical and ultrasonic modulation of isolated neurons and of the parkinsonian cortex-basal ganglia-thalamus loop.Approach. Continuous-wave or pulsed electrical and ultrasonic neuromodulation is applied to isolated Otsuka plateau-potential generating subthalamic nucleus (STN) and Pospischil regular, fast and low-threshold spiking cortical cells in a temporally alternating or simultaneous manner. Similar combinations of electrical/ultrasonic waveforms are applied to a parkinsonian biophysical cortex-basal ganglia-thalamus neuronal network. Ultrasound-neuron interaction is modelled respectively for isolated neurons and the neuronal network with the NICE and SONIC implementations of the bilayer sonophore underlying mechanism. Reduction inα-βspectral energy is used as a proxy to express improvement in Parkinson's disease by insonication and electrostimulation.Main results. Simultaneous electro-acoustic stimulation achieves a given level of neuronal activity at lower intensities compared to the separate stimulation modalities. Conversely, temporally alternating stimulation with50 Hzelectrical and ultrasound pulses is capable of eliciting100 HzSTN firing rates. Furthermore, combination of ultrasound with hyperpolarizing currents can alter cortical cell relative spiking regimes. In the parkinsonian neuronal network, continuous-wave and pulsed ultrasound reduce pathological oscillations by different mechanisms. High-frequency pulsed separated electrical and ultrasonic deep brain stimulation (DBS) reduce pathologicalα-βpower by entraining STN-neurons. In contrast, continuous-wave ultrasound reduces pathological oscillations by silencing the STN. Compared to the separated stimulation modalities, temporally simultaneous or alternating electro-acoustic stimulation can achieve higher reductions inα-βpower for the same safety contraints on electrical/ultrasonic intensity.Significance. Focused ultrasound has the potential of becoming a non-invasive alternative of conventional DBS for the treatment of Parkinson's disease. Here, we elaborate on proposed benefits of combined electro-acoustic stimulation in terms of improved dynamic range, efficiency, spatial resolution, and neuronal selectivity.

Nonlinear transient amplification in recurrent neural networks with short-term plasticity

To rapidly process information, neural circuits have to amplify specific activity patterns transiently. How the brain performs this nonlinear operation remains elusive. Hebbian assemblies are one possibility whereby strong recurrent excitatory connections boost neuronal activity. However, such Hebbian amplification is often associated with dynamical slowing of network dynamics, non-transient attractor states, and pathological run-away activity. Feedback inhibition can alleviate these effects but typically linearizes responses and reduces amplification gain. Here, we study nonlinear transient amplification (NTA), a plausible alternative mechanism that reconciles strong recurrent excitation with rapid amplification while avoiding the above issues. NTA has two distinct temporal phases. Initially, positive feedback excitation selectively amplifies inputs that exceed a critical threshold. Subsequently, short-term plasticity quenches the run-away dynamics into an inhibition-stabilized network state. By characterizing NTA in supralinear network models, we establish that the resulting onset transients are stimulus selective and well-suited for speedy information processing. Further, we find that excitatory-inhibitory co-tuning widens the parameter regime in which NTA is possible in the absence of persistent activity. In summary, NTA provides a parsimonious explanation for how excitatory-inhibitory co-tuning and short-term plasticity collaborate in recurrent networks to achieve transient amplification.

Low frequency novel interictal EEG biomarker for localizing seizures and predicting outcomes

Localizing hyperexcitable brain tissue to treat focal seizures remains challenging. We want to identify the seizure onset zone from interictal EEG biomarkers. We hypothesize that a combination of interictal EEG biomarkers, including a novel low frequency marker, can predict mesial temporal involvement and can assist in prognosis related to surgical resections. Interictal direct current wide bandwidth invasive EEG recordings from 83 patients implanted with 5111 electrodes were retrospectively studied. Logistic regression was used to classify electrodes and patient outcomes. A feed-forward neural network was implemented to understand putative mechanisms. Interictal infraslow frequency EEG activity was decreased for seizure onset zone electrodes while faster frequencies such as delta (2-4 Hz) and beta-gamma (20-50 Hz) activity were increased. These spectral changes comprised a novel interictal EEG biomarker that was significantly increased for mesial temporal seizure onset zone electrodes compared to non-seizure onset zone electrodes. Interictal EEG biomarkers correctly classified mesial temporal seizure onset zone electrodes with a specificity of 87% and positive predictive value of 80%. These interictal EEG biomarkers also correctly classified patient outcomes after surgical resection with a specificity of 91% and positive predictive value of 87%. Interictal infraslow EEG activity is decreased near the seizure onset zone while higher frequency power is increased, which may suggest distinct underlying physiologic mechanisms. Narrowband interictal EEG power bands provide information about the seizure onset zone and can help predict mesial temporal involvement in seizure onset. Narrowband interictal EEG power bands may be less useful for predictions related to non-mesial temporal electrodes. Together with interictal epileptiform discharges and high-frequency oscillations, these interictal biomarkers may provide prognostic information prior to surgical resection. Computational modelling suggests changes in neural adaptation may be related to the observed low frequency power changes.

A model of naturalistic decision making in preference tests

Decisions as to whether to continue with an ongoing activity or to switch to an alternative are a constant in an animal’s natural world, and in particular underlie foraging behavior and performance in food preference tests. Stimuli experienced by the animal both impact the choice and are themselves impacted by the choice, in a dynamic back and forth. Here, we present model neural circuits, based on spiking neurons, in which the choice to switch away from ongoing behavior instantiates this back and forth, arising as a state transition in neural activity. We analyze two classes of circuit, which differ in whether state transitions result from a loss of hedonic input from the stimulus (an “entice to stay” model) or from aversive stimulus-input (a “repel to leave” model). In both classes of model, we find that the mean time spent sampling a stimulus decreases with increasing value of the alternative stimulus, a fact that we linked to the inclusion of depressing synapses in our model. The competitive interaction is much greater in “entice to stay” model networks, which has qualitative features of the marginal value theorem, and thereby provides a framework for optimal foraging behavior. We offer suggestions as to how our models could be discriminatively tested through the analysis of electrophysiological and behavioral data.

Many decisions are of the ilk of whether to continue sampling a stimulus or to switch to an alternative, a key feature of foraging behavior. We produce two classes of model for such stay-switch decisions, which differ in how decisions to switch stimuli can arise. In an “entice-to-stay” model, a reduction in the necessary positive stimulus input causes switching decisions. In a “repel-to-leave” model, a rise in aversive stimulus input produces a switch decision. We find that in tasks where the sampling of one stimulus follows another, adaptive biological processes arising from a highly hedonic stimulus can reduce the time spent at the following stimulus, by up to ten-fold in the “entice-to-stay” models. Along with potentially observable behavioral differences that could distinguish the classes of networks, we also found signatures in neural activity, such as oscillation of neural firing rates and a rapid change in rates preceding the time of choice to leave a stimulus. In summary, our model findings lead to testable predictions and suggest a neural circuit-based framework for explaining foraging choices.

Linear-nonlinear cascades capture synaptic dynamics

Short-term synaptic dynamics differ markedly across connections and strongly regulate how action potentials communicate information. To model the range of synaptic dynamics observed in experiments, we have developed a flexible mathematical framework based on a linear-nonlinear operation. This model can capture various experimentally observed features of synaptic dynamics and different types of heteroskedasticity. Despite its conceptual simplicity, we show that it is more adaptable than previous models. Combined with a standard maximum likelihood approach, synaptic dynamics can be accurately and efficiently characterized using naturalistic stimulation patterns. These results make explicit that synaptic processing bears algorithmic similarities with information processing in convolutional neural networks.

Synaptic Iontronic Devices for Brain-Mimicking Functions: Fundamentals and Applications

Inspired by the information transmission mechanism in the central nervous systems of life, synapse-mimicking devices have been designed and fabricated for the purpose of breaking the bottleneck of von Neumann architecture and realizing the construction of effective hardware-based artificial intelligence. In this case, synaptic iontronic devices, dealing with current information with ions instead of electrons, have attracted enormous scientific interests owing to their unique characteristics provided by ions, such as the designability of charge carriers and the diversity of chemical regulation. Herein, the basic conception, working mechanism, performance metrics, and advanced applications of synaptic iontronic devices based on three-terminal transistors and two-terminal memristors are systematically reviewed and comprehensively discussed. This Review provides a prospect on how to realize artificial synaptic functions based on the regulation of ions and raises a series of further challenges unsolved in this area.

Temporal interference stimulation targets deep brain regions by modulating neural oscillations

BACKGROUND

Temporal interference (TI) stimulation of the brain generates amplitude-modulated electric fields oscillating in the kHz range with the goal of non-invasive targeted deep brain stimulation. Yet, the current intensities required in human (sensitivity) to modulate deep brain activity and if superficial brain region are spared (selectivity) at these intensities remains unclear.

OBJECTIVE

We developed an experimentally constrained theory for TI sensitivity to kHz electric field given the attenuation by membrane low-pass filtering property, and for TI selectivity to deep structures given the distribution of modulated and unmodulated electric fields in brain.

METHODS

The electric field threshold to modulate carbachol-induced gamma oscillations in rat hippocampal slices was determined for unmodulated 0.05-2 kHz sine waveforms, and 5 Hz amplitude-modulated waveforms with 0.1-2 kHz carrier frequencies. The neuronal effects are replicated with a computational network model to explore the underlying mechanisms, and then coupled to a validated current-flow model of the human head.

RESULTS

Amplitude-modulated electric fields are stronger in deep brain regions, while unmodulated electric fields are maximal at the cortical regions. Both experiment and model confirmed the hypothesis that spatial selectivity of temporal interference stimulation depends on the phasic modulation of neural oscillations only in deep brain regions. Adaptation mechanism (e.g. GABA_b) enhanced sensitivity to amplitude modulated waveform in contrast to unmodulated kHz and produced selectivity in modulating gamma oscillation (i.e. Higher gamma modulation in amplitude modulated vs unmodulated kHz stimulation). Selection of carrier frequency strongly affected sensitivity to amplitude modulation stimulation. Amplitude modulated stimulation with 100 Hz carrier frequency required ∼5 V/m (corresponding to ∼13 mA at the scalp surface), whereas, 1 kHz carrier frequency ∼60 V/m (∼160 mA) and 2 kHz carrier frequency ∼80 V/m (∼220 mA) to significantly modulate gamma oscillation. Sensitivity is increased (scalp current required decreased) for theoretical neuronal membranes with faster time constants.

CONCLUSION

The TI sensitivity (current required at the scalp) depends on the neuronal membrane time-constant (e.g. axons) approaching the kHz carrier frequency. TI selectivity is governed by network adaption (e.g. GABA_b) that is faster than the amplitude-modulation frequency. Thus, we show neuronal and network oscillations time-constants determine the scalp current required and the selectivity achievable with TI in humans.

Short-term depression and long-term plasticity together tune sensitive range of synaptic plasticity

Synaptic efficacy is subjected to activity-dependent changes on short- and long time scales. While short-term changes decay over minutes, long-term modifications last from hours up to a lifetime and are thought to constitute the basis of learning and memory. Both plasticity mechanisms have been studied extensively but how their interaction shapes synaptic dynamics is little known. To investigate how both short- and long-term plasticity together control the induction of synaptic depression and potentiation, we used numerical simulations and mathematical analysis of a calcium-based model, where pre- and postsynaptic activity induces calcium transients driving synaptic long-term plasticity. We found that the model implementing known synaptic short-term dynamics in the calcium transients can be successfully fitted to long-term plasticity data obtained in visual- and somatosensory cortex. Interestingly, the impact of spike-timing and firing rate changes on plasticity occurs in the prevalent firing rate range, which is different in both cortical areas considered here. Our findings suggest that short- and long-term plasticity are together tuned to adapt plasticity to area-specific activity statistics such as firing rates.

Synaptic long-term plasticity, the long-lasting change in efficacy of connections between neurons, is believed to underlie learning and memory. Synapses furthermore change their efficacy reversibly in an activity-dependent manner on the subsecond time scale, referred to as short-term plasticity. It is not known how both synaptic plasticity mechanisms—long- and short-term—interact during activity epochs. To address this question, we used a biologically-inspired plasticity model in which calcium drives changes in synaptic efficacy. We applied the model to plasticity data from visual- and somatosensory cortex and found that synaptic changes occur in very different firing rate ranges, which correspond to the prevalent firing rates in both structures. Our results suggest that short- and long-term plasticity act in a well concerted fashion.

Attractor-state itinerancy in neural circuits with synaptic depression

Neural populations with strong excitatory recurrent connections can support bistable states in their mean firing rates. Multiple fixed points in a network of such bistable units can be used to model memory retrieval and pattern separation. The stability of fixed points may change on a slower timescale than that of the dynamics due to short-term synaptic depression, leading to transitions between quasi-stable point attractor states in a sequence that depends on the history of stimuli. To better understand these behaviors, we study a minimal model, which characterizes multiple fixed points and transitions between them in response to stimuli with diverse time- and amplitude-dependencies. The interplay between the fast dynamics of firing rate and synaptic responses and the slower timescale of synaptic depression makes the neural activity sensitive to the amplitude and duration of square-pulse stimuli in a nontrivial, history-dependent manner. Weak cross-couplings further deform the basins of attraction for different fixed points into intricate shapes. We find that while short-term synaptic depression can reduce the total number of stable fixed points in a network, it tends to strongly increase the number of fixed points visited upon repetitions of fixed stimuli. Our analysis provides a natural explanation for the system's rich responses to stimuli of different durations and amplitudes while demonstrating the encoding capability of bistable neural populations for dynamical features of incoming stimuli.

Cholinergic and Adenosinergic Modulation of Synaptic Release

In this review we will discuss the effect of two neuromodulatory transmitters, acetylcholine (ACh) and adenosine, on the synaptic release probability and short-term synaptic plasticity. ACh and adenosine differ fundamentally in the way they are released into the extracellular space. ACh is released mostly from synaptic terminals and axonal bouton of cholinergic neurons in the basal forebrain (BF). Its mode of action on synaptic release probability is complex because it activate both ligand-gated ion channels, so-called nicotinic ACh receptors and G-protein coupled muscarinic ACh receptors. In contrast, adenosine is released from both neurons and glia via nucleoside transporters or diffusion over the cell membrane in a non-vesicular, non-synaptic fashion; its receptors are exclusively G-protein coupled receptors. We show that ACh and adenosine effects are highly specific for an identified synaptic connection and depend mostly on the presynaptic but also on the postsynaptic receptor type and discuss the functional implications of these differences.

Short term depression, presynaptic inhibition and local neuron diversity play key functional roles in the insect antennal lobe

As the oldest, but least understood sensory system in evolution, the olfactory system represents one of the most challenging research targets in sensory neurobiology. Although a large number of computational models of the olfactory system have been proposed, they do not account for the diversity in physiology, connectivity of local neurons, and several recent discoveries in the insect antennal lobe, a major olfactory organ in insects. Recent studies revealed that the response of some projection neurons were reduced by application of a GABA antagonist, and that insects are sensitive to odor pulse frequency. To account for these observations, we propose a spiking neural circuit model of the insect antennal lobe. Based on recent anatomical and physiological studies, we included three sub-types of local neurons as well as synaptic short-term depression (STD) in the model and showed that the interaction between STD and local neurons resulted in frequency-sensitive responses. We further discovered that the unexpected response of the projection neurons to the GABA antagonist is the result of complex interactions between STD and presynaptic inhibition, which is required for enhancing sensitivity to odor stimuli. Finally, we found that odor discrimination is improved if the innervation of the local neurons in the glomeruli follows a specific pattern. Our findings suggest that STD, presynaptic inhibition and diverse physiology and connectivity of local neurons are not independent properties, but they interact to play key roles in the function of antennal lobes.

Distinct subtypes of inhibitory interneurons differentially promote the propagation of rate and temporal codes in the feedforward neural network

Sensory information is believed to be encoded in neuronal spikes using two different neural codes, the rate code (spike firing rate) and the temporal code (precisely-timed spikes). Since the sensory cortex has a highly hierarchical feedforward structure, sensory information-carrying neural codes should reliably propagate across the feedforward network (FFN) of the cortex. Experimental evidence suggests that inhibitory interneurons, such as the parvalbumin-positive (PV) and somatostatin-positive (SST) interneurons, that have distinctively different electrophysiological and synaptic properties, modulate the neural codes during sensory information processing in the cortex. However, how PV and SST interneurons impact on the neural code propagation in the cortical FFN is unknown. We address this question by building a five-layer FFN model consisting of a physiologically realistic Hodgkin-Huxley-type models of excitatory neurons and PV/SST interneurons at different ratios. In response to different firing rate inputs (20-80 Hz), a higher ratio of PV over SST interneurons promoted a reliable propagation of all ranges of firing rate inputs. In contrast, in response to a range of precisely-timed spikes in the form of pulse-packets [with a different number of spikes (α, 40-400 spikes) and degree of dispersion (σ, 0-20 ms)], a higher ratio of SST over PV interneurons promoted a reliable propagation of pulse-packets. Our simulation results show that PV and SST interneurons differentially promote a reliable propagation of the rate and temporal codes, respectively, indicating that the dynamic recruitment of PV and SST interneurons may play critical roles in a reliable propagation of sensory information-carrying neural codes in the cortical FFN.

Direct current stimulation boosts hebbian plasticity in vitro

BACKGROUND

There is evidence that transcranial direct current stimulation (tDCS) can improve learning performance. Arguably, this effect is related to long term potentiation (LTP), but the precise biophysical mechanisms remain unknown.

HYPOTHESIS

We propose that direct current stimulation (DCS) causes small changes in postsynaptic membrane potential during ongoing endogenous synaptic activity. The altered voltage dynamics in the postsynaptic neuron then modify synaptic strength via the machinery of endogenous voltage-dependent Hebbian plasticity. This hypothesis predicts that DCS should exhibit Hebbian properties, namely pathway specificity and associativity.

METHODS

We studied the effects of DCS applied during the induction of LTP in the CA1 region of rat hippocampal slices and using a biophysical computational model.

RESULTS

DCS enhanced LTP, but only at synapses that were undergoing plasticity, confirming that DCS respects Hebbian pathway specificity. When different synaptic pathways cooperated to produce LTP, DCS enhanced this cooperation, boosting Hebbian associativity. Further slice experiments and computer simulations support a model where polarization of postsynaptic pyramidal neurons drives these plasticity effects through endogenous Hebbian mechanisms. The model is able to reconcile several experimental results by capturing the complex interaction between the induced electric field, neuron morphology, and endogenous neural activity.

CONCLUSIONS

These results suggest that tDCS can enhance associative learning. We propose that clinical tDCS should be applied during tasks that induce Hebbian plasticity to harness this phenomenon, and that the effects should be task specific through their interaction with endogenous plasticity mechanisms. Models that incorporate brain state and plasticity mechanisms may help to improve prediction of tDCS outcomes.

Serotonin, norepinephrine, and acetylcholine differentially affect astrocytic potassium clearance to modulate somatosensory signaling in male mice

Changes in extracellular potassium ([K⁺ ]_e ) modulate neuronal networks via changes in membrane potential, voltage-gated channel activity, and alteration to transmission at the synapse. Given the limited extracellular space in the central nervous system, potassium clearance is crucial. As activity-induced potassium transients are rapidly managed by astrocytic Kir4.1 and astrocyte-specific Na⁺ /K⁺ -ATPase, any neurotransmitter/neuromodulator that can regulate their function may have indirect influence on network activity. Neuromodulators differentially affect cortical/thalamic networks to align sensory processing with differing behavioral states. Given serotonin (5HT), norepinephrine (NE), and acetylcholine (ACh) differentially affect spike frequency adaptation and signal fidelity ("signal-to-noise") in somatosensory cortex, we hypothesize that [K⁺ ]_e may be differentially regulated by the different neuromodulators to exert their individual effects on network function. This study aimed to compare effects of individually applied 5HT, NE, and ACh on regulating [K⁺ ]_e in connection to effects on cortical-evoked response amplitude and adaptation in male mice. Using extracellular field and K⁺ ion-selective recordings of somatosensory stimulation, we found that differential effects of 5HT, NE, and ACh on [K⁺ ]_e regulation mirrored differential effects on amplitude and adaptation. 5HT effects on transient K⁺ recovery, adaptation, and field post-synaptic potential amplitude were disrupted by barium (200 µM), whereas NE and ACh effects were disrupted by ouabain (1 µM) or iodoacetate (100 µM). Considering the impact [K⁺ ]_e can have on many network functions; it seems highly efficient that neuromodulators regulate [K⁺ ]_e to exert their many effects. This study provides functional significance for astrocyte-mediated buffering of [K⁺ ]_e in neuromodulator-mediated shaping of cortical network activity.

Disparities in Short-Term Depression Among Prefrontal Cortex Synapses Sustain Persistent Activity in a Balanced Network

Persistent activity of cue-representing neurons in the prefrontal cortex (PFC) is regarded as a neural basis for working memory. The contribution of short-term synaptic plasticity (STP) at different types of synapses comprising the cortical network to persistent activity, however, remains unclear. Characterizing STP at synapses of the rat PFC layer 5 network, we found that PFC synapses exhibit distinct STP patterns according to presynaptic and postsynaptic identities. Excitatory postsynaptic currents (EPSCs) from corticopontine (Cpn) neurons were well sustained throughout continued activity, with stronger depression at synapses onto fast-spiking interneurons than those onto pyramidal cells. Inhibitory postsynaptic currents (IPSCs) were sustained at a weaker level compared with EPSC from Cpn synapses. Computational modeling of a balanced network incorporating empirically observed STP revealed that little depression at recurrent excitatory synapses, combined with stronger depression at other synapses, could provide the PFC with a unique synaptic mechanism for the generation and maintenance of persistent activity.

Short-term depression shapes information transmission in a constitutively active GABAergic synapse

Short-term depression is a low-pass filter of synaptic information, reducing synaptic information transfer at high presynaptic firing frequencies. Consequently, during elevated presynaptic firing, little information passes to the postsynaptic neuron. However, many neurons fire at relatively high frequencies all the time. Does depression silence their synapses? We tested this apparent contradiction in the indirect pathway of the basal ganglia. Using numerical modeling and whole-cell recordings from single entopeduncular nucleus (EP) neurons in rat brain slices, we investigated how different firing rates of globus pallidus (GP) neurons affect information transmission to the EP. Whole-cell recordings showed significant variability in steady-state depression, which decreased as stimulation frequency increased. Modeling predicted that this variability would translate into different postsynaptic noise levels during constitutive presynaptic activity. Our simulations further predicted that individual GP-EP synapses mediate gain control. However, when we consider the integration of multiple inputs, the broad range of GP firing rates would enable different modes of information transmission. Finally, we predict that changes in dopamine levels can shift the action of GP neurons from rate coding to gain modulation. Our results thus demonstrate how short-term depression shapes information transmission in the basal ganglia in particular and via GABAergic synapses in general.

Cellular and synaptic phenotypes lead to disrupted information processing in Fmr1-KO mouse layer 4 barrel cortex

Sensory hypersensitivity is a common and debilitating feature of neurodevelopmental disorders such as Fragile X Syndrome (FXS). How developmental changes in neuronal function culminate in network dysfunction that underlies sensory hypersensitivities is unknown. By systematically studying cellular and synaptic properties of layer 4 neurons combined with cellular and network simulations, we explored how the array of phenotypes in Fmr1-knockout (KO) mice produce circuit pathology during development. We show that many of the cellular and synaptic pathologies in Fmr1-KO mice are antagonistic, mitigating circuit dysfunction, and hence may be compensatory to the primary pathology. Overall, the layer 4 network in the Fmr1-KO exhibits significant alterations in spike output in response to thalamocortical input and distorted sensory encoding. This developmental loss of layer 4 sensory encoding precision would contribute to subsequent developmental alterations in layer 4-to-layer 2/3 connectivity and plasticity observed in Fmr1-KO mice, and circuit dysfunction underlying sensory hypersensitivity.

Estimating the Readily-Releasable Vesicle Pool Size at Synaptic Connections in the Neocortex

Previous studies based on the 'Quantal Model' for synaptic transmission suggest that neurotransmitter release is mediated by a single release site at individual synaptic contacts in the neocortex. However, recent studies seem to contradict this hypothesis and indicate that multi-vesicular release (MVR) could better explain the synaptic response variability observed in vitro. In this study we present a novel method to estimate the number of release sites per synapse, also known as the size of the readily releasable pool (N_RRP), from paired whole-cell recordings of connections between layer 5 thick tufted pyramidal cell (L5_TTPC) in the juvenile rat somatosensory cortex. Our approach extends the work of Loebel et al. (2009) by leveraging a recently published data-driven biophysical model of neocortical tissue. Using this approach, we estimated N_RRP to be between two to three for synaptic connections between L5_TTPCs. To constrain N_RRP values for other connections in the microcircuit, we developed and validated a generalization approach using published data on the coefficient of variation (CV) of the amplitudes of post-synaptic potentials (PSPs) from literature and comparing them against in silico experiments. Our study predicts that transmitter release at synaptic connections in the neocortex could be mediated by MVR and provides a data-driven approach to constrain the MVR model parameters in the microcircuit.

Model-Based Inference of Synaptic Transmission

Synaptic computation is believed to underlie many forms of animal behavior. A correct identification of synaptic transmission properties is thus crucial for a better understanding of how the brain processes information, stores memories and learns. Recently, a number of new statistical methods for inferring synaptic transmission parameters have been introduced. Here we review and contrast these developments, with a focus on methods aimed at inferring both synaptic release statistics and synaptic dynamics. Furthermore, based on recent proposals we discuss how such methods can be applied to data across different levels of investigation: from intracellular paired experiments to in vivo network-wide recordings. Overall, these developments open the window to reliably estimating synaptic parameters in behaving animals.

Unsupervised Temporal Contiguity Experience Does Not Break the Invariance of Orientation Selectivity Across Spatial Frequency

The images projected onto the retina can vary widely for a single object. Despite these transformations primates can quickly and reliably recognize objects. At the neural level, transformation tolerance in monkey inferotemporal cortex is affected by the temporal contiguity statistics of the visual input. Here we investigated whether temporal contiguity learning also influences the basic feature detectors in lower levels of the visual hierarchy, in particular the independent coding of orientation and spatial frequency (SF) in primary visual cortex. Eight male Long Evans rats were repeatedly exposed to a temporal transition between two gratings that changed in SF and had either the same (control SF) or a different (swap SF) orientation. Electrophysiological evidence showed that the responses of single neurons during this exposure were sensitive to the change in orientation. Nevertheless, the tolerance of orientation selectivity for changes in SF was unaffected by the temporal contiguity manipulation, as observed in 239 single neurons isolated pre-exposure and 234 post-exposure. Temporal contiguity learning did not affect orientation selectivity in V1. The basic filter mechanisms that characterize V1 processing seem unaffected by temporal contiguity manipulations.

Discrimination of Motion Direction in a Robot Using a Phenomenological Model of Synaptic Plasticity

Recognizing and tracking the direction of moving stimuli is crucial to the control of much animal behaviour. In this study, we examine whether a bio-inspired model of synaptic plasticity implemented in a robotic agent may allow the discrimination of motion direction of real-world stimuli. Starting with a well-established model of short-term synaptic plasticity (STP), we develop a microcircuit motif of spiking neurons capable of exhibiting preferential and nonpreferential responses to changes in the direction of an orientation stimulus in motion. While the robotic agent processes sensory inputs, the STP mechanism introduces direction-dependent changes in the synaptic connections of the microcircuit, resulting in a population of units that exhibit a typical cortical response property observed in primary visual cortex (V1), namely, direction selectivity. Visually evoked responses from the model are then compared to those observed in multielectrode recordings from V1 in anesthetized macaque monkeys, while sinusoidal gratings are displayed on a screen. Overall, the model highlights the role of STP as a complementary mechanism in explaining the direction selectivity and applies these insights in a physical robot as a method for validating important response characteristics observed in experimental data from V1, namely, direction selectivity.

Effect of adenosine on short-term synaptic plasticity in mouse piriform cortex in vitro: adenosine acts as a high-pass filter

We examined the effect of adenosine and of adenosine A1 receptor blockage on short-term synaptic plasticity in slices of adult mouse anterior piriform cortex maintained in vitro in an in vivo-like ACSF. Extracellular recording of postsynaptic responses was performed in layer 1a while repeated electrical stimulation (5-pulse-trains, frequency between 3.125 and 100 Hz) was applied to the lateral olfactory tract. Our stimulation protocol was aimed at covering the frequency range of oscillatory activities observed in the olfactory bulb in vivo. In control condition, postsynaptic response amplitude showed a large enhancement for stimulation frequencies in the beta and gamma frequency range. A phenomenological model of short-term synaptic plasticity fitted to the data suggests that this frequency-dependent enhancement can be explained by the interplay between a short-term facilitation mechanism and two short-term depression mechanisms, with fast and slow recovery time constants. In the presence of adenosine, response amplitude evoked by low-frequency stimulation decreased in a dose-dependent manner (IC50  = 70 μmol/L). Yet short-term plasticity became more dominated by facilitation and less influenced by depression. Both changes compensated for the initial decrease in response amplitude in a way that depended on stimulation frequency: compensation was strongest at high frequency, up to restoring response amplitudes to values similar to those measured in control condition. The model suggested that the main effects of adenosine were to decrease neurotransmitter release probability and to attenuate short-term depression mechanisms. Overall, these results suggest that adenosine does not merely inhibit neuronal activity but acts in a more subtle, frequency-dependent manner.

The Virtual Electrode Recording Tool for EXtracellular Potentials (VERTEX) Version 2.0: Modelling in vitro electrical stimulation of brain tissue

Neuronal circuits can be modelled in detail allowing us to predict the effects of stimulation on individual neurons. Electrical stimulation of neuronal circuits in vitro and in vivo excites a range of neurons within the tissue and measurements of neural activity, e.g the local field potential (LFP), are again an aggregate of a large pool of cells. The previous version of our Virtual Electrode Recording Tool for EXtracellular Potentials (VERTEX) allowed for the simulation of the LFP generated by a patch of brain tissue. Here, we extend VERTEX to simulate the effect of electrical stimulation through a focal electric field. We observe both direct changes in neural activity and changes in synaptic plasticity. Testing our software in a model of a rat neocortical slice, we determine the currents contributing to the LFP, the effects of paired pulse stimulation to induce short term plasticity (STP), and the effect of theta burst stimulation (TBS) to induce long term potentiation (LTP).

A Single-Cell Level and Connectome-Derived Computational Model of the Drosophila Brain

Computer simulations play an important role in testing hypotheses, integrating knowledge, and providing predictions of neural circuit functions. While considerable effort has been dedicated into simulating primate or rodent brains, the fruit fly (Drosophila melanogaster) is becoming a promising model animal in computational neuroscience for its small brain size, complex cognitive behavior, and abundancy of data available from genes to circuits. Moreover, several Drosophila connectome projects have generated a large number of neuronal images that account for a significant portion of the brain, making a systematic investigation of the whole brain circuit possible. Supported by FlyCircuit (http://www.flycircuit.tw), one of the largest Drosophila neuron image databases, we began a long-term project with the goal to construct a whole-brain spiking network model of the Drosophila brain. In this paper, we report the outcome of the first phase of the project. We developed the Flysim platform, which (1) identifies the polarity of each neuron arbor, (2) predicts connections between neurons, (3) translates morphology data from the database into physiology parameters for computational modeling, (4) reconstructs a brain-wide network model, which consists of 20,089 neurons and 1,044,020 synapses, and (5) performs computer simulations of the resting state. We compared the reconstructed brain network with a randomized brain network by shuffling the connections of each neuron. We found that the reconstructed brain can be easily stabilized by implementing synaptic short-term depression, while the randomized one exhibited seizure-like firing activity under the same treatment. Furthermore, the reconstructed Drosophila brain was structurally and dynamically more diverse than the randomized one and exhibited both Poisson-like and patterned firing activities. Despite being at its early stage of development, this single-cell level brain model allows us to study some of the fundamental properties of neural networks including network balance, critical behavior, long-term stability, and plasticity.

A multi-scale layer-resolved spiking network model of resting-state dynamics in macaque visual cortical areas

Cortical activity has distinct features across scales, from the spiking statistics of individual cells to global resting-state networks. We here describe the first full-density multi-area spiking network model of cortex, using macaque visual cortex as a test system. The model represents each area by a microcircuit with area-specific architecture and features layer- and population-resolved connectivity between areas. Simulations reveal a structured asynchronous irregular ground state. In a metastable regime, the network reproduces spiking statistics from electrophysiological recordings and cortico-cortical interaction patterns in fMRI functional connectivity under resting-state conditions. Stable inter-area propagation is supported by cortico-cortical synapses that are moderately strong onto excitatory neurons and stronger onto inhibitory neurons. Causal interactions depend on both cortical structure and the dynamical state of populations. Activity propagates mainly in the feedback direction, similar to experimental results associated with visual imagery and sleep. The model unifies local and large-scale accounts of cortex, and clarifies how the detailed connectivity of cortex shapes its dynamics on multiple scales. Based on our simulations, we hypothesize that in the spontaneous condition the brain operates in a metastable regime where cortico-cortical projections target excitatory and inhibitory populations in a balanced manner that produces substantial inter-area interactions while maintaining global stability.

The mammalian cortex fulfills its complex tasks by operating on multiple temporal and spatial scales from single cells to entire areas comprising millions of cells. These multi-scale dynamics are supported by specific network structures at all levels of organization. Since models of cortex hitherto tend to concentrate on a single scale, little is known about how cortical structure shapes the multi-scale dynamics of the network. We here present dynamical simulations of a multi-area network model at neuronal and synaptic resolution with population-specific connectivity based on extensive experimental data which accounts for a wide range of dynamical phenomena. Our model elucidates relationships between local and global scales in cortex and provides a platform for future studies of cortical function.

Two Forms of Synaptic Depression Produced by Differential Neuromodulation of Presynaptic Calcium Channels

Neuromodulators are important regulators of synaptic transmission throughout the brain. At the presynaptic terminal, neuromodulation of calcium channels (CaVs) can affect transmission not only by changing neurotransmitter release probability, but also by shaping short-term plasticity (STP). Indeed, changes in STP are often considered a requirement for defining a presynaptic site of action. Nevertheless, some synapses exhibit non-canonical forms of neuromodulation, where release probability is altered without a corresponding change in STP. Here, we identify biophysical mechanisms whereby both canonical and non-canonical presynaptic neuromodulation can occur at the same synapse. At a subset of glutamatergic terminals in prefrontal cortex, GABAB and D1/D5 dopamine receptors suppress release probability with and without canonical increases in short-term facilitation by modulating different aspects of presynaptic CaV function. These findings establish a framework whereby signaling from multiple neuromodulators can converge on presynaptic CaVs to differentially tune release dynamics at the same synapse.

Architecture of the Mouse Brain Synaptome

Synapses are found in vast numbers in the brain and contain complex proteomes. We developed genetic labeling and imaging methods to examine synaptic proteins in individual excitatory synapses across all regions of the mouse brain. Synapse catalogs were generated from the molecular and morphological features of a billion synapses. Each synapse subtype showed a unique anatomical distribution, and each brain region showed a distinct signature of synapse subtypes. Whole-brain synaptome cartography revealed spatial architecture from dendritic to global systems levels and previously unknown anatomical features. Synaptome mapping of circuits showed correspondence between synapse diversity and structural and functional connectomes. Behaviorally relevant patterns of neuronal activity trigger spatiotemporal postsynaptic responses sensitive to the structure of synaptome maps. Areas controlling higher cognitive function contain the greatest synapse diversity, and mutations causing cognitive disorders reorganized synaptome maps. Synaptome technology and resources have wide-ranging application in studies of the normal and diseased brain.

A Distance-Dependent Distribution of Presynaptic Boutons Tunes Frequency-Dependent Dendritic Integration

How presynaptic inputs and neurotransmitter release dynamics are distributed along a dendritic tree is not well established. Here, we show that presynaptic boutons that form onto basal dendrites of CA1 pyramidal neurons display a decrease in active zone (AZ) size with distance from the soma, resulting in a distance-dependent increase in short-term facilitation. Our findings suggest that the spatial distribution of short-term facilitation serves to compensate for the electrotonic attenuation of subthreshold distal inputs during repeated stimulation and fine-tunes the preferred input frequency of dendritic domains.

•

Presynaptic inputs decrease in size with distance along CA1 basal dendrites

•

Release probability decreases with distance along basal dendrites

•

Short-term facilitation increases with distance along basal dendrites

•

Increased synaptic facilitation offsets passive decay and boosts supralinear events

Layer-Dependent Short-Term Synaptic Plasticity Between Excitatory Neurons in the C2 Barrel Column of Mouse Primary Somatosensory Cortex

Neurons process information through spatiotemporal integration of synaptic input. Synaptic transmission between any given pair of neurons is typically a dynamic process with presynaptic action potentials (APs) evoking depressing or facilitating postsynaptic potentials when presynaptic APs occur within hundreds of milliseconds of each other. In order to understand neocortical function, it is therefore important to investigate such short-term synaptic plasticity at synapses between different types of neocortical neurons. Here, we examine short-term synaptic dynamics between excitatory neurons in different layers of the mouse C2 barrel column through in vitro whole-cell recordings. We find layer-dependent short-term plasticity, with depression being dominant at many synaptic connections. Interestingly, however, presynaptic layer 2 neurons predominantly give rise to facilitating excitatory synaptic output at short interspike intervals of 10 and 30 ms. Previous studies have found prominent burst firing of excitatory neurons in supragranular layers of awake mice. The facilitation we observed in the synaptic output of layer 2 may, therefore, be functionally relevant, possibly serving to enhance the postsynaptic impact of burst firing.

NeuroMatic: An Integrated Open-Source Software Toolkit for Acquisition, Analysis and Simulation of Electrophysiological Data

Acquisition, analysis and simulation of electrophysiological properties of the nervous system require multiple software packages. This makes it difficult to conserve experimental metadata and track the analysis performed. It also complicates certain experimental approaches such as online analysis. To address this, we developed NeuroMatic, an open-source software toolkit that performs data acquisition (episodic, continuous and triggered recordings), data analysis (spike rasters, spontaneous event detection, curve fitting, stationarity) and simulations (stochastic synaptic transmission, synaptic short-term plasticity, integrate-and-fire and Hodgkin-Huxley-like single-compartment models). The merging of a wide range of tools into a single package facilitates a more integrated style of research, from the development of online analysis functions during data acquisition, to the simulation of synaptic conductance trains during dynamic-clamp experiments. Moreover, NeuroMatic has the advantage of working within Igor Pro, a platform-independent environment that includes an extensive library of built-in functions, a history window for reviewing the user's workflow and the ability to produce publication-quality graphics. Since its original release, NeuroMatic has been used in a wide range of scientific studies and its user base has grown considerably. NeuroMatic version 3.0 can be found at http://www.neuromatic.thinkrandom.com and https://github.com/SilverLabUCL/NeuroMatic.