Analysis of oscillatory patterns in the human sleep EEG using a novel detection algorithm

Transient oscillation dynamics during sleep provide a robust basis for electroencephalographic phenotyping and biomarker identification

Transient oscillatory events in the sleep electroencephalogram represent short-term coordinated network activity. Of particular importance, sleep spindles are transient oscillatory events associated with memory consolidation, which are altered in aging and in several psychiatric and neurodegenerative disorders. Spindle identification, however, currently contains implicit assumptions derived from what waveforms were historically easiest to discern by eye, and has recently been shown to select only a high-amplitude subset of transient events. Moreover, spindle activity is typically averaged across a sleep stage, collapsing continuous dynamics into discrete states. What information can be gained by expanding our view of transient oscillatory events and their dynamics? In this paper, we develop a novel approach to electroencephalographic phenotyping, characterizing a generalized class of transient time-frequency events across a wide frequency range using continuous dynamics. We demonstrate that the complex temporal evolution of transient events during sleep is highly stereotyped when viewed as a function of slow oscillation power (an objective, continuous metric of depth-of-sleep) and phase (a correlate of cortical up/down states). This two-fold power-phase representation has large intersubject variability-even within healthy controls-yet strong night-to-night stability for individuals, suggesting a robust basis for phenotyping. As a clinical application, we then analyze patients with schizophrenia, confirming established spindle (12-15 Hz) deficits as well as identifying novel differences in transient non-rapid eye movement events in low-alpha (7-10 Hz) and theta (4-6 Hz) ranges. Overall, these results offer an expanded view of transient activity, describing a broad class of events with properties varying continuously across spatial, temporal, and phase-coupling dimensions.

Sleep spindles comprise a subset of a broader class of electroencephalogram events

STUDY OBJECTIVES

Sleep spindles are defined based on expert observations of waveform features in the electroencephalogram (EEG) traces. This is a potentially limiting characterization, as transient oscillatory bursts like spindles are easily obscured in the time domain by higher amplitude activity at other frequencies or by noise. It is therefore highly plausible that many relevant events are missed by current approaches based on traditionally defined spindles. Given their oscillatory structure, we reexamine spindle activity from first principles, using time-frequency activity in comparison to scored spindles.

METHODS

Using multitaper spectral analysis, we observe clear time-frequency peaks in the sigma (10-16 Hz) range (TFσ peaks). While nearly every scored spindle coincides with a TFσ peak, numerous similar TFσ peaks remain undetected. We therefore perform statistical analyses of spindles and TFσ peaks using manual and automated detection methods, comparing event cooccurrence, morphological similarities, and night-to-night consistency across multiple datasets.

RESULTS

On average, TFσ peaks have more than three times the rate of spindles (mean rate: 9.8 vs. 3.1 events/minute). Moreover, spindles subsample the most prominent TFσ peaks with otherwise identical spectral morphology. We further demonstrate that detected TFσ peaks have stronger night-to-night rate stability (ρ = 0.98) than spindles (ρ = 0.67), while covarying with spindle rates across subjects (ρ = 0.72).

CONCLUSIONS

These results provide compelling evidence that traditionally defined spindles constitute a subset of a more generalized class of EEG events. TFσ peaks are therefore a more complete representation of the underlying phenomenon, providing a more consistent and robust basis for future experiments and analyses.

Estimation of Time-Varying Spectral Peaks and Decomposition of EEG Spectrograms

Detection of spectral peaks and estimation of their properties, including frequency and amplitude, are fundamental to many applications of signal processing. Electroencephalography (EEG) of sleep, in particular, displays characteristic oscillations that change continuously throughout the night. Capturing these dynamics is essential to understanding the sleep process and characterizing the heterogeneity observed across individuals. Most sleep EEG analyses rely on either time-averaged spectra or bandpassed amplitude/power. Unfortunately, these approaches obscure the time-variability of peak properties, require specification of a priori criteria, and cannot distinguish power from nearby oscillations. More sophisticated approaches, using various spectral models, have been proposed to better estimate oscillatory properties, but these too have limitations. We present an improved approach to spectrogram decomposition, tracking time-varying parameterized peak functions and dynamically estimating their parameters using a modified form of the iterated extended Kalman filter (IEKF) that incorporates discrete On/Off-switching of peak combinations and a sampling step to draw the initial reference trajectory. We evaluate this approach on two types of simulated examples-one nearly within the model class and one outside. We find excellent performance, in terms of spectral fits and accuracy of estimated states, for both simulation types. We then apply the approach to real EEG data of sleep onset, obtaining quality spectral estimates with estimated peak combinations closely matching the expert-scored sleep stages. This approach offers not only the ability to estimate time-varying parameters of spectral peaks but, moving forward, the potential to estimate the governing dynamics and analyze their variability across nights, subjects, and clinical groups.

