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.

Emergent effects of synaptic connectivity on the dynamics of global and local slow waves in a large-scale thalamocortical network model of the human brain.

Slow-wave sleep (SWS), characterized by slow oscillations (SOs, <1Hz) of alternating active and silent states in the thalamocortical network, is a primary brain state during Non-Rapid Eye Movement (NREM) sleep. In the last two decades, the traditional view of SWS as a global and uniform whole-brain state has been challenged by a growing body of evidence indicating that SO can be local and can coexist with wake-like activity. However, the mechanisms by which global and local SOs arise from micro-scale neuronal dynamics and network connectivity remain poorly understood. We developed a multi-scale, biophysically realistic human whole-brain thalamocortical network model capable of transitioning between the awake state and SWS, and we investigated the role of connectivity in the spatio-temporal dynamics of sleep SO. We found that the overall strength and a relative balance between long and short-range synaptic connections determined the network state. Importantly, for a range of synaptic strengths, the model demonstrated complex mixed SO states, where periods of synchronized global slow-wave activity were intermittent with the periods of asynchronous local slow-waves. An increase in the overall synaptic strength led to synchronized global SO, while a decrease in synaptic connectivity produced only local slow-waves that would not propagate beyond local areas. These results were compared to human data to validate probable models of biophysically realistic SO. The model producing mixed states provided the best match to the spatial coherence profile and the functional connectivity estimated from human subjects. These findings shed light on how the spatio-temporal properties of SO emerge from local and global cortical connectivity and provide a framework for further exploring the mechanisms and functions of SWS in health and disease.

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.

Stimulation Augments Spike Sequence Replay and Memory Consolidation during Slow-Wave Sleep.

Newly acquired memory traces are spontaneously reactivated during slow-wave sleep (SWS), leading to the consolidation of recent memories. Empirical studies found that sensory stimulation during SWS can selectively enhance memory consolidation with the effect depending on the phase of stimulation. In this new study, we aimed to understand the mechanisms behind the role of sensory stimulation on memory consolidation using computational models implementing effects of neuromodulators to simulate transitions between awake and SWS sleep, and synaptic plasticity to allow the change of synaptic connections due to the training in awake or replay during sleep. We found that when closed-loop stimulation was applied during the Down states of sleep slow oscillation, particularly right before the transition from Down to Up state, it significantly affected the spatiotemporal pattern of the slow waves and maximized memory replay. In contrast, when the stimulation was presented during the Up states, it did not have a significant impact on the slow waves or memory performance after sleep. For multiple memories trained in awake, presenting stimulation cues associated with specific memory trace could selectively augment replay and enhance consolidation of that memory and interfere with consolidation of the others (particularly weak) memories. Our study proposes a synaptic-level mechanism of how memory consolidation is affected by sensory stimulation during sleep.SIGNIFICANCE STATEMENT Stimulation, such as training-associated cues or auditory stimulation, during sleep can augment consolidation of the newly encoded memories. In this study, we used a computational model of the thalamocortical system to describe the mechanisms behind the role of stimulation in memory consolidation during slow-wave sleep. Our study suggests that stimulation preferentially strengthens memory traces when delivered at a specific phase of the slow oscillation, just before the Down to Up state transition when it makes the largest impact on the spatiotemporal pattern of sleep slow waves. In the presence of multiple memories, presenting sensory cues during sleep could selectively strengthen selected memories. Our study proposes a synaptic-level mechanism of how memory consolidation is affected by sensory stimulation during sleep.

Replay in Deep Learning: Current Approaches and Missing Biological Elements.

Replay is the reactivation of one or more neural patterns that are similar to the activation patterns experienced during past waking experiences. Replay was first observed in biological neural networks during sleep, and it is now thought to play a critical role in memory formation, retrieval, and consolidation. Replay-like mechanisms have been incorporated in deep artificial neural networks that learn over time to avoid catastrophic forgetting of previous knowledge. Replay algorithms have been successfully used in a wide range of deep learning methods within supervised, unsupervised, and reinforcement learning paradigms. In this letter, we provide the first comprehensive comparison between replay in the mammalian brain and replay in artificial neural networks. We identify multiple aspects of biological replay that are missing in deep learning systems and hypothesize how they could be used to improve artificial neural networks.

