Sleep and circadian rhythms modeling: From hypothalamic regulatory networks to cortical dynamics and behavior

Sleep and circadian rhythms are regulated by dynamic physiologic processes that operate across multiple spatial and temporal scales. These include, but are not limited to, genetic oscillators, clearance of waste products from the brain, dynamic interplay among brain regions, and propagation of local dynamics across the cortex. The combination of these processes, modulated by environmental cues, such as light-dark cycles and work schedules, represents a complex multiscale system that regulates sleep-wake cycles and brain dynamics. Physiology-based mathematical models have successfully explained the mechanisms underpinning dynamics at specific scales and are a useful tool to investigate interactions across multiple scales. They can help answer questions such as how do electroencephalographic (EEG) features relate to subthalamic neuron activity? Or how are local cortical dynamics regulated by the homeostatic and circadian mechanisms? In this chapter, we review two types of models that are well-positioned to consider such interactions. Part I of the chapter focuses on the subthalamic sleep regulatory networks and a model of arousal dynamics capable of predicting sleep, circadian rhythms, and cognitive outputs. Part II presents a model of corticothalamic circuits, capable of predicting spatial and temporal EEG features. We then discuss existing approaches and unsolved challenges in developing unified multiscale models.

The human claustrum tracks slow waves during sleep

Slow waves are a distinguishing feature of non-rapid-eye-movement (NREM) sleep, an evolutionarily conserved process critical for brain function. Non-human studies suggest that the claustrum, a small subcortical nucleus, coordinates slow waves. We show that, in contrast to neurons from other brain regions, claustrum neurons in the human brain increase their spiking activity and track slow waves during NREM sleep, suggesting that the claustrum plays a role in coordinating human sleep architecture.

Oscillatory and nonoscillatory sleep electroencephalographic biomarkers of the epileptic network

OBJECTIVE

In addition to the oscillatory brain activity, the nonoscillatory (scale-free) components of the background electroencephalogram (EEG) may provide further information about the complexity of the underlying neuronal network. As epilepsy is considered a network disease, such scale-free metrics might help to delineate the epileptic network. Here, we performed an analysis of the sleep oscillatory (spindle, slow wave, and rhythmic spectral power) and nonoscillatory (H exponent) intracranial EEG using multiple interictal features to estimate whether and how they deviate from normalcy in 38 adults with drug-resistant epilepsy.

METHODS

To quantify intracranial EEG abnormalities within and outside the seizure onset areas, patients' values were adjusted based on normative maps derived from the open-access Montreal Neurological Institute open iEEG Atlas. In a subset of 29 patients who underwent resective surgery, we estimated the predictive value of these features to identify the epileptogenic zone in those with a good postsurgical outcome.

RESULTS

We found that distinct sleep oscillatory and nonoscillatory metrics behave differently across the epileptic network, with the strongest differences observed for (1) a reduction in spindle activity (spindle rates and rhythmic sigma power in the 10-16 Hz band), (2) a higher rhythmic gamma power (30-80 Hz), and (3) a higher H exponent (steeper 1/f slope). As expected, epileptic spikes were also highest in the seizure onset areas. Furthermore, in surgical patients, the H exponent achieved the highest performance (balanced accuracy of .76) for classifying resected versus nonresected channels in good outcome patients.

SIGNIFICANCE

This work suggests that nonoscillatory components of the intracranial EEG signal could serve as promising interictal sleep candidates of epileptogenicity in patients with drug-resistant epilepsy. Our findings further advance the understanding of epilepsy as a disease, whereby absence or loss of sleep physiology may provide information complementary to pathological epileptic processes.

Neural Dynamics Associated with Biological Variation in Normal Human Brain Regions

The processes involved in encoding and decoding signals in the human brain are a continually studied topic, as neuronal information flow involves complex nonlinear dynamics. This study examines awake human intracranial electroencephalography (iEEG) data from normal brain regions to explore how biological sex influences these dynamics. The iEEG data were analyzed using permutation entropy and statistical complexity in the time domain and power spectrum calculations in the frequency domain. The Bandt and Pompe method was used to assess time series causality by associating probability distributions based on ordinal patterns with the signals. Due to the invasive nature of data acquisition, the study encountered limitations such as small sample sizes and potential sources of error. Nevertheless, the high spatial resolution of iEEG allows detailed analysis and comparison of specific brain regions. The results reveal differences between sexes in brain regions, observed through power spectrum, entropy, and complexity analyses. Significant differences were found in the left supramarginal gyrus, posterior cingulate, supplementary motor cortex, middle temporal gyrus, and right superior temporal gyrus. This study emphasizes the importance of considering sex as a biological variable in brain dynamics research, which is essential for improving the diagnosis and treatment of neurological and psychiatric disorders.

Rényi entropy-complexity causality space: a novel neurocomputational tool for detecting scale-free features in EEG/iEEG data

Scale-free brain activity, linked with learning, the integration of different time scales, and the formation of mental models, is correlated with a metastable cognitive basis. The spectral slope, a key aspect of scale-free dynamics, was proposed as a potential indicator to distinguish between different sleep stages. Studies suggest that brain networks maintain a consistent scale-free structure across wakefulness, anesthesia, and recovery. Although differences in anesthetic sensitivity between the sexes are recognized, these variations are not evident in clinical electroencephalographic recordings of the cortex. Recently, changes in the slope of the power law exponent of neural activity were found to correlate with changes in Rényi entropy, an extended concept of Shannon's information entropy. These findings establish quantifiers as a promising tool for the study of scale-free dynamics in the brain. Our study presents a novel visual representation called the Rényi entropy-complexity causality space, which encapsulates complexity, permutation entropy, and the Rényi parameter q. The main goal of this study is to define this space for classical dynamical systems within theoretical bounds. In addition, the study aims to investigate how well different time series mimicking scale-free activity can be discriminated. Finally, this tool is used to detect dynamic features in intracranial electroencephalography (iEEG) signals. To achieve these goals, the study implementse the Bandt and Pompe method for ordinal patterns. In this process, each signal is associated with a probability distribution, and the causal measures of Rényi entropy and complexity are computed based on the parameter q. This method is a valuable tool for analyzing simulated time series. It effectively distinguishes elements of correlated noise and provides a straightforward means of examining differences in behaviors, characteristics, and classifications. For the iEEG experimental data, the REM state showed a greater number of significant sex-based differences, while the supramarginal gyrus region showed the most variation across different modes and analyzes. Exploring scale-free brain activity with this framework could provide valuable insights into cognition and neurological disorders. The results may have implications for understanding differences in brain function between the sexes and their possible relevance to neurological disorders.

Spatiotemporal EEG dynamics of the sleep onset process in preadolescence

BACKGROUND

During preadolescence the sleep electroencephalography undergoes massive qualitative and quantitative modifications. Despite these relevant age-related peculiarities, the specific EEG pattern of the wake-sleep transition in preadolescence has not been exhaustively described.

METHODS

The aim of the present study is to characterize regional and temporal electrophysiological features of the sleep onset (SO) process in a group of 23 preadolescents (9-14 years) and to compare the topographical pattern of slow wave activity and delta/beta ratio of preadolescents with the EEG pattern of young adults.

RESULTS

Results showed in preadolescence the same dynamics known for adults, but with peculiarities in the delta and beta activity, likely associated with developmental cerebral modifications: the delta power showed a widespread increase during the SO with central maxima, and the lower bins of the beta activity showed a power increase after SO. Compared to adults, preadolescents during the SO exhibited higher delta power only in the slowest bins of the band: before SO slow delta activity was higher in prefrontal, frontal and occipital areas in preadolescents, and, after SO the younger group had higher slow delta activity in occipital areas. In preadolescents delta/beta ratio was higher in more posterior areas both before and after the wake-sleep transition and, after SO, preadolescents showed also a lower delta/beta ratio in frontal areas, compared to adults.

CONCLUSION

Results point to a general higher homeostatic drive for the developing areas, consistently with plastic-related maturational modifications, that physiologically occur during preadolescence.

Alpha anteriorization and theta posteriorization during deep sleep

Brain states (wake, sleep, general anesthesia, etc.) are profoundly associated with the spatiotemporal dynamics of brain oscillations. Previous studies showed that the EEG alpha power shifted from the occipital cortex to the frontal cortex (alpha anteriorization) after being induced into a state of general anesthesia via propofol. The sleep research literature suggests that slow waves and sleep spindles are generated locally and propagated gradually to different brain regions. Since sleep and general anesthesia are conceptualized under the same framework of consciousness, the present study examines whether alpha anteriorization similarly occurs during sleep and how the EEG power in other frequency bands changes during different sleep stages. The results from the analysis of three polysomnography datasets of 234 participants show consistent alpha anteriorization during the sleep stages N2 and N3, beta anteriorization during stage REM, and theta posteriorization during stages N2 and N3. Although it is known that the neural circuits responsible for sleep are not exactly the same for general anesthesia, the findings of alpha anteriorization in this study suggest that, at macro level, the circuits for alpha oscillations are organized in the similar cortical areas. The spatial shifts of EEG power in different frequency bands during sleep may offer meaningful neurophysiological markers for the level of consciousness.

Sleep waves in a large-scale corticothalamic model constrained by activities intrinsic to neocortical networks and single thalamic neurons

AIM

Many biophysical and non-biophysical models have been able to reproduce the corticothalamic activities underlying different EEG sleep rhythms but none of them included the known ability of neocortical networks and single thalamic neurons to generate some of these waves intrinsically.

METHODS

We built a large-scale corticothalamic model with a high fidelity in anatomical connectivity consisting of a single cortical column and first- and higher-order thalamic nuclei. The model is constrained by different neocortical excitatory and inhibitory neuronal populations eliciting slow (<1 Hz) oscillations and by thalamic neurons generating sleep waves when isolated from the neocortex.

RESULTS

Our model faithfully reproduces all EEG sleep waves and the transition from a desynchronized EEG to spindles, slow (<1 Hz) oscillations, and delta waves by progressively increasing neuronal membrane hyperpolarization as it occurs in the intact brain. Moreover, our model shows that slow (<1 Hz) waves most often start in a small assembly of thalamocortical neurons though they can also originate in cortical layer 5. Moreover, the input of thalamocortical neurons increases the frequency of EEG slow (<1 Hz) waves compared to those generated by isolated cortical networks.

CONCLUSION

Our simulations challenge current mechanistic understanding of the temporal dynamics of sleep wave generation and suggest testable predictions.

Scalp and hippocampal sleep correlates of memory function in drug-resistant temporal lobe epilepsy

Seminal animal studies demonstrated the role of sleep oscillations such as cortical slow waves, thalamocortical spindles, and hippocampal ripples in memory consolidation. In humans, whether ripples are involved in sleep-related memory processes is less clear. Here, we explored the interactions between sleep oscillations (measured as traits) and general episodic memory abilities in 26 adults with drug-resistant temporal lobe epilepsy who performed scalp-intracranial electroencephalographic recordings and neuropsychological testing, including two analogous hippocampal-dependent verbal and nonverbal memory tasks. We explored the relationships between hemispheric scalp (spindles, slow waves) and hippocampal physiological and pathological oscillations (spindles, slow waves, ripples, and epileptic spikes) and material-specific memory function. To differentiate physiological from pathological ripples, we used multiple unbiased data-driven clustering approaches. At the individual level, we found material-specific cerebral lateralization effects (left-verbal memory, right-nonverbal memory) for all scalp spindles (rs > 0.51, ps < 0.01) and fast spindles (rs > 0.61, ps < 0.002). Hippocampal epileptic spikes and short pathological ripples, but not physiological oscillations, were negatively (rs > -0.59, ps < 0.01) associated with verbal learning and retention scores, with left lateralizing and antero-posterior effects. However, data-driven clustering failed to separate the ripple events into defined clusters. Correlation analyses with the resulting clusters revealed no meaningful or significant associations with the memory scores. Our results corroborate the role of scalp spindles in memory processes in patients with drug-resistant temporal lobe epilepsy. Yet, physiological and pathological ripples were not separable when using data-driven clustering, and thus our findings do not provide support for a role of sleep ripples as trait-like characteristics of general memory abilities in epilepsy.

The human claustrum tracks slow waves during sleep

Slow waves are a distinguishing feature of non-rapid-eye-movement (NREM) sleep, an evolutionarily conserved process critical for brain function. Non-human studies posit that the claustrum, a small subcortical nucleus, coordinates slow waves. We recorded claustrum neurons in humans during sleep. In contrast to neurons from other brain regions, claustrum neurons increased their activity and tracked slow waves during NREM sleep suggesting that the claustrum plays a role in human sleep architecture.

Interpretation of the Intracranial Electroencephalogram of the Human Hippocampus

Understanding and discriminating the normal and abnormal elements of the intracranial electroencephalogram (iEEG) is essential in decision-making for epilepsy surgery. The hippocampus is widely acknowledged as a key structure in decision-making processes for surgical treatment in temporal lobe epilepsy and epilepsies that involve the mesial temporal structures. This review will provide a summary of the current state of our knowledge and understanding regarding normal and abnormal features of the iEEG of the human hippocampus.

Fractal Neurodynamics

The neuronal ongoing electrical activity in the brain network, the neurodynamics, reflects the structure and functionality of generating neuronal pools. The activity of neurons due to their excitatory and inhibitory projections is associated with specific brain functions. Here, the purpose was to investigate if the local ongoing electrical activity exhibits its characteristic spectral and fractal features in wakefulness and sleep across and within subjects. Moreover, we aimed to show that measures typical of complex systems catch physiological features missed by linear spectral analyses. For this study, we concentrated on the evaluation of the power spectral density (PSD) and Higuchi fractal dimension (HFD) measures. Relevant clinical impact of the specific features of neurodynamics identification stands primarily in the potential of classifying cortical parcels according to their neurodynamics as well as enhancing the effectiveness of neuromodulation interventions to cure symptoms secondary to neuronal activity unbalances.

Wake slow waves in focal human epilepsy impact network activity and cognition

Slow waves of neuronal activity are a fundamental component of sleep that are proposed to have homeostatic and restorative functions. Despite this, their interaction with pathology is unclear and there is only indirect evidence of their presence during wakefulness. Using intracortical recordings from the temporal lobe of 25 patients with epilepsy, we demonstrate the existence of local wake slow waves (LoWS) with key features of sleep slow waves, including a down-state of neuronal firing. Consistent with a reduction in neuronal activity, LoWS were associated with slowed cognitive processing. However, we also found that LoWS showed signatures of a homeostatic relationship with interictal epileptiform discharges (IEDs): exhibiting progressive adaptation during the build-up of network excitability before an IED and reducing the impact of subsequent IEDs on network excitability. We therefore propose an epilepsy homeostasis hypothesis: that slow waves in epilepsy reduce aberrant activity at the price of transient cognitive impairment.

Differential cortical network engagement during states of un/consciousness in humans

What happens in the human brain when we are unconscious? Despite substantial work, we are still unsure which brain regions are involved and how they are impacted when consciousness is disrupted. Using intracranial recordings and direct electrical stimulation, we mapped global, network, and regional involvement during wake vs. arousable unconsciousness (sleep) vs. non-arousable unconsciousness (propofol-induced general anesthesia). Information integration and complex processing we`re reduced, while variability increased in any type of unconscious state. These changes were more pronounced during anesthesia than sleep and involved different cortical engagement. During sleep, changes were mostly uniformly distributed across the brain, whereas during anesthesia, the prefrontal cortex was the most disrupted, suggesting that the lack of arousability during anesthesia results not from just altered overall physiology but from a disconnection between the prefrontal and other brain areas. These findings provide direct evidence for different neural dynamics during loss of consciousness compared with loss of arousability.

Regional variability in intracerebral properties of NREM to REM sleep transitions in humans

Transitions between wake and sleep states show a progressive pattern underpinned by local sleep regulation. In contrast, little evidence is available on non-rapid eye movement (NREM) to rapid eye movement (REM) sleep boundaries, considered as mainly reflecting subcortical regulation. Using polysomnography (PSG) combined with stereoelectroencephalography (SEEG) in humans undergoing epilepsy presurgical evaluation, we explored the dynamics of NREM-to-REM transitions. PSG was used to visually score transitions and identify REM sleep features. SEEG-based local transitions were determined automatically with a machine learning algorithm using features validated for automatic intra-cranial sleep scoring (10.5281/zenodo.7410501). We analyzed 2988 channel-transitions from 29 patients. The average transition time from all intracerebral channels to the first visually marked REM sleep epoch was 8 s ± 1 min 58 s, with a great heterogeneity between brain areas. Transitions were observed first in the lateral occipital cortex, preceding scalp transition by 1 min 57 s ± 2 min 14 s (d = -0.83), and close to the first sawtooth wave marker. Regions with late transitions were the inferior frontal and orbital gyri (1 min 1 s ± 2 min 1 s, d = 0.43, and 1 min 1 s ± 2 min 5 s, d = 0.43, after scalp transition). Intracranial transitions were earlier than scalp transitions as the night advanced (last sleep cycle, d = -0.81). We show a reproducible gradual pattern of REM sleep initiation, suggesting the involvement of cortical mechanisms of regulation. This provides clues for understanding oneiric experiences occurring at the NREM/REM boundary.

Pallidal activities during sleep and sleep decoding in dystonia, Huntington's, and Parkinson's disease

BACKGROUND

Sleep disturbances are highly prevalent in movement disorders, potentially due to the malfunctioning of basal ganglia structures. Pallidal deep brain stimulation (DBS) has been widely used for multiple movement disorders and been reported to improve sleep. We aimed to investigate the oscillatory pattern of pallidum during sleep and explore whether pallidal activities can be utilized to differentiate sleep stages, which could pave the way for sleep-aware adaptive DBS.

METHODS

We directly recorded over 500 h of pallidal local field potentials during sleep from 39 subjects with movement disorders (20 dystonia, 8 Huntington's disease, and 11 Parkinson's disease). Pallidal spectrum and cortical-pallidal coherence were computed and compared across sleep stages. Machine learning approaches were utilized to build sleep decoders for different diseases to classify sleep stages through pallidal oscillatory features. Decoding accuracy was further associated with the spatial localization of the pallidum.

RESULTS

Pallidal power spectra and cortical-pallidal coherence were significantly modulated by sleep-stage transitions in three movement disorders. Differences in sleep-related activities between diseases were identified in non-rapid eye movement (NREM) and REM sleep. Machine learning models using pallidal oscillatory features can decode sleep-wake states with over 90% accuracy. Decoding accuracies were higher in recording sites within the internus-pallidum than the external-pallidum, and can be precited using structural (P < 0.0001) and functional (P < 0.0001) whole-brain neuroimaging connectomics.

CONCLUSION

Our findings revealed strong sleep-stage dependent distinctions in pallidal oscillations in multiple movement disorders. Pallidal oscillatory features were sufficient for sleep stage decoding. These data may facilitate the development of adaptive DBS systems targeting sleep problems that have broad translational prospects.

Chaotic dynamics of the Hénon map and neuronal input-output: A comparison with neurophysiological data

In this study, the Hénon map was analyzed using quantifiers from information theory in order to compare its dynamics to experimental data from brain regions known to exhibit chaotic behavior. The goal was to investigate the potential of the Hénon map as a model for replicating chaotic brain dynamics in the treatment of Parkinson's and epilepsy patients. The dynamic properties of the Hénon map were compared with data from the subthalamic nucleus, the medial frontal cortex, and a q-DG model of neuronal input-output with easy numerical implementation to simulate the local behavior of a population. Using information theory tools, Shannon entropy, statistical complexity, and Fisher's information were analyzed, taking into account the causality of the time series. For this purpose, different windows over the time series were considered. The findings revealed that neither the Hénon map nor the q-DG model could perfectly replicate the dynamics of the brain regions studied. However, with careful consideration of the parameters, scales, and sampling used, they were able to model some characteristics of neural activity. According to these results, normal neural dynamics in the subthalamic nucleus region may present a more complex spectrum within the complexity-entropy causality plane that cannot be represented by chaotic models alone. The dynamic behavior observed in these systems using these tools is highly dependent on the studied temporal scale. As the size of the sample studied increases, the dynamics of the Hénon map become increasingly different from those of biological and artificial neural systems.

Local neurodynamics as a signature of cortical areas: new insights from sleep

Sleep crucial for the animal survival is accompanied by huge changes in neuronal electrical activity over time, the neurodynamics. Here, drawing on intracranial stereo-electroencephalographic (sEEG) recordings from the Montreal Neurological Institute (MNI), we analyzed local neurodynamics in the waking state at rest and during the N2, N3, and rapid eye movement (REM) sleep phases. Higuchi fractal dimension (HFD)-a measure of signal complexity-was studied as a feature of the local neurodynamics of the primary motor (M1), somatosensory (S1), and auditory (A1) cortices. The key working hypothesis, that the relationships between local neurodynamics preserve in all sleep phases despite the neurodynamics complexity reduces in sleep compared with wakefulness, was supported by the results. In fact, while HFD awake > REM > N2 > N3 (P < 0.001 consistently), HFD in M1 > S1 > A1 in awake and all sleep stages (P < 0.05 consistently). Also power spectral density was studied for consistency with previous investigations. Meaningfully, we found a local specificity of neurodynamics, well quantified by the fractal dimension, expressed in wakefulness and during sleep. We reinforce the idea that neurodynamic may become a new criterion for cortical parcellation, prospectively improving the understanding and ability of compensatory interventions for behavioral disorders.

A Subpopulation of Spikes Predicts Successful Epilepsy Surgery Outcome

OBJECTIVE

Epileptic spikes are the traditional interictal electroencephalographic (EEG) biomarker for epilepsy. Given their low specificity for identifying the epileptogenic zone (EZ), they are given only moderate attention in presurgical evaluation. This study aims to demonstrate that it is possible to identify specific spike features in intracranial EEG that optimally define the EZ and predict surgical outcome.

METHODS

We analyzed spike features on stereo-EEG segments from 83 operated patients from 2 epilepsy centers (37 Engel IA) in wakefulness, non-rapid eye movement sleep, and rapid eye movement sleep. After automated spike detection, we investigated 135 spike features based on rate, morphology, propagation, and energy to determine the best feature or feature combination to discriminate the EZ in seizure-free and non-seizure-free patients by applying 4-fold cross-validation.

RESULTS

The rate of spikes with preceding gamma activity in wakefulness performed better for surgical outcome classification (4-fold area under receiver operating characteristics curve [AUC] = 0.755 ± 0.07) than the seizure onset zone, the current gold standard (AUC = 0.563 ± 0.05, p = 0.015) and the ripple rate, an emerging seizure-independent biomarker (AUC = 0.537 ± 0.07, p = 0.006). Channels with a spike-gamma rate exceeding 1.9/min had an 80% probability of being in the EZ. Combining features did not improve the results.

INTERPRETATION

Resection of brain regions with high spike-gamma rates in wakefulness is associated with a high probability of achieving seizure freedom. This rate could be applied to determine the minimal number of spiking channels requiring resection. In addition to quantitative analysis, this feature is easily accessible to visual analysis, which could aid clinicians during presurgical evaluation. ANN NEUROL 2023;93:522-535.

Enhanced phase-amplitude coupling of human electrocorticography selectively in the posterior cortical region during rapid eye movement sleep

The spatiotemporal dynamics of interaction between slow (delta or infraslow) waves and fast (gamma) activities during wakefulness and sleep are yet to be elucidated in human electrocorticography (ECoG). We evaluated phase-amplitude coupling (PAC), which reflects neuronal coding in information processing, using ECoG in 11 patients with intractable focal epilepsy. PAC was observed between slow waves of 0.5-0.6 Hz and gamma activities, not only during light sleep and slow-wave sleep (SWS) but even during wakefulness and rapid eye movement (REM) sleep. While PAC was high over a large region during SWS, it was stronger in the posterior cortical region around the temporoparietal junction than in the frontal cortical region during REM sleep. PAC tended to be higher in the posterior cortical region than in the frontal cortical region even during wakefulness. Our findings suggest that the posterior cortical region has a functional role in REM sleep and may contribute to the maintenance of the dreaming experience.

Thalamocortical coherence and causality in different sleep stages using deep brain stimulation recordings

Previous research has shown an interplay between the thalamus and cerebral cortex during NREM sleep in humans, however the directionality of the thalamocortical synchronization is as yet unknown. In this study thalamocortical connectivity during different NREM sleep stages using sleep scalp electroencephalograms and local field potentials from the left and right anterior thalamus was measured in three epilepsy patients implanted with deep brain stimulation electrodes. Connectivity was assessed as debiased weighted phase lag index and granger causality between the thalamus and cortex for the NREM sleep stages N1, N2 and N3. Results showed connectivity was most prominently directed from cortex to thalamus. Moreover, connectivity varied in strength between the different sleep stages, but barely in direction or frequency. These results imply relatively stable thalamocortical connectivity during NREM sleep directed from the cortex to the thalamus.

A consensus statement on detection of hippocampal sharp wave ripples and differentiation from other fast oscillations

Decades of rodent research have established the role of hippocampal sharp wave ripples (SPW-Rs) in consolidating and guiding experience. More recently, intracranial recordings in humans have suggested their role in episodic and semantic memory. Yet, common standards for recording, detection, and reporting do not exist. Here, we outline the methodological challenges involved in detecting ripple events and offer practical recommendations to improve separation from other high-frequency oscillations. We argue that shared experimental, detection, and reporting standards will provide a solid foundation for future translational discovery.

Spike Timing-Dependent Plasticity with Enhanced Long-Term Depression Leads to an Increase of Statistical Complexity

Synaptic plasticity is characterized by remodeling of existing synapses caused by strengthening and/or weakening of connections. This is represented by long-term potentiation (LTP) and long-term depression (LTD). The occurrence of a presynaptic spike (or action potential) followed by a temporally nearby postsynaptic spike induces LTP; conversely, if the postsynaptic spike precedes the presynaptic spike, it induces LTD. This form of synaptic plasticity induction depends on the order and timing of the pre- and postsynaptic action potential, and has been termed spike time-dependent plasticity (STDP). After an epileptic seizure, LTD plays an important role as a depressor of synapses, which may lead to their complete disappearance together with that of their neighboring connections until days after the event. Added to the fact that after an epileptic seizure the network seeks to regulate the excess activity through two key mechanisms: depressed connections and neuronal death (eliminating excitatory neurons from the network), LTD becomes of great interest in our study. To investigate this phenomenon, we develop a biologically plausible model that privileges LTD at the triplet level while maintaining the pairwise structure in the STPD and study how network dynamics are affected as neuronal damage increases. We find that the statistical complexity is significantly higher for the network where LTD presented both types of interactions. While in the case where the STPD is defined with purely pairwise interactions an increase is observed as damage becomes higher for both Shannon Entropy and Fisher information.

Interictal epileptiform discharges show distinct spatiotemporal and morphological patterns across wake and sleep

Presurgical evaluation of mesial temporal and neocortical focal pharmacoresistant epilepsy patients using intracranial EEG recordings has led to the generation of extensive data on interictal epileptiform discharges, located within or remotely from seizure onset zones. In this study, we used this data to investigate how interictal epileptiform discharges are modulated and how their spatial distribution changes during wake and sleep and analysed the relationship between these discharge events and seizure onset zones. Preoperative evaluation data from 11 adult patients with focal pharmacoresistant epilepsy were extracted from the Epilepsiae database. Interictal epileptiform discharges were automatically detected during wakefulness and over several hours of continuous seizure-free sleep (total duration of EEG recordings:106.7 h; mean per patient: 9.7 h), and analysed across four brain areas (mesial temporal, lateral neocortical, basal cortical and the temporal pole). Sleep stages were classified manually from scalp EEG. Discharge events were characterized according to their rate and morphology (amplitude, sharpness and duration). Eight patients had a seizure onset zone over mesial areas and three patients over lateral neocortical areas. Overall, discharge rates varied across brain areas during wakefulness and sleep [wake/sleep stages × brain areas interaction; Wald χ 2(df = 6) = 31.1, P < 0.0001]. N2-N3 non-rapid eye movement sleep increased interictal epileptiform discharges in mesial areas compared with wakefulness and rapid eye movement sleep (P < 0.0001), and to other areas (P < 0.0001 for all comparisons). This mesial pattern was observed both within and outside of seizure onset zones. During wakefulness, the rate of interictal epileptiform discharges was significantly higher than during N2-N3 non-rapid eye movement sleep (P = 0.04), and rapid eye movement sleep (P = 0.01) in lateral neocortical areas (referred to as lateral neocortical pattern), a finding that was more pronounced in seizures onset zones (P = 0.004). The morphological characteristics of the discharge events were modulated during wakefulness and sleep stages across brain areas. The effect of seizure onset zones on discharge morphology was conditioned by brain area and was particularly marked in temporal pole areas. Our analysis of discharge patterns in relation to cerebral localization, vigilance state and the anatomical affiliation of seizure onset zones revealed the global and local aspects of the complex relationship between interictal discharges, sleep and seizure onset zones. This novel approach may lead to a better understanding of cognitive decline and responses to therapy, as well as to adaptation of surgical interventions for epileptic patients.

Focal epilepsy disrupts spindle structure and function

Sleep spindles are the hallmark of N2 sleep and are attributed a key role in cognition. Little is known about the impact of epilepsy on sleep oscillations underlying sleep-related functions. This study assessed changes in the global spindle rate in patients with epilepsy, analysed the distribution of spindles in relation to the epileptic focus, and performed correlations with neurocognitive function. Twenty-one patients with drug-resistant focal epilepsy (12 females; mean age 32.6 ± 10.7 years [mean ± SD]) and 12 healthy controls (3 females; 24.5 ± 3.3 years) underwent combined whole-night high-density electroencephalography and polysomnography. Global spindle rates during N2 were lower in epilepsy patients compared to controls (mean = 5.78/min ± 0.72 vs. 6.49/min ± 0.71, p = 0.02, d =  − 0.70). Within epilepsy patients, spindle rates were lower in the region of the epileptic focus compared to the contralateral region (median = 4.77/min [range 2.53–6.18] vs. 5.26/min [2.53–6.56], p = 0.02, rank biserial correlation RC =  − 0.57). This decrease was driven by fast spindles (12–16 Hz) (1.50/min [0.62–4.08] vs. 1.65/min [0.51–4.28], p = 0.002, RC =  − 0.76). The focal reduction in spindles was negatively correlated with two scales of attention (r =  − 0.54, p = 0.01; r =  − 0.51, p = 0.025). Patients with focal epilepsy show a reduction in global and local spindle rates dependent on the region of the epileptic focus. This may play a role in impaired cognitive functioning. Future work will show if the local reduction in spindles can be used as potential marker of the epileptic focus.

Waveform detection by deep learning reveals multi-area spindles that are selectively modulated by memory load

Sleep is generally considered to be a state of large-scale synchrony across thalamus and neocortex; however, recent work has challenged this idea by reporting isolated sleep rhythms such as slow oscillations and spindles. What is the spatial scale of sleep rhythms? To answer this question, we adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves in high-noise settings for analysis of neural recordings in sleep. We then studied sleep spindles in non-human primate electrocorticography (ECoG), human electroencephalogram (EEG), and clinical intracranial electroencephalogram (iEEG) recordings in the human. Within each recording type, we find widespread spindles occur much more frequently than previously reported. We then analyzed the spatiotemporal patterns of these large-scale, multi-area spindles and, in the EEG recordings, how spindle patterns change following a visual memory task. Our results reveal a potential role for widespread, multi-area spindles in consolidation of memories in networks widely distributed across primate cortex.

The brain processes memories as we sleep, generating rhythms of electrical activity called ‘sleep spindles’. Sleep spindles were long thought to be a state where the entire brain was fully synchronized by this rhythm. This was based on EEG recordings, short for electroencephalogram, a technique that uses electrodes on the scalp to measure electrical activity in the outermost layer of the brain, the cortex. But more recent intracranial recordings of people undergoing brain surgery have challenged this idea and suggested that sleep spindles may not be a state of global brain synchronization, but rather localised to specific areas.

Mofrad et al. sought to clarify the extent to which spindles co-occur at multiple sites in the brain, which could shed light on how networks of neurons coordinate memory storage during sleep. To analyse highly variable brain wave recordings, Mofrad et al. adapted deep learning algorithms initially developed for detecting earthquakes and gravitational waves. The resulting algorithm, designed to more sensitively detect spindles amongst other brain activity, was then applied to a range of sleep recordings from humans and macaque monkeys.

The analyses revealed that widespread and complex patterns of spindle rhythms, spanning multiple areas in the cortex of the brain, actually appear much more frequently than previously thought. This finding was consistent across all the recordings analysed, even recordings under the skull, which provide the clearest window into brain circuits. Further analyses found that these multi-area spindles occurred more often in sleep after people had completed tasks that required holding many visual scenes in memory, as opposed to control conditions with fewer visual scenes.

In summary, Mofrad et al. show that neuroscientists had previously not appreciated the complex and dynamic patterns in this sleep rhythm. These patterns in sleep spindles may be able to adapt based on the demands needed for memory storage, and this will be the subject of future work. Moreover, the findings support the idea that sleep spindles help coordinate the consolidation of memories in brain circuits that stretch across the cortex. Understanding this mechanism may provide insights into how memory falters in aging and sleep-related diseases, such as Alzheimer’s disease. Lastly, the algorithm developed by Mofrad et al. stands to be a useful tool for analysing other rhythmic waveforms in noisy recordings.

EEG Activation Does Not Differ in Simple and Complex Episodes of Disorders of Arousal: A Spectral Analysis Study

Purpose

Disorders of arousal (DoA) are characterized by incomplete awakening from NREM sleep, with the admixture of both deep sleep and wake EEG activity. Previous observations suggested that changes in EEG activity could be detected in the seconds preceding DoA episodes. The aims of this work were to characterize the topography of EEG spectral changes prior to DoA episodes and to investigate whether or not behavioral complexity could be predicted by changes in EEG immediately preceding behavioral onsets.

Patients and Methods

We collected 103 consecutive video-polysomnographic recordings of 53 DoA adult patients and classified all episodes as simple, rising and complex arousal movements. For each episode, a 5-second window preceding its motor onset ("pre-event") and a 60-second window from 2 to 3 minutes before the episodes ("baseline") were compared. Subsequently, a between-group comparison was performed for the pre-event of simpler versus the more complex episodes.

Results

Spectral analysis over 325 DoA episodes showed an absolute significant increase prior to DoA episodes in all frequency bands excluding sigma, which displayed the opposite effect. In normalized maps, the increase was relatively higher over the central/anterior areas for both slow and fast frequency bands. No significant differences emerged from the comparison between simpler and more complex episodes.

Conclusion

Taken together, these results show that deep sleep and wake-like EEG rhythms coexist over overlapping areas before DoA episodes, suggesting an alteration of local sleep mechanisms. Episodes of different complexity are preceded by a similar EEG activation, implying that they possibly share a similar pathophysiology.

SleepSEEG: Automatic sleep scoring using intracranial EEG recordings only

OBJECTIVE

To perform automatic sleep scoring based only on intracranial EEG, without the need for scalp electroencephalography (EEG), electrooculography (EOG) and electromyography (EMG), in order to study sleep, epilepsy, and their interaction.

APPROACH

Data from 33 adult patients was used for development and training of the automatic scoring algorithm using both oscillatory and non-oscillatory spectral features. The first step consisted in unsupervised clustering of channels based on feature variability. For each cluster the classification was done in two steps, a multiclass tree followed by binary classification trees to distinguish the more challenging stage N1. The test data consisted in 11 patients, in whom the classification was done independently for each channel and then combined to get a single stage per epoch.

MAIN RESULTS

An overall agreement of 78% was observed in the test set between the sleep scoring of the algorithm and two human experts scoring based on scalp EEG, EOG and EMG. Balanced sensitivity and specificity were obtained for the different sleep stages. The performance was excellent for stages W, N2, and N3, and good for stage R, but with high variability across patients. The performance for the challenging stage N1 was poor, but at a similar level as for published algorithms based on scalp EEG. High confidence epochs in different stages (other than N1) can be identified with median per patient specificity >80%.

SIGNIFICANCE

The automatic algorithm can perform sleep scoring of long term recordings of patients with intracranial electrodes undergoing presurgical evaluation in the absence of scalp EEG, EOG and EMG, which are normally required to define sleep stages but are difficult to use in the context of intracerebral studies. It also constitutes a valuable tool to generate hypotheses regarding local aspects of sleep, and will be significant for sleep evaluation in clinical epileptology and neuroscience research.

Region-specific complexity of the intracranial EEG in the sleeping human brain

As the brain is a complex system with occurrence of self-similarity at different levels, a dedicated analysis of the complexity of brain signals is of interest to elucidate the functional role of various brain regions across the various stages of vigilance. We exploited intracranial electroencephalogram data from 38 cortical regions using the Higuchi fractal dimension (HFD) as measure to assess brain complexity, on a dataset of 1772 electrode locations. HFD values depended on sleep stage and topography. HFD increased with higher levels of vigilance, being highest during wakefulness in the frontal lobe. HFD did not change from wake to stage N2 in temporo-occipital regions. The transverse temporal gyrus was the only area in which the HFD did not differ between any two vigilance stages. Interestingly, HFD of wakefulness and stage R were different mainly in the precentral gyrus, possibly reflecting motor inhibition in stage R. The fusiform and parahippocampal gyri were the only areas showing no difference between wakefulness and N2. Stages R and N2 were similar only for the postcentral gyrus. Topographical analysis of brain complexity revealed that sleep stages are clearly differentiated in fronto-central brain regions, but that temporo-occipital regions sleep differently.

Correcting for physiological ripples improves epileptic focus identification and outcome prediction

Objective

The integration of high‐frequency oscillations (HFOs; ripples [80–250 Hz], fast ripples [250–500 Hz]) in epilepsy evaluation is hampered by physiological HFOs, which cannot be reliably differentiated from pathological HFOs. We evaluated whether defining abnormal HFO rates by statistical comparison to region‐specific physiological HFO rates observed in the healthy brain improves identification of the epileptic focus and surgical outcome prediction.

Methods

We detected HFOs in 151 consecutive patients who underwent stereo‐electroencephalography and subsequent resective epilepsy surgery at two tertiary epilepsy centers. We compared how HFOs identified the resection cavity and predicted seizure‐free outcome using two thresholds from the literature (HFO rate > 1/min; 50% of the total number of a patient's HFOs) and three thresholds based on normative rates from the Montreal Neurological Institute Open iEEG Atlas (https://mni‐open‐ieegatlas.research.mcgill.ca/): global Atlas threshold, regional Atlas threshold, and regional + 10% threshold after regional Atlas correction.

Results

Using ripples, the regional + 10% threshold performed best for focus identification (77.3% accuracy, 27% sensitivity, 97.1% specificity, 80.6% positive predictive value [PPV], 78.2% negative predictive value [NPV]) and outcome prediction (69.5% accuracy, 58.6% sensitivity, 76.3% specificity, 60.7% PPV, 74.7% NPV). This was an improvement for focus identification (+1.1% accuracy, +17.0% PPV; p < .001) and outcome prediction (+12.0% sensitivity, +1.0% PPV; p = .05) compared to the 50% threshold. The improvement was particularly marked for foci in cortex, where physiological ripples are frequent (outcome: +35.3% sensitivity, +5.3% PPV; p = .014). In these cases, the regional + 10% threshold outperformed fast ripple rate > 1/min (+3.6% accuracy, +26.5% sensitivity, +21.6% PPV; p < .001) and seizure onset zone (+13.5% accuracy, +29.4% sensitivity, +17.0% PPV; p < .05–.01) for outcome prediction. Normalization did not improve the performance of fast ripples.

Significance

Defining abnormal HFO rates by statistical comparison to rates in healthy tissue overcomes an important weakness in the clinical use of ripples. It improves focus identification and outcome prediction compared to standard HFO measures, increasing their clinical applicability.

Hippocampal spindles and barques are normal intracranial electroencephalographic entities

OBJECTIVE

To assess whether hippocampal spindles and barques are markers of epileptogenicity.

METHODS

Focal epilepsy patients that underwent stereo-electroencephalography implantation with at least one electrode in their hippocampus were selected (n = 75). The occurrence of spindles and barques in the hippocampus was evaluated in each patient. We created pairs of pathologic and pathology-free groups according to two sets of criteria: 1. Non-invasive diagnostic criteria (patients grouped according to focal epilepsy classification). 2. Intracranial neurophysiological criteria (patient's hippocampi grouped according to their seizure onset involvement).

RESULTS

Hippocampal spindles and barques appear equally often in both pathologic and pathology-free groups, both for non-invasive (P_spindles = 0.73; P_barques = 0.46) and intracranial criteria (P_spindles = 0.08; P_barques = 0.26). In Engel Class I patients, spindles occurred with similar incidence both within the non-invasive (P = 0.67) and the intracranial criteria group (P = 0.20). Barques were significantly more frequent in extra-temporal lobe epilepsy defined by either non-invasive (P = 0.01) or intracranial (P = 0.01) criteria.

CONCLUSIONS

Both spindles and barques are normal entities of the hippocampal intracranial electroencephalogram. The presence of barques may also signify lack of epileptogenic properties in the hippocampus.

SIGNIFICANCE

Understanding that hippocampal spindles and barques do not reflect epileptogenicity is critical for correct interpretation of epilepsy surgery evaluations and appropriate surgical treatment selection.

A study on the relationship between non-epileptic fast (40 - 200 Hz) oscillations in scalp EEG and development in children

OBJECTIVE

Physiological gamma and ripple activities may be linked to neurocognitive functions. This study investigated the relationship between development and non-epileptic, probably physiological, fast (40-200 Hz) oscillations (FOs) including gamma (40 - 80 Hz) and ripple (80 - 200 Hz) oscillations in scalp EEG in children with neurodevelopmental disorders.

METHODS

Participants were 124 children with autism spectrum disorder (ASD) and/or attention deficit/hyperactivity disorder (ADHD). Gamma and ripple oscillations were explored from 60-second-long sleep EEG data in each subject using a semi-automatic detection tool supplemented with visual confirmation and time-frequency analysis.

RESULTS

Gamma and ripple oscillations were detected in 25 (20.2%) and 22 (17.7%) children, respectively. The observation of one or more occurrence(s) of ripple oscillations, but not gamma oscillations, was significantly related to lower age at EEG recording (odds ratio, OR: 0.727 [95% confidence interval, CI: 0.568-0.929]), higher intelligence/developmental quotient (OR: 1.041, 95% CI: 1.002-1.082), and lack of a diagnosis with ADHD (OR: 0.191, 95% CI: 0.039 - 0.937) according to a binominal logistic regression analysis that included diagnosis with ASD, sex, history of perinatal complications, history of febrile seizures, and use of a sedative agent for the EEG recording as the other non-significant parameters. Diagnostic group was not related to frequency or power of spectral peaks of FOs.

CONCLUSION

The production of non-epileptic scalp ripples was confirmed to be associated with brain development and function/dysfunction in childhood. Further investigation is necessary to interpret all of the information on higher brain functions that may be embedded in scalp FOs.

Dynamics of hippocampus and orbitofrontal cortex activity during arousing reactions from sleep: An intracranial electroencephalographic study

Sleep is punctuated by transient elevations of vigilance level called arousals or awakenings depending on their durations. Understanding the dynamics of brain activity modifications during these transitional phases could help to better understand the changes in cognitive functions according to vigilance states. In this study, we investigated the activity of memory‐related areas (hippocampus and orbitofrontal cortex) during short (3 s to 2 min) arousing reactions detected from thalamic activity, using intracranial recordings in four drug‐resistant epilepsy patients. The average power of the signal between 0.5 and 128 Hz was compared across four time windows: 10 s of preceding sleep, the first part and the end of the arousal/awakening, and 10 s of wakefulness. We observed that (a) in most frequency bands, the spectral power during hippocampal arousal/awakenings is intermediate between wakefulness and sleep whereas frontal cortex shows an early increase in low and fast activities during non‐rapid‐eye‐movement (NREM) sleep arousals/awakenings; (b) this pattern depends on the preceding sleep stage with fewer modifications for REM than for non‐REM sleep arousal/awakenings, potentially reflecting the EEG similarities between REM sleep and wakefulness; (c) a greater activation at the arousing reaction onset in the prefrontal cortex predicts longer arousals/awakenings. Our findings suggest that hippocampus and prefrontal arousals/awakenings are progressive phenomena modulated by sleep stage, and, in the neocortex, by the intensity of the early activation. This pattern of activity could underlie the link between sleep stage, arousal/awakening duration and restoration of memory abilities including dream recall.

The distinctive sleep pattern of the human calcarine cortex: a stereo-electroencephalographic study

STUDY OBJECTIVES

The aim of this study was to describe the spontaneous electroencephalographic (EEG) features of sleep in the human calcarine cortex, comparing them with the well-established pattern of the parietal cortex.

METHODS

We analyzed presurgical intracerebral EEG activity in calcarine and parietal cortices during non-rapid eye movement (NREM) and rapid eye movement (REM) sleep in seven patients with drug-resistant focal epilepsy. The time course of the EEG spectral power and NREM vs REM differences was assessed. Sleep spindles were automatically detected. To assess homeostatic dynamics, we considered the first vs second half of the night ratio in the delta frequency range (0.5-4 Hz) and the rise rate of delta activity during the first sleep cycle.

RESULTS

While the parietal area showed the classically described NREM and REM sleep hallmarks, the calcarine cortex exhibited a distinctive pattern characterized by: (1) the absence of sleep spindles; (2) a large similarity between EEG power spectra of NREM and REM; and (3) reduced signs of homeostatic dynamics, with a decreased delta ratio between the first and the second half of the night, a reduced rise rate of delta activity during the first NREM sleep cycle, and lack of correlation between these measures.

CONCLUSIONS

Besides describing for the first time the peculiar sleep EEG pattern in the human calcarine cortex, our findings provide evidence that different cortical areas may exhibit specific sleep EEG pattern, supporting the view of sleep as a local process and promoting the idea that the functional role of sleep EEG features should be considered at a regional level.

Focal limbic sources create the large slow oscillations of the EEG in human deep sleep

BACKGROUND

Initial observations with the human electroencephalogram (EEG) have interpreted slow oscillations (SOs) of the EEG during deep sleep (N3) as reflecting widespread surface-negative traveling waves that originate in frontal regions and propagate across the neocortex. However, mapping SOs with a high-density array shows the simultaneous appearance of posterior positive voltage fields in the EEG at the time of the frontal-negative fields, with the typical inversion point (apparent source) around the temporal lobe.

METHODS

Overnight 256-channel EEG recordings were gathered from 10 healthy young adults. Individual head conductivity models were created using each participant's own structural MRI. Source localization of SOs during N3 was then performed.

RESULTS

Electrical source localization models confirmed that these large waves were created by focal discharges within the ventral limbic cortex, including medial temporal and caudal orbitofrontal cortex.

CONCLUSIONS

Although the functional neurophysiology of deep sleep involves interactions between limbic and neocortical networks, the large EEG deflections of deep sleep are not created by distributed traveling waves in lateral neocortex but instead by relatively focal limbic discharges.

Functional Connectivity in Dorsolateral Frontal Cortex: Intracranial Electroencephalogram Study

Motivation: Mechanisms underlying the variation in the appearance of electroencephalogram (EEG) over human head are not well characterized. We hypothesized that spatial variation of the EEG, being ultimately linked to variations in cortical neurobiology, was dependent on cortical connectivity patterns. Specifically, we explored the relationship of resting-state functional connectivity derived from intracranial EEG (iEEG) data in seven (N = 7) human epilepsy patients with the intrinsic dynamic variability of the local iEEG. We asked whether primary and association brain areas over the lateral frontal lobe-due to their sharply different connectivity patterns-were thus dissociable in "EEG space." Methods: Functional connectivity between pairs of subdural grid electrodes was averaged to yield an electrode connectivity (EC) whose time-average yielded mean electrode connectivity (mEC), compared with that electrode's time-averaged sample entropy (SE; mean electrode sample entropy, mESE). Results: We found that mEC and mESE were generally in inverse proportion to each other. Extreme values of mEC and mESE occurred over the Rolandic region and were part of a more general rostrocaudal gradient observed in all patients, with larger (smaller) values of mEC (mESE) occurring anteriorly. Conclusions: Brain networks influence brain dynamics. Over the lateral frontal lobe, mEC and mESE demonstrate a rostrocaudal topography, consistent with current notions regarding the structural and functional parcellation of the human frontal lobe. Our findings distinguish the frontal association cortex from primary sensorimotor cortex, effectively "diagnosing" Rolandic iEEG independent of the classical mu rhythm associated with the latter brain region. Impact statement Electroencephalographic rhythms (electroencephalogram [EEG]) exhibit well-recognized spatial variation over the brain surface. How such variation pertains to the biology of the cortex is poorly understood. Here we identify a novel relationship between sample entropy of the local EEG and the connectivity of that local cortical region to the rest of the brain. Due to the differing connectivities of primary and association motor areas, our methods identify new differences in the EEG arising from those respective brain areas. Our work demonstrates that aspects of brain dynamics (i.e., EEG entropy) may be understood in terms of brain architecture (i.e., functional connectivity) and vice versa.

Neurophysiological brain mapping of human sleep-wake states

OBJECTIVE

We recently proposed a spectrum-based model of the awake intracranial electroencephalogram (iEEG) (Kalamangalam et al., 2020), based on a publicly-available normative database (Frauscher et al., 2018). The latter has been expanded to include data from non-rapid eye movement (NREM) and rapid eye movement (REM) sleep (von Ellenrieder et al., 2020), and the present work extends our methods to those data.

METHODS

Normalized amplitude spectra on semi-logarithmic axes from all four arousal states (wake, N2, N3 and REM) were averaged region-wise and fitted to a multi-component Gaussian distribution. A reduced model comprising five key parameters per brain region was color-coded on to cortical surface models.

RESULTS

The lognormal Gaussian mixture model described the iEEG accurately from all brain regions, in all sleep-wake states. There was smooth variation in model parameters as sleep and wake states yielded to each other. Specific observations unrelated to the model were that the primary cortical areas of vision, motor function and audition, in addition to the hippocampus, did not participate in the 'awakening' of the cortex during REM sleep.

CONCLUSIONS

Despite the significant differences in the appearance of the time-domain EEG in wakefulness and sleep, the iEEG in all arousal states was successfully described by a parametric spectral model of low dimension.

SIGNIFICANCE

Spectral variation in the iEEG is continuous in space (across different cortical regions) and time (stage of circadian cycle), arguing for a 'continuum' hypothesis in the generative processes of sleep and wakefulness in human brain.

What do intracerebral electrodes measure?

OBJECTIVE

Gain insight and improve our interpretation of measurements from intracerebral electrodes. Determine if interpretation of intracerebral EEG is dependent on electrode characteristics.

METHODS

We use intracerebral EEG measurements differing only in the recording electrodes (Dixi or homemade electrodes), and numerical simulations to determine the spatial sensitivity of intracerebral electrodes and its dependence on several parameters.

RESULTS

There is a difference in the high frequency (>20 Hz) power depending on the electrode type, which cannot be explained by the different contact sizes or distance between contacts. Simulations show that the width of the gap between electrode and brain and the extent of the generators have an effect on sensitivity, while other parameters are less important.

CONCLUSIONS

The sensitivity of intracerebral electrodes is not affected in an important way by the dimensions of the contacts, but depends on the extent of generators. The unusual insertion technique of homemade electrodes resulting in a large gap between functional brain and electrodes, explains the observed signal difference.

SIGNIFICANCE

Numerical simulation is a useful tool in the choice or design of intracerebral electrodes, and increases our understanding of their measurements. The interpretation of intracerebral EEG is not affected by differences between typical commercially available electrodes.

Technical Aspects of SEEG and Its Interpretation in the Delineation of the Epileptogenic Zone

Stereo-electroencephalography (SEEG) has gained global popularity in recent years. In Japan, a country in which invasive studies using subdural electrodes (SDEs) have been the mainstream, SEEG has been approved for insurance coverage in 2020 and is expected to gain in popularity. Some concepts supporting SEEG methodology are fundamentally different from that of SDE studies. Clinicians interested in utilizing SEEG in their practice should be aware of those aspects in which they differ. Success in utilizing the SEEG methodology relies heavily on the construction of an a priori hypothesis regarding the putative seizure onset zone (SOZ) and propagation. This article covers the technical and theoretical aspects of SEEG, including the surgical techniques and precautions, hypothesis construction, and the interpretation of the recording, all with the aim of providing an introductory guide to SEEG.

Rapid Eye Movement Sleep Sawtooth Waves Are Associated with Widespread Cortical Activations

Sawtooth waves (STW) are bursts of frontocentral slow oscillations recorded in the scalp electroencephalogram (EEG) during rapid eye movement (REM) sleep. Little is known about their cortical generators and functional significance. Stereo-EEG performed for presurgical epilepsy evaluation offers the unique possibility to study neurophysiology in situ in the human brain. We investigated intracranial correlates of scalp-detected STW in 26 patients (14 women) undergoing combined stereo-EEG/polysomnography. We visually marked STW segments in scalp EEG and selected stereo-EEG channels exhibiting normal activity for intracranial analyses. Channels were grouped in 30 brain regions. The spectral power in each channel and frequency band was computed during STW and non-STW control segments. Ripples (80-250 Hz) were automatically detected during STW and control segments. The spectral power in the different frequency bands and the ripple rates were then compared between STW and control segments in each brain region. An increase in 2-4 Hz power during STW segments was found in all brain regions, except the occipital lobe, with large effect sizes in the parietotemporal junction, the lateral and orbital frontal cortex, the anterior insula, and mesiotemporal structures. A widespread increase in high-frequency activity, including ripples, was observed concomitantly, involving the sensorimotor cortex, associative areas, and limbic structures. This distribution showed a high spatiotemporal heterogeneity. Our results suggest that STW are associated with widely distributed, but locally regulated REM sleep slow oscillations. By driving fast activities, STW may orchestrate synchronized reactivations of multifocal activities, allowing tagging of complex representations necessary for REM sleep-dependent memory consolidation.SIGNIFICANCE STATEMENT Sawtooth waves (STW) present as scalp electroencephalographic (EEG) bursts of slow waves contrasting with the low-voltage fast desynchronized activity of REM sleep. Little is known about their cortical origin and function. Using combined stereo-EEG/polysomnography possible only in the human brain during presurgical epilepsy evaluation, we explored the intracranial correlates of STW. We found that a large set of regions in the parietal, frontal, and insular cortices shows increases in 2-4 Hz power during scalp EEG STW, that STW are associated with a strong and widespread increase in high frequencies, and that these slow and fast activities exhibit a high spatiotemporal heterogeneity. These electrophysiological properties suggest that STW may be involved in cognitive processes during REM sleep.

Brain-scale cortico-cortical functional connectivity in the delta-theta band is a robust signature of conscious states: an intracranial and scalp EEG study

Long-range cortico-cortical functional connectivity has long been theorized to be necessary for conscious states. In the present work, we estimate long-range cortical connectivity in a series of intracranial and scalp EEG recordings experiments. In the two first experiments intracranial-EEG (iEEG) was recorded during four distinct states within the same individuals: conscious wakefulness (CW), rapid-eye-movement sleep (REM), stable periods of slow-wave sleep (SWS) and deep propofol anaesthesia (PA). We estimated functional connectivity using the following two methods: weighted Symbolic-Mutual-Information (wSMI) and phase-locked value (PLV). Our results showed that long-range functional connectivity in the delta-theta frequency band specifically discriminated CW and REM from SWS and PA. In the third experiment, we generalized this original finding on a large cohort of brain-injured patients. FC in the delta-theta band was significantly higher in patients being in a minimally conscious state (MCS) than in those being in a vegetative state (or unresponsive wakefulness syndrome). Taken together the present results suggest that FC of cortical activity in this slow frequency band is a new and robust signature of conscious states.

Interpretation of the Intracranial Stereoelectroencephalography Signal

The intracranial electroencephalogram (iEEG) is essential in decision making for epilepsy surgery. Although localization of epileptogenic brain regions by means of iEEG has been the gold standard for surgical decision-making for more than 70 years, established guidelines for what constitutes genuine iEEG epileptic activity and what is normal brain activity are not available. This review provides a summary of the current state of knowledge and understanding on normal iEEG entities and variants, the effects of sleep on regional and lobar iEEG, iEEG patterns of interictal and ictal epileptic activity and their relation to well-described epileptogenic pathologies and surgical outcome.

The human K-complex: Insights from combined scalp-intracranial EEG recordings

Sleep spindles and K-complexes (KCs) are a hallmark of N2 sleep. While the functional significance of spindles is comparatively well investigated, there is still ongoing debate about the role of the KC: it is unclear whether it is a cortical response to an arousing stimulus (either external or internal) or whether it has sleep-promoting properties. Invasive intracranial EEG recordings from individuals with drug-resistant epilepsy offer a unique opportunity to study in-situ human brain physiology. To better understand the function of the KC, we aimed to (i) investigate the intracranial correlates of spontaneous scalp KCs, and (ii) compare the intracranial activity of scalp KCs associated or not with arousals. Whole-night recordings from adults with drug-resistant focal epilepsy who underwent combined intracranial-scalp EEG for pre-surgical evaluation at the Montreal Neurological Institute between 2010 and 2018 were selected. KCs were visually marked in the scalp and categorized according to the presence of microarousals: (i) Pre-microarousal KCs; (ii) KCs during an ongoing microarousal; and (iii) KCs without microarousal. Power in different spectral bands was computed to compare physiological intracranial EEG activity at the time of scalp KCs relative to the background, as well as to compare microarousal subcategories. A total of 1198 scalp KCs selected from 40 subjects were analyzed, resulting in 32,504 intracranial KC segments across 992 channels. Forty-seven percent of KCs were without microarousal, 30% were pre-microarousal, and 23% occurred during microarousals. All scalp KCs were accompanied by widespread cortical increases in delta band power (0.3-4 ​Hz) relative to the background: the highest percentages were observed in the parietal (60-65%) and frontal cortices (52-58%). Compared to KCs without microarousal, pre-microarousal KCs were accompanied by increases (66%) in beta band power (16-30 ​Hz) in the motor cortex, which was present before the peak of the KC. In addition, spatial distribution of spectral power changes following each KC without microarousal revealed that certain brain regions were associated with increases in delta power (25-62%) or decreases in alpha/beta power (11-24%), suggesting a sleep-promoting pattern, whereas others were accompanied by increases of higher frequencies (12-27%), suggesting an arousal-related pattern. This study shows that KCs can be generated across widespread cortical areas. Interestingly, the motor cortex shows awake-like EEG activity before the onset of KCs followed by microarousals. Our findings also highlight region-specific sleep- or arousal-promoting responses following KCs, suggesting a dual role for the human KC.