Cholinergic modulation supports dynamic switching of resting state networks through selective DMN suppression.

Brain activity during the resting state is widely used to examine brain organization, cognition and alterations in disease states. While it is known that neuromodulation and the state of alertness impact resting-state activity, neural mechanisms behind such modulation of resting-state activity are unknown. In this work, we used a computational model to demonstrate that change in excitability and recurrent connections, due to cholinergic modulation, impacts resting-state activity. The results of such modulation in the model match closely with experimental work on direct cholinergic modulation of Default Mode Network (DMN) in rodents. We further extended our study to the human connectome derived from diffusion-weighted MRI. In human resting-state simulations, an increase in cholinergic input resulted in a brain-wide reduction of functional connectivity. Furthermore, selective cholinergic modulation of DMN closely captured experimentally observed transitions between the baseline resting state and states with suppressed DMN fluctuations associated with attention to external tasks. Our study thus provides insight into potential neural mechanisms for the effects of cholinergic neuromodulation on resting-state activity and its dynamics.

Sleep prevents catastrophic forgetting in spiking neural networks by forming a joint synaptic weight representation.

Artificial neural networks overwrite previously learned tasks when trained sequentially, a phenomenon known as catastrophic forgetting. In contrast, the brain learns continuously, and typically learns best when new training is interleaved with periods of sleep for memory consolidation. Here we used spiking network to study mechanisms behind catastrophic forgetting and the role of sleep in preventing it. The network could be trained to learn a complex foraging task but exhibited catastrophic forgetting when trained sequentially on different tasks. In synaptic weight space, new task training moved the synaptic weight configuration away from the manifold representing old task leading to forgetting. Interleaving new task training with periods of off-line reactivation, mimicking biological sleep, mitigated catastrophic forgetting by constraining the network synaptic weight state to the previously learned manifold, while allowing the weight configuration to converge towards the intersection of the manifolds representing old and new tasks. The study reveals a possible strategy of synaptic weights dynamics the brain applies during sleep to prevent forgetting and optimize learning.

Bidirectional Interaction of Hippocampal Ripples and Cortical Slow Waves Leads to Coordinated Spiking Activity During NREM Sleep.

The dialogue between cortex and hippocampus is known to be crucial for sleep-dependent memory consolidation. During slow wave sleep, memory replay depends on slow oscillation (SO) and spindles in the (neo)cortex and sharp wave-ripples (SWRs) in the hippocampus. The mechanisms underlying interaction of these rhythms are poorly understood. We examined the interaction between cortical SO and hippocampal SWRs in a model of the hippocampo-cortico-thalamic network and compared the results with human intracranial recordings during sleep. We observed that ripple occurrence peaked following the onset of an Up-state of SO and that cortical input to hippocampus was crucial to maintain this relationship. A small fraction of ripples occurred during the Down-state and controlled initiation of the next Up-state. We observed that the effect of ripple depends on its precise timing, which supports the idea that ripples occurring at different phases of SO might serve different functions, particularly in the context of encoding the new and reactivation of the old memories during memory consolidation. The study revealed complex bidirectional interaction of SWRs and SO in which early hippocampal ripples influence transitions to Up-state, while cortical Up-states control occurrence of the later ripples, which in turn influence transition to Down-state.

Multi-layer network utilizing rewarded spike time dependent plasticity to learn a foraging task.

Neural networks with a single plastic layer employing reward modulated spike time dependent plasticity (STDP) are capable of learning simple foraging tasks. Here we demonstrate advanced pattern discrimination and continuous learning in a network of spiking neurons with multiple plastic layers. The network utilized both reward modulated and non-reward modulated STDP and implemented multiple mechanisms for homeostatic regulation of synaptic efficacy, including heterosynaptic plasticity, gain control, output balancing, activity normalization of rewarded STDP and hard limits on synaptic strength. We found that addition of a hidden layer of neurons employing non-rewarded STDP created neurons that responded to the specific combinations of inputs and thus performed basic classification of the input patterns. When combined with a following layer of neurons implementing rewarded STDP, the network was able to learn, despite the absence of labeled training data, discrimination between rewarding patterns and the patterns designated as punishing. Synaptic noise allowed for trial-and-error learning that helped to identify the goal-oriented strategies which were effective in task solving. The study predicts a critical set of properties of the spiking neuronal network with STDP that was sufficient to solve a complex foraging task involving pattern classification and decision making.

Classification of odorants across layers in locust olfactory pathway.

