Modeling Sleep and Wakefulness in the Thalamocortical System

Hypothalamic feedforward inhibition of thalamocortical network controls arousal and consciousness

During non-rapid eye movement (NREM) sleep, synchronous synaptic activity in the thalamocortical network generates predominantly low-frequency oscillations (<4 Hz) that are modulated by inhibitory inputs from the thalamic reticular nucleus (TRN). Whether TRN cells integrate sleep-wake signals from subcortical circuits remains unclear. We found that GABA neurons from the lateral hypothalamus (LHGABA) exert a strong inhibitory control over TRN GABA neurons (TRNGABA). We found that optogenetic activation of this circuit recapitulated state-dependent changes of TRN neuron activity in behaving mice and induced rapid arousal during NREM, but not REM, sleep. During deep anesthesia, activation of this circuit induced sustained cortical arousal. In contrast, optogenetic silencing of LHGABA-TRNGABA transmission increased the duration of NREM sleep and amplitude of delta (1-4 Hz) oscillations. Collectively, these results demonstrate that TRN cells integrate subcortical arousal inputs selectively during NREM sleep and may participate in sleep intensity.

Computing the Local Field Potential (LFP) from Integrate-and-Fire Network Models

Leaky integrate-and-fire (LIF) network models are commonly used to study how the spiking dynamics of neural networks changes with stimuli, tasks or dynamic network states. However, neurophysiological studies in vivo often rather measure the mass activity of neuronal microcircuits with the local field potential (LFP). Given that LFPs are generated by spatially separated currents across the neuronal membrane, they cannot be computed directly from quantities defined in models of point-like LIF neurons. Here, we explore the best approximation for predicting the LFP based on standard output from point-neuron LIF networks. To search for this best "LFP proxy", we compared LFP predictions from candidate proxies based on LIF network output (e.g, firing rates, membrane potentials, synaptic currents) with "ground-truth" LFP obtained when the LIF network synaptic input currents were injected into an analogous three-dimensional (3D) network model of multi-compartmental neurons with realistic morphology, spatial distributions of somata and synapses. We found that a specific fixed linear combination of the LIF synaptic currents provided an accurate LFP proxy, accounting for most of the variance of the LFP time course observed in the 3D network for all recording locations. This proxy performed well over a broad set of conditions, including substantial variations of the neuronal morphologies. Our results provide a simple formula for estimating the time course of the LFP from LIF network simulations in cases where a single pyramidal population dominates the LFP generation, and thereby facilitate quantitative comparison between computational models and experimental LFP recordings in vivo.

Plasticity of brain wave network interactions and evolution across physiologic states

Neural plasticity transcends a range of spatio-temporal scales and serves as the basis of various brain activities and physiologic functions. At the microscopic level, it enables the emergence of brain waves with complex temporal dynamics. At the macroscopic level, presence and dominance of specific brain waves is associated with important brain functions. The role of neural plasticity at different levels in generating distinct brain rhythms and how brain rhythms communicate with each other across brain areas to generate physiologic states and functions remains not understood. Here we perform an empirical exploration of neural plasticity at the level of brain wave network interactions representing dynamical communications within and between different brain areas in the frequency domain. We introduce the concept of time delay stability (TDS) to quantify coordinated bursts in the activity of brain waves, and we employ a system-wide Network Physiology integrative approach to probe the network of coordinated brain wave activations and its evolution across physiologic states. We find an association between network structure and physiologic states. We uncover a hierarchical reorganization in the brain wave networks in response to changes in physiologic state, indicating new aspects of neural plasticity at the integrated level. Globally, we find that the entire brain network undergoes a pronounced transition from low connectivity in Deep Sleep and REM to high connectivity in Light Sleep and Wake. In contrast, we find that locally, different brain areas exhibit different network dynamics of brain wave interactions to achieve differentiation in function during different sleep stages. Moreover, our analyses indicate that plasticity also emerges in frequency-specific networks, which represent interactions across brain locations mediated through a specific frequency band. Comparing frequency-specific networks within the same physiologic state we find very different degree of network connectivity and link strength, while at the same time each frequency-specific network is characterized by a different signature pattern of sleep-stage stratification, reflecting a remarkable flexibility in response to change in physiologic state. These new aspects of neural plasticity demonstrate that in addition to dominant brain waves, the network of brain wave interactions is a previously unrecognized hallmark of physiologic state and function.

Shaping Neuronal Network Activity by Presynaptic Mechanisms

Neuronal microcircuits generate oscillatory activity, which has been linked to basic functions such as sleep, learning and sensorimotor gating. Although synaptic release processes are well known for their ability to shape the interaction between neurons in microcircuits, most computational models do not simulate the synaptic transmission process directly and hence cannot explain how changes in synaptic parameters alter neuronal network activity. In this paper, we present a novel neuronal network model that incorporates presynaptic release mechanisms, such as vesicle pool dynamics and calcium-dependent release probability, to model the spontaneous activity of neuronal networks. The model, which is based on modified leaky integrate-and-fire neurons, generates spontaneous network activity patterns, which are similar to experimental data and robust under changes in the model's primary gain parameters such as excitatory postsynaptic potential and connectivity ratio. Furthermore, it reliably recreates experimental findings and provides mechanistic explanations for data obtained from microelectrode array recordings, such as network burst termination and the effects of pharmacological and genetic manipulations. The model demonstrates how elevated asynchronous release, but not spontaneous release, synchronizes neuronal network activity and reveals that asynchronous release enhances utilization of the recycling vesicle pool to induce the network effect. The model further predicts a positive correlation between vesicle priming at the single-neuron level and burst frequency at the network level; this prediction is supported by experimental findings. Thus, the model is utilized to reveal how synaptic release processes at the neuronal level govern activity patterns and synchronization at the network level.

The activity of neuronal networks underlies basic neural functions such as sleep, learning and sensorimotor gating. Computational models of neuronal networks have been developed to capture the complexity of the network activity and predict how neuronal networks generate spontaneous activity. However, most computational models do not simulate the intricate synaptic release process that governs the interaction between neurons and has been shown to significantly impact neuronal network activity and animal behavior, learning and memory. Our paper demonstrates the importance of simulating the elaborate synaptic release process to understand how neuronal networks generate spontaneous activity and respond to manipulations of the release process. The model provides mechanistic explanations and predictions for experimental pharmacological and genetic manipulations. Thus, the model presents a novel computational platform to understand how mechanistic changes in the synaptic release process modulate network oscillatory activity that might impact basic neural functions.

Scalability of Asynchronous Networks Is Limited by One-to-One Mapping between Effective Connectivity and Correlations

Network models are routinely downscaled compared to nature in terms of numbers of nodes or edges because of a lack of computational resources, often without explicit mention of the limitations this entails. While reliable methods have long existed to adjust parameters such that the first-order statistics of network dynamics are conserved, here we show that limitations already arise if also second-order statistics are to be maintained. The temporal structure of pairwise averaged correlations in the activity of recurrent networks is determined by the effective population-level connectivity. We first show that in general the converse is also true and explicitly mention degenerate cases when this one-to-one relationship does not hold. The one-to-one correspondence between effective connectivity and the temporal structure of pairwise averaged correlations implies that network scalings should preserve the effective connectivity if pairwise averaged correlations are to be held constant. Changes in effective connectivity can even push a network from a linearly stable to an unstable, oscillatory regime and vice versa. On this basis, we derive conditions for the preservation of both mean population-averaged activities and pairwise averaged correlations under a change in numbers of neurons or synapses in the asynchronous regime typical of cortical networks. We find that mean activities and correlation structure can be maintained by an appropriate scaling of the synaptic weights, but only over a range of numbers of synapses that is limited by the variance of external inputs to the network. Our results therefore show that the reducibility of asynchronous networks is fundamentally limited.

