A continuum model for the dynamics of the phase transition from slow-wave sleep to REM sleep

Previous studies have shown that activated cortical states (awake and rapid eye-movement (REM) sleep), are associated with increased cholinergic input into the cerebral cortex. However, the mechanisms that underlie the detailed dynamics of the cortical transition from slow-wave to REM sleep have not been quantitatively modeled. How does the sequence of abrupt changes in the cortical dynamics (as detected in the electrocorticogram) result from the more gradual change in subcortical cholinergic input? We compare the output from a continuum model of cortical neuronal dynamics with experimentally-derived rat electrocorticogram data. The output from the computer model was consistent with experimental observations. In slow-wave sleep, 0.5–2-Hz oscillations arise from the cortex jumping between “up” and “down” states on the stationary-state manifold. As cholinergic input increases, the upper state undergoes a bifurcation to an 8-Hz oscillation. The coexistence of both oscillations is similar to that found in the intermediate stage of sleep of the rat. Further cholinergic input moves the trajectory to a point where the lower part of the manifold in not available, and thus the slow oscillation abruptly ceases (REM sleep). The model provides a natural basis to explain neuromodulator-induced changes in cortical activity, and indicates that a cortical phase change, rather than a brainstem “flip-flop”, may describe the transition from slow-wave sleep to REM

Phase transitions in single neurons and neural populations: Critical slowing, anesthesia, and sleep cycles

The firing of an action potential by a biological neuron represents a dramatic transition from small-scale linear stochastics (subthreshold voltage fluctuations) to gross-scale nonlinear dynamics (birth of a 1-ms voltage spike). In populations of neurons we see similar, but slower, switch-like there-and-back transitions between low-firing background states and high-firing activated states. These state transitions are controlled by varying levels of input current (single neuron), varying amounts of GABAergic drug (anesthesia), or varying concentrations of neuromodulators and neurotransmitters (natural sleep), and all occur within a milieu of unrelenting biological noise. By tracking the altering responsiveness of the excitable membrane to noisy stimulus, we can infer how close the neuronal system (single unit or entire population) is to switching threshold. We can quantify this “nearness to switching” in terms of the altering eigenvalue structure: the dominant eigenvalue approaches zero, leading to a growth in correlated, low-frequency power, with exaggerated responsiveness to small perturbations, the responses becoming larger and slower as the neural population approaches its critical point–-this is critical slowing. In this chapter we discuss phase-transition predictions for both single-neuron and neural-population models, comparing theory with laboratory and clinical measurement

Interacting Turing-Hopf Instabilities Drive Symmetry-Breaking Transitions in a Mean-Field Model of the Cortex: A Mechanism for the Slow Oscillation

Electrical recordings of brain activity during the transition from wake to anesthetic coma show temporal and spectral alterations that are correlated with gross changes in the underlying brain state. Entry into anesthetic unconsciousness is signposted by the emergence of large, slow oscillations of electrical activity (≲1  Hz) similar to the slow waves observed in natural sleep. Here we present a two-dimensional mean-field model of the cortex in which slow spatiotemporal oscillations arise spontaneously through a Turing (spatial) symmetry-breaking bifurcation that is modulated by a Hopf (temporal) instability. In our model, populations of neurons are densely interlinked by chemical synapses, and by interneuronal gap junctions represented as an inhibitory diffusive coupling. To demonstrate cortical behavior over a wide range of distinct brain states, we explore model dynamics in the vicinity of a general-anesthetic-induced transition from “wake” to “coma.” In this region, the system is poised at a codimension-2 point where competing Turing and Hopf instabilities coexist. We model anesthesia as a moderate reduction in inhibitory diffusion, paired with an increase in inhibitory postsynaptic response, producing a coma state that is characterized by emergent low-frequency oscillations whose dynamics is chaotic in time and space. The effect of long-range axonal white-matter connectivity is probed with the inclusion of a single idealized point-to-point connection. We find that the additional excitation from the long-range connection can provoke seizurelike bursts of cortical activity when inhibitory diffusion is weak, but has little impact on an active cortex. Our proposed dynamic mechanism for the origin of anesthetic slow waves complements—and contrasts with—conventional explanations that require cyclic modulation of ion-channel conductances. We postulate that a similar bifurcation mechanism might underpin the slow waves of natural sleep and comment on the possible consequences of chaotic dynamics for memory processing and learning

Quantifying sleep architecture dynamics and individual differences using big data and Bayesian networks

The pattern of sleep stages across a night (sleep architecture) is influenced by biological, behavioral, and clinical variables. However, traditional measures of sleep architecture such as stage proportions, fail to capture sleep dynamics. Here we quantify the impact of individual differences on the dynamics of sleep architecture and determine which factors or set of factors best predict the next sleep stage from current stage information. We investigated the influence of age, sex, body mass index, time of day, and sleep time on static (e.g. minutes in stage, sleep efficiency) and dynamic measures of sleep architecture (e.g. transition probabilities and stage duration distributions) using a large dataset of 3202 nights from a non-clinical population. Multi-level regressions show that sex effects duration of all Non-Rapid Eye Movement (NREM) stages, and age has a curvilinear relationship for Wake After Sleep Onset (WASO) and slow wave sleep (SWS) minutes. Bayesian network modeling reveals sleep architecture depends on time of day, total sleep time, age and sex, but not BMI. Older adults, and particularly males, have shorter bouts (more fragmentation) of Stage 2, SWS, and they transition less frequently to these stages. Additionally, we showed that the next sleep stage and its duration can be optimally predicted by the prior 2 stages and age. Our results demonstrate the potential benefit of big data and Bayesian network approaches in quantifying static and dynamic architecture of normal sleep

Deep residual networks for automatic sleep stage classification of raw polysomnographic waveforms

We have developed an automatic sleep stage classification algorithm based on deep residual neural networks and raw polysomnogram signals. Briefly, the raw data is passed through 50 convolutional layers before subsequent classification into one of five sleep stages. Three model configurations were trained on 1850 polysomnogram recordings and subsequently tested on 230 independent recordings. Our best performing model yielded an accuracy of 84.1% and a Cohen's kappa of 0.746, improving on previous reported results by other groups also using only raw polysomnogram data. Most errors were made on non-REM stage 1 and 3 decisions, errors likely resulting from the definition of these stages. Further testing on independent cohorts is needed to verify performance for clinical use

Interacting Turing-Hopf Instabilities Drive Symmetry-Breaking Transitions in a Mean-Field Model of the Cortex: A Mechanism for the Slow Oscillation

Electrical recordings of brain activity during the transition from wake to anesthetic coma show temporal and spectral alterations that are correlated with gross changes in the underlying brain state. Entry into anesthetic unconsciousness is signposted by the emergence of large, slow oscillations of electrical activity (≲1  Hz) similar to the slow waves observed in natural sleep. Here we present a two-dimensional mean-field model of the cortex in which slow spatiotemporal oscillations arise spontaneously through a Turing (spatial) symmetry-breaking bifurcation that is modulated by a Hopf (temporal) instability. In our model, populations of neurons are densely interlinked by chemical synapses, and by interneuronal gap junctions represented as an inhibitory diffusive coupling. To demonstrate cortical behavior over a wide range of distinct brain states, we explore model dynamics in the vicinity of a general-anesthetic-induced transition from “wake” to “coma.” In this region, the system is poised at a codimension-2 point where competing Turing and Hopf instabilities coexist. We model anesthesia as a moderate reduction in inhibitory diffusion, paired with an increase in inhibitory postsynaptic response, producing a coma state that is characterized by emergent low-frequency oscillations whose dynamics is chaotic in time and space. The effect of long-range axonal white-matter connectivity is probed with the inclusion of a single idealized point-to-point connection. We find that the additional excitation from the long-range connection can provoke seizurelike bursts of cortical activity when inhibitory diffusion is weak, but has little impact on an active cortex. Our proposed dynamic mechanism for the origin of anesthetic slow waves complements—and contrasts with—conventional explanations that require cyclic modulation of ion-channel conductances. We postulate that a similar bifurcation mechanism might underpin the slow waves of natural sleep and comment on the possible consequences of chaotic dynamics for memory processing and learning