Origin of slow spontaneous resting-state neuronal fluctuations in brain networks.

Resting- or baseline-state low-frequency (0.01-0.2 Hz) brain activity is observed in fMRI, EEG, and local field potential recordings. These fluctuations were found to be correlated across brain regions and are thought to reflect neuronal activity fluctuations between functionally connected areas of the brain. However, the origin of these infra-slow resting-state fluctuations remains unknown. Here, using a detailed computational model of the brain network, we show that spontaneous infra-slow (<0.05 Hz) activity could originate due to the ion concentration dynamics. The computational model implemented dynamics for intra- and extracellular K+ and Na+ and intracellular Cl- ions, Na+/K+ exchange pump, and KCC2 cotransporter. In the network model simulating resting awake-like brain state, we observed infra-slow fluctuations in the extracellular K+ concentration, Na+/K+ pump activation, firing rate of neurons, and local field potentials. Holding K+ concentration constant prevented generation of the infra-slow fluctuations. The amplitude and peak frequency of this activity were modulated by the Na+/K+ pump, AMPA/GABA synaptic currents, and glial properties. Further, in a large-scale network with long-range connections based on CoCoMac connectivity data, the infra-slow fluctuations became synchronized among remote clusters similar to the resting-state activity observed in vivo. Overall, our study proposes that ion concentration dynamics mediated by neuronal and glial activity may contribute to the generation of very slow spontaneous fluctuations of brain activity that are reported as the resting-state fluctuations in fMRI and EEG recordings.

Differential Covariance: A New Method to Estimate Functional Connectivity in fMRI.

Measuring functional connectivity from fMRI recordings is important in understanding processing in cortical networks. However, because the brain's connection pattern is complex, currently used methods are prone to producing false functional connections. We introduce differential covariance analysis, a new method that uses derivatives of the signal for estimating functional connectivity. We generated neural activities from dynamical causal modeling and a neural network of Hodgkin-Huxley neurons and then converted them to hemodynamic signals using the forward balloon model. The simulated fMRI signals, together with the ground-truth connectivity pattern, were used to benchmark our method with other commonly used methods. Differential covariance achieved better results in complex network simulations. This new method opens an alternative way to estimate functional connectivity.

Timing between Cortical Slow Oscillations and Heart Rate Bursts during Sleep Predicts Temporal Processing Speed, but Not Offline Consolidation.

Central and autonomic nervous system activities are coupled during sleep. Cortical slow oscillations (SOs; <1 Hz) coincide with brief bursts in heart rate (HR), but the functional consequence of this coupling in cognition remains elusive. We measured SO-HR temporal coupling (i.e., the peak-to-peak interval between downstate of SO event and HR burst) during a daytime nap and asked whether this SO-HR timing measure was associated with temporal processing speed and learning on a texture discrimination task by testing participants before and after a nap. The coherence of SO-HR events during sleep strongly correlated with an individual's temporal processing speed in the morning and evening test sessions, but not with their change in performance after the nap (i.e., consolidation). We confirmed this result in two additional experimental visits and also discovered that this association was visit-specific, indicating a state (not trait) marker. Thus, we introduce a novel physiological index that may be a useful marker of state-dependent processing speed of an individual.

Ionic and synaptic mechanisms of seizure generation and epileptogenesis.