Automatic detection of microsleep episodes with feature-based machine learning

STUDY OBJECTIVES

Microsleep episodes (MSEs) are brief episodes of sleep, mostly defined to be shorter than 15 s. In the electroencephalogram (EEG), MSEs are mainly characterized by a slowing in frequency. The identification of early signs of sleepiness and sleep (e.g. MSEs) is of considerable clinical and practical relevance. Under laboratory conditions, the maintenance of wakefulness test (MWT) is often used for assessing vigilance.

METHODS

We analyzed MWT recordings of 76 patients referred to the Sleep-Wake-Epilepsy-Center. MSEs were scored by experts defined by the occurrence of theta dominance on ≥1 occipital derivation lasting 1-15 s, whereas the eyes were at least 80% closed. We calculated spectrograms using an autoregressive model of order 16 of 1 s epochs moved in 200 ms steps in order to visualize oscillatory activity and derived seven features per derivation: power in delta, theta, alpha and beta bands, ratio theta/(alpha + beta), quantified eye movements, and median frequency. Three algorithms were used for MSE classification: support vector machine (SVM), random forest (RF), and an artificial neural network (long short-term memory [LSTM] network). Data of 53 patients were used for the training of the classifiers, and 23 for testing.

RESULTS

MSEs were identified with a high performance (sensitivity, specificity, precision, accuracy, and Cohen's kappa coefficient). Training revealed that delta power and the ratio theta/(alpha + beta) were most relevant features for the RF classifier and eye movements for the LSTM network.

CONCLUSIONS

The automatic detection of MSEs was successful for our EEG-based definition of MSEs, with good performance of all algorithms applied.

Sleep Spindle Detection Using RUSBoost and Synchrosqueezed Wavelet Transform

Sleep spindles are important electroencephalographic (EEG) waveforms in sleep medicine; however, it is burdensome even for experts to detect spindles, so automatic spindle detection methodologies have been investigated. Conventional methods utilize waveforms template matching or machine learning for detecting spindles. In the former approach, it is necessary to tune thresholds for individual adaptation, while the latter approach has the problem of imbalanced data because the amount of sleep spindles is small compared with the entire EEG data. The present work proposes a sleep spindle detection method that combines wavelet synchrosqueezed transform (SST) and random under-sampling boosting (RUSBoost). SST is a time-frequency analysis method suitable for extracting features of spindle waveforms. RUSBoost is a framework for coping with the imbalanced data problem. The proposed SST-RUS can deal with the imbalanced data in spindle detection and does not require threshold tuning because RUSBoost uses majority voting of weak classifiers for discrimination. The performance of SST-RUS was validated using an open-access database called the Montreal archives of sleep studies cohort 1 (MASS-C1), which showed an F-measure of 0.70 with a sensitivity of 76.9% and a positive predictive value of 61.2%. The proposed method can reduce the burden of PSG scoring.

Sleep Spindles: Mechanisms and Functions

Sleep spindles are burstlike signals in the electroencephalogram (EEG) of the sleeping mammalian brain and electrical surface correlates of neuronal oscillations in thalamus. As one of the most inheritable sleep EEG signatures, sleep spindles probably reflect the strength and malleability of thalamocortical circuits that underlie individual cognitive profiles. We review the characteristics, organization, regulation, and origins of sleep spindles and their implication in non-rapid-eye-movement sleep (NREMS) and its functions, focusing on human and rodent. Spatially, sleep spindle-related neuronal activity appears on scales ranging from small thalamic circuits to functional cortical areas, and generates a cortical state favoring intracortical plasticity while limiting cortical output. Temporally, sleep spindles are discrete events, part of a continuous power band, and elements grouped on an infraslow time scale over which NREMS alternates between continuity and fragility. We synthesize diverse and seemingly unlinked functions of sleep spindles for sleep architecture, sensory processing, synaptic plasticity, memory formation, and cognitive abilities into a unifying sleep spindle concept, according to which sleep spindles 1) generate neural conditions of large-scale functional connectivity and plasticity that outlast their appearance as discrete EEG events, 2) appear preferentially in thalamic circuits engaged in learning and attention-based experience during wakefulness, and 3) enable a selective reactivation and routing of wake-instated neuronal traces between brain areas such as hippocampus and cortex. Their fine spatiotemporal organization reflects NREMS as a physiological state coordinated over brain and body and may indicate, if not anticipate and ultimately differentiate, pathologies in sleep and neurodevelopmental, -degenerative, and -psychiatric conditions.