Prolonged wakefulness alters neuronal responsiveness to local electrical stimulation of the neocortex in awake rats

Prolonged wakefulness or a lack of sleep lead to cognitive deficits, but little is known about the underlying cellular mechanisms. We recently found that sleep deprivation affects spontaneous neuronal activity in the neocortex of sleeping and awake rats. While it is well known that synaptic responses are modulated by ongoing cortical activity, it remains unclear whether prolonged waking affects responsiveness of cortical neurons to incoming stimuli. By applying local electrical microstimulation to the frontal area of the neocortex, we found that after a 4 h period of waking the initial neuronal response in the contralateral frontal cortex was stronger and more synchronous, and was followed by a more profound inhibition of neuronal spiking as compared with the control condition. These changes in evoked activity suggest increased neuronal excitability and indicate that, after staying awake, cortical neurons become transiently bistable. We propose that some of the detrimental effects of sleep deprivation may be a result of altered neuronal responsiveness to incoming intrinsic and extrinsic inputs.

Thalamocortical Sensorimotor Circuit Damage Associated with Disorders of Consciousness for Diffuse Axonal Injury Patients

The relationship of structural and functional brain damage and disorders of consciousness (DOC) for diffuse axonal injury (DAI) is still not fully explored. We employed diffusion tensor imaging (DTI) and resting-state fMRI (RS-fMRI) to examine the changes of resting activations and white matter (WM) integrity for DAI with DOC. WM damages were observed in the body and genu of the corpus callosum, right external capsule (EC) and superior corona radiate (SCR), left superior cerebellar peduncle (SCP) and posterior thalamic radiation (PTR). The RS-fMRI revealed augmented amplitude of low-frequency fluctuation (ALFF) in the anterior cingulate cortex, hippocampus, insula, amygdala and putamen, and reduced ALFF in the precuneus, thalamus, pre-central and post-central gyri. Correlation analysis identified positive associations between the Glasgow Coma Scale (GCS) and activation of the precuneus and between GCS and DTI measurements in the left PTR and SCP, but a negative correlation was found between GCS and activation of the thalamus. Cross modality association analyses indicated that activations of the amygdala and postcentral gyrus were correlated with DTI measurements of the right EC and left PTR respectively. These results implicate that the WM damages in thalamocortical sensorimotor circuit and aberrant brain activity responding to self-awareness and sensation are critical factors to DOC, which expand the current understanding of the neural mechanisms underlying DAI.

A distinct class of slow (~0.2-2 Hz) intrinsically bursting layer 5 pyramidal neurons determines UP/DOWN state dynamics in the neocortex

During sleep and anesthesia, neocortical neurons exhibit rhythmic UP/DOWN membrane potential states. Although UP states are maintained by synaptic activity, the mechanisms that underlie the initiation and robust rhythmicity of UP states are unknown. Using a physiologically validated model of UP/DOWN state generation in mouse neocortical slices whereby the cholinergic tone present in vivo is reinstated, we show that the regular initiation of UP states is driven by an electrophysiologically distinct subset of morphologically identified layer 5 neurons, which exhibit intrinsic rhythmic low-frequency burst firing at ~0.2-2 Hz. This low-frequency bursting is resistant to block of glutamatergic and GABAergic transmission but is absent when slices are maintained in a low Ca(2+) medium (an alternative, widely used model of cortical UP/DOWN states), thus explaining the lack of rhythmic UP states and abnormally prolonged DOWN states in this condition. We also characterized the activity of various other pyramidal and nonpyramidal neurons during UP/DOWN states and found that an electrophysiologically distinct subset of layer 5 regular spiking pyramidal neurons fires earlier during the onset of network oscillations compared with all other types of neurons recorded. This study, therefore, identifies an important role for cell-type-specific neuronal activity in driving neocortical UP states.

Looking for a precursor of spontaneous Sleep Slow Oscillations in human sleep: The role of the sigma activity

Sleep Slow Oscillations (SSOs), paradigmatic EEG markers of cortical bistability (alternation between cellular downstates and upstates), and sleep spindles, paradigmatic EEG markers of thalamic rhythm, are two hallmarks of sleeping brain. Selective thalamic lesions are reportedly associated to reductions of spindle activity and its spectrum ~14 Hz (sigma), and to alterations of SSO features. This apparent, parallel behavior suggests that thalamo-cortical entrainment favors cortical bistability. Here we investigate temporally-causal associations between thalamic sigma activity and shape, topology, and dynamics of SSOs. We recorded sleep EEG and studied whether spatio-temporal variability of SSO amplitude, negative slope (synchronization in downstate falling) and detection rate are driven by cortical-sigma-activity expression (12-18Hz), in 3 consecutive 1s-EEG-epochs preceding each SSO event (Baselines). We analyzed: (i) spatial variability, comparing maps of baseline sigma power and of SSO features, averaged over the first sleep cycle; (ii) event-by-event shape variability, computing for each electrode correlations between baseline sigma power and amplitude/slope of related SSOs; (iii) event-by-event spreading variability, comparing baseline sigma power in electrodes showing an SSO event with the homologous ones, spared by the event. The scalp distribution of baseline sigma power mirrored those of SSO amplitude and slope; event-by-event variability in baseline sigma power was associated with that in SSO amplitude in fronto-central areas; within each SSO event, electrodes involved in cortical bistability presented higher baseline sigma activity than those free of SSO. In conclusion, spatio-temporal variability of thalamocortical entrainment, measured by background sigma activity, is a reliable estimate of the cortical proneness to bistability.

Mapping human brain networks with cortico-cortical evoked potentials

The cerebral cortex forms a sheet of neurons organized into a network of interconnected modules that is highly expanded in humans and presumably enables our most refined sensory and cognitive abilities. The links of this network form a fundamental aspect of its organization, and a great deal of research is focusing on understanding how information flows within and between different regions. However, an often-overlooked element of this connectivity regards a causal, hierarchical structure of regions, whereby certain nodes of the cortical network may exert greater influence over the others. While this is difficult to ascertain non-invasively, patients undergoing invasive electrode monitoring for epilepsy provide a unique window into this aspect of cortical organization. In this review, we highlight the potential for cortico-cortical evoked potential (CCEP) mapping to directly measure neuronal propagation across large-scale brain networks with spatio-temporal resolution that is superior to traditional neuroimaging methods. We first introduce effective connectivity and discuss the mechanisms underlying CCEP generation. Next, we highlight how CCEP mapping has begun to provide insight into the neural basis of non-invasive imaging signals. Finally, we present a novel approach to perturbing and measuring brain network function during cognitive processing. The direct measurement of CCEPs in response to electrical stimulation represents a potentially powerful clinical and basic science tool for probing the large-scale networks of the human cerebral cortex.

Sleep smart-optimizing sleep for declarative learning and memory

The last decade has witnessed a spurt of new publications documenting sleep's essential contribution to the brains ability to form lasting memories. For the declarative memory domain, slow wave sleep (the deepest sleep stage) has the greatest beneficial effect on the consolidation of memories acquired during preceding wakefulness. The finding that newly encoded memories become reactivated during subsequent sleep fostered the idea that reactivation leads to the strengthening and transformation of the memory trace. According to the active system consolidation account, trace reactivation leads to the redistribution of the transient memory representations from the hippocampus to the long-lasting knowledge networks of the cortex. Apart from consolidating previously learned information, sleep also facilitates the encoding of new memories after sleep, which probably relies on the renormalization of synaptic weights during sleep as suggested by the synaptic homeostasis theory. During wakefulness overshooting potentiation causes an imbalance in synaptic weights that is countered by synaptic downscaling during subsequent sleep. This review briefly introduces the basic concepts and central findings of the research on sleep and memory, and discusses implications of this lab-based work for everyday applications to make the best possible use of sleep's beneficial effect on learning and memory.