The biophysical mechanisms underlying epileptogenesis and the generation of seizures remain to be better understood. Among many factors triggering epileptogenesis are traumatic brain injury breaking normal synaptic homeostasis and genetic mutations disrupting ionic concentration homeostasis. Impairments in these mechanisms, as seen in various brain diseases, may push the brain network to a pathological state characterized by increased susceptibility to unprovoked seizures. Here, we review recent computational studies exploring the roles of ionic concentration dynamics in the generation, maintenance, and termination of seizures. We further discuss how ionic and synaptic homeostatic mechanisms may give rise to conditions which prime brain networks to exhibit recurrent spontaneous seizures and epilepsy.

Coupling of autonomic and central events during sleep benefits declarative memory consolidation.

While anatomical pathways between forebrain cognitive and brainstem autonomic nervous centers are well-defined, autonomic-central interactions during sleep and their contribution to waking performance are not understood. Here, we analyzed simultaneous central activity via electroencephalography (EEG) and autonomic heart beat-to-beat intervals (RR intervals) from electrocardiography (ECG) during wake and daytime sleep. We identified bursts of ECG activity that lasted 4-5 s and predominated in non-rapid-eye-movement sleep (NREM). Using event-based analysis of NREM sleep, we found an increase in delta (0.5-4 Hz) and sigma (12-15 Hz) power and an elevated density of slow oscillations (0.5-1 Hz) about 5 s prior to peak of the heart rate burst, as well as a surge in vagal activity, assessed by high-frequency (HF) component of RR intervals. Using regression framework, we show that these Autonomic/Central Events (ACE) positively predicted post-nap improvement in a declarative memory task after controlling for the effects of spindles and slow oscillations from sleep periods without ACE. No such relation was found between memory performance and a control nap. Additionally, NREM ACE negatively correlated with REM sleep and learning in a non-declarative memory task. These results provide the first evidence that coordinated autonomic and central events play a significant role in declarative memory consolidation.

Differential roles of sleep spindles and sleep slow oscillations in memory consolidation.

Sleep plays an important role in the consolidation of recent memories. However, the cellular and synaptic mechanisms of consolidation during sleep remain poorly understood. In this study, using a realistic computational model of the thalamocortical network, we tested the role of Non-Rapid Eye Movement (NREM) sleep in memory consolidation. We found that sleep spindles (the hallmark of N2 stage sleep) and slow oscillations (the hallmark of N3 stage sleep) both promote replay of the spike sequences learned in the awake state and replay was localized at the trained network locations. Memory performance improved after a period of NREM sleep but not after the same time period in awake. When multiple memories were trained, the local nature of the spike sequence replay during spindles allowed replay of the distinct memory traces independently, while slow oscillations promoted competition that could prevent replay of the weak memories in a presence of the stronger memory traces. This could lead to extinction of the weak memories unless when sleep spindles (N2 sleep) preceded slow oscillations (N3 sleep), as observed during the natural sleep cycle. Our study presents a mechanistic explanation for the role of sleep rhythms in memory consolidation and proposes a testable hypothesis how the natural structure of sleep stages provides an optimal environment to consolidate memories.

Thalamocortical and intracortical laminar connectivity determines sleep spindle properties.

Sleep spindles are brief oscillatory events during non-rapid eye movement (NREM) sleep. Spindle density and synchronization properties are different in MEG versus EEG recordings in humans and also vary with learning performance, suggesting spindle involvement in memory consolidation. Here, using computational models, we identified network mechanisms that may explain differences in spindle properties across cortical structures. First, we report that differences in spindle occurrence between MEG and EEG data may arise from the contrasting properties of the core and matrix thalamocortical systems. The matrix system, projecting superficially, has wider thalamocortical fanout compared to the core system, which projects to middle layers, and requires the recruitment of a larger population of neurons to initiate a spindle. This property was sufficient to explain lower spindle density and higher spatial synchrony of spindles in the superficial cortical layers, as observed in the EEG signal. In contrast, spindles in the core system occurred more frequently but less synchronously, as observed in the MEG recordings. Furthermore, consistent with human recordings, in the model, spindles occurred independently in the core system but the matrix system spindles commonly co-occurred with core spindles. We also found that the intracortical excitatory connections from layer III/IV to layer V promote spindle propagation from the core to the matrix system, leading to widespread spindle activity. Our study predicts that plasticity of intra- and inter-cortical connectivity can potentially be a mechanism for increased spindle density as has been observed during learning.