Oscillatory patterns in the electroencephalogram at sleep onset

Falling asleep is a gradually unfolding process. We investigated the role of various oscillatory activities including sleep spindles and alpha and delta oscillations at sleep onset (SO) by automatically detecting oscillatory events. We used two datasets of healthy young males, eight with four baseline recordings, and eight with a baseline and recovery sleep after 40 h of sustained wakefulness. We analyzed the 2-min interval before SO (stage 2) and the five consecutive 2-min intervals after SO. The incidence of delta/theta events reached its maximum in the first 2-min episode after SO, while the frequency of them was continuously decreasing from stage 1 onwards, continuing over SO and further into deeper sleep. Interestingly, this decrease of the frequencies of the oscillations were not affected by increased sleep pressure, in contrast to the incidence which increased. We observed an increasing number of alpha events after SO, predominantly frontally, with their prevalence varying strongly across individuals. Sleep spindles started to occur after SO, with first an increasing then a decreasing incidence and a continuous decrease in their frequency. Again, the frequency of the spindles was not altered after sleep deprivation. Oscillatory events revealed derivation dependent aspects. However, these regional aspects were not specific of the process of SO but rather reflect a general sleep related phenomenon. No individual traits of SO features (incidence and frequency of oscillations) and their dynamics were observed. Delta/theta events are important features for the analysis of SO in addition to slow waves.

Brain dynamics during the sleep onset transition: An EEG source localization study

EEG source localization is an essential tool to reveal the cortical sources underlying brain oscillatory activity. We applied LORETA, a technique of EEG source localization, to identify the principal brain areas involved in the process of falling asleep (sleep onset, SO). We localized the contributing brain areas of activity in the classical frequency bands and tracked their temporal evolution (in 2-min intervals from 2 min prior to SO up to 10 min after SO) during a baseline night and subsequent recovery sleep after total sleep deprivation of 40 h.

Delta activity (0.5–5 Hz) gradually increased both in baseline and recovery sleep, starting in frontal areas and finally involving the entire cortex. This increase was steeper in the recovery condition. The evolution of sigma activity (12–16 Hz) resembled an inverted U-shape in both conditions and the activity was most salient in the parietal cortex. In recovery, sigma activity reached its maximum faster than in baseline, but attained lower levels. Theta activity (5–8 Hz) increased with time in large parts of the occipital lobe (baseline and recovery) and in recovery involved additionally frontal areas. Changes in alpha activity (8–12 Hz) at sleep onset involved large areas of the cortex, whereas activity in the beta range (16–24 Hz) was restricted to small cortical areas. The dynamics in recovery could be considered as a “fast-forward version” of the one in baseline.

Our results confirm that the process of falling asleep is neither spatially nor temporally a uniform process and that different brain areas might be falling asleep at a different speed potentially reflecting use dependent aspects of sleep regulation.

•

LORETA is a valuable tool to reveal cortical sources of brain activity at sleep onset.

•

Spectral bands had location dependent dynamics; brain areas fell asleep asynchronously

•

BA 11 was the most relevant brain region associated with delta activity.

•

Spindle dynamics resembled an inverted U-shape.

•

During recovery from sleep deprivation capacity for spindle generation was reduced.

Discovery of key whole-brain transitions and dynamics during human wakefulness and non-REM sleep

The modern understanding of sleep is based on the classification of sleep into stages defined by their electroencephalography (EEG) signatures, but the underlying brain dynamics remain unclear. Here we aimed to move significantly beyond the current state-of-the-art description of sleep, and in particular to characterise the spatiotemporal complexity of whole-brain networks and state transitions during sleep. In order to obtain the most unbiased estimate of how whole-brain network states evolve through the human sleep cycle, we used a Markovian data-driven analysis of continuous neuroimaging data from 57 healthy participants falling asleep during simultaneous functional magnetic resonance imaging (fMRI) and EEG. This Hidden Markov Model (HMM) facilitated discovery of the dynamic choreography between different whole-brain networks across the wake-non-REM sleep cycle. Notably, our results reveal key trajectories to switch within and between EEG-based sleep stages, while highlighting the heterogeneities of stage N1 sleep and wakefulness before and after sleep.