Neocortical focus: experimental view

All brain normal or pathological activities occur in one of the states of vigilance: wake, slow-wave sleep, or REM sleep. Neocortical seizures preferentially occur during slow-wave sleep. We provide a description of neuronal behavior and mechanisms mediating such a behavior within neocortex taking place in natural states of vigilance as well as during seizures pointing to similarities and differences exhibited during sleep and seizures. A concept of epileptic focus is described using a model of cortical undercut, because in that model, the borders of the focus are well defined. In this model, as in other models of acquired epilepsy, the main factor altering excitability is deafferentation, which upregulates neuronal excitability that promotes generation of seizures. Periods of disfacilitation recorded during slow-wave sleep further upregulate neuronal excitability. It appears that the state of neurons and neuronal network in the epileptic focus produced by deafferentation are such that seizures cannot be generated there. Instead, seizures always start around the perimeter of the undercut cortex. Therefore, we define these areas as the seizure focus. In this zone, neuronal connectivity and excitability are moderately enhanced, lowering the threshold for seizure generation.

Stochastic transitions into silence cause noise correlations in cortical circuits

The spiking activity of cortical neurons is highly variable. This variability is generally correlated among nearby neurons, an effect commonly interpreted to reflect the coactivation of neurons due to anatomically shared inputs. Recent findings, however, indicate that correlations can be dynamically modulated, suggesting that the underlying mechanisms are not well understood. Here, we investigate the hypothesis that correlations are dominated by neuronal coinactivation: the occurrence of brief silent periods during which all neurons in the local network stop firing. We recorded spiking activity from large populations of neurons in the auditory cortex of anesthetized rats across different brain states. During spontaneous activity, the reduction of correlation accompanying brain state desynchronization was largely explained by a decrease in the density of the silent periods. The presentation of a stimulus caused an initial drop of correlations followed by a rebound, a time course that was mimicked by the instantaneous silence density. We built a rate network model with fluctuation-driven transitions between a silent and an active attractor and assumed that neurons fired Poisson spike trains with a rate following the model dynamics. Variations of the network external input altered the transition rate into the silent attractor and reproduced the relation between correlation and silence density found in the data, both in spontaneous and evoked conditions. This suggests that the observed changes in correlation, occurring gradually with brain state variations or abruptly with sensory stimulation, are due to changes in the likeliness of the microcircuit to transiently cease firing.

Slow-wave oscillations in a corticothalamic model of sleep and wake

A physiologically-based corticothalamic neural field model is used to study slow wave oscillations including cortical UP and DOWN states in deep sleep by extending it to incorporate bursting dynamics of neurons in the thalamic reticular nucleus. The interplay of local bursting dynamics and network interactions produces the cortical UP and DOWN states of slow wave sleep while preserving previously verified model predictions in the wake state. Results show that EEG spectral features in wake and sleep are reproduced. The bursting is subthreshold but acts to intensify the amplitude of oscillations in slow wave sleep with deep UP/DOWN oscillations on the cortex emerging naturally. Furthermore, there is a continuous cycle between the two regimes, rather than a flip-flop between discrete states.

Contribution of parvalbumin and somatostatin-expressing GABAergic neurons to slow oscillations and the balance in beta-gamma oscillations across cortical layers

Cortical interneurons are classified into several subtypes that contribute to cortical oscillatory activity. Parvalbumin (PV)-expressing cells, a type of inhibitory interneuron, are involved in the gamma oscillations of local field potentials (LFPs). Under ketamine-xylazine anesthesia or sleep, mammalian cortical circuits exhibit slow oscillations in which the active-up state and silent-down state alternate at ~1 Hz. The up state is composed of various high-frequency oscillations, including gamma oscillations. However, it is unclear how PV cells and somatostatin (SOM) cells contribute to the slow oscillations and the high-frequency oscillations nested in the up state. To address these questions, we used mice lacking glutamate decarboxylase 67, primarily in PV cells (PV-GAD67 mice) or in SOM cells (SOM-GAD67 mice). We then compared LFPs between PV-GAD67 mice and SOM-GAD67 mice. PV cells target the proximal regions of pyramidal cells, whereas SOM cells are dendrite-preferring interneurons. We found that the up state was shortened in duration in the PV-GAD67 mice, but tended to be longer in SOM-GAD67 mice. Firing rate tended to increase in PV-GAD67 mice, but tended to decrease in SOM-GAD67 mice. We also found that delta oscillations tended to increase in SOM-GAD67 mice, but tended to decrease in PV-GAD67 mice. Current source density and wavelet analyses were performed to determine the depth profiles of various high-frequency oscillations. High gamma and ripple (60–200 Hz) power decreased in the neocortical upper layers specifically in PV-GAD67 mice, but not in SOM-GAD67. In addition, beta power (15–30 Hz) increased in the deep layers, specifically in PV-GAD67 mice. These results suggest that PV cells play important roles in persistence of the up state and in the balance between gamma and beta bands across cortical layers, whereas SOM and PV cells may make an asymmetric contribution to regulate up-state and delta oscillations.

Enhancement of sleep slow waves: underlying mechanisms and practical consequences

Even modest sleep restriction, especially the loss of sleep slow wave activity (SWA), is invariably associated with slower electroencephalogram (EEG) activity during wake, the occurrence of local sleep in an otherwise awake brain, and impaired performance due to cognitive and memory deficits. Recent studies not only confirm the beneficial role of sleep in memory consolidation, but also point to a specific role for sleep slow waves. Thus, the implementation of methods to enhance sleep slow waves without unwanted arousals or lightening of sleep could have significant practical implications. Here we first review the evidence that it is possible to enhance sleep slow waves in humans using transcranial direct-current stimulation (tDCS) and transcranial magnetic stimulation. Since these methods are currently impractical and their safety is questionable, especially for chronic long-term exposure, we then discuss novel data suggesting that it is possible to enhance slow waves using sensory stimuli. We consider the physiology of the K-complex (KC), a peripheral evoked slow wave, and show that, among different sensory modalities, acoustic stimulation is the most effective in increasing the magnitude of slow waves, likely through the activation of non-lemniscal ascending pathways to the thalamo-cortical system. In addition, we discuss how intensity and frequency of the acoustic stimuli, as well as exact timing and pattern of stimulation, affect sleep enhancement. Finally, we discuss automated algorithms that read the EEG and, in real-time, adjust the stimulation parameters in a closed-loop manner to obtain an increase in sleep slow waves and avoid undesirable arousals. In conclusion, while discussing the mechanisms that underlie the generation of sleep slow waves, we review the converging evidence showing that acoustic stimulation is safe and represents an ideal tool for slow wave sleep (SWS) enhancement.

Adaptation to prolonged neuromodulation in cortical cultures: an invariable return to network synchrony

Background

Prolonged neuromodulatory regimes, such as those critically involved in promoting arousal and suppressing sleep-associated synchronous activity patterns, might be expected to trigger adaptation processes and, consequently, a decline in neuromodulator-driven effects. This possibility, however, has rarely been addressed.

Results