Computational model of brain-stem circuit for state-dependent control of hypoglossal motoneurons.

In patients with obstructive sleep apnea (OSA), the pharyngeal muscles become relaxed during sleep, which leads to a partial or complete closure of upper airway. Experimental studies suggest that withdrawal of noradrenergic and serotonergic drives importantly contributes to depression of hypoglossal motoneurons and, therefore, may contribute to OSA pathophysiology; however, specific cellular and synaptic mechanisms remain unknown. In this new study, we developed a biophysical network model to test the hypothesis that, to explain experimental observations, the neuronal network for monoaminergic control of excitability of hypoglossal motoneurons needs to include excitatory and inhibitory perihypoglossal interneurons that mediate noradrenergic and serotonergic drives to hypoglossal motoneurons. In the model, the state-dependent activation of the hypoglossal motoneurons was in qualitative agreement with in vivo data during simulated rapid eye movement (REM) and non-REM sleep. The model was applied to test the mechanisms of action of noradrenergic and serotonergic drugs during REM sleep as observed in vivo. We conclude that the proposed minimal neuronal circuit is sufficient to explain in vivo data and supports the hypothesis that perihypoglossal interneurons may mediate state-dependent monoaminergic drive to hypoglossal motoneurons. The population of the hypothesized perihypoglossal interneurons may serve as novel targets for pharmacological treatment of OSA. NEW & NOTEWORTHY In vivo studies suggest that during rapid eye movement sleep, withdrawal of noradrenergic and serotonergic drives critically contributes to depression of hypoglossal motoneurons (HMs), which innervate the tongue muscles. By means of a biophysical model, which is consistent with a broad range of empirical data, we demonstrate that the neuronal network controlling the excitability of HMs needs to include excitatory and inhibitory interneurons that mediate noradrenergic and serotonergic drives to HMs.

Role of KCC2-dependent potassium efflux in 4-Aminopyridine-induced Epileptiform synchronization.

A balance between excitation and inhibition is necessary to maintain stable brain network dynamics. Traditionally, seizure activity is believed to arise from the breakdown of this delicate balance in favor of excitation with loss of inhibition. Surprisingly, recent experimental evidence suggests that this conventional view may be limited, and that inhibition plays a prominent role in the development of epileptiform synchronization. Here, we explored the role of the KCC2 co-transporter in the onset of inhibitory network-induced seizures. Our experiments in acute mouse brain slices, of either sex, revealed that optogenetic stimulation of either parvalbumin- or somatostatin-expressing interneurons induced ictal discharges in rodent entorhinal cortex during 4-aminopyridine application. These data point to a proconvulsive role of GABAA receptor signaling that is independent of the inhibitory input location (i.e., dendritic vs. somatic). We developed a biophysically realistic network model implementing dynamics of ion concentrations to explore the mechanisms leading to inhibitory network-induced seizures. In agreement with experimental results, we found that stimulation of the inhibitory interneurons induced seizure-like activity in a network with reduced potassium A-current. Our model predicts that interneuron stimulation triggered an increase of interneuron firing, which was accompanied by an increase in the intracellular chloride concentration and a subsequent KCC2-dependent gradual accumulation of the extracellular potassium promoting epileptiform ictal activity. When the KCC2 activity was reduced, stimulation of the interneurons was no longer able to induce ictal events. Overall, our study provides evidence for a proconvulsive role of GABAA receptor signaling that depends on the involvement of the KCC2 co-transporter.

Synaptic Mechanisms of Memory Consolidation during Sleep Slow Oscillations.