Sleep is composed of a number of different stages, each associated with a different pattern of brain activity. Here, using a data-driven Hidden Markov Model (HMM) of fMRI data, the authors discover a more complex set of neural activity states underlying the conventional stages of non-REM sleep.

Absent sleep EEG spindle activity in GluA1 (Gria1) knockout mice: relevance to neuropsychiatric disorders

Sleep EEG spindles have been implicated in attention, sensory processing, synaptic plasticity and memory consolidation. In humans, deficits in sleep spindles have been reported in a wide range of neurological and psychiatric disorders, including schizophrenia. Genome-wide association studies have suggested a link between schizophrenia and genes associated with synaptic plasticity, including the Gria1 gene which codes for the GluA1 subunit of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor. Gria1^-/- mice exhibit a phenotype relevant for neuropsychiatric disorders, including reduced synaptic plasticity and, at the behavioural level, attentional deficits leading to aberrant salience. In this study we report a striking reduction of EEG power density including the spindle-frequency range (10-15 Hz) during sleep in Gria1^-/- mice. The reduction of spindle-activity in Gria1^-/- mice was accompanied by longer REM sleep episodes, increased EEG slow-wave activity in the occipital derivation during baseline sleep, and a reduced rate of decline of EEG slow wave activity (0.5-4 Hz) during NREM sleep after sleep deprivation. These data provide a novel link between glutamatergic dysfunction and sleep abnormalities in a schizophrenia-relevant mouse model.

Sleep- and Wake-Like States in Small Networks In Vivo and In Vitro

Wakefulness and sleep are highly complex and heterogeneous processes, involving multiple neurotransmitter systems and a sophisticated interplay between global and local networks of neurons and non-neuronal cells. Macroscopic approaches applied at the level of the whole organism, view sleep as a global behaviour and allow for investigation into aspects such as the effects of insufficient or disrupted sleep on cognitive function, metabolism, thermoregulation and sensory processing. While significant progress has been achieved using such large-scale approaches, the inherent complexity of sleep-wake regulation has necessitated the development of methods which tackle specific aspects of sleep in isolation. One way this may be achieved is by investigating specific cellular or molecular phenomena in the whole organism in situ, either during spontaneous or induced sleep-wake states. This approach has greatly advanced our knowledge about the electrophysiology and pharmacology of ion channels, specific receptors, intracellular pathways and the small networks implicated in the control and regulation of the sleep-wake cycle. Importantly though, there are a variety of external and internal factors that influence global behavioural states which are difficult to control for using these approaches. For this reason, over the last few decades, ex vivo experimental models have become increasingly popular and have greatly advanced our understanding of many fundamental aspects of sleep, including the neuroanatomy and neurochemistry of sleep states, sleep regulation, the origin and dynamics of specific sleep oscillations, network homeostasis as well as the functional roles of sleep. This chapter will focus on the use of small neuronal networks as experimental models and will highlight the most significant and novel insights these approaches have provided.

Developmental Changes in Sleep Oscillations during Early Childhood

Although quantitative analysis of the sleep electroencephalogram (EEG) has uncovered important aspects of brain activity during sleep in adolescents and adults, similar findings from preschool-age children remain scarce. This study utilized our time-frequency method to examine sleep oscillations as characteristic features of human sleep EEG. Data were collected from a longitudinal sample of young children (n = 8; 3 males) at ages 2, 3, and 5 years. Following sleep stage scoring, we detected and characterized oscillatory events across age and examined how their features corresponded to spectral changes in the sleep EEG. Results indicated a developmental decrease in the incidence of delta and theta oscillations. Spindle oscillations, however, were almost absent at 2 years but pronounced at 5 years. All oscillatory event changes were stronger during light sleep than slow-wave sleep. Large interindividual differences in sleep oscillations and their characteristics (e.g., "ultrafast" spindle-like oscillations, theta oscillation incidence/frequency) also existed. Changes in delta and spindle oscillations across early childhood may indicate early maturation of the thalamocortical system. Our analytic approach holds promise for revealing novel types of sleep oscillatory events that are specific to periods of rapid normal development across the lifespan and during other times of aberrant changes in neurobehavioral function.

A study of problems encountered in Granger causality analysis from a neuroscience perspective