Using networks of cultured cortical neurons, acetylcholine microinjections and a novel closed-loop ‘synchrony-clamp’ system, we found that acetylcholine pulses strongly suppressed network synchrony. Over the course of many hours, however, synchrony invariably reemerged, even when feedback was used to compensate for declining cholinergic efficacy. Network synchrony also reemerged following its initial suppression by noradrenaline, but this did not occlude the suppression of synchrony or its gradual reemergence following subsequent cholinergic input. Importantly, cholinergic efficacy could be restored and preserved over extended time scales by periodically withdrawing cholinergic input.

Conclusions

These findings indicate that the capacity of neuromodulators to suppress network synchrony is constrained by slow-acting, reactive processes. A multiplicity of neuromodulators and ultimately neuromodulator withdrawal periods might thus be necessary to cope with an inevitable reemergence of network synchrony.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0083-3) contains supplementary material, which is available to authorized users.

A cortical circuit for gain control by behavioral state

The brain's response to sensory input is strikingly modulated by behavioral state. Notably, the visual response of mouse primary visual cortex (V1) is enhanced by locomotion, a tractable and accessible example of a time-locked change in cortical state. The neural circuits that transmit behavioral state to sensory cortex to produce this modulation are unknown. In vivo calcium imaging of behaving animals revealed that locomotion activates vasoactive intestinal peptide (VIP)-positive neurons in mouse V1 independent of visual stimulation and largely through nicotinic inputs from basal forebrain. Optogenetic activation of VIP neurons increased V1 visual responses in stationary awake mice, artificially mimicking the effect of locomotion, and photolytic damage of VIP neurons abolished the enhancement of V1 responses by locomotion. These findings establish a cortical circuit for the enhancement of visual response by locomotion and provide a potential common circuit for the modulation of sensory processing by behavioral state.

Engineering a thalamo-cortico-thalamic circuit on SpiNNaker: a preliminary study toward modeling sleep and wakefulness

We present a preliminary study of a thalamo-cortico-thalamic (TCT) implementation on SpiNNaker (Spiking Neural Network architecture), a brain inspired hardware platform designed to incorporate the inherent biological properties of parallelism, fault tolerance and energy efficiency. These attributes make SpiNNaker an ideal platform for simulating biologically plausible computational models. Our focus in this work is to design a TCT framework that can be simulated on SpiNNaker to mimic dynamical behavior similar to Electroencephalogram (EEG) time and power-spectra signatures in sleep-wake transition. The scale of the model is minimized for simplicity in this proof-of-concept study; thus the total number of spiking neurons is ≈1000 and represents a "mini-column" of the thalamocortical tissue. All data on model structure, synaptic layout and parameters is inspired from previous studies and abstracted at a level that is appropriate to the aims of the current study as well as computationally suitable for model simulation on a small 4-chip SpiNNaker system. The initial results from selective deletion of synaptic connectivity parameters in the model show similarity with EEG power spectra characteristics of sleep and wakefulness. These observations provide a positive perspective and a basis for future implementation of a very large scale biologically plausible model of thalamo-cortico-thalamic interactivity-the essential brain circuit that regulates the biological sleep-wake cycle and associated EEG rhythms.

CyNEST: a maintainable Cython-based interface for the NEST simulator

NEST is a simulator for large-scale networks of spiking point neuron models (Gewaltig and Diesmann, 2007). Originally, simulations were controlled via the Simulation Language Interpreter (SLI), a built-in scripting facility implementing a language derived from PostScript (Adobe Systems, Inc., 1999). The introduction of PyNEST (Eppler et al., 2008), the Python interface for NEST, enabled users to control simulations using Python. As the majority of NEST users found PyNEST easier to use and to combine with other applications, it immediately displaced SLI as the default NEST interface. However, developing and maintaining PyNEST has become increasingly difficult over time. This is partly because adding new features requires writing low-level C++ code intermixed with calls to the Python/C API, which is unrewarding. Moreover, the Python/C API evolves with each new version of Python, which results in a proliferation of version-dependent code branches. In this contribution we present the re-implementation of PyNEST in the Cython language, a superset of Python that additionally supports the declaration of C/C++ types for variables and class attributes, and provides a convenient foreign function interface (FFI) for invoking C/C++ routines (Behnel et al., 2011). Code generation via Cython allows the production of smaller and more maintainable bindings, including increased compatibility with all supported Python releases without additional burden for NEST developers. Furthermore, this novel approach opens up the possibility to support alternative implementations of the Python language at no cost given a functional Cython back-end for the corresponding implementation, and also enables cross-compilation of Python bindings for embedded systems and supercomputers alike.

Modeling the dynamical effects of anesthesia on brain circuits

General anesthesia is a neurophysiological state that consists of unconsciousness, amnesia, analgesia, and immobility along with maintenance of physiological stability. General anesthesia has been used in the United States for more than 167 years. Now, using systems neuroscience paradigms how anesthetics act in the brain and central nervous system to create the states of general anesthesia is being understood. Propofol is one of the most widely used and the most widely studied anesthetics. When administered for general anesthesia or sedation, the electroencephalogram (EEG) under propofol shows highly structured, rhythmic activity that is strongly associated with changes in the patient's level of arousal. These highly structured oscillations lend themselves readily to mathematical descriptions using dynamical systems models. We review recent model descriptions of the commonly observed EEG patterns associated with propofol: paradoxical excitation, strong frontal alpha oscillations, anteriorization and burst suppression. Our analysis suggests that propofol's actions at GABAergic networks in the cortex, thalamus and brainstem induce profound brain dynamics that are one of the likely mechanisms through which this anesthetic induces altered arousal states from sedation to unconsciousness. Because these dynamical effects are readily observed in the EEG, the mathematical descriptions of how propofol's EEG signatures relate to its mechanisms of action in neural circuits provide anesthesiologists with a neurophysiologically based approach to monitoring the brain states of patients receiving anesthesia care.

Quantifying cortical EEG responses to TMS in (un)consciousness

We normally assess another individual's level of consciousness based on her or his ability to interact with the surrounding environment and communicate. Usually, if we observe purposeful behavior, appropriate responses to sensory inputs, and, above all, appropriate answers to questions, we can be reasonably sure that the person is conscious. However, we know that consciousness can be entirely within the brain, even in the absence of any interaction with the external world; this happens almost every night, while we dream. Yet, to this day, we lack an objective, dependable measure of the level of consciousness that is independent of processing sensory inputs and producing appropriate motor outputs. Theoretically, consciousness is thought to require the joint presence of functional integration and functional differentiation, otherwise defined as brain complexity. Here we review a series of recent studies in which Transcranial Magnetic Stimulation combined with electroencephalography (TMS/EEG) has been employed to quantify brain complexity in wakefulness and during physiological (sleep), pharmacological (anesthesia) and pathological (brain injury) loss of consciousness. These studies invariably show that the complexity of the cortical response to TMS collapses when consciousness is lost during deep sleep, anesthesia and vegetative state following severe brain injury, while it recovers when consciousness resurges in wakefulness, during dreaming, in the minimally conscious state or locked-in syndrome. The present paper will also focus on how this approach may contribute to unveiling the pathophysiology of disorders of consciousness affecting brain-injured patients. Finally, we will underline some crucial methodological aspects concerning TMS/EEG measurements of brain complexity.

Sleep and synaptic homeostasis