Sleep is critical for regulation of synaptic efficacy, memories, and learning. However, the underlying mechanisms of how sleep rhythms contribute to consolidating memories acquired during wakefulness remain unclear. Here we studied the role of slow oscillations, 0.2-1 Hz rhythmic transitions between Up and Down states during stage 3/4 sleep, on dynamics of synaptic connectivity in the thalamocortical network model implementing spike-timing-dependent synaptic plasticity. We found that the spatiotemporal pattern of Up-state propagation determines the changes of synaptic strengths between neurons. Furthermore, an external input, mimicking hippocampal ripples, delivered to the cortical network results in input-specific changes of synaptic weights, which persisted after stimulation was removed. These synaptic changes promoted replay of specific firing sequences of the cortical neurons. Our study proposes a neuronal mechanism on how an interaction between hippocampal input, such as mediated by sharp wave-ripple events, cortical slow oscillations, and synaptic plasticity, may lead to consolidation of memories through preferential replay of cortical cell spike sequences during slow-wave sleep. SIGNIFICANCE STATEMENT: Sleep is critical for memory and learning. Replay during sleep of temporally ordered spike sequences related to a recent experience was proposed to be a neuronal substrate of memory consolidation. However, specific mechanisms of replay or how spike sequence replay leads to synaptic changes that underlie memory consolidation are still poorly understood. Here we used a detailed computational model of the thalamocortical system to report that interaction between slow cortical oscillations and synaptic plasticity during deep sleep can underlie mapping hippocampal memory traces to persistent cortical representation. This study provided, for the first time, a mechanistic explanation of how slow-wave sleep may promote consolidation of recent memory events.

Differential Covariance: A New Class of Methods to Estimate Sparse Connectivity from Neural Recordings.

With our ability to record more neurons simultaneously, making sense of these data is a challenge. Functional connectivity is one popular way to study the relationship of multiple neural signals. Correlation-based methods are a set of currently well-used techniques for functional connectivity estimation. However, due to explaining away and unobserved common inputs (Stevenson, Rebesco, Miller, & Körding, 2008 ), they produce spurious connections. The general linear model (GLM), which models spike trains as Poisson processes (Okatan, Wilson, & Brown, 2005 ; Truccolo, Eden, Fellows, Donoghue, & Brown, 2005 ; Pillow et al., 2008 ), avoids these confounds. We develop here a new class of methods by using differential signals based on simulated intracellular voltage recordings. It is equivalent to a regularized AR(2) model. We also expand the method to simulated local field potential recordings and calcium imaging. In all of our simulated data, the differential covariance-based methods achieved performance better than or similar to the GLM method and required fewer data samples. This new class of methods provides alternative ways to analyze neural signals.

Hippocampal CA1 Ripples as Inhibitory Transients.

Memories are stored and consolidated as a result of a dialogue between the hippocampus and cortex during sleep. Neurons active during behavior reactivate in both structures during sleep, in conjunction with characteristic brain oscillations that may form the neural substrate of memory consolidation. In the hippocampus, replay occurs within sharp wave-ripples: short bouts of high-frequency activity in area CA1 caused by excitatory activation from area CA3. In this work, we develop a computational model of ripple generation, motivated by in vivo rat data showing that ripples have a broad frequency distribution, exponential inter-arrival times and yet highly non-variable durations. Our study predicts that ripples are not persistent oscillations but result from a transient network behavior, induced by input from CA3, in which the high frequency synchronous firing of perisomatic interneurons does not depend on the time scale of synaptic inhibition. We found that noise-induced loss of synchrony among CA1 interneurons dynamically constrains individual ripple duration. Our study proposes a novel mechanism of hippocampal ripple generation consistent with a broad range of experimental data, and highlights the role of noise in regulating the duration of input-driven oscillatory spiking in an inhibitory network.

Modeling of Age-Dependent Epileptogenesis by Differential Homeostatic Synaptic Scaling.

