Large-Scale Cortical Dynamics of Sleep Slow Waves

Slow waves constitute the main signature of sleep in the electroencephalogram (EEG). They reflect alternating periods of neuronal hyperpolarization and depolarization in cortical networks. While recent findings have demonstrated their functional role in shaping and strengthening neuronal networks, a large-scale characterization of these two processes remains elusive in the human brain. In this study, by using simultaneous scalp EEG and intracranial recordings in 10 epileptic subjects, we examined the dynamics of hyperpolarization and depolarization waves over a large extent of the human cortex. We report that both hyperpolarization and depolarization processes can occur with two different characteristic time durations which are consistent across all subjects. For both hyperpolarization and depolarization waves, their average speed over the cortex was estimated to be approximately 1 m/s. Finally, we characterized their propagation pathways by studying the preferential trajectories between most involved intracranial contacts. For both waves, although single events could begin in almost all investigated sites across the entire cortex, we found that the majority of the preferential starting locations were located in frontal regions of the brain while they had a tendency to end in posterior and temporal regions

Robust Off- and Online Separation of Intracellularly Recorded Up and Down Cortical States

BACKGROUND: The neuronal cortical network generates slow (<1 Hz) spontaneous rhythmic activity that emerges from the recurrent connectivity. This activity occurs during slow wave sleep or anesthesia and also in cortical slices, consisting of alternating up (active, depolarized) and down (silent, hyperpolarized) states. The search for the underlying mechanisms and the possibility of analyzing network dynamics in vitro has been subject of numerous studies. This exposes the need for a detailed quantitative analysis of the membrane fluctuating behavior and computerized tools to automatically characterize the occurrence of up and down states. METHODOLOGY/PRINCIPAL FINDINGS: Intracellular recordings from different areas of the cerebral cortex were obtained from both in vitro and in vivo preparations during slow oscillations. A method that separates up and down states recorded intracellularly is defined and analyzed here. The method exploits the crossover of moving averages, such that transitions between up and down membrane regimes can be anticipated based on recent and past voltage dynamics. We demonstrate experimentally the utility and performance of this method both offline and online, the online use allowing to trigger stimulation or other events in the desired period of the rhythm. This technique is compared with a histogram-based approach that separates the states by establishing one or two discriminating membrane potential levels. The robustness of the method presented here is tested on data that departs from highly regular alternating up and down states. CONCLUSIONS/SIGNIFICANCE: We define a simple method to detect cortical states that can be applied in real time for offline processing of large amounts of recorded data on conventional computers. Also, the online detection of up and down states will facilitate the study of cortical dynamics. An open-source MATLAB toolbox, and Spike 2-compatible version are made freely available

The Neuronal Transition Probability (NTP) Model for the Dynamic Progression of Non-REM Sleep EEG: The Role of the Suprachiasmatic Nucleus

Little attention has gone into linking to its neuronal substrates the dynamic structure of non-rapid-eye-movement (NREM) sleep, defined as the pattern of time-course power in all frequency bands across an entire episode. Using the spectral power time-courses in the sleep electroencephalogram (EEG), we showed in the typical first episode, several moves towards-and-away from deep sleep, each having an identical pattern linking the major frequency bands beta, sigma and delta. The neuronal transition probability model (NTP) – in fitting the data well – successfully explained the pattern as resulting from stochastic transitions of the firing-rates of the thalamically-projecting brainstem-activating neurons, alternating between two steady dynamic-states (towards-and-away from deep sleep) each initiated by a so-far unidentified flip-flop. The aims here are to identify this flip-flop and to demonstrate that the model fits well all NREM episodes, not just the first. Using published data on suprachiasmatic nucleus (SCN) activity we show that the SCN has the information required to provide a threshold-triggered flip-flop for timing the towards-and-away alternations, information provided by sleep-relevant feedback to the SCN. NTP then determines the pattern of spectral power within each dynamic-state. NTP was fitted to individual NREM episodes 1–4, using data from 30 healthy subjects aged 20–30 years, and the quality of fit for each NREM measured. We show that the model fits well all NREM episodes and the best-fit probability-set is found to be effectively the same in fitting all subject data. The significant model-data agreement, the constant probability parameter and the proposed role of the SCN add considerable strength to the model. With it we link for the first time findings at cellular level and detailed time-course data at EEG level, to give a coherent picture of NREM dynamics over the entire night and over hierarchic brain levels all the way from the SCN to the EEG