In the last decades a substantial knowledge about sleep mechanisms has been accumulated. However, the function of sleep still remains elusive. The difficulty with unraveling sleep's function may arise from the lack of understanding of how the multitude of processes associated with waking and sleep-from gene expression and single neuron activity to the whole brain dynamics and behavior-functionally and mechanistically relate to each other. Therefore, novel conceptual frameworks, which integrate and take into account the variety of phenomena occurring during waking and sleep at different levels, will likely lead to advances in our understanding of the function of sleep, above and beyond what merely descriptive or correlative approaches can provide. One such framework, the synaptic homeostasis hypothesis, focuses on wake- and sleep-dependent changes in synaptic strength. The core claim of this hypothesis is that learning and experience during wakefulness are associated with a net increase in synaptic strength. In turn, the proposed function of sleep is to provide synaptic renormalization, which has important implications with respect to energy needs, intracranial space, metabolic supplies, and, importantly, enables further plastic changes. In this article we review the empirical evidence for this hypothesis, which was obtained at several levels-from gene expression and cellular excitability to structural synaptic modifications and behavioral outcomes. We conclude that although the mechanisms behind the proposed role of sleep in synaptic homeostasis are undoubtedly complex, this conceptual framework offers a unique opportunity to provide mechanistic and functional explanation for many previously disparate observations, and define future research strategies.

Brain repair after stroke--a novel neurological model

Following stroke, patients are commonly left with debilitating motor and speech impairments. This article reviews the state of the art in neurological repair for stroke and proposes a new model for the future. We suggest that stroke treatment--from the time of the ictus itself to living with the consequences--must be fundamentally neurological, from limiting the extent of injury at the outset, to repairing the consequent damage. Our model links brain and behaviour by targeting brain circuits, and we illustrate the model though action observation treatment, which aims to enhance brain network connectivity. The model is based on the assumptions that the mechanisms of neural repair inherently involve cellular and circuit plasticity, that brain plasticity is a synaptic phenomenon that is largely stimulus-dependent, and that brain repair required both physical and behavioural interventions that are tailored to reorganize specific brain circuits. We review current approaches to brain repair after stroke and present our new model, and discuss the biological foundations, rationales, and data to support our novel approach to upper-extremity and language rehabilitation. We believe that by enhancing plasticity at the level of brain network interactions, this neurological model for brain repair could ultimately lead to a cure for stroke.

How delays matter in an oscillatory whole-brain spiking-neuron network model for MEG alpha-rhythms at rest

In recent years the study of the intrinsic brain dynamics in a relaxed awake state in the absence of any specific task has gained increasing attention, as spontaneous neural activity has been found to be highly structured at a large scale. This so called resting-state activity has been found to be comprised by nonrandom spatiotemporal patterns and fluctuations, and several Resting-State Networks (RSN) have been found in BOLD-fMRI as well as in MEG signal power envelope correlations. The underlying anatomical connectivity structure between areas of the brain has been identified as being a key to the observed functional network connectivity, but the mechanisms behind this are still underdetermined. Theoretical large-scale brain models for fMRI data have corroborated the importance of the connectome in shaping network dynamics, while the importance of delays and noise differ between studies and depend on the models' specific dynamics. In the current study, we present a spiking neuron network model that is able to produce noisy, distributed alpha-oscillations, matching the power peak in the spectrum of group resting-state MEG recordings. We studied how well the model captured the inter-node correlation structure of the alpha-band power envelopes for different delays between brain areas, and found that the model performs best for propagation delays inside the physiological range (5-10 m/s). Delays also shift the transition from noisy to bursting oscillations to higher global coupling values in the model. Thus, in contrast to the asynchronous fMRI state, delays are important to consider in the presence of oscillation.

Plastic Changes Following Imitation-Based Speech and Language Therapy for Aphasia

Background. Objective measurement of plastic brain changes induced by a novel rehabilitative approach is a key requirement for validating its biological rationale linking the potential therapeutic gains to the changes in brain physiology. Objective. Based on an emerging notion linking cortical plastic changes to EEG sleep slow-wave activity (SWA) regulation, we aimed to assess the acute plastic changes induced by an imitation-based speech therapy in individuals with aphasia by comparing sleep SWA changes before and after therapy. Methods. A total of 13 left-hemispheric stroke patients underwent language assessment with the Western Aphasia Battery (WAB) before and after 2 consecutive high-density (hd) EEG sleep recordings interleaved by a daytime session of imitation-based speech therapy (Intensive Mouth Imitation and Talking for Aphasia Therapeutic Effects [IMITATE]). This protocol is thought to stimulate bilateral connections between the inferior parietal lobule and the ventral premotor areas. Results. A single exposure to IMITATE resulted in increases in local EEG SWA during subsequent sleep over the same regions predicted by the therapeutic rationale, particularly over the right hemisphere (unaffected by the lesion). Furthermore, changes in SWA over the left-precentral areas predicted changes in WAB repetition scores in our group, supporting the role of perilesional areas in predicting positive functional responses. Conclusions. Our results suggest that SWA changes occurring in brain areas activated during imitation-based aphasia therapy may reflect the acute plastic changes induced by this intervention. Further testing will be needed to evaluate SWA as a non-invasive assessment of changes induced by the therapy and as a predictor of positive long-term clinical outcome.

Conductance-based refractory density model of primary visual cortex

A layered continual population model of primary visual cortex has been constructed, which reproduces a set of experimental data, including postsynaptic responses of single neurons on extracellular electric stimulation and spatially distributed activity patterns in response to visual stimulation. In the model, synaptically interacting excitatory and inhibitory neuronal populations are described by a conductance-based refractory density approach. Populations of two-compartment excitatory and inhibitory neurons in cortical layers 2/3 and 4 are distributed in the 2-d cortical space and connected by AMPA, NMDA and GABA type synapses. The external connections are pinwheel-like, according to the orientation of a stimulus. Intracortical connections are isotropic local and patchy between neurons with similar orientations. The model proposes better temporal resolution and more detailed elaboration than conventional mean-field models. In comparison to large network simulations, it excludes a posteriori statistical data manipulation and provides better computational efficiency and minimal parametrization.

Modeling resting-state functional networks when the cortex falls asleep: local and global changes

The transition from wakefulness to sleep represents the most conspicuous change in behavior and the level of consciousness occurring in the healthy brain. It is accompanied by similarly conspicuous changes in neural dynamics, traditionally exemplified by the change from "desynchronized" electroencephalogram activity in wake to globally synchronized slow wave activity of early sleep. However, unit and local field recordings indicate that the transition is more gradual than it might appear: On one hand, local slow waves already appear during wake; on the other hand, slow sleep waves are only rarely global. Studies with functional magnetic resonance imaging also reveal changes in resting-state functional connectivity (FC) between wake and slow wave sleep. However, it remains unclear how resting-state networks may change during this transition period. Here, we employ large-scale modeling of the human cortico-cortical anatomical connectivity to evaluate changes in resting-state FC when the model "falls asleep" due to the progressive decrease in arousal-promoting neuromodulation. When cholinergic neuromodulation is parametrically decreased, local slow waves appear, while the overall organization of resting-state networks does not change. Furthermore, we show that these local slow waves are structured macroscopically in networks that resemble the resting-state networks. In contrast, when the neuromodulator decrease further to very low levels, slow waves become global and resting-state networks merge into a single undifferentiated, broadly synchronized network.

Is sleep essential for neural plasticity in humans, and how does it affect motor and cognitive recovery?

There is a general consensus that sleep is strictly linked to memory, learning, and, in general, to the mechanisms of neural plasticity, and that this link may directly affect recovery processes. In fact, a coherent pattern of empirical findings points to beneficial effect of sleep on learning and plastic processes, and changes in synaptic plasticity during wakefulness induce coherent modifications in EEG slow wave cortical topography during subsequent sleep. However, the specific nature of the relation between sleep and synaptic plasticity is not clear yet. We reported findings in line with two models conflicting with respect to the underlying mechanisms, that is, the “synaptic homeostasis hypothesis” and the “consolidation” hypothesis, and some recent results that may reconcile them. Independently from the specific mechanisms involved, sleep loss is associated with detrimental effects on plastic processes at a molecular and electrophysiological level. Finally, we reviewed growing evidence supporting the notion that plasticity-dependent recovery could be improved managing sleep quality, while monitoring EEG during sleep may help to explain how specific rehabilitative paradigms work. We conclude that a better understanding of the sleep-plasticity link could be crucial from a rehabilitative point of view.