Olfactory processing takes place across multiple layers of neurons from the transduction of odorants in the periphery, to odor quality processing, learning, and decision making in higher olfactory structures. In insects, projection neurons (PNs) in the antennal lobe send odor information to the Kenyon cells (KCs) of the mushroom bodies and lateral horn neurons (LHNs). To examine the odor information content in different structures of the insect brain, antennal lobe, mushroom bodies and lateral horn, we designed a model of the olfactory network based on electrophysiological recordings made in vivo in the locust. We found that populations of all types (PNs, LHNs, and KCs) had lower odor classification error rates than individual cells of any given type. This improvement was quantitatively different from that observed using uniform populations of identical neurons compared with spatially structured population of neurons tuned to different odor features. This result, therefore, reflects an emergent network property. Odor classification improved with increasing stimulus duration: for similar odorants, KC and LHN ensembles reached optimal discrimination within the first 300-500 ms of the odor response. Performance improvement with time was much greater for a population of cells than for individual neurons. We conclude that, for PNs, LHNs, and KCs, ensemble responses are always much more informative than single-cell responses, despite the accumulation of noise along with odor information.

Feed-Forward versus Feedback Inhibition in a Basic Olfactory Circuit.

Inhibitory interneurons play critical roles in shaping the firing patterns of principal neurons in many brain systems. Despite difference in the anatomy or functions of neuronal circuits containing inhibition, two basic motifs repeatedly emerge: feed-forward and feedback. In the locust, it was proposed that a subset of lateral horn interneurons (LHNs), provide feed-forward inhibition onto Kenyon cells (KCs) to maintain their sparse firing--a property critical for olfactory learning and memory. But recently it was established that a single inhibitory cell, the giant GABAergic neuron (GGN), is the main and perhaps sole source of inhibition in the mushroom body, and that inhibition from this cell is mediated by a feedback (FB) loop including KCs and the GGN. To clarify basic differences in the effects of feedback vs. feed-forward inhibition in circuit dynamics we here use a model of the locust olfactory system. We found both inhibitory motifs were able to maintain sparse KCs responses and provide optimal odor discrimination. However, we further found that only FB inhibition could create a phase response consistent with data recorded in vivo. These findings describe general rules for feed-forward versus feedback inhibition and suggest GGN is potentially capable of providing the primary source of inhibition to the KCs. A better understanding of how inhibitory motifs impact post-synaptic neuronal activity could be used to reveal unknown inhibitory structures within biological networks.

Python/NEURON code for simulating biophysically realistic thalamocortical dynamics during sleep.

Understanding the function of sleep and its associated neural rhythms is an important goal in neuroscience. While many theoretical models of neural dynamics during sleep exist, few include the effects of neuromodulators on sleep oscillations and describe transitions between sleep and wake states or different sleep stages. Here, we started with a C++-based thalamocortical network model that describes characteristic thalamic and cortical oscillations specific to sleep. This model, which includes a biophysically realistic description of intrinsic and synaptic channels, allows for testing the effects of different neuromodulators, intrinsic cell properties, and synaptic connectivity on neural dynamics during sleep. We present a complete reimplementation of this previously-published sleep model in the standardized NEURON/Python framework, making it more accessible to the wider scientific community.

Tackling the challenges of group network inference from intracranial EEG data.

Introduction: Intracranial EEG (iEEG) data is a powerful way to map brain function, characterized by high temporal and spatial resolution, allowing the study of interactions among neuronal populations that orchestrate cognitive processing. However, the statistical inference and analysis of brain networks using iEEG data faces many challenges related to its sparse brain coverage, and its inhomogeneity across patients. Methods: We review these challenges and develop a methodological pipeline for estimation of network structure not obtainable from any single patient, illustrated on the inference of the interaction among visual streams using a dataset of 27 human iEEG recordings from a visual experiment employing visual scene stimuli. 100 ms sliding window and multiple band-pass filtered signals are used to provide temporal and spectral resolution. For the connectivity analysis we showcase two connectivity measures reflecting different types of interaction between regions of interest (ROI): Phase Locking Value as a symmetric measure of synchrony, and Directed Transfer Function-asymmetric measure describing causal interaction. For each two channels, initial uncorrected significance testing at p < 0.05 for every time-frequency point is carried out by comparison of the data-derived connectivity to a baseline surrogate-based null distribution, providing a binary time-frequency connectivity map. For each ROI pair, a connectivity density map is obtained by averaging across all pairs of channels spanning them, effectively agglomerating data across relevant channels and subjects. Finally, the difference of the mean map value after and before the stimulation is compared to the same statistic in surrogate data to assess link significance. Results: The analysis confirmed the function of the parieto-medial temporal pathway, mediating visuospatial information between dorsal and ventral visual streams during visual scene analysis. Moreover, we observed the anterior hippocampal connectivity with more posterior areas in the medial temporal lobe, and found the reciprocal information flow between early processing areas and medial place area. Discussion: To summarize, we developed an approach for estimating network connectivity, dealing with the challenge of sparse individual coverage of intracranial EEG electrodes. Its application provided new insights into the interaction between the dorsal and ventral visual streams, one of the iconic dualities in human cognition.

