Deterministic and Stochastic Components of Cortical Down States: Dynamics and Modulation.

Slow oscillations are an emergent activity of the cerebral cortex network consisting of alternating periods of activity (Up states) and silence (Down states). Up states are periods of persistent cortical activity that share properties with that of underlying wakefulness. However, the occurrence of Down states is almost invariably associated with unconsciousness, both in animal models and clinical studies. Down states have been attributed relevant functions, such as being a resetting mechanism or breaking causal interactions between cortical areas. But what do Down states consist of? Here, we explored in detail the network dynamics (e.g., synchronization and phase) during these silent periods in vivo (male mice), in vitro (ferrets, either sex), and in silico, investigating various experimental conditions that modulate them: anesthesia levels, excitability (electric fields), and excitation/inhibition balance. We identified metastability as two complementary phases composing such quiescence states: a highly synchronized "deterministic" period followed by a low-synchronization "stochastic" period. The balance between these two phases determines the dynamical properties of the resulting rhythm, as well as the responsiveness to incoming inputs or refractoriness. We propose detailed Up and Down state cycle dynamics that bridge cortical properties emerging at the mesoscale with their underlying mechanisms at the microscale, providing a key to understanding unconscious states.SIGNIFICANCE STATEMENT The cerebral cortex expresses slow oscillations consisting of Up (active) and Down (silent) states. Such activity emerges not only in slow wave sleep, but also under anesthesia and in brain lesions. Down states functionally disconnect the network, and are associated with unconsciousness. Based on a large collection of data, novel data analysis approaches and computational modeling, we thoroughly investigate the nature of Down states. We identify two phases: a highly synchronized "deterministic" period, followed by a low-synchronization "stochastic" period. The balance between these two phases determines the dynamic properties of the resulting rhythm and responsiveness to incoming inputs. This finding reconciles different theories of slow rhythm generation and provides clues about how the brain switches from conscious to unconscious brain states.

Impact of GABAA and GABAB Inhibition on Cortical Dynamics and Perturbational Complexity during Synchronous and Desynchronized States.

Quantitative estimations of spatiotemporal complexity of cortical activity patterns are used in the clinic as a measure of consciousness levels, but the cortical mechanisms involved are not fully understood. We used a version of the perturbational complexity index (PCI) adapted to multisite recordings from the ferret (either sex) cerebral cortex in vitro (sPCI) to investigate the role of GABAergic inhibition in cortical complexity. We studied two dynamical states: slow-wave activity (synchronous state) and desynchronized activity, that express low and high causal complexity respectively. Progressive blockade of GABAergic inhibition during both regimes revealed its impact on the emergent cortical activity and on sPCI. Gradual GABAA receptor blockade resulted in higher synchronization, being able to drive the network from a desynchronized to a synchronous state, with a progressive decrease of complexity (sPCI). Blocking GABAB receptors also resulted in a reduced sPCI, in particular when in a synchronous, slow wave state. Our findings demonstrate that physiological levels of inhibition contribute to the generation of dynamical richness and spatiotemporal complexity. However, if inhibition is diminished or enhanced, cortical complexity decreases. Using a computational model, we explored a larger parameter space in this relationship and demonstrate a link between excitatory/inhibitory balance and the complexity expressed by the cortical network.SIGNIFICANCE STATEMENT The spatiotemporal complexity of the activity expressed by the cerebral cortex is a highly revealing feature of the underlying network's state. Complexity varies with physiological brain states: it is higher during awake than during sleep states. But it also informs about pathologic states: in disorders of consciousness, complexity is lower in an unresponsive wakefulness syndrome than in a minimally conscious state. What are the network parameters that modulate complexity? Here we investigate how inhibition, mediated by either GABAA or GABAA receptors, influences cortical complexity. And we do this departing from two extreme functional states: a highly synchronous, slow-wave state, and a desynchronized one that mimics wakefulness. We find that there is an optimal level of inhibition in which complexity is highest.

Modulation of cortical slow oscillations and complexity across anesthesia levels.

The ability of different groups of cortical neurons to engage in causal interactions that are at once differentiated and integrated results in complex dynamic patterns. Complexity is low during periods of unconsciousness (deep sleep, anesthesia, unresponsive wakefulness syndrome) in which the brain tends to generate a stereotypical pattern consisting of alternating active and silent periods of neural activity-slow oscillations- and is high during wakefulness. But how is cortical complexity built up? Is it a continuum? An open question is whether cortical complexity can vary within the same brain state. Here we recorded with 32-channel multielectrode arrays from the cortical surface of the mouse and used both spontaneous dynamics (wave propagation entropy and functional complexity) and a perturbational approach (a variation of the perturbation complexity index) to measure complexity at different anesthesia levels. Variations in anesthesia level within the bistable regime of slow oscillations (0.1-1.5 Hz) resulted in a modulation of the slow oscillation frequency. Both perturbational and spontaneous complexity increased with decreasing anesthesia levels, in correlation with the decrease in coherence of the underlying network. Changes in complexity level are related to, but not dependent on, changes in excitability. We conclude that cortical complexity can vary within a single brain state dominated by slow oscillations, building up to the higher complexity associated with consciousness.