A spontaneous state of weakly correlated synaptic excitation and inhibition in visual cortex

Cortical spontaneous activity reflects an animal's behavioral state and affects neural responses to sensory stimuli. The correlation between excitatory and inhibitory synaptic input to single neurons is a key parameter in models of cortical circuitry. Recent measurements demonstrated highly correlated synaptic excitation and inhibition during spontaneous "up-and-down" states, during which excitation accounted for approximately 80% of inhibitory variance (Shu et al., 2003; Haider et al., 2006). Here we report in vivo whole-cell estimates of the correlation between excitation and inhibition in the rat visual cortex under pentobarbital anesthesia, during which up-and-down states are absent. Excitation and inhibition are weakly correlated, relative to the up-and-down state: excitation accounts for less than 40% of inhibitory variance. Although these correlations are lower than when the circuit cycles between up-and-down states, both behaviors may arise from the same circuitry. Our observations provide evidence that different correlational patterns of excitation and inhibition underlie different cortical states.

Forebrain HCN1 channels contribute to hypnotic actions of ketamine

BACKGROUND

Ketamine is a commonly used anesthetic, but the mechanistic basis for its clinically relevant actions remains to be determined. The authors previously showed that HCN1 channels are inhibited by ketamine and demonstrated that global HCN1 knockout mice are twofold less sensitive to hypnotic actions of ketamine. Although that work identified HCN1 channels as a viable molecular target for ketamine, it did not determine the relevant neural substrate.

METHODS

To localize the brain region responsible for HCN1-mediated hypnotic actions of ketamine, the authors used a conditional knockout strategy to delete HCN1 channels selectively in excitatory cells of the mouse forebrain. A combination of molecular, immunohistochemical, and cellular electrophysiologic approaches was used to verify conditional HCN1 deletion; a loss-of-righting reflex assay served to ascertain effects of forebrain HCN1 channel ablation on hypnotic actions of ketamine.

RESULTS

In conditional knockout mice, HCN1 channels were selectively deleted in cortex and hippocampus, with expression retained in cerebellum. In cortical pyramidal neurons from forebrain-selective HCN1 knockout mice, effects of ketamine on HCN1-dependent membrane properties were absent; notably, ketamine was unable to evoke membrane hyperpolarization or enhance synaptic inputs. Finally, the EC50 for ketamine-induced loss-of-righting reflex was shifted to significantly higher concentrations (by approximately 31%).

CONCLUSIONS

These data indicate that forebrain principal cells represent a relevant neural substrate for HCN1-mediated hypnotic actions of ketamine. The authors suggest that ketamine inhibition of HCN1 shifts cortical neuron electroresponsive properties to contribute to ketamine-induced hypnosis.

Disrupted directed connectivity along the cingulate cortex determines vigilance after sleep deprivation

The cingulate cortex is regarded as the backbone of structural and functional connectivity of the brain. While its functional connectivity has been intensively studied, little is known about its effective connectivity, its modulation by behavioral states, and its involvement in cognitive performance. Given the previously reported effects on cingulate functional connectivity, we investigated how eye-closure and sleep deprivation changed cingulate effective connectivity, estimated from resting-state high-density electroencephalography (EEG) using a novel method to calculate Granger Causality directly in source space. Effective connectivity along the cingulate cortex was dominant in the forward direction. Eyes-open connectivity in the forward direction was greater compared to eyes-closed, in well-rested participants. The difference between eyes-open and eyes-closed connectivity was attenuated and no longer significant after sleep deprivation. Individual variability in the forward connectivity after sleep deprivation predicted subsequent task performance, such that those subjects who showed a greater increase in forward connectivity between the eyes-open and the eyes-closed periods also performed better on a sustained attention task. Effective connectivity in the opposite, backward, direction was not affected by whether the eyes were open or closed or by sleep deprivation. These findings indicate that the effective connectivity from posterior to anterior cingulate regions is enhanced when a well-rested subject has his eyes open compared to when they are closed. Sleep deprivation impairs this directed information flow, proportional to its deleterious effect on vigilance. Therefore, sleep may play a role in the maintenance of waking effective connectivity.

Cortical information flow in Parkinson's disease: a composite network/field model

The basal ganglia play a crucial role in the execution of movements, as demonstrated by the severe motor deficits that accompany Parkinson's disease (PD). Since motor commands originate in the cortex, an important question is how the basal ganglia influence cortical information flow, and how this influence becomes pathological in PD. To explore this, we developed a composite neuronal network/neural field model. The network model consisted of 4950 spiking neurons, divided into 15 excitatory and inhibitory cell populations in the thalamus and cortex. The field model consisted of the cortex, thalamus, striatum, subthalamic nucleus, and globus pallidus. Both models have been separately validated in previous work. Three field models were used: one with basal ganglia parameters based on data from healthy individuals, one based on data from individuals with PD, and one purely thalamocortical model. Spikes generated by these field models were then used to drive the network model. Compared to the network driven by the healthy model, the PD-driven network had lower firing rates, a shift in spectral power toward lower frequencies, and higher probability of bursting; each of these findings is consistent with empirical data on PD. In the healthy model, we found strong Granger causality between cortical layers in the beta and low gamma frequency bands, but this causality was largely absent in the PD model. In particular, the reduction in Granger causality from the main “input” layer of the cortex (layer 4) to the main “output” layer (layer 5) was pronounced. This may account for symptoms of PD that seem to reflect deficits in information flow, such as bradykinesia. In general, these results demonstrate that the brain's large-scale oscillatory environment, represented here by the field model, strongly influences the information processing that occurs within its subnetworks. Hence, it may be preferable to drive spiking network models with physiologically realistic inputs rather than pure white noise.

DOUBLE-LAYERED MODELS CAN EXPLAIN MACRO AND MICRO STRUCTURE OF HUMAN SLEEP

The model simulates the activity of three neural populations using a Lotka–Volterra predator–prey system and, based on neuro-anatomical and neuro-physiological recent findings, assumes that a functional thalamo-cortical gate should be crossed by 'queuing' thalamic signals and that a sleep promoting substance acts as a modulator. The resultant activity accounts for the sleep stage transitions. In accordance with sleep cycles timing, the model proves to be able to reproduce the clustering and randomness of those peculiar transient synchronized EEG patterns (TSEP) described in normal human sleep and supposed to be related to the dynamic building up of NREM sleep until its stabilization against perturbations.

DOC2B and Munc13-1 differentially regulate neuronal network activity

Alterations in the levels of synaptic proteins affect synaptic transmission and synaptic plasticity. However, the precise effects on neuronal network activity are still enigmatic. Here, we utilized microelectrode array (MEA) to elucidate how manipulation of the presynaptic release process affects the activity of neuronal networks. By combining pharmacological tools and genetic manipulation of synaptic proteins, we show that overexpression of DOC2B and Munc13-1, proteins known to promote vesicular maturation and release, elicits opposite effects on the activity of the neuronal network. Although both cause an increase in the overall number of spikes, the distribution of spikes is different. While DOC2B enhances, Munc13-1 reduces the firing rate within bursts of spikes throughout the network; however, Munc13-1 increases the rate of network bursts. DOC2B's effects were mimicked by Strontium that elevates asynchronous release but not by a DOC2B mutant that enhances spontaneous release rate. This suggests for the first time that increased asynchronous release on the single-neuron level promotes bursting activity in the network level. This innovative study demonstrates the complementary role of the network level in explaining the physiological relevance of the cellular activity of presynaptic proteins and the transformation of synaptic release manipulation from the neuron to the network level.