Cortical Layer 1 and Layer 2/3 Astrocytes Exhibit Distinct Calcium Dynamics In Vivo

Cumulative evidence supports bidirectional interactions between astrocytes and neurons, suggesting glial involvement of neuronal information processing in the brain. Cytosolic calcium (Ca2+) concentration is important for astrocytes as Ca2+ surges co-occur with gliotransmission and neurotransmitter reception. Cerebral cortex is organized in layers which are characterized by distinct cytoarchitecture. We asked if astrocyte-dominant layer 1 (L1) of the somatosensory cortex was different from layer 2/3 (L2/3) in spontaneous astrocytic Ca2+ activity and if it was influenced by background neural activity. Using a two-photon laser scanning microscope, we compared spontaneous Ca2+ activity of astrocytic somata and processes in L1 and L2/3 of anesthetized mature rat somatosensory cortex. We also assessed the contribution of background neural activity to the spontaneous astrocytic Ca2+ dynamics by investigating two distinct EEG states (“synchronized” vs. “de-synchronized” states). We found that astrocytes in L1 had nearly twice higher Ca2+ activity than L2/3. Furthermore, Ca2+ fluctuations of processes within an astrocyte were independent in L1 while those in L2/3 were synchronous. Pharmacological blockades of metabotropic receptors for glutamate, ATP, and acetylcholine, as well as suppression of action potentials did not have a significant effect on the spontaneous somatic Ca2+ activity. These results suggest that spontaneous astrocytic Ca2+ surges occurred in large part intrinsically, rather than neural activity-driven. Our findings propose a new functional segregation of layer 1 and 2/3 that is defined by autonomous astrocytic activity

Characterization of K-Complexes and Slow Wave Activity in a Neural Mass Model

NREM sleep is characterized by two hallmarks, namely K-complexes (KCs) during sleep stage N2 and cortical slow oscillations (SOs) during sleep stage N3. While the underlying dynamics on the neuronal level is well known and can be easily measured, the resulting behavior on the macroscopic population level remains unclear. On the basis of an extended neural mass model of the cortex, we suggest a new interpretation of the mechanisms responsible for the generation of KCs and SOs. As the cortex transitions from wake to deep sleep, in our model it approaches an oscillatory regime via a Hopf bifurcation. Importantly, there is a canard phenomenon arising from a homoclinic bifurcation, whose orbit determines the shape of large amplitude SOs. A KC corresponds to a single excursion along the homoclinic orbit, while SOs are noise-driven oscillations around a stable focus. The model generates both time series and spectra that strikingly resemble real electroencephalogram data and points out possible differences between the different stages of natural sleep

Activity of cortical and thalamic neurons during the slow (<1 Hz) rhythm in the mouse in vivo