Human brain structural connectivity matrices-ready for modelling.

The human brain represents a complex computational system, the function and structure of which may be measured using various neuroimaging techniques focusing on separate properties of the brain tissue and activity. We capture the organization of white matter fibers acquired by diffusion-weighted imaging using probabilistic diffusion tractography. By segmenting the results of tractography into larger anatomical units, it is possible to draw inferences about the structural relationships between these parts of the system. This pipeline results in a structural connectivity matrix, which contains an estimate of connection strength among all regions. However, raw data processing is complex, computationally intensive, and requires expert quality control, which may be discouraging for researchers with less experience in the field. We thus provide brain structural connectivity matrices in a form ready for modelling and analysis and thus usable by a wide community of scientists. The presented dataset contains brain structural connectivity matrices together with the underlying raw diffusion and structural data, as well as basic demographic data of 88 healthy subjects.

Responses of single neurons and neuronal ensembles in frog first- and second-order olfactory neurons.

A major challenge in sensory neuroscience is to elucidate the coding and processing of stimulus representations in successive populations of neurons. Here we recorded the spiking activity of receptor neurons (RNs) and mitral/tufted cells (MCs) in the frog olfactory epithelium and olfactory bulb respectively, in response to four odorants applied at precisely controlled concentrations. We compared how RN responses are translated in MCs. We examined the time course of the instantaneous firing frequency before and after stimulation in neuron ensembles and the dependency on odorant concentration of the number of action potentials fired in a preselected 5-s time window (dose-response curves) in both single neurons and neuron ensembles. In RNs and MCs, the dose-response curves typically increase then decrease and are well described by alpha functions. We established the main quantitative properties of these curves, including the distributions of concentrations at threshold and maximum responses. We showed that the main transformations occurring in the transition from RNs to MCs is the lowering of the firing threshold and a large decrease in the total number of spikes fired. We also found that the number of action potentials fired by recorded neurons and hence their energy consumption is independent of odorant concentration, and that this is a consequence of their time- and concentration-dependent activities. This article is part of a Special Issue entitled Neural Coding 2012.

Stochastic interpolation model of the medial superior olive neural circuit.

This article presents a stochastic model of binaural hearing in the medial superior olive (MSO) circuit. This model is a variant of the slope encoding models. First, a general framework is developed describing the elementary neural operations realized on spike trains in individual parts of the circuit and how the neurons converging onto the MSO are connected. Random delay, coincidence detection of spikes, divergence and convergence of spike trains are operations implemented by the following modules: spike generator, jitter generator, and coincidence detector. Subsequent processing of spike trains computes the sound azimuth in the circuit. The circuit parameters that influence efficiency of slope encoding are studied. In order to measure the overall circuit performance the concept of an ideal observer is used instead of a detailed model of higher relays in the auditory pathway. This makes it possible to bridge the gap between psychophysical observations in humans and recordings taken of small rodents. Most of the results are obtained through numerical simulations of the model.

Effect of stimulation on the input parameters of stochastic leaky integrate-and-fire neuronal model.

The Ornstein-Uhlenbeck neuronal model is specified by two types of parameters. One type corresponds to the properties of the neuronal membrane, whereas the second type (local average rate of the membrane depolarization and its variability) corresponds to the input of the neuron. In this article, we estimate the parameters of the second type from an intracellular record during neuronal firing caused by stimulation (audio signal). We compare the obtained estimates with those from the spontaneous part of the record. As predicted from the model construction, the values of the input parameters are larger for the periods when neuron is stimulated than for the spontaneous ones. Finally, the firing regimen of the model is checked. It is confirmed that the neuron is in the suprathreshold regimen during the stimulation.

The parameters of the stochastic leaky integrate-and-fire neuronal model.

Five parameters of one of the most common neuronal models, the diffusion leaky integrate-and-fire model, also known as the Ornstein-Uhlenbeck neuronal model, were estimated on the basis of intracellular recording. These parameters can be classified into two categories. Three of them (the membrane time constant, the resting potential and the firing threshold) characterize the neuron itself. The remaining two characterize the neuronal input. The intracellular data were collected during spontaneous firing, which in this case is characterized by a Poisson process of interspike intervals. Two methods for the estimation were applied, the regression method and the maximum-likelihood method. Both methods permit to estimate the input parameters and the membrane time constant in a short time window (a single interspike interval). We found that, at least in our example, the regression method gave more consistent results than the maximum-likelihood method. The estimates of the input parameters show the asymptotical normality, which can be further used for statistical testing, under the condition that the data are collected in different experimental situations. The model neuron, as deduced from the determined parameters, works in a subthreshold regimen. This result was confirmed by both applied methods. The subthreshold regimen for this model is characterized by the Poissonian firing. This is in a complete agreement with the observed interspike interval data.