Granger causality methods were developed to analyze the flow of information between time series. These methods have become more widely applied in neuroscience. Frequency-domain causality measures, such as those of Geweke, as well as multivariate methods, have particular appeal in neuroscience due to the prevalence of oscillatory phenomena and highly multivariate experimental recordings. Despite its widespread application in many fields, there are ongoing concerns regarding the applicability of Granger causality methods in neuroscience. When are these methods appropriate? How reliably do they recover the system structure underlying the observed data? What do frequency-domain causality measures tell us about the functional properties of oscillatory neural systems? In this paper, we analyze fundamental properties of Granger-Geweke (GG) causality, both computational and conceptual. Specifically, we show that (i) GG causality estimates can be either severely biased or of high variance, both leading to spurious results; (ii) even if estimated correctly, GG causality estimates alone are not interpretable without examining the component behaviors of the system model; and (iii) GG causality ignores critical components of a system's dynamics. Based on this analysis, we find that the notion of causality quantified is incompatible with the objectives of many neuroscience investigations, leading to highly counterintuitive and potentially misleading results. Through the analysis of these problems, we provide important conceptual clarification of GG causality, with implications for other related causality approaches and for the role of causality analyses in neuroscience as a whole.

Wake-sleep transition as a noisy bifurcation

A recent physiologically based model of the ascending arousal system is used to analyze the dynamics near the transition from wake to sleep, which corresponds to a saddle-node bifurcation at a critical point. A normal form is derived by approximating the dynamics by those of a particle in a parabolic potential well with dissipation. This mechanical analog is used to calculate the power spectrum of fluctuations in response to a white noise drive, and the scalings of fluctuation variance and spectral width are derived versus distance from the critical point. The predicted scalings are quantitatively confirmed by numerical simulations, which show that the variance increases and the spectrum undergoes critical slowing, both in accord with theory. These signals can thus serve as potential precursors to indicate imminent wake-sleep transition, with potential application to safety-critical occupations in transport, air-traffic control, medicine, and heavy industry.

Sleep Spindles as an Electrographic Element: Description and Automatic Detection Methods

Sleep spindle is a peculiar oscillatory brain pattern which has been associated with a number of sleep (isolation from exteroceptive stimuli, memory consolidation) and individual characteristics (intellectual quotient). Oddly enough, the definition of a spindle is both incomplete and restrictive. In consequence, there is no consensus about how to detect spindles. Visual scoring is cumbersome and user dependent. To analyze spindle activity in a more robust way, automatic sleep spindle detection methods are essential. Various algorithms were developed, depending on individual research interest, which hampers direct comparisons and meta-analyses. In this review, sleep spindle is first defined physically and topographically. From this general description, we tentatively extract the main characteristics to be detected and analyzed. A nonexhaustive list of automatic spindle detection methods is provided along with a description of their main processing principles. Finally, we propose a technique to assess the detection methods in a robust and comparable way.

Using a quadratic parameter sinusoid model to characterize the structure of EEG sleep spindles

Sleep spindles are essentially non-stationary signals that display time and frequency-varying characteristics within their envelope, which makes it difficult to accurately identify its instantaneous frequency and amplitude. To allow a better parameterization of the structure of spindle, we propose modeling spindles using a Quadratic Parameter Sinusoid (QPS). The QPS is well suited to model spindle activity as it utilizes a quadratic representation to capture the inherent duration and frequency variations within spindles. The effectiveness of our proposed model and estimation technique was quantitatively evaluated in parameter determination experiments using simulated spindle-like signals and real spindles in the presence of background EEG. We used the QPS parameters to predict the energy and frequency of spindles with a mean accuracy of 92.34 and 97.73% respectively. We also show that the QPS parameters provide a quantification of the amplitude and frequency variations occurring within sleep spindles that can be observed visually and related to their characteristic "waxing and waning" shape. We analyze the variations in the parameters values to present how they can be used to understand the inter- and intra-participant variations in spindle structure. Finally, we present a comparison of the QPS parameters of spindles and non-spindles, which shows a substantial difference in parameter values between the two classes.

Combining time-frequency and spatial information for the detection of sleep spindles

EEG sleep spindles are short (0.5-2.0 s) bursts of activity in the 11-16 Hz band occurring during non-rapid eye movement (NREM) sleep. This sporadic activity is thought to play a role in memory consolidation, brain plasticity, and protection of sleep integrity. Many automatic detectors have been proposed to assist or replace experts for sleep spindle scoring. However, these algorithms usually detect too many events making it difficult to achieve a good tradeoff between sensitivity (Se) and false detection rate (FDr). In this work, we propose a semi-automatic detector comprising a sensitivity phase based on well-established criteria followed by a specificity phase using spatial and spectral criteria. In the sensitivity phase, selected events are those which amplitude in the 10-16 Hz band and spectral ratio characteristics both reject a null hypothesis (p < 0.1) stating that the considered event is not a spindle. This null hypothesis is constructed from events occurring during rapid eye movement (REM) sleep epochs. In the specificity phase, a hierarchical clustering of the selected candidates is done based on events' frequency and spatial position along the anterior-posterior axis. Only events from the classes grouping most (at least 80%) spindles scored by an expert are kept. We obtain Se = 93.2% and FDr = 93.0% in the first phase and Se = 85.4% and FDr = 86.2% in the second phase. For these two phases, Matthew's correlation coefficients are respectively 0.228 and 0.324. Results suggest that spindles are defined by specific spatio-spectral properties and that automatic detection methods can be improved by considering these features.