High-resolution mapping of infraslow cortical brain activity enabled by graphene microtransistors.

Recording infraslow brain signals (<0.1 Hz) with microelectrodes is severely hampered by current microelectrode materials, primarily due to limitations resulting from voltage drift and high electrode impedance. Hence, most recording systems include high-pass filters that solve saturation issues but come hand in hand with loss of physiological and pathological information. In this work, we use flexible epicortical and intracortical arrays of graphene solution-gated field-effect transistors (gSGFETs) to map cortical spreading depression in rats and demonstrate that gSGFETs are able to record, with high fidelity, infraslow signals together with signals in the typical local field potential bandwidth. The wide recording bandwidth results from the direct field-effect coupling of the active transistor, in contrast to standard passive electrodes, as well as from the electrochemical inertness of graphene. Taking advantage of such functionality, we envision broad applications of gSGFET technology for monitoring infraslow brain activity both in research and in the clinic.

Quantifying the magnitude of pharyngeal obstruction during sleep using airflow shape.

RATIONALE AND OBJECTIVES: Non-invasive quantification of the severity of pharyngeal airflow obstruction would enable recognition of obstructive versus central manifestation of sleep apnoea, and identification of symptomatic individuals with severe airflow obstruction despite a low apnoea-hypopnoea index (AHI). Here we provide a novel method that uses simple airflow-versus-time ("shape") features from individual breaths on an overnight sleep study to automatically and non-invasively quantify the severity of airflow obstruction without oesophageal catheterisation. METHODS: 41 individuals with suspected/diagnosed obstructive sleep apnoea (AHI range 0-91 events·h-1) underwent overnight polysomnography with gold-standard measures of airflow (oronasal pneumotach: "flow") and ventilatory drive (calibrated intraoesophageal diaphragm electromyogram: "drive"). Obstruction severity was defined as a continuous variable (flow:drive ratio). Multivariable regression used airflow shape features (inspiratory/expiratory timing, flatness, scooping, fluttering) to estimate flow:drive ratio in 136 264 breaths (performance based on leave-one-patient-out cross-validation). Analysis was repeated using simultaneous nasal pressure recordings in a subset (n=17). RESULTS: Gold-standard obstruction severity (flow:drive ratio) varied widely across individuals independently of AHI. A multivariable model (25 features) estimated obstruction severity breath-by-breath (R2=0.58 versus gold-standard, p<0.00001; mean absolute error 22%) and the median obstruction severity across individual patients (R2=0.69, p<0.00001; error 10%). Similar performance was achieved using nasal pressure. CONCLUSIONS: The severity of pharyngeal obstruction can be quantified non-invasively using readily available airflow shape information. Our work overcomes a major hurdle necessary for the recognition and phenotyping of patients with obstructive sleep disordered breathing.

The temporal asymmetry of cortical dynamics as a signature of brain states.

The brain is a complex non-equilibrium system capable of expressing many different dynamics as well as the transitions between them. We hypothesized that the level of non-equilibrium can serve as a signature of a given brain state, which was quantified using the arrow of time (the level of irreversibility). Using this thermodynamic framework, the irreversibility of emergent cortical activity was quantified from local field potential recordings in male Lister-hooded rats at different anesthesia levels and during the sleep-wake cycle. This measure was carried out on five distinct brain states: slow-wave sleep, awake, deep anesthesia-slow waves, light anesthesia-slow waves, and microarousals. Low levels of irreversibility were associated with synchronous activity found both in deep anesthesia and slow-wave sleep states, suggesting that slow waves were the state closest to the thermodynamic equilibrium (maximum symmetry), thus requiring minimum energy. Higher levels of irreversibility were found when brain dynamics became more asynchronous, for example, in wakefulness. These changes were also reflected in the hierarchy of cortical dynamics across different cortical areas. The neural dynamics associated with different brain states were characterized by different degrees of irreversibility and hierarchy, also acting as markers of brain state transitions. This could open new routes to monitoring, controlling, and even changing brain states in health and disease.

Energy-Based Hierarchical Clustering of Cortical Slow Waves in Multi-Electrode Recordings.