Relationships between Electroencephalographic Spectral Peaks Across Frequency Bands

The degree to which electroencephalographic spectral peaks are independent, and the relationships between their frequencies have been debated. A novel fitting method was used to determine peak parameters in the range 2-35 Hz from a large sample of eyes-closed spectra, and their interrelationships were investigated. Findings were compared with a mean-field model of thalamocortical activity, which predicts near-harmonic relationships between peaks. The subject set consisted of 1424 healthy subjects from the Brain Resource International Database. Peaks in the theta range occurred on average near half the alpha peak frequency, while peaks in the beta range tended to occur near twice and three times the alpha peak frequency on an individual-subject basis. Moreover, for the majority of subjects, alpha peak frequencies were significantly positively correlated with frequencies of peaks in the theta and low and high beta ranges. Such a harmonic progression agrees semiquantitatively with theoretical predictions from the mean-field model. These findings indicate a common or analogous source for different rhythms, and help to define appropriate individual frequency bands for peak identification.

Population rate dynamics and multineuron firing patterns in sensory cortex

Cortical circuits encode sensory stimuli through the firing of neuronal ensembles, and also produce spontaneous population patterns in the absence of sensory drive. This population activity is often characterized experimentally by the distribution of multineuron "words" (binary firing vectors), and a match between spontaneous and evoked word distributions has been suggested to reflect learning of a probabilistic model of the sensory world. We analyzed multineuron word distributions in sensory cortex of anesthetized rats and cats, and found that they are dominated by fluctuations in population firing rate rather than precise interactions between individual units. Furthermore, cortical word distributions change when brain state shifts, and similar behavior is seen in simulated networks with fixed, random connectivity. Our results suggest that similarity or dissimilarity in multineuron word distributions could primarily reflect similarity or dissimilarity in population firing rate dynamics, and not necessarily the precise interactions between neurons that would indicate learning of sensory features.

Influence of slow oscillation on hippocampal activity and ripples through cortico-hippocampal synaptic interactions, analyzed by a cortical-CA3-CA1 network model

Hippocampal sharp wave-ripple complexes (SWRs) involve the synchronous discharge of thousands of cells throughout the CA3-CA1-subiculum-entorhinal cortex axis. Their strong transient output affects cortical targets, rendering SWRs a possible means for memory transfer from the hippocampus to the neocortex for long-term storage. Neurophysiological observations of hippocampal activity modulation by the cortical slow oscillation (SO) during deep sleep and anesthesia, and correlations between ripples and UP states, support the role of SWRs in memory consolidation through a cortico-hippocampal feedback loop. We couple a cortical network exhibiting SO with a hippocampal CA3-CA1 computational network model exhibiting SWRs, in order to model such cortico-hippocampal correlations and uncover important parameters and coupling mechanisms controlling them. The cortical oscillatory output entrains the CA3 network via connections representing the mossy fiber input, and the CA1 network via the temporoammonic pathway (TA). The spiking activity in CA3 and CA1 is shown to depend on the excitation-to-inhibition ratio, induced by combining the two hippocampal inputs, with mossy fiber input controlling the UP-state correlation of CA3 population bursts and corresponding SWRs, whereas the temporoammonic input affects the overall CA1 spiking activity. Ripple characteristics and pyramidal spiking participation to SWRs are shaped by the strength of the Schaffer collateral drive. A set of in vivo recordings from the rat hippocampus confirms a model-predicted segregation of pyramidal cells into subgroups according to the SO state where they preferentially fire and their response to SWRs. These groups can potentially play distinct functional roles in the replay of spike sequences.

Local experience-dependent changes in the wake EEG after prolonged wakefulness

STUDY OBJECTIVES

Prolonged wakefulness leads to a progressive increase in sleep pressure, reflected in a global increase in slow wave activity (SWA, 0.5-4.5 Hz) in the sleep electroencephalogram (EEG). A global increase in wake theta activity (5-9 Hz) also occurs. Recently, it was shown that prolonged wakefulness in rodents leads to signs of "local sleep" in an otherwise awake brain, accompanied by a slow/theta wave (2-6 Hz) in the local EEG that occurs at different times in different cortical areas. Compelling evidence in animals and humans also indicates that sleep is locally regulated by the amount of experience-dependent plasticity. Here, we asked whether the extended practice of tasks that involve specific brain circuits results in increased occurrence of local intermittent theta waves in the human EEG, above and beyond the global EEG changes previously described.

DESIGN

Participants recorded with high-density EEG completed 2 experiments during which they stayed awake ≥ 24 h practicing a language task (audiobook listening [AB]) or a visuomotor task (driving simulator [DS]).

SETTING

Sleep laboratory.

PATIENTS OR PARTICIPANTS

16 healthy participants (7 females).

INTERVENTIONS

Two extended wake periods.

MEASUREMENTS AND RESULTS

Both conditions resulted in global increases in resting wake EEG theta power at the end of 24 h of wake, accompanied by increased sleepiness. Moreover, wake theta power as well as the occurrence and amplitude of theta waves showed regional, task-dependent changes, increasing more over left frontal derivations in AB, and over posterior parietal regions in DS. These local changes in wake theta power correlated with similar local changes in sleep low frequencies including SWA.

CONCLUSIONS

Extended experience-dependent plasticity of specific circuits results in a local increase of the wake theta EEG power in those regions, followed by more intense sleep, as reflected by SWA, over the same areas.

The cell-type specific cortical microcircuit: relating structure and activity in a full-scale spiking network model

In the past decade, the cell-type specific connectivity and activity of local cortical networks have been characterized experimentally to some detail. In parallel, modeling has been established as a tool to relate network structure to activity dynamics. While available comprehensive connectivity maps ( Thomson, West, et al. 2002; Binzegger et al. 2004) have been used in various computational studies, prominent features of the simulated activity such as the spontaneous firing rates do not match the experimental findings. Here, we analyze the properties of these maps to compile an integrated connectivity map, which additionally incorporates insights on the specific selection of target types. Based on this integrated map, we build a full-scale spiking network model of the local cortical microcircuit. The simulated spontaneous activity is asynchronous irregular and cell-type specific firing rates are in agreement with in vivo recordings in awake animals, including the low rate of layer 2/3 excitatory cells. The interplay of excitation and inhibition captures the flow of activity through cortical layers after transient thalamic stimulation. In conclusion, the integration of a large body of the available connectivity data enables us to expose the dynamical consequences of the cortical microcircuitry.

Creating, documenting and sharing network models

As computational neuroscience matures, many simulation environments are available that are useful for neuronal network modeling. However, methods for successfully documenting models for publication and for exchanging models and model components among these projects are still under development. Here we briefly review existing software and applications for network model creation, documentation and exchange. Then we discuss a few of the larger issues facing the field of computational neuroscience regarding network modeling and suggest solutions to some of these problems, concentrating in particular on standardized network model terminology, notation, and descriptions and explicit documentation of model scaling. We hope this will enable and encourage computational neuroscientists to share their models more systematically in the future.

Synaptic depression and slow oscillatory activity in a biophysical network model of the cerebral cortex