During NREM sleep and under certain types of anaesthesia, the mammalian brain exhibits a distinctive slow (<1 Hz) rhythm. At the cellular level, this rhythm correlates with so-called UP and DOWN membrane potential states. In the neocortex, these UP and DOWN states correspond to periods of intense network activity and widespread neuronal silence, respectively, whereas in thalamocortical (TC) neurons, UP/DOWN states take on a more stereotypical oscillatory form, with UP states commencing with a low-threshold Ca2+ potential (LTCP). Whilst these properties are now well recognised for neurons in cats and rats, whether or not they are also shared by neurons in the mouse is not fully known. To address this issue, we obtained intracellular recordings from neocortical and TC neurons during the slow (<1 Hz) rhythm in anaesthetised mice. We show that UP/DOWN states in this species are broadly similar to those observed in cats and rats, with UP states in neocortical neurons being characterised by a combination of action potential output and intense synaptic activity, whereas UP states in TC neurons always commence with an LTCP. In some neocortical and TC neurons, we observed ‘spikelets’ during UP states, supporting the possible presence of electrical coupling. Lastly, we show that, upon tonic depolarisation, UP/DOWN states in TC neurons are replaced by rhythmic high-threshold bursting at ~5 Hz, as predicted by in vitro studies. Thus, UP/DOWN state generation appears to be an elemental and conserved process in mammals that underlies the slow (<1 Hz) rhythm in several species, including humans

Morning and Evening-Type Differences in Slow Waves during NREM Sleep Reveal Both Trait and State-Dependent Phenotypes

Brain recovery after prolonged wakefulness is characterized by increased density, amplitude and slope of slow waves (SW, <4 Hz) during non-rapid eye movement (NREM) sleep. These SW comprise a negative phase, during which cortical neurons are mostly silent, and a positive phase, in which most neurons fire intensively. Previous work showed, using EEG spectral analysis as an index of cortical synchrony, that Morning-types (M-types) present faster dynamics of sleep pressure than Evening-types (E-types). We thus hypothesized that single SW properties will also show larger changes in M-types than in E-types in response to increased sleep pressure. SW density (number per minute) and characteristics (amplitude, slope between negative and positive peaks, frequency and duration of negative and positive phases) were compared between chronotypes for a baseline sleep episode (BL) and for recovery sleep (REC) after two nights of sleep fragmentation. While SW density did not differ between chronotypes, M-types showed higher SW amplitude and steeper slope than E-types, especially during REC. SW properties were also averaged for 3 NREM sleep periods selected for their decreasing level of sleep pressure (first cycle of REC [REC1], first cycle of BL [BL1] and fourth cycle of BL [BL4]). Slope was significantly steeper in M-types than in E-types in REC1 and BL1. SW frequency was consistently higher and duration of positive and negative phases constantly shorter in M-types than in E-types. Our data reveal that specific properties of cortical synchrony during sleep differ between M-types and E-types, although chronotypes show a similar capacity to generate SW. These differences may involve 1) stable trait characteristics independent of sleep pressure (i.e., frequency and durations) likely linked to the length of silent and burst-firing phases of individual neurons, and 2) specific responses to increased sleep pressure (i.e., slope and amplitude) expected to depend on the synchrony between neurons

Extensive astrocyte synchronization advances neuronal coupling in slow wave activity in vivo

Slow wave activity (SWA) is a characteristic brain oscillation in sleep and quiet wakefulness. Although the cell types contributing to SWA genesis are not yet identified, the principal role of neurons in the emergence of this essential cognitive mechanism has not been questioned. To address the possibility of astrocytic involvement in SWA, we used a transgenic rat line expressing a calcium sensitive fluorescent protein in both astrocytes and interneurons and simultaneously imaged astrocytic and neuronal activity in vivo. Here we demonstrate, for the first time, that the astrocyte network display synchronized recurrent activity in vivo coupled to UP states measured by field recording and neuronal calcium imaging. Furthermore, we present evidence that extensive synchronization of the astrocytic network precedes the spatial build-up of neuronal synchronization. The earlier extensive recruitment of astrocytes in the synchronized activity is reinforced by the observation that neurons surrounded by active astrocytes are more likely to join SWA, suggesting causality. Further supporting this notion, we demonstrate that blockade of astrocytic gap junctional communication or inhibition of astrocytic Ca2+ transients reduces the ratio of both astrocytes and neurons involved in SWA. These in vivo findings conclusively suggest a causal role of the astrocytic syncytium in SWA generation