Homeostatic synaptic plasticity (HSP) has been implicated in the development of hyperexcitability and epileptic seizures following traumatic brain injury (TBI). Our in vivo experimental studies in cats revealed that the severity of TBI-mediated epileptogenesis depends on the age of the animal. To characterize mechanisms of these differences, we studied the properties of the TBI-induced epileptogenesis in a biophysically realistic cortical network model with dynamic ion concentrations. After deafferentation, which was induced by dissection of the afferent inputs, there was a reduction of the network activity and upregulation of excitatory connections leading to spontaneous spike-and-wave type seizures. When axonal sprouting was implemented, the seizure threshold increased in the model of young but not the older animals, which had slower or unidirectional homeostatic processes. Our study suggests that age-related changes in the HSP mechanisms are sufficient to explain the difference in the likelihood of seizure onset in young versus older animals. Significance statement: Traumatic brain injury (TBI) is one of the leading causes of intractable epilepsy. Likelihood of developing epilepsy and seizures following severe brain trauma has been shown to increase with age. Specific mechanisms of TBI-related epileptogenesis and how these mechanisms are affected by age remain to be understood. We test a hypothesis that the failure of homeostatic synaptic regulation, a slow negative feedback mechanism that maintains neural activity within a physiological range through activity-dependent modulation of synaptic strength, in older animals may augment TBI-induced epileptogenesis. Our results provide new insight into understanding this debilitating disorder and may lead to novel avenues for the development of effective treatments of TBI-induced epilepsy.

Electrogenic properties of the Na⁺/K⁺ ATPase control transitions between normal and pathological brain states.

Ionic concentrations fluctuate significantly during epileptic seizures. In this study, using a combination of in vitro electrophysiology, computer modeling, and dynamical systems analysis, we demonstrate that changes in the potassium and sodium intra- and extracellular ion concentrations ([K(+)] and [Na(+)], respectively) during seizure affect the neuron dynamics by modulating the outward Na(+)/K(+) pump current. First, we show that an increase of the outward Na(+)/K(+) pump current mediates termination of seizures when there is a progressive increase in the intracellular [Na(+)]. Second, we show that the Na(+)/K(+) pump current is crucial in maintaining the stability of the physiological network state; a reduction of this current leads to the onset of seizures via a positive-feedback loop. We then present a novel dynamical mechanism for bursting in neurons with a reduced Na(+)/K(+) pump. Overall, our study demonstrates the profound role of the current mediated by Na(+)/K(+) ATPase on the stability of neuronal dynamics that was previously unknown.

Emergent effects of synaptic connectivity on the dynamics of global and local slow waves in a large-scale thalamocortical network model of the human brain.

Slow-wave sleep (SWS), characterized by slow oscillations (SO, <1Hz) of alternating active and silent states in the thalamocortical network, is a primary brain state during Non-Rapid Eye Movement (NREM) sleep. In the last two decades, the traditional view of SWS as a global and uniform whole-brain state has been challenged by a growing body of evidence indicating that SO can be local and can coexist with wake-like activity. However, the understanding of how global and local SO emerges from micro-scale neuron dynamics and network connectivity remains unclear. We developed a multi-scale, biophysically realistic human whole-brain thalamocortical network model capable of transitioning between the awake state and slow-wave sleep, and we investigated the role of connectivity in the spatio-temporal dynamics of sleep SO. We found that the overall strength and a relative balance between long and short-range synaptic connections determined the network state. Importantly, for a range of synaptic strengths, the model demonstrated complex mixed SO states, where periods of synchronized global slow-wave activity were intermittent with the periods of asynchronous local slow-waves. Increase of the overall synaptic strength led to synchronized global SO, while decrease of synaptic connectivity produced only local slow-waves that would not propagate beyond local area. These results were compared to human data to validate probable models of biophysically realistic SO. The model producing mixed states provided the best match to the spatial coherence profile and the functional connectivity estimated from human subjects. These findings shed light on how the spatio-temporal properties of SO emerge from local and global cortical connectivity and provide a framework for further exploring the mechanisms and functions of SWS in health and disease.