Large-scale cortical dynamics of sleep slow waves

Cerebral Gray Matter May Not Explain Sleep Slow-Wave Characteristics after Severe Brain Injury

Sleep slow waves are the hallmark of deeper non-rapid eye movement sleep. It is generally assumed that gray matter properties predict slow-wave density, morphology, and spectral power in healthy adults. Here, we tested the association between gray matter volume (GMV) and slow-wave characteristics in 27 patients with moderate-to-severe traumatic brain injury (TBI, 32.0 ± 12.2 years old, eight women) and compared that with 32 healthy controls (29.2 ± 11.5 years old, nine women). Participants underwent overnight polysomnography and cerebral MRI with a 3 Tesla scanner. A whole-brain voxel-wise analysis was performed to compare GMV between groups. Slow-wave density, morphology, and spectral power (0.4-6 Hz) were computed, and GMV was extracted from the thalamus, cingulate, insula, precuneus, and orbitofrontal cortex to test the relationship between slow waves and gray matter in regions implicated in the generation and/or propagation of slow waves. Compared with controls, TBI patients had significantly lower frontal and temporal GMV and exhibited a subtle decrease in slow-wave frequency. Moreover, higher GMV in the orbitofrontal cortex, insula, cingulate cortex, and precuneus was associated with higher slow-wave frequency and slope, but only in healthy controls. Higher orbitofrontal GMV was also associated with higher slow-wave density in healthy participants. While we observed the expected associations between GMV and slow-wave characteristics in healthy controls, no such associations were observed in the TBI group despite lower GMV. This finding challenges the presumed role of GMV in slow-wave generation and morphology.

A modular and adaptable analysis pipeline to compare slow cerebral rhythms across heterogeneous datasets

Neuroscience is moving toward a more integrative discipline where understanding brain function requires consolidating the accumulated evidence seen across experiments, species, and measurement techniques. A remaining challenge on that path is integrating such heterogeneous data into analysis workflows such that consistent and comparable conclusions can be distilled as an experimental basis for models and theories. Here, we propose a solution in the context of slow-wave activity (<1 Hz), which occurs during unconscious brain states like sleep and general anesthesia and is observed across diverse experimental approaches. We address the issue of integrating and comparing heterogeneous data by conceptualizing a general pipeline design that is adaptable to a variety of inputs and applications. Furthermore, we present the Collaborative Brain Wave Analysis Pipeline (Cobrawap) as a concrete, reusable software implementation to perform broad, detailed, and rigorous comparisons of slow-wave characteristics across multiple, openly available electrocorticography (ECoG) and calcium imaging datasets.

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.

Local slow-wave activity over the right prefrontal cortex reveals individual risk preferences

In everyday life, we have to make decisions under varying degrees of risk. Even though previous research has shown that the manipulation of sleep affects risky decision-making, it remains unknown whether individual, temporally stable neural sleep characteristics relate to individual differences in risk preferences. Here, we collected sleep data under normal conditions in fifty-four healthy adults using a portable high-density EEG at participants' home. Whole-brain corrected for multiple testing, we found that lower slow-wave activity (SWA, an indicator of sleep depth) in a cluster of electrodes over the right prefrontal cortex (PFC) is associated with higher individual risk propensity. Importantly, the association between local sleep depth and risk preferences remained significant when controlling for total sleep time and for time spent in deep sleep, i.e., sleep stages N2 and N3. Moreover, the association between risk preferences and SWA over the right PFC was very similar in all sleep cycles. Because the right PFC plays a central role in cognitive control functions, we speculate that local sleep depth in this area, as reflected by SWA, might serve as a dispositional indicator of self-regulatory ability, which in turn reflects risk preferences.

Intracranial recordings reveal ubiquitous in-phase and in-antiphase functional connectivity between homotopic brain regions in humans

Whether neuronal populations exhibit zero-lag (in-phase or in-antiphase) functional connectivity is a fundamental question when conceptualizing communication between cell assemblies. It also has profound implications on how we assess such interactions. Given that the brain is a delayed network due to the finite conduction velocity of the electrical impulses traveling across its fibers, the existence of long-distance zero-lag functional connectivity may be considered improbable. However, in this study, using human intracranial recordings we demonstrate that most interhemispheric connectivity between homotopic cerebral regions is zero-lagged and that this type of connectivity is ubiquitous. Volume conduction can be safely discarded as a confounding factor since it is known to drop almost completely within short interelectrode distances (<20 mm) in intracranial recordings. This finding should guide future electrophysiological connectivity studies and highlight the importance of considering the role of zero-lag connectivity in our understanding of communication between cell assemblies.

Nonlinear analysis of periodic waves in a neural field model

Various types of brain activity, including motor, visual, and language, are accompanied by the propagation of periodic waves of electric potential in the cortex, possibly providing the synchronization of the epicenters involved in these activities. One example is cortical electrical activity propagating during sleep and described as traveling waves [Massimini et al., J. Neurosci. 24, 6862-6870 (2004)]. These waves modulate cortical excitability as they progress. Clinically related examples include cortical spreading depression in which a wave of depolarization propagates not only in migraine but also in stroke, hemorrhage, or traumatic brain injury [Whalen et al., Sci. Rep. 8, 1-9 (2018)]. Here, we consider the possible role of epicenters and explore a neural field model with two nonlinear integrodifferential equations for the distributions of activating and inhibiting signals. It is studied with symmetric connectivity functions characterizing signal exchange between two populations of neurons, excitatory and inhibitory. Bifurcation analysis is used to investigate the emergence of periodic traveling waves and of standing oscillations from the stationary, spatially homogeneous solutions, and the stability of these solutions. Both types of solutions can be started by local oscillations indicating a possible role of epicenters in the initiation of wave propagation.

Characterizing Sleep Spindles in Sheep

Sleep spindles are distinctive transient patterns of brain activity that typically occur during non-rapid eye movement (NREM) sleep in humans and other mammals. Thought to be important for the consolidation of learning, they may also be useful for indicating the progression of aging and neurodegenerative diseases. The aim of this study was to characterize sleep spindles in sheep (Ovis aries). We recorded electroencephalographs wirelessly from six sheep over a continuous period containing 2 nights and a day. We detected and characterized spindles using an automated algorithm. We found that sheep sleep spindles fell within the classical range seen in humans (10–16 Hz), but we did not see a further separation into fast and slow bands. Spindles were detected predominantly during NREM sleep. Spindle characteristics (frequency, duration, density, topography) varied between individuals, but were similar within individuals between nights. Spindles that occurred during NREM sleep in daytime were indistinguishable from those found during NREM sleep at night. Surprisingly, we also detected numerous spindle-like events during unequivocal periods of wake during the day. These events were mainly local (detected at single sites), and their characteristics differed from spindles detected during sleep. These “wake spindles” are likely to be events that are commonly categorized as “spontaneous alpha activity” during wake. We speculate that wake and sleep spindles are generated via different mechanisms, and that wake spindles play a role in cognitive processes that occur during the daytime.

Complex Dynamics of Propagating Waves in a Two-Dimensional Neural Field

Propagating waves with complex dynamics have been widely observed in neural population activity. To understand their formation mechanisms, we investigate a type of two-dimensional neural field model by systematically varying its recurrent excitatory and inhibitory inputs. We show that the neural field model exhibits a rich repertoire of dynamical activity states when the relevant strength of excitation and inhibition is increased, ranging from localized rotating and traveling waves to global waves. Particularly, near the transition between stable states of rotating and traveling waves, the model exhibits a bistable state; that is, both the rotating and the traveling waves can exist, and the inclusion of noise can induce spontaneous transitions between them. Furthermore, we demonstrate that when there are multiple propagating waves, they exhibit rich collective propagation dynamics with variable propagating speeds and trajectories. We use techniques from time series analysis such detrended fluctuation analysis to characterize the effect of the strength of excitation and inhibition on these collective dynamics, which range from purely random motion to motion with long-range spatiotemporal correlations. These results provide insights into the possible contribution of excitation and inhibition toward a range of previously observed spatiotemporal wave phenomena.

Metastable brain waves

Traveling patterns of neuronal activity-brain waves-have been observed across a breadth of neuronal recordings, states of awareness, and species, but their emergence in the human brain lacks a firm understanding. Here we analyze the complex nonlinear dynamics that emerge from modeling large-scale spontaneous neural activity on a whole-brain network derived from human tractography. We find a rich array of three-dimensional wave patterns, including traveling waves, spiral waves, sources, and sinks. These patterns are metastable, such that multiple spatiotemporal wave patterns are visited in sequence. Transitions between states correspond to reconfigurations of underlying phase flows, characterized by nonlinear instabilities. These metastable dynamics accord with empirical data from multiple imaging modalities, including electrical waves in cortical tissue, sequential spatiotemporal patterns in resting-state MEG data, and large-scale waves in human electrocorticography. By moving the study of functional networks from a spatially static to an inherently dynamic (wave-like) frame, our work unifies apparently diverse phenomena across functional neuroimaging modalities and makes specific predictions for further experimentation.

Increased cortical involvement and synchronization during CAP A1 slow waves

Slow waves recorded with EEG in NREM sleep are indicative of the strength and spatial extent of synchronized firing in neuronal assemblies of the cerebral cortex. Slow waves often appear in the A1 part of the cyclic alternating patterns (CAP), which correlate with a number of behavioral and biological parameters, but their physiological significance is not adequately known. We automatically detected slow waves from the scalp recordings of 37 healthy patients, visually identified CAP A1 events and compared slow waves during CAP A1 with those during NCAP. For each slow wave, we computed the amplitude, slopes, frequency, synchronization (synchronization likelihood) between specific cortical areas, as well as the location of origin and scalp propagation of individual waves. CAP A1 slow waves were characterized by greater spatial extent and amplitude, steeper slopes and greater cortical synchronization, but a similar prominence in frontal areas and similar propagation patterns to other areas on the scalp. Our results indicate that CAP A1 represents a period of highly synchronous neuronal firing over large areas of the cortical mantle. This feature may contribute to the role CAP A1 plays in both normal synaptic homeostasis and in the generation of epileptiform phenomena in epileptic patients.

Desynchronization of slow oscillations in the basal ganglia during natural sleep

Slow oscillations of neuronal activity alternating between firing and silence are a hallmark of slow-wave sleep (SWS). These oscillations reflect the default activity present in all mammalian species, and are ubiquitous to anesthesia, brain slice preparations, and neuronal cultures. In all these cases, neuronal firing is highly synchronous within local circuits, suggesting that oscillation-synchronization coupling may be a governing principle of sleep physiology regardless of anatomical connectivity. To investigate whether this principle applies to overall brain organization, we recorded the activity of individual neurons from basal ganglia (BG) structures and the thalamocortical (TC) network over 70 full nights of natural sleep in two vervet monkeys. During SWS, BG neurons manifested slow oscillations (∼0.5 Hz) in firing rate that were as prominent as in the TC network. However, in sharp contrast to any neural substrate explored thus far, the slow oscillations in all BG structures were completely desynchronized between individual neurons. Furthermore, whereas in the TC network single-cell spiking was locked to slow oscillations in the local field potential (LFP), the BG LFP exhibited only weak slow oscillatory activity and failed to entrain nearby cells. We thus show that synchrony is not inherent to slow oscillations, and propose that the BG desynchronization of slow oscillations could stem from its unique anatomy and functional connectivity. Finally, we posit that BG slow-oscillation desynchronization may further the reemergence of slow-oscillation traveling waves from multiple independent origins in the frontal cortex, thus significantly contributing to normal SWS.

State-Dependent Modulation of Slow Wave Motifs towards Awakening

Slow cortical waves that propagate across the cerebral cortex forming large-scale spatiotemporal propagation patterns are a hallmark of non-REM sleep and anesthesia, but also occur during resting wakefulness. To investigate how the spatial temporal properties of slow waves change with the depth of anesthetic, we optically imaged population voltage transients generated by mouse layer 2/3 pyramidal neurons across one or two cortical hemispheres dorsally with a genetically encoded voltage indicator (GEVI). From deep barbiturate anesthesia to light barbiturate sedation, depolarizing wave events recruiting at least 50% of the imaged cortical area consistently appeared as a conserved repertoire of distinct wave motifs. Toward awakening, the incidence of individual motifs changed systematically (the motif propagating from visual to motor areas increased while that from somatosensory to visual areas decreased) and both local and global cortical dynamics accelerated. These findings highlight that functional endogenous interactions between distant cortical areas are not only constrained by anatomical connectivity, but can also be modulated by the brain state.

Behavioral state classification in epileptic brain using intracranial electrophysiology

OBJECTIVE

Automated behavioral state classification can benefit next generation implantable epilepsy devices. In this study we explored the feasibility of automated awake (AW) and slow wave sleep (SWS) classification using wide bandwidth intracranial EEG (iEEG) in patients undergoing evaluation for epilepsy surgery.

APPROACH

Data from seven patients (age [Formula: see text], 4 women) who underwent intracranial depth electrode implantation for iEEG monitoring were included. Spectral power features (0.1-600 Hz) spanning several frequency bands from a single electrode were used to train and test a support vector machine classifier.

MAIN RESULTS

Classification accuracy of 97.8  ±  0.3% (normal tissue) and 89.4  ±  0.8% (epileptic tissue) across seven subjects using multiple spectral power features from a single electrode was achieved. Spectral power features from electrodes placed in normal temporal neocortex were found to be more useful (accuracy 90.8  ±  0.8%) for sleep-wake state classification than electrodes located in normal hippocampus (87.1  ±  1.6%). Spectral power in high frequency band features (Ripple (80-250 Hz), Fast Ripple (250-600 Hz)) showed comparable performance for AW and SWS classification as the best performing Berger bands (Alpha, Beta, low Gamma) with accuracy  ⩾90% using a single electrode contact and single spectral feature.

SIGNIFICANCE

Automated classification of wake and SWS should prove useful for future implantable epilepsy devices with limited computational power, memory, and number of electrodes. Applications include quantifying patient sleep patterns and behavioral state dependent detection, prediction, and electrical stimulation therapies.

Large Scale Cortical Functional Networks Associated with Slow-Wave and Spindle-Burst-Related Spontaneous Activity

Cortical sensory systems are active with rich patterns of activity during sleep and under light anesthesia. Remarkably, this activity shares many characteristics with those present when the awake brain responds to sensory stimuli. We review two specific forms of such activity: slow-wave activity (SWA) in the adult brain and spindle bursts in developing brain. SWA is composed of 0.5-4 Hz resting potential fluctuations. Although these fluctuations synchronize wide regions of cortex, recent large-scale imaging has shown spatial details of their distribution that reflect underlying cortical structural projections and networks. These networks are regulated, as prior awake experiences alter both the spatial and temporal features of SWA in subsequent sleep. Activity patterns of the immature brain, however, are very different from those of the adult. SWA is absent, and the dominant pattern is spindle bursts, intermittent high frequency oscillations superimposed on slower depolarizations within sensory cortices. These bursts are driven by intrinsic brain activity, which act to generate peripheral inputs, for example via limb twitches. They are present within developing sensory cortex before they are mature enough to exhibit directed movements and respond to external stimuli. Like in the adult, these patterns resemble those evoked by sensory stimulation when awake. It is suggested that spindle-burst activity is generated purposefully by the developing nervous system as a proxy for true external stimuli. While the sleep-related functions of both slow-wave and spindle-burst activity may not be entirely clear, they reflect robust regulated phenomena which can engage select wide-spread cortical circuits. These circuits are similar to those activated during sensory processing and volitional events. We highlight these two patterns of brain activity because both are prominent and well-studied forms of spontaneous activity that will yield valuable insights into brain function in the coming years.

Optimizing the Detection of Wakeful and Sleep-Like States for Future Electrocorticographic Brain Computer Interface Applications

Previous studies suggest stable and robust control of a brain-computer interface (BCI) can be achieved using electrocorticography (ECoG). Translation of this technology from the laboratory to the real world requires additional methods that allow users operate their ECoG-based BCI autonomously. In such an environment, users must be able to perform all tasks currently performed by the experimenter, including manually switching the BCI system on/off. Although a simple task, it can be challenging for target users (e.g., individuals with tetraplegia) due to severe motor disability. In this study, we present an automated and practical strategy to switch a BCI system on or off based on the cognitive state of the user. Using a logistic regression, we built probabilistic models that utilized sub-dural ECoG signals from humans to estimate in pseudo real-time whether a person is awake or in a sleep-like state, and subsequently, whether to turn a BCI system on or off. Furthermore, we constrained these models to identify the optimal anatomical and spectral parameters for delineating states. Other methods exist to differentiate wake and sleep states using ECoG, but none account for practical requirements of BCI application, such as minimizing the size of an ECoG implant and predicting states in real time. Our results demonstrate that, across 4 individuals, wakeful and sleep-like states can be classified with over 80% accuracy (up to 92%) in pseudo real-time using high gamma (70-110 Hz) band limited power from only 5 electrodes (platinum discs with a diameter of 2.3 mm) located above the precentral and posterior superior temporal gyrus.

Distribution, Amplitude, Incidence, Co-Occurrence, and Propagation of Human K-Complexes in Focal Transcortical Recordings

K-complexes (KCs) are thought to play a key role in sleep homeostasis and memory consolidation; however, their generation and propagation remain unclear. The commonly held view from scalp EEG findings is that KCs are primarily generated in medial frontal cortex and propagate parietally, whereas an electrocorticography (ECOG) study suggested dorsolateral prefrontal generators and an absence of KCs in many areas. In order to resolve these differing views, we used unambiguously focal bipolar depth electrode recordings in patients with intractable epilepsy to investigate spatiotemporal relationships of human KCs. KCs were marked manually on each channel, and local generation was confirmed with decreased gamma power. In most cases (76%), KCs occurred in a single location, and rarely (1%) in all locations. However, if automatically detected KC-like phenomena were included, only 15% occurred in a single location, and 27% occurred in all recorded locations. Locally generated KCs were found in all sampled areas, including cingulate, ventral temporal, and occipital cortices. Surprisingly, KCs were smallest and occurred least frequently in anterior prefrontal channels. When KCs occur on two channels, their peak order is consistent in only 13% of cases, usually from prefrontal to lateral temporal. Overall, the anterior-posterior separation of electrode pairs explained only 2% of the variance in their latencies. KCs in stages 2 and 3 had similar characteristics. These results open a novel view where KCs overall are universal cortical phenomena, but each KC may variably involve small or large cortical regions and spread in variable directions, allowing flexible and heterogeneous contributions to sleep homeostasis and memory consolidation.

Neuronal avalanches differ from wakefulness to deep sleep--evidence from intracranial depth recordings in humans

Neuronal activity differs between wakefulness and sleep states. In contrast, an attractor state, called self-organized critical (SOC), was proposed to govern brain dynamics because it allows for optimal information coding. But is the human brain SOC for each vigilance state despite the variations in neuronal dynamics? We characterized neuronal avalanches--spatiotemporal waves of enhanced activity--from dense intracranial depth recordings in humans. We showed that avalanche distributions closely follow a power law--the hallmark feature of SOC--for each vigilance state. However, avalanches clearly differ with vigilance states: slow wave sleep (SWS) shows large avalanches, wakefulness intermediate, and rapid eye movement (REM) sleep small ones. Our SOC model, together with the data, suggested first that the differences are mediated by global but tiny changes in synaptic strength, and second, that the changes with vigilance states reflect small deviations from criticality to the subcritical regime, implying that the human brain does not operate at criticality proper but close to SOC. Independent of criticality, the analysis confirms that SWS shows increased correlations between cortical areas, and reveals that REM sleep shows more fragmented cortical dynamics.

Individual differences in white matter diffusion affect sleep oscillations

The characteristic oscillations of the sleeping brain, spindles and slow waves, show trait-like, within-subject stability and a remarkable interindividual variability that correlates with functionally relevant measures such as memory performance and intelligence. Yet, the mechanisms underlying these interindividual differences are largely unknown. Spindles and slow waves are affected by the recent history of learning and neuronal activation, indicating sensitivity to changes in synaptic strength and thus to the connectivity of the neuronal network. Because the structural backbone of this network is formed by white matter tracts, we hypothesized that individual differences in spindles and slow waves depend on the white matter microstructure across a distributed network. We recorded both diffusion-weighted magnetic resonance images and whole-night, high-density electroencephalography and investigated whether individual differences in sleep spindle and slow wave parameters were associated with diffusion tensor imaging metrics; white matter fractional anisotropy and axial diffusivity were quantified using tract-based spatial statistics. Individuals with higher spindle power had higher axial diffusivity in the forceps minor, the anterior corpus callosum, fascicles in the temporal lobe, and the tracts within and surrounding the thalamus. Individuals with a steeper rising slope of the slow wave had higher axial diffusivity in the temporal fascicle and frontally located white matter tracts (forceps minor, anterior corpus callosum). These results indicate that the profiles of sleep oscillations reflect not only the dynamics of the neuronal network at the synaptic level, but also the localized microstructural properties of its structural backbone, the white matter tracts.