NREM and REM Sleep: Complementary Roles in Recovery after Wakefulness

The overall function of sleep is hypothesized to provide "recovery" after preceding waking activities, thereby ensuring optimal functioning during subsequent wakefulness. However, the functional significance of the temporal dynamics of sleep, manifested in the slow homeostatic process and the alternation between non-rapid eye movement (NREM) and REM sleep remains unclear. We propose that NREM and REM sleep have distinct and complementary contributions to the overall function of sleep. Specifically, we suggest that cortical slow oscillations, occurring within specific functionally interconnected neuronal networks during NREM sleep, enable information processing, synaptic plasticity, and prophylactic cellular maintenance ("recovery process"). In turn, periodic excursions into an activated brain state-REM sleep-appear to be ideally placed to perform "selection" of brain networks, which have benefited from the process of "recovery," based on their offline performance. Such two-stage modus operandi of the sleep process would ensure that its functions are fulfilled according to the current need and in the shortest time possible. Our hypothesis accounts for the overall architecture of normal sleep and opens up new perspectives for understanding pathological conditions associated with abnormal sleep patterns.

Effect of prolonged wakefulness on electroencephalographic oscillatory activity during sleep

The human sleep electroencephalogram (EEG) is characterized by the occurrence of distinct oscillatory events such as delta waves, sleep spindles and alpha activity. We applied a previously proposed algorithm for the detection of such events and investigated their incidence and frequency in baseline and recovery sleep after 40 h of sustained wakefulness in 27 healthy young subjects. The changes in oscillatory events induced by sleep deprivation were compared to the corresponding spectral changes. Both approaches revealed, on average, an increase in low frequency activity and a decrease in spindle activity after sleep deprivation. However, the increase of oscillatory events in the delta range and decrease in the sigma range occurred in a more restricted frequency range compared to spectral changes. The mean relative power spectra showed a significant increase in theta and alpha activity after sleep deprivation while, on average, the event analysis showed only a weak effect in the theta band. The reason for this discrepancy is that the spectral analysis does not distinguish between diffuse activity and clearly visible temporally localized oscillations, while the event analysis would detect only the latter. Additionally, only a few individuals clearly showed activity in the theta or alpha frequency bands. Conversely, event analysis revealed that some individuals showed an increased rate of sleep spindles after sleep deprivation, a fact that was not evident in the relative power spectra due to a decrease in background activity. The two methods complement each other and facilitate the interpretation of distinct changes induced by prolonged wakefulness in sleep EEG.

Detecting alpha spindle events in EEG time series using adaptive autoregressive models

Background

Rhythmic oscillatory activity is widely observed during a variety of subject behaviors and is believed to play a central role in information processing and control. A classic example of rhythmic activity is alpha spindles, which consist of short (0.5-2 s) bursts of high frequency alpha activity. Recent research has shown that alpha spindles in the parietal/occipital area are statistically related to fatigue and drowsiness. These spindles constitute sharp changes in the underlying statistical properties of the signal. Our hypothesis is that change point detection models can be used to identify the onset and duration of spindles in EEG. In this work we develop an algorithm that accurately identifies sudden bursts of narrowband oscillatory activity in EEG using techniques derived from change point analysis. Our motivating example is detection of alpha spindles in the parietal/occipital areas of the brain. Our goal is to develop an algorithm that can be applied to any type of rhythmic oscillatory activity of interest for accurate online detection.

Methods

In this work we propose modeling the alpha band EEG time series using discounted autoregressive (DAR) modeling. The DAR model uses a discounting rate to weigh points measured further in the past less heavily than points more recently observed. This model is used together with predictive loss scoring to identify periods of EEG data that are statistically significant.

Results

Our algorithm accurately captures changes in the statistical properties of the alpha frequency band. These statistical changes are highly correlated with alpha spindle occurrences and form a reliable measure for detecting alpha spindles in EEG. We achieve approximately 95% accuracy in detecting alpha spindles, with timing precision to within approximately 150 ms, for two datasets from an experiment of prolonged simulated driving, as well as in simulated EEG. Sensitivity and specificity values are above 0.9, and in many cases are above 0.95, for our analysis.

Conclusion

Modeling the alpha band EEG using discounted AR models provides an efficient method for detecting oscillatory alpha activity in EEG. The method is based on statistical principles and can generally be applied to detect rhythmic activity in any frequency band or brain region.

Extracting more information from EEG recordings for a better description of sleep

We are introducing and validating an EEG data-based model of the sleep process with an arbitrary number of different sleep stages and a high time resolution allowing modeling of sleep microstructure. In contrast to the standard practice of sleep staging, defined by scoring rules, we describe sleep via posterior probabilities of a finite number of states, not necessarily reflecting the traditional sleep stages. To test the proposed probabilistic sleep model (PSM) for validity, we correlate statistics derived from the state posteriors with the results of psychometric tests, physiological variables and questionnaires collected before and after sleep. Considering short, in this study 3s long, data window the PSM allows describing the sleep process on finer time scale in comparison to the traditional sleep staging based on 20 or 30s long data segments visual inspection. By combining sleep states and using two measures derived from the posterior curves we show that the average absolute correlations between the measures and subjective and objective sleep quality measures are considerably higher when compared with the analogous measures derived from hypnograms based on sleep staging. In most cases these differences are significant. The results obtained with the PSM support its wider use in sleep process modeling research and these results also suggest that EEG signals contain more information about sleep than what sleep profiles based on discrete stages can reveal. Therefore the standardized scoring of sleep may not be sufficient to reveal important sleep changes related to subjective and objective sleep quality indexes. The proposed PSM represents a promising alternative.

Evolution of genuine cross-correlation strength of focal onset seizures

To quantify the evolution of genuine zero-lag cross-correlations of focal onset seizures, we apply a recently introduced multivariate measure to broad band and to narrow-band EEG data. For frequency components below 12.5 Hz, the strength of genuine cross-correlations decreases significantly during the seizure and the immediate postseizure period, while higher frequency bands show a tendency of elevated cross-correlations during the same period. We conclude that in terms of genuine zero-lag cross-correlations, the electrical brain activity as assessed by scalp electrodes shows a significant spatial fragmentation, which might promote seizure offset.

Increased Alpha (8–12 Hz) Activity during Slow Wave Sleep as a Marker for the Transition from Implicit Knowledge to Explicit Insight

The number reduction task (NRT) allows us to study the transition from implicit knowledge of hidden task regularities to explicit insight into these regularities. To identify sleep-associated neurophysiological indicators of this restructuring of knowledge representations, we measured frequency-specific power of EEG while participants slept during the night between two sessions of the NRT. Alpha (8–12 Hz) EEG power during slow wave sleep (SWS) emerged as a specific marker of the transformation of presleep implicit knowledge to postsleep explicit knowledge (ExK). Beta power during SWS was increased whenever ExK was attained after sleep, irrespective of presleep knowledge. No such EEG predictors of insight were found during Sleep Stage 2 and rapid eye movement sleep. These results support the view that it is neuronal memory reprocessing during sleep, in particular during SWS, that lays the foundations for restructuring those task-related representations in the brain that are necessary for promoting the gain of ExK.

The multiple time scales of sleep dynamics as a challenge for modelling the sleeping brain

A particular property of the sleeping brain is that it exhibits dynamics on very different time scales ranging from the typical sleep oscillations such as sleep spindles and slow waves that can be observed in electroencephalogram (EEG) segments of several seconds duration over the transitions between the different sleep stages on a time scale of minutes to the dynamical processes involved in sleep regulation with typical time constants in the range of hours. There is an increasing body of work on mathematical and computational models addressing these different dynamics, however, usually considering only processes on a single time scale. In this paper, we review and present a new analysis of the dynamics of human sleep EEG at the different time scales and relate the findings to recent modelling efforts pointing out both the achievements and remaining challenges.

Validating an automated sleep spindle detection algorithm using an individualized approach

The goal of the current investigation was to develop a systematic method to validate the accuracy of an automated method of sleep spindle detection that takes into consideration individual differences in spindle amplitude. The benchmarking approach used here could be employed more generally to validate automated spindle scoring from other detection algorithms. In a sample of Stage 2 sleep from 10 healthy young subjects, spindles were identified both manually and automatically. The minimum amplitude threshold used by the Prana (PhiTools, Strasbourg, France) software spindle detection algorithm to identify a spindle was subject-specific and determined based upon each subject's mean peak spindle amplitude. Overall sensitivity and specificity values were 98.96 and 88.49%, respectively, when compared to manual scoring. Selecting individual amplitude thresholds for spindle detection based on systematic benchmarking data may validate automated spindle detection methods and improve reproducibility of experimental results. Given that interindividual differences are accounted for, we feel that automatic spindle detection provides an accurate and efficient alternative approach for detecting sleep spindles.

Functional connectivity between motor cortex and globus pallidus in human non-REM sleep

Recent evidence suggests that the motor system undergoes very specific modulation in its functional state during the different sleep stages. Here we test the hypothesis that changes in the functional organization of the motor system involve both cortical and subcortical levels and that these distributed changes are interrelated in defined frequency bands. To this end we evaluated functional connectivity between motor and non-motor cortical sites (fronto-central, parieto-occipital) and the globus pallidus (GP) in human non-REM sleep in seven patients undergoing deep brain stimulation (DBS) for dystonia using a variety of spectral measures (power, coherence, partial coherence and directed transfer function (DTF)). We found significant coherence between GP and fronto-central cortex as well as between GP and parieto-occipital cortex in circumscribed frequency bands that correlated with sleep specific oscillations in 'light sleep' (N2) and 'slow-wave sleep' (N3). These sleep specific oscillations were also reflected in significant coherence between the two cortical sites corroborating previous studies. Importantly, we found two different physiological activities represented within the broad band of significant coherence between 9.5 and 17 Hz. One component occurred in the frequency range of sleep spindles (12.5-17 Hz) and was maximal in the coherence between fronto-central and parieto-occipital cortex as well as between GP and both cortical sites during N2. This component was still present between fronto-central and parieto-occipital cortex in N3. Functional connectivity in this frequency band may be due to a common input to both GP and cortex. The second component consisted of a spectral peak over 9.5-12.5 Hz. Coherence was elevated in this band for all topographical constellations in both N2 and N3, but especially between GP and fronto-central cortex. The DTF suggested that the 9.5-12.5 Hz activity consisted of a preferential drive from GP to the fronto-central cortex in N2, whereas in N3 the DTF between GP and fronto-central cortex was symmetrical. Partial coherence supported distinctive patterns for the 9.5-12.5 and 12.5 and 17 Hz component, so that only coherence in the 9.5-12.5 Hz band was reduced when the effects of GP were removed from the coherence between the two cortical sites. The data suggest that activities in the GP and fronto-central cortex are functionally connected over 9.5-12.5 Hz, possibly as a specific signature of the motor system in human non-REM sleep. This finding is pertinent to the longstanding debate about the nature of alpha-delta sleep as a physiological or pathological feature of non-REM sleep.

Analysis of the temporal organization of sleep spindles in the human sleep EEG using a phenomenological modeling approach

The sleep electroencephalogram (EEG) is characterized by typical oscillatory patterns such as sleep spindles and slow waves. Recently, we proposed a method to detect and analyze these patterns using linear autoregressive models for short (approximately 1 s) data segments. We analyzed the temporal organization of sleep spindles and discuss to what extent the observed interevent intervals correspond to properties of stationary stochastic processes and whether additional slow processes, such as slow oscillations, have to be assumed. We have found evidence for such an additional slow process, most pronounced in sleep stage 2.

Detection of Active and Silent States in Neocortical Neurons from the Field Potential Signal during Slow-Wave Sleep

Oscillations of the local field potentials (LFPs) or electroencephalogram (EEG) at frequencies below 1 Hz are a hallmark of the slow-wave sleep. However, the timing of the underlying cellular events, which is an alternation of active and silent states of thalamocortical network, can be assessed only approximately from the phase of slow waves. Is it possible to detect, using the LFP or EEG, the timing of each episode of cellular activity or silence? With simultaneous recordings of the LFP and intracellular activity of 2-3 neocortical cells, we show that high-gamma-range (20-100 Hz) components in the LFP have significantly higher power when cortical cells are in active states as compared with silent-state periods. Exploiting this difference we have developed a new method, which uses the LFP signal to detect episodes of activity and silence of neocortical neurons. The method allows robust, reliable, and precise detection of timing of each episode of activity and silence of the neocortical network. It works with both surface and depth EEG, and its performance is affected little by the EEG prefiltering during recording. These results open new perspectives for studying differential operation of neural networks during periods of activity and silence, which rapidly alternate on the subsecond scale.