Short-term synaptic depression (STD) is a form of synaptic plasticity that has a large impact on network computations. Experimental results suggest that STD is modulated by cortical activity, decreasing with activity in the network and increasing during silent states. Here, we explored different activity-modulation protocols in a biophysical network model for which the model displayed less STD when the network was active than when it was silent, in agreement with experimental results. Furthermore, we studied how trains of synaptic potentials had lesser decay during periods of activity (UP states) than during silent periods (DOWN states), providing new experimental predictions. We next tackled the inverse question of what is the impact of modifying STD parameters on the emergent activity of the network, a question difficult to answer experimentally. We found that synaptic depression of cortical connections had a critical role to determine the regime of rhythmic cortical activity. While low STD resulted in an emergent rhythmic activity with short UP states and long DOWN states, increasing STD resulted in longer and more frequent UP states interleaved with short silent periods. A still higher synaptic depression set the network into a non-oscillatory firing regime where DOWN states no longer occurred. The speed of propagation of UP states along the network was not found to be modulated by STD during the oscillatory regime; it remained relatively stable over a range of values of STD. Overall, we found that the mutual interactions between synaptic depression and ongoing network activity are critical to determine the mechanisms that modulate cortical emergent patterns.

Astrocyte regulation of sleep circuits: experimental and modeling perspectives

Integrated within neural circuits, astrocytes have recently been shown to modulate brain rhythms thought to mediate sleep function. Experimental evidence suggests that local impact of astrocytes on single synapses translates into global modulation of neuronal networks and behavior. We discuss these findings in the context of current conceptual models of sleep generation and function, each of which have historically focused on neural mechanisms. We highlight the implications and the challenges introduced by these results from a conceptual and computational perspective. We further provide modeling directions on how these data might extend our knowledge of astrocytic properties and sleep function. Given our evolving understanding of how local cellular activities during sleep lead to functional outcomes for the brain, further mechanistic and theoretical understanding of astrocytic contribution to these dynamics will undoubtedly be of great basic and translational benefit.

Diversity and noise effects in a model of homeostatic regulation of the sleep-wake cycle

Recent advances in sleep neurobiology have allowed development of physiologically based mathematical models of sleep regulation that account for the neuronal dynamics responsible for the regulation of sleep-wake cycles and allow detailed examination of the underlying mechanisms. Neuronal systems in general, and those involved in sleep regulation in particular, are noisy and heterogeneous by their nature. It has been shown in various systems that certain levels of noise and diversity can significantly improve signal encoding. However, these phenomena, especially the effects of diversity, are rarely considered in the models of sleep regulation. The present paper is focused on a neuron-based physiologically motivated model of sleep-wake cycles that proposes a novel mechanism of the homeostatic regulation of sleep based on the dynamics of a wake-promoting neuropeptide orexin. Here this model is generalized by the introduction of intrinsic diversity and noise in the orexin-producing neurons, in order to study the effect of their presence on the sleep-wake cycle. A simple quantitative measure of the quality of a sleep-wake cycle is introduced and used to systematically study the generalized model for different levels of noise and diversity. The model is shown to exhibit a clear diversity-induced resonance: that is, the best wake-sleep cycle turns out to correspond to an intermediate level of diversity at the synapses of the orexin-producing neurons. On the other hand, only a mild evidence of stochastic resonance is found, when the level of noise is varied. These results show that disorder, especially in the form of quenched diversity, can be a key-element for an efficient or optimal functioning of the homeostatic regulation of the sleep-wake cycle. Furthermore, this study provides an example of a constructive role of diversity in a neuronal system that can be extended beyond the system studied here.

Information coding in a laminar computational model of cat primary visual cortex

Neural populations across cortical layers perform different computational tasks. However, it is not known whether information in different layers is encoded using a common neural code or whether it depends on the specific layer. Here we studied the laminar distribution of information in a large-scale computational model of cat primary visual cortex. We analyzed the amount of information about the input stimulus conveyed by the different representations of the cortical responses. In particular, we compared the information encoded in four possible neural codes: (1) the information carried by the firing rate of individual neurons; (2) the information carried by spike patterns within a time window; (3) the rate-and-phase information carried by the firing rate labelled by the phase of the Local Field Potentials (LFP); (4) the pattern-and-phase information carried by the spike patterns tagged with the LFP phase. We found that there is substantially more information in the rate-and-phase code compared with the firing rate alone for low LFP frequency bands (less than 30 Hz). When comparing how information is encoded across layers, we found that the extra information contained in a rate-and-phase code may reach 90 % in Layer 4, while in other layers it reaches only 60 %, compared to the information carried by the firing rate alone. These results suggest that information processing in primary sensory cortices could rely on different coding strategies across different layers.

Long-term relationships between cholinergic tone, synchronous bursting and synaptic remodeling

Cholinergic neuromodulation plays key roles in the regulation of neuronal excitability, network activity, arousal, and behavior. On longer time scales, cholinergic systems play essential roles in cortical development, maturation, and plasticity. Presumably, these processes are associated with substantial synaptic remodeling, yet to date, long-term relationships between cholinergic tone and synaptic remodeling remain largely unknown. Here we used automated microscopy combined with multielectrode array recordings to study long-term relationships between cholinergic tone, excitatory synapse remodeling, and network activity characteristics in networks of cortical neurons grown on multielectrode array substrates. Experimental elevations of cholinergic tone led to the abrupt suppression of episodic synchronous bursting activity (but not of general activity), followed by a gradual growth of excitatory synapses over hours. Subsequent blockage of cholinergic receptors led to an immediate restoration of synchronous bursting and the gradual reversal of synaptic growth. Neither synaptic growth nor downsizing was governed by multiplicative scaling rules. Instead, these occurred in a subset of synapses, irrespective of initial synaptic size. Synaptic growth seemed to depend on intrinsic network activity, but not on the degree to which bursting was suppressed. Intriguingly, sustained elevations of cholinergic tone were associated with a gradual recovery of synchronous bursting but not with a reversal of synaptic growth. These findings show that cholinergic tone can strongly affect synaptic remodeling and synchronous bursting activity, but do not support a strict coupling between the two. Finally, the reemergence of synchronous bursting in the presence of elevated cholinergic tone indicates that the capacity of cholinergic neuromodulation to indefinitely suppress synchronous bursting might be inherently limited.

Control of sleep and wakefulness

This review summarizes the brain mechanisms controlling sleep and wakefulness. Wakefulness promoting systems cause low-voltage, fast activity in the electroencephalogram (EEG). Multiple interacting neurotransmitter systems in the brain stem, hypothalamus, and basal forebrain converge onto common effector systems in the thalamus and cortex. Sleep results from the inhibition of wake-promoting systems by homeostatic sleep factors such as adenosine and nitric oxide and GABAergic neurons in the preoptic area of the hypothalamus, resulting in large-amplitude, slow EEG oscillations. Local, activity-dependent factors modulate the amplitude and frequency of cortical slow oscillations. Non-rapid-eye-movement (NREM) sleep results in conservation of brain energy and facilitates memory consolidation through the modulation of synaptic weights. Rapid-eye-movement (REM) sleep results from the interaction of brain stem cholinergic, aminergic, and GABAergic neurons which control the activity of glutamatergic reticular formation neurons leading to REM sleep phenomena such as muscle atonia, REMs, dreaming, and cortical activation. Strong activation of limbic regions during REM sleep suggests a role in regulation of emotion. Genetic studies suggest that brain mechanisms controlling waking and NREM sleep are strongly conserved throughout evolution, underscoring their enormous importance for brain function. Sleep disruption interferes with the normal restorative functions of NREM and REM sleep, resulting in disruptions of breathing and cardiovascular function, changes in emotional reactivity, and cognitive impairments in attention, memory, and decision making.