The recent development of novel multi-electrode recording technologies has revealed the existence of traveling patterns of cortical activity in many species and under different states of awareness. Among these, slow activation waves occurring under sleep and anesthesia have been widely investigated as they provide unique insights into network features such as excitability, connectivity, structure, and dynamics of the cerebral cortex. Such characterization is usually based on clustering methods which are constrained by a priori assumptions as to the number of clusters to be used or rely on wave-by-wave pattern reconstruction. Here, we introduce a new computational tool based on modal analysis of fluid flows which is robustly applied to multivariate electrophysiological data from cortical networks, namely the Energy-based Hierarchical Waves Clustering method (EHWC). EHWC is composed of three main steps: (1) detecting the occurrence of global waves; (2) reducing the data dimensionality via singular value decomposition; (3) clustering hierarchically the singled-out waves. The analysis does not require the single-channel contribution to the waves, which is a typical bottleneck in this kind of analysis due to the unavoidable intrinsic variability of locally recorded activity. For testing and validation, here we used in vivo extracellular recordings from mice cortex under three different levels of anesthesia. As a result, we found slow waves with an increasing number of propagation modes as the anesthesia level decreases, giving an estimate of the increasing complexity of network dynamics. This and other wave's features replicate and extend the findings from previous literature, paving the way to extend the same approach to non-invasive electrophysiological recordings like EEG and fMRI used clinically for the characterization of brain dynamics and clinical stratification in brain lesions.

Altered Neocortical Dynamics in a Mouse Model of Williams-Beuren Syndrome.

Williams-Beuren syndrome (WBS) is a rare neurodevelopmental disorder characterized by moderate intellectual disability and learning difficulties alongside behavioral abnormalities such as hypersociability. Several structural and functional brain alterations are characteristic of this syndrome, as well as disturbed sleep and sleeping patterns. However, the detailed physiological mechanisms underlying WBS are mostly unknown. Here, we characterized the cortical dynamics in a mouse model of WBS previously reported to replicate most of the behavioral alterations described in humans. We recorded the laminar local field potential generated in the frontal cortex during deep anesthesia and characterized the properties of the emergent slow oscillation activity. Moreover, we performed micro-electrocorticogram recordings using multielectrode arrays covering the cortical surface of one hemisphere. We found significant differences between the cortical emergent activity and functional connectivity between wild-type mice and WBS model mice. Slow oscillations displayed Up states with diminished firing rate and lower high-frequency content in the gamma range. Lower firing rates were also recorded in the awake WBS animals while performing a marble burying task and could be associated with the decreased spine density and thus synaptic connectivity in this cortical area. We also found an overall increase in functional connectivity between brain areas, reflected in lower clustering and abnormally high integration, especially in the gamma range. These results expand previous findings in humans, suggesting that the cognitive deficits characterizing WBS might be associated with reduced excitability, plus an imbalance in the capacity to functionally integrate and segregate information.

Validating an Algorithm for Automatic Scoring of Inspiratory Flow Limitation Within a Range of Recording Settings.

Inspiratory Flow Limitation (IFL) is a phenomenon associated with narrowing of the upper airway, preventing an increase in inspiratory airflow despite an elevation in intrathoracic pressure. It has been shown that quantification of IFL might complement information provided by standard indices such as the apnea-hypopnea index (AHI) in characterizing sleep disordered breathing and identifying subclinical disease. Defining guidelines for visual scoring of IFL has been of increasing interest, and automated methods are desirable to avoid inter-scorer variability and allow analysis of large datasets. In addition, as recording instrumentation and practices may vary across hospitals and laboratories, it is useful to assess the influence of the recording parameters on the accuracy of the automated classification. We employed nasal pressure signals recorded as part of polysomnography (PSG) studies in 7 patients. Two experts independently classified approximately 2000 breaths per subject as IFL or non-IFL, and we used the consensus scoring as the gold standard. For each breath, we derived features indicative of the shape and frequency content of the signals and used them to train and validate a Support Vector Machine (SVM) to distinguish IFL from non-IFL breaths. We also assessed the effect of signal filtering (down-sampling and baseline-removal) on classification performance. The performance of the classifier was excellent (accuracy ~93%) for the raw signals (collected at 125 Hz with no filtering), and decreased for increasing high-pass cut-off frequencies (fc = [0.05, 0.1, 0.15, 0.2] Hz) down to 84% for fc= 0.2 Hz and for decreasing sampling rate (fs = [20, 50, 75, 100] Hz) down to ~85% for fs=20 Hz. Loss of performance was minimized when the classifier was re-trained using data with matched filtering characteristics (accuracy > 89%). We can conclude that the SVM feature-based algorithm provides a reliable and efficient tool for breath-by-breath classification.