Rotating waves during human sleep spindles organize global patterns of activity that repeat precisely through the night

Brain states as wave-like motifs

There is ample evidence of wave-like activity in the brain at multiple scales and levels. This emerging literature supports the broader adoption of a wave perspective of brain activity. Specifically, a brain state can be described as a set of recurring, sequential patterns of propagating brain activity, namely a wave. We examine a collective body of experimental work investigating wave-like properties. Based on these works, we consider brain states as waves using a scale-agnostic framework across time and space. Emphasis is placed on the sequentiality and periodicity associated with brain activity. We conclude by discussing the implications, prospects, and experimental opportunities of this framework.

The direction of theta and alpha travelling waves modulates human memory processing

To support a range of behaviours, the brain must flexibly coordinate neural activity across widespread brain regions. One potential mechanism for this coordination is a travelling wave, in which a neural oscillation propagates across the brain while organizing the order and timing of activity across regions. Although travelling waves are present across the brain in various species, their potential functional relevance has remained unknown. Here, using rare direct human brain recordings, we demonstrate a distinct functional role for travelling waves of theta- and alpha-band (2-13 Hz) oscillations in the cortex. Travelling waves propagate in different directions during separate cognitive processes. In episodic memory, travelling waves tended to propagate in a posterior-to-anterior direction during successful memory encoding and in an anterior-to-posterior direction during recall. Because travelling waves of oscillations correspond to local neuronal spiking, these patterns indicate that rhythmic pulses of activity move across the brain in different directions for separate behaviours. More broadly, our results suggest a fundamental role for travelling waves and oscillations in dynamically coordinating neural connectivity, by flexibly organizing the timing and directionality of network interactions across the cortex to support cognition and behaviour.

Human connectome topology directs cortical traveling waves and shapes frequency gradients

Traveling waves and neural oscillation frequency gradients are pervasive in the human cortex. While the direction of traveling waves has been linked to brain function and dysfunction, the factors that determine this direction remain elusive. We hypothesized that structural connectivity instrength gradients - defined as the gradually varying sum of incoming connection strengths across the cortex - could shape both traveling wave direction and frequency gradients. We confirm the presence of instrength gradients in the human connectome across diverse cohorts and parcellations. Using a cortical network model, we demonstrate how these instrength gradients direct traveling waves and shape frequency gradients. Our model fits resting-state MEG functional connectivity best in a regime where instrength-directed traveling waves and frequency gradients emerge. We further show how structural subnetworks of the human connectome generate opposing wave directions and frequency gradients observed in the alpha and beta bands. Our findings suggest that structural connectivity instrength gradients affect both traveling wave direction and frequency gradients.

An electroencephalogram microdisplay to visualize neuronal activity on the brain surface

Functional mapping during brain surgery is applied to define brain areas that control critical functions and cannot be removed. Currently, these procedures rely on verbal interactions between the neurosurgeon and electrophysiologist, which can be time-consuming. In addition, the electrode grids that are used to measure brain activity and to identify the boundaries of pathological versus functional brain regions have low resolution and limited conformity to the brain surface. Here, we present the development of an intracranial electroencephalogram (iEEG)-microdisplay that consists of freestanding arrays of 2048 GaN light-emitting diodes laminated on the back of micro-electrocorticography electrode grids. With a series of proof-of-concept experiments in rats and pigs, we demonstrate that these iEEG-microdisplays allowed us to perform real-time iEEG recordings and display cortical activities by spatially corresponding light patterns on the surface of the brain in the surgical field. Furthermore, iEEG-microdisplays allowed us to identify and display cortical landmarks and pathological activities from rat and pig models. Using a dual-color iEEG-microdisplay, we demonstrated coregistration of the functional cortical boundaries with one color and displayed the evolution of electrical potentials associated with epileptiform activity with another color. The iEEG-microdisplay holds promise to facilitate monitoring of pathological brain activity in clinical settings.

Neural assemblies coordinated by cortical waves are associated with waking and hallucinatory brain states

The relationship between sensory stimuli and perceptions is brain-state dependent: in wakefulness, suprathreshold stimuli evoke perceptions; under anesthesia, perceptions are abolished; and during dreaming and in dissociated states, percepts are internally generated. Here, we exploit this state dependence to identify brain activity associated with internally generated or stimulus-evoked perceptions. In awake mice, visual stimuli phase reset spontaneous cortical waves to elicit 3-6 Hz feedback traveling waves. These stimulus-evoked waves traverse the cortex and entrain visual and parietal neurons. Under anesthesia as well as during ketamine-induced dissociation, visual stimuli do not disrupt spontaneous waves. Uniquely, in the dissociated state, spontaneous waves traverse the cortex caudally and entrain visual and parietal neurons, akin to stimulus-evoked waves in wakefulness. Thus, coordinated neuronal assemblies orchestrated by traveling cortical waves emerge in states in which perception can manifest. The awake state is privileged in that this coordination is reliably elicited by external visual stimuli.

Sequentially activated discrete modules appear as traveling waves in neuronal measurements with limited spatiotemporal sampling

Numerous studies have identified traveling waves in the cortex and suggested they play important roles in brain processing. These waves are most often measured using macroscopic methods that are unable to assess the local spiking activity underlying wave dynamics. Here, we investigated the possibility that waves may not be traveling at the single neuron scale. We first show that sequentially activating two discrete brain areas can appear as traveling waves in EEG simulations. We next reproduce these results using an analytical model of two sequentially activated regions. Using this model, we were able to generate wave-like activity with variable directions, velocities, and spatial patterns, and to map the discriminability limits between traveling waves and modular sequential activations. Finally, we investigated the link between field potentials and single neuron excitability using large-scale measurements from turtle cortex ex vivo. We found that while field potentials exhibit wave-like dynamics, the underlying spiking activity was better described by consecutively activated spatially adjacent groups of neurons. Taken together, this study suggests caution when interpreting phase delay measurements as continuously propagating wavefronts in two different spatial scales. A careful distinction between modular and wave excitability profiles across scales will be critical for understanding the nature of cortical computations.

Intraoperative application and early experience with novel high-resolution, high-channel-count thin-film electrodes for human microelectrocorticography

OBJECTIVE

The study objective was to evaluate intraoperative experience with newly developed high-spatial-resolution microelectrode grids composed of poly(3,4-ethylenedioxythiophene) with polystyrene sulfonate (PEDOT:PSS), and those composed of platinum nanorods (PtNRs).

METHODS

A cohort of patients who underwent craniotomy for pathological tissue resection and who had high-spatial-resolution microelectrode grids placed intraoperatively were evaluated. Patient demographic and baseline clinical variables as well as relevant microelectrode grid characteristic data were collected. The primary and secondary outcome measures of interest were successful microelectrode grid utilization with usable resting-state or task-related data, and grid-related adverse intraoperative events and/or grid dysfunction.

RESULTS

Included in the analysis were 89 cases of patients who underwent a craniotomy for resection of neoplasms (n = 58) or epileptogenic tissue (n = 31). These cases accounted for 94 grids: 58 PEDOT:PSS and 36 PtNR grids. Of these 94 grids, 86 were functional and used successfully to obtain cortical recordings from 82 patients. The mean cortical grid recording duration was 15.3 ± 1.15 minutes. Most recordings in patients were obtained during experimental tasks (n = 52, 58.4%), involving language and sensorimotor testing paradigms, or were obtained passively during resting state (n = 32, 36.0%). There were no intraoperative adverse events related to grid placement. However, there were instances of PtNR grid dysfunction (n = 8) related to damage incurred by suboptimal preoperative sterilization (n = 7) and improper handling (n = 1); intraoperative recordings were not performed. Vaporized peroxide sterilization was the most optimal sterilization method for PtNR grids, providing a significantly greater number of usable channels poststerilization than did steam-based sterilization techniques (median 905.0 [IQR 650.8-935.5] vs 356.0 [IQR 18.0-597.8], p = 0.0031).

CONCLUSIONS

High-spatial-resolution microelectrode grids can be readily incorporated into appropriately selected craniotomy cases for clinical and research purposes. Grids are reliable when preoperative handling and sterilization considerations are accounted for. Future investigations should compare the diagnostic utility of these high-resolution grids to commercially available counterparts and assess whether diagnostic discrepancies relate to clinical outcomes.

Frequency Analysis of EEG Microstate Sequences in Wakefulness and NREM Sleep

The majority of EEG microstate analyses concern wakefulness, and the existing sleep studies have focused on changes in spatial microstate properties and on microstate transitions between adjacent time points, the shortest available time scale. We present a more extensive time series analysis of unsmoothed EEG microstate sequences in wakefulness and non-REM sleep stages across many time scales. Very short time scales are assessed with Markov tests, intermediate time scales by the entropy rate and long time scales by a spectral analysis which identifies characteristic microstate frequencies. During the descent from wakefulness to sleep stage N3, we find that the increasing mean microstate duration is a gradual phenomenon explained by a continuous slowing of microstate dynamics as described by the relaxation time of the transition probability matrix. The finite entropy rate, which considers longer microstate histories, shows that microstate sequences become more predictable (less random) with decreasing vigilance level. Accordingly, the Markov property is absent in wakefulness but in sleep stage N3, 10/19 subjects have microstate sequences compatible with a second-order Markov process. A spectral microstate analysis is performed by comparing the time-lagged mutual information coefficients of microstate sequences with the autocorrelation function of the underlying EEG. We find periodic microstate behavior in all vigilance states, linked to alpha frequencies in wakefulness, theta activity in N1, sleep spindle frequencies in N2, and in the delta frequency band in N3. In summary, we show that EEG microstates are a dynamic phenomenon with oscillatory properties that slow down in sleep and are coupled to specific EEG frequencies across several sleep stages.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

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.

Identifying sources of human interictal discharges with travelling wave and white matter propagation

Interictal epileptiform discharges have been shown to propagate from focal epileptogenic sources as travelling waves or through more rapid white matter conduction. We hypothesize that both modes of propagation are necessary to explain interictal discharge timing delays. We propose a method that, for the first time, incorporates both propagation modes to identify unique potential sources of interictal activity. We retrospectively analysed 38 focal epilepsy patients who underwent intracranial EEG recordings and diffusion-weighted imaging for epilepsy surgery evaluation. Interictal discharges were detected and localized to the most likely source based on relative delays in time of arrival across electrodes, incorporating travelling waves and white matter propagation. We assessed the influence of white matter propagation on distance of spread, timing and clinical interpretation of interictal activity. To evaluate accuracy, we compared our source localization results to earliest spiking regions to predict seizure outcomes. White matter propagation helps to explain the timing delays observed in interictal discharge sequences, underlying rapid and distant propagation. Sources identified based on differences in time of receipt of interictal discharges are often distinct from the leading electrode location. Receipt of activity propagating rapidly via white matter can occur earlier than more local activity propagating via slower cortical travelling waves. In our cohort, our source localization approach was more accurate in predicting seizure outcomes than the leading electrode location. Inclusion of white matter in addition to travelling wave propagation in our model of discharge spread did not improve overall accuracy but allowed for identification of unique and at times distant potential sources of activity, particularly in patients with persistent postoperative seizures. Since distant white matter propagation can occur more rapidly than local travelling wave propagation, combined modes of propagation within an interictal discharge sequence can decouple the commonly assumed relationship between spike timing and distance from the source. Our findings thus highlight the clinical importance of recognizing the presence of dual modes of propagation during interictal discharges, as this may be a cause of clinical mislocalization.

How far neuroscience is from understanding brains

The cellular biology of brains is relatively well-understood, but neuroscientists have not yet generated a theory explaining how brains work. Explanations of how neurons collectively operate to produce what brains can do are tentative and incomplete. Without prior assumptions about the brain mechanisms, I attempt here to identify major obstacles to progress in neuroscientific understanding of brains and central nervous systems. Most of the obstacles to our understanding are conceptual. Neuroscience lacks concepts and models rooted in experimental results explaining how neurons interact at all scales. The cerebral cortex is thought to control awake activities, which contrasts with recent experimental results. There is ambiguity distinguishing task-related brain activities from spontaneous activities and organized intrinsic activities. Brains are regarded as driven by external and internal stimuli in contrast to their considerable autonomy. Experimental results are explained by sensory inputs, behavior, and psychological concepts. Time and space are regarded as mutually independent variables for spiking, post-synaptic events, and other measured variables, in contrast to experimental results. Dynamical systems theory and models describing evolution of variables with time as the independent variable are insufficient to account for central nervous system activities. Spatial dynamics may be a practical solution. The general hypothesis that measurements of changes in fundamental brain variables, action potentials, transmitter releases, post-synaptic transmembrane currents, etc., propagating in central nervous systems reveal how they work, carries no additional assumptions. Combinations of current techniques could reveal many aspects of spatial dynamics of spiking, post-synaptic processing, and plasticity in insects and rodents to start with. But problems defining baseline and reference conditions hinder interpretations of the results. Furthermore, the facts that pooling and averaging of data destroy their underlying dynamics imply that single-trial designs and statistics are necessary.

Beta traveling waves in monkey frontal and parietal areas encode recent reward history

Brain function depends on neural communication, but the mechanisms of this communication are not well understood. Recent studies suggest that one form of neural communication is through traveling waves (TWs)-patterns of neural oscillations that propagate within and between brain areas. We show that TWs are robust in microarray recordings in frontal and parietal cortex and encode recent reward history. Two adult male monkeys made saccades to obtain probabilistic rewards and were sensitive to the (statistically irrelevant) reward on the previous trial. TWs in frontal and parietal areas were stronger in trials that followed a prior reward versus a lack of reward and, in the frontal lobe, correlated with the monkeys' behavioral sensitivity to the prior reward. The findings suggest that neural communication mediated by TWs within the frontal and parietal lobes contribute to maintaining information about recent reward history and mediating the impact of this history on the monkeys' expectations.

Sleep and vigilance states: Embracing spatiotemporal dynamics

The classic view of sleep and vigilance states is a global stationary perspective driven by the interaction between neuromodulators and thalamocortical systems. However, recent data are challenging this view by demonstrating that vigilance states are highly dynamic and regionally complex. Spatially, sleep- and wake-like states often co-occur across distinct brain regions, as in unihemispheric sleep, local sleep in wakefulness, and during development. Temporally, dynamic switching prevails around state transitions, during extended wakefulness, and in fragmented sleep. This knowledge, together with methods monitoring brain activity across multiple regions simultaneously at millisecond resolution with cell-type specificity, is rapidly shifting how we consider vigilance states. A new perspective incorporating multiple spatial and temporal scales may have important implications for considering the governing neuromodulatory mechanisms, the functional roles of vigilance states, and their behavioral manifestations. A modular and dynamic view highlights novel avenues for finer spatiotemporal interventions to improve sleep function.

Interacting spiral wave patterns underlie complex brain dynamics and are related to cognitive processing

The large-scale activity of the human brain exhibits rich and complex patterns, but the spatiotemporal dynamics of these patterns and their functional roles in cognition remain unclear. Here by characterizing moment-by-moment fluctuations of human cortical functional magnetic resonance imaging signals, we show that spiral-like, rotational wave patterns (brain spirals) are widespread during both resting and cognitive task states. These brain spirals propagate across the cortex while rotating around their phase singularity centres, giving rise to spatiotemporal activity dynamics with non-stationary features. The properties of these brain spirals, such as their rotational directions and locations, are task relevant and can be used to classify different cognitive tasks. We also demonstrate that multiple, interacting brain spirals are involved in coordinating the correlated activations and de-activations of distributed functional regions; this mechanism enables flexible reconfiguration of task-driven activity flow between bottom-up and top-down directions during cognitive processing. Our findings suggest that brain spirals organize complex spatiotemporal dynamics of the human brain and have functional correlates to cognitive processing.

A CMOS-based highly scalable flexible neural electrode interface

Perception, thoughts, and actions are encoded by the coordinated activity of large neuronal populations spread over large areas. However, existing electrophysiological devices are limited by their scalability in capturing this cortex-wide activity. Here, we developed an electrode connector based on an ultra-conformable thin-film electrode array that self-assembles onto silicon microelectrode arrays enabling multithousand channel counts at a millimeter scale. The interconnects are formed using microfabricated electrode pads suspended by thin support arms, termed Flex2Chip. Capillary-assisted assembly drives the pads to deform toward the chip surface, and van der Waals forces maintain this deformation, establishing Ohmic contact. Flex2Chip arrays successfully measured extracellular action potentials ex vivo and resolved micrometer scale seizure propagation trajectories in epileptic mice. We find that seizure dynamics in absence epilepsy in the Scn8a+/- model do not have constant propagation trajectories.

Cortical activity emerges in region-specific patterns during early brain development

The development of precise neural circuits in the brain requires spontaneous patterns of neural activity prior to functional maturation. In the rodent cerebral cortex, patchwork and wave patterns of activity develop in somatosensory and visual regions, respectively, and are present at birth. However, whether such activity patterns occur in noneutherian mammals, as well as when and how they arise during development, remain open questions relevant for understanding brain formation in health and disease. Since the onset of patterned cortical activity is challenging to study prenatally in eutherians, here we offer an approach in a minimally invasive manner using marsupial dunnarts, whose cortex forms postnatally. We discovered similar patchwork and travelling waves in the dunnart somatosensory and visual cortices at stage 27 (equivalent to newborn mice) and examined earlier stages of development to determine the onset of these patterns and how they first emerge. We observed that these patterns of activity emerge in a region-specific and sequential manner, becoming evident as early as stage 24 in somatosensory and stage 25 in visual cortices (equivalent to embryonic day 16 and 17, respectively, in mice), as cortical layers establish and thalamic axons innervate the cortex. In addition to sculpting synaptic connections of existing circuits, evolutionarily conserved patterns of neural activity could therefore help regulate other early events in cortical development.

Cytoelectric Coupling: Electric fields sculpt neural activity and "tune" the brain's infrastructure

We propose and present converging evidence for the Cytoelectric Coupling Hypothesis: Electric fields generated by neurons are causal down to the level of the cytoskeleton. This could be achieved via electrodiffusion and mechanotransduction and exchanges between electrical, potential and chemical energy. Ephaptic coupling organizes neural activity, forming neural ensembles at the macroscale level. This information propagates to the neuron level, affecting spiking, and down to molecular level to stabilize the cytoskeleton, "tuning" it to process information more efficiently.

Characterization of spatiotemporal dynamics in EEG data during picture naming with optical flow patterns

In this study, we investigate the spatiotemporal dynamics of the neural oscillations by analyzing the electric potential that arises from neural activity. We identify two types of dynamics based on the frequency and phase of oscillations: standing waves or as out-of-phase and modulated waves, which represent a combination of standing and moving waves. To characterize these dynamics, we use optical flow patterns such as sources, sinks, spirals and saddles. We compare analytical and numerical solutions with real EEG data acquired during a picture-naming task. Analytical approximation of standing waves helps us to establish some properties of pattern location and number. Specifically, sources and sinks are mainly located in the same location, while saddles are positioned between them. The number of saddles correlates with the sum of all the other patterns. These properties are confirmed in both the simulated and real EEG data. In particular, source and sink clusters in the EEG data overlap with each other with median percentages around 60%, and hence have high spatial correlation, while source/sink clusters overlap with saddle clusters in less than 1%, and have different locations. Our statistical analysis showed that saddles account for about 45% of all patterns, while the remaining patterns are present in similar proportions.

Subsets of cortico-cortical evoked potentials propagate as traveling waves

Emerging evidence suggests that the temporal dynamics of cortico-cortical evoked potentials (CCEPs) may be used to characterize the patterns of information flow between and within brain networks. At present, however, the spatiotemporal dynamics of CCEP propagation cortically and subcortically are incompletely understood. We hypothesized that CCEPs propagate as an evoked traveling wave emanating from the site of stimulation. To elicit CCEPs, we applied single-pulse stimulation to stereoelectroencephalography (SEEG) electrodes implanted in 21 adult patients with intractable epilepsy. For each robust CCEP, we measured the timing of the maximal descent in evoked local field potentials and broadband high-gamma power (70-150 Hz) envelopes relative to the distance between the recording and stimulation contacts using three different metrics (i.e., Euclidean distance, path length, geodesic distance), representing direct, subcortical, and transcortical propagation, respectively. Many evoked responses to single-pulse electrical stimulation appear to propagate as traveling waves (~17-30%), even in the sparsely sampled, three-dimensional SEEG space. These results provide new insights into the spatiotemporal dynamics of CCEP propagation.

Constructing 2D maps of human spinal cord activity and isolating the functional midline with high-density microelectrode arrays

Intraoperative neuromonitoring (IONM) is a widely used practice in spine surgery for early detection and minimization of neurological injury. IONM is most commonly conducted by indirectly recording motor and somatosensory evoked potentials from either muscles or the scalp, which requires large-amplitude electrical stimulation and provides limited spatiotemporal information. IONM may inform of inadvertent events during neurosurgery after they occur, but it does not guide safe surgical procedures when the anatomy of the diseased spinal cord is distorted. To overcome these limitations and to increase our understanding of human spinal cord neurophysiology, we applied a microelectrode array with hundreds of channels to the exposed spinal cord during surgery and resolved spatiotemporal dynamics with high definition. We used this method to construct two-dimensional maps of responsive channels and define with submillimeter precision the electrophysiological midline of the spinal cord. The high sensitivity of our microelectrode array allowed us to record both epidural and subdural responses at stimulation currents that are well below those used clinically and to resolve postoperative evoked potentials when IONM could not. Together, these advances highlight the potential of our microelectrode arrays to capture previously unexplored spinal cord neural activity and its spatiotemporal dynamics at high resolution, offering better electrophysiological markers that can transform IONM.

Spatiotemporal patterns of sleep spindle activity in human anterior thalamus and cortex

Sleep spindles (8 - 16 Hz) are transient electrophysiological events during non-rapid eye movement sleep. While sleep spindles are routinely observed in the cortex using scalp electroencephalography (EEG), recordings of their thalamic counterparts have not been widely studied in humans. Based on a few existing studies, it has been hypothesized that spindles occur as largely local phenomena. We investigated intra-thalamic and thalamocortical spindle co-occurrence, which may underlie thalamocortical communication. We obtained scalp EEG and thalamic recordings from 7 patients that received bilateral deep brain stimulation (DBS) electrodes to the anterior thalamus for the treatment of drug resistant focal epilepsy. Spindles were categorized into subtypes based on their main frequency (i.e., slow (10±2 Hz) or fast (14±2 Hz)) and their level of thalamic involvement (spanning one channel, or spreading uni- or bilaterally within the thalamus). For the first time, we contrasted observed spindle patterns with permuted data to estimate random spindle co-occurrence. We found that multichannel spindle patterns were systematically coordinated at the thalamic and thalamocortical level. Importantly, distinct topographical patterns of thalamocortical spindle overlap were associated with slow and fast subtypes of spindles. These observations provide further evidence for coordinated spindle activity in thalamocortical networks.

Precision dynamical mapping using topological data analysis reveals a hub-like transition state at rest

In the absence of external stimuli, neural activity continuously evolves from one configuration to another. Whether these transitions or explorations follow some underlying arrangement or lack a predictable ordered plan remains to be determined. Here, using fMRI data from highly sampled individuals (~5 hours of resting-state data per individual), we aimed to reveal the rules that govern transitions in brain activity at rest. Our Topological Data Analysis based Mapper approach characterized a highly visited transition state of the brain that acts as a switch between different neural configurations to organize the spontaneous brain activity. Further, while the transition state was characterized by a uniform representation of canonical resting-state networks (RSNs), the periphery of the landscape was dominated by a subject-specific combination of RSNs. Altogether, we revealed rules or principles that organize spontaneous brain activity using a precision dynamics approach.

The space-time profiles of sleep spindles and their coordination with slow oscillations on the electrode manifold

Sleep spindles are important for sleep quality and cognitive functions, with their coordination with slow oscillations (SOs) potentially organizing cross-region reactivation of memory traces. Here, we describe the organization of spindles on the electrode manifold and their relation to SOs. We analyzed the sleep night EEG of 34 subjects and detected spindles and SOs separately at each electrode. We compared spindle properties (frequency, duration, and amplitude) in slow wave sleep (SWS) and Stage 2 sleep (S2); and in spindles that coordinate with SOs or are uncoupled. We identified different topographical spindle types using clustering analysis that grouped together spindles co-detected across electrodes within a short delay (±300 ms). We then analyzed the properties of spindles of each type, and coordination to SOs. We found that SWS spindles are shorter than S2 spindles, and spindles at frontal electrodes have higher frequencies in S2 compared to SWS. Furthermore, S2 spindles closely following an SO (about 10% of all spindles) show faster frequency, shorter duration, and larger amplitude than uncoupled ones. Clustering identified Global, Local, Posterior, Frontal-Right and Left spindle types. At centro-parietal locations, Posterior spindles show faster frequencies compared to other types. Furthermore, the infrequent SO-spindle complexes are preferentially recruiting Global SO waves coupled with fast Posterior spindles. Our results suggest a non-uniform participation of spindles to complexes, especially evident in S2. This suggests the possibility that different mechanisms could initiate an SO-spindle complex compared to SOs and spindles separately. This has implications for understanding the role of SOs-spindle complexes in memory reactivation.

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.

Spontaneous Spiking Is Governed by Broadband Fluctuations

Populations of cortical neurons generate rhythmic fluctuations in their ongoing spontaneous activity. These fluctuations can be seen in the local field potential (LFP), which reflects summed return currents from synaptic activity in the local population near a recording electrode. The LFP is spectrally broad, and many researchers view this breadth as containing many narrowband oscillatory components that may have distinct functional roles. This view is supported by the observation that the phase of narrowband oscillations is often correlated with cortical excitability and can relate to the timing of spiking activity and the fidelity of sensory evoked responses. Accordingly, researchers commonly tune in to these channels by narrowband filtering the LFP. Alternatively, neural activity may be fundamentally broadband and composed of transient, nonstationary rhythms that are difficult to approximate as oscillations. In this view, the instantaneous state of the broad ensemble relates directly to the excitability of the local population with no particular allegiance to any frequency band. To test between these alternatives, we asked whether the spiking activity of neocortical neurons in marmoset of either sex is better aligned with the phase of the LFP within narrow frequency bands or with a broadband measure. We find that the phase of broadband LFP fluctuations provides a better predictor of spike timing than the phase after filtering in narrow bands. These results challenge the view of the neocortex as a system composed of narrowband oscillators and supports a view in which neural activity fluctuations are intrinsically broadband.SIGNIFICANCE STATEMENT Research into the dynamical state of neural populations often attributes unique significance to the state of narrowband oscillatory components. However, rhythmic fluctuations in cortical activity are nonstationary and broad spectrum. We find that the timing of spontaneous spiking activity is better captured by the state of broadband fluctuations over any latent oscillatory component. These results suggest narrowband interpretations of rhythmic population activity may be limited, and broader representations may provide higher fidelity in describing moment-to-moment fluctuations in cortical activity.

Emerging principles of spacetime in brains: Meeting report on spatial neurodynamics

How do neurons and networks of neurons interact spatially? Here, we overview recent discoveries revealing how spatial dynamics of spiking and postsynaptic activity efficiently expose and explain fundamental brain and brainstem mechanisms behind detection, perception, learning, and behavior.

Human Spindle Variability

In humans, sleep spindles are 10- to 16-Hz oscillations lasting approximately 0.5-2 s. Spindles, along with cortical slow oscillations, may facilitate memory consolidation by enabling synaptic plasticity. Early recordings of spindles at the scalp found anterior channels had overall slower frequency than central-posterior channels. This robust, topographical finding led to dichotomizing spindles as "slow" versus "fast," modeled as two distinct spindle generators in frontal versus posterior cortex. Using a large dataset of intracranial stereoelectroencephalographic (sEEG) recordings from 20 patients (13 female, 7 male) and 365 bipolar recordings, we show that the difference in spindle frequency between frontal and parietal channels is comparable to the variability in spindle frequency within the course of individual spindles, across different spindles recorded by a given site, and across sites within a given region. Thus, fast and slow spindles only capture average differences that obscure a much larger underlying overlap in frequency. Furthermore, differences in mean frequency are only one of several ways that spindles differ. For example, compared with parietal, frontal spindles are smaller, tend to occur after parietal when both are engaged, and show a larger decrease in frequency within-spindles. However, frontal and parietal spindles are similar in being longer, less variable, and more widespread than occipital, temporal, and Rolandic spindles. These characteristics are accentuated in spindles which are highly phase-locked to posterior hippocampal spindles. We propose that rather than a strict parietal-fast/frontal-slow dichotomy, spindles differ continuously and quasi-independently in multiple dimensions, with variability due about equally to within-spindle, within-region, and between-region factors.SIGNIFICANCE STATEMENT Sleep spindles are 10- to 16-Hz neural oscillations generated by cortico-thalamic circuits that promote memory consolidation. Spindles are often dichotomized into slow-anterior and fast-posterior categories for cognitive and clinical studies. Here, we show that the anterior-posterior difference in spindle frequency is comparable to that observed between different cycles of individual spindles, between spindles from a given site, or from different sites within a region. Further, we show that spindles vary on other dimensions such as duration, amplitude, spread, primacy and consistency, and that these multiple dimensions vary continuously and largely independently across cortical regions. These findings suggest that multiple continuous variables rather than a strict frequency dichotomy may be more useful biomarkers for memory consolidation or psychiatric disorders.

Spontaneous neuronal oscillations in the human insula are hierarchically organized traveling waves

The insula plays a fundamental role in a wide range of adaptive human behaviors, but its electrophysiological dynamics are poorly understood. Here, we used human intracranial electroencephalographic recordings to investigate the electrophysiological properties and hierarchical organization of spontaneous neuronal oscillations within the insula. We analyzed the neuronal oscillations of the insula directly and found that rhythms in the theta and beta frequency oscillations are widespread and spontaneously present. These oscillations are largely organized along the anterior–posterior (AP) axis of the insula. Both the left and right insula showed anterior-­to-posterior decreasing gradients for the power of oscillations in the beta frequency band. The left insula also showed a posterior-to-anterior decreasing frequency gradient and an anterior-to-posterior decreasing power gradient in the theta frequency band. In addition to measuring the power of these oscillations, we also examined the phase of these signals across simultaneous recording channels and found that the insula oscillations in the theta and beta bands are traveling waves. The strength of the traveling waves in each frequency was positively correlated with the amplitude of each oscillation. However, the theta and beta traveling waves were uncoupled to each other in terms of phase and amplitude, which suggested that insular traveling waves in the theta and beta bands operate independently. Our findings provide new insights into the spatiotemporal dynamics and hierarchical organization of neuronal oscillations within the insula, which, given its rich connectivity with widespread cortical regions, indicates that oscillations and traveling waves have an important role in intrainsular and interinsular communications.

Nine-day continuous recording of EEG and 2-hour of high-density EEG under chronic sleep restriction in mice

This work provides an EEG dataset collected from nine mice during the sleep deprivation (SD) paradigm for the sleep science community. It includes 9-day of continuous recording of the frontal and parietal EEG, accelerometer, and 2-hour of high-density EEG (HD-EEG) under SD and SD-free conditions. Eighteen hours of SD were conducted on 5 consecutive days. The HD-EEG data were saved in the EEGLAB format and stored as the brain imaging data structure (BIDS). These datasets can be used to (i) compare mouse HD-EEG to human HD-EEG, (ii) track oscillatory activities of the sleep EEG (e.g., slow waves, spindles) across the cortical regions under different conditions of sleep pressure, and (iii) investigate the cortical traveling waves in the mouse brain. We also provided Python code for basic analyses of this dataset, including the detection of slow waves and sleep spindles. We hope that our dataset will reveal hidden activities during sleep and lead to a better understanding of the functions and mechanisms of sleep.

Brain-wide interactions during hippocampal sharp wave ripples

During periods of disengagement from the environment, transient population bursts, known as sharp wave ripples (SPW-Rs), occur sporadically. While numerous experiments have characterized the bidirectional relationship between SPW-Rs and activity in chosen brain areas, the topographic relationship between different segments of the hippocampus and brain-wide target areas has not been studied at high temporal and spatial resolution. Yet, such knowledge is necessary to infer the direction of communication. We analyzed two publicly available datasets with simultaneous high-density silicon probe recordings from across the mouse forebrain. We found that SPW-Rs coincide with a transient brain-wide increase in functional connectivity. In addition, we show that the diversity in SPW-R features, such as their incidence, magnitude, and intrahippocampal topography in the septotemporal axis, are correlated with slower excitability fluctuations in cortical and subcortical areas. Further, variations in SPW-R features correlated with the timing, sign, and magnitude of downstream responses with large-amplitude SPW-Rs followed by transient silence in extrahippocampal structures. Our findings expand on previous results and demonstrate that the activity patterns in extrahippocampal structures depend both on the intrahippocampal topographic origin and magnitude of hippocampal SPW-Rs.

Event-Related Phase Synchronization Propagates Rapidly across Human Ventral Visual Cortex

Visual inputs to early visual cortex integrate with semantic, linguistic and memory inputs in higher visual cortex, in a manner that is rapid and accurate, and enables complex computations such as face recognition and word reading. This implies the existence of fundamental organizational principles that enable such efficiency. To elaborate on this, we performed intracranial recordings in 82 individuals while they performed tasks of varying visual and cognitive complexity. We discovered that visual inputs induce highly organized posterior-to-anterior propagating patterns of phase modulation across the ventral occipitotemporal cortex. At individual electrodes there was a stereotyped temporal pattern of phase progression following both stimulus onset and offset, consistent across trials and tasks. The phase of low frequency activity in anterior regions was predicted by the prior phase in posterior cortical regions. This spatiotemporal propagation of phase likely serves as a feed-forward organizational influence enabling the integration of information across the ventral visual stream. This phase modulation manifests as the early components of the event related potential; one of the most commonly used measures in human electrophysiology. These findings illuminate fundamental organizational principles of the higher order visual system that enable the rapid recognition and characterization of a variety of inputs.

Propofol Anesthesia Alters Cortical Traveling Waves

Oscillatory dynamics in cortex seem to organize into traveling waves that serve a variety of functions. Recent studies show that propofol, a widely used anesthetic, dramatically alters cortical oscillations by increasing slow-delta oscillatory power and coherence. It is not known how this affects traveling waves. We compared traveling waves across the cortex of non-human primates before, during, and after propofol-induced loss of consciousness (LOC). After LOC, traveling waves in the slow-delta (∼1 Hz) range increased, grew more organized, and traveled in different directions relative to the awake state. Higher frequency (8-30 Hz) traveling waves, by contrast, decreased, lost structure, and switched to directions where the slow-delta waves were less frequent. The results suggest that LOC may be due, in part, to increases in the strength and direction of slow-delta traveling waves that, in turn, alter and disrupt traveling waves in the higher frequencies associated with cognition.

Slow waves form expanding, memory-rich mesostates steered by local excitability in fading anesthesia

In the arousal process, the brain restores its integrative activity from the synchronized state of slow wave activity (SWA). The mechanisms underpinning this state transition remain, however, to be elucidated. Here we simultaneously probed neuronal assemblies throughout the whole cortex with micro-electrocorticographic recordings in mice. We investigated the progressive shaping of propagating SWA at different levels of isoflurane. We found a form of memory of the wavefront shapes at deep anesthesia, tightly alternating posterior-anterior-posterior patterns. At low isoflurane, metastable patterns propagated in more directions, reflecting an increased complexity. The wandering across these mesostates progressively increased its randomness, as predicted by simulations of a network of spiking neurons, and confirmed in our experimental data. The complexity increase is explained by the elevated excitability of local assemblies with no modifications of the network connectivity. These results shed new light on the functional reorganization of the cortical network as anesthesia fades out.

•

Complexity of isoflurane-induced slow waves reliably determines anesthesia level

•

In deep anesthesia, the propagation strictly alternates between front-back-front patterns

•

In light anesthesia, there is a continuum of directions and faster propagation

•

Local excitability underpins the cortical reorganization in fading anesthesia

Cortical neural dynamics unveil the rhythm of natural visual behavior in marmosets

Numerous studies have shown that the visual system consists of functionally distinct ventral and dorsal streams; however, its exact spatial-temporal dynamics during natural visual behavior remain to be investigated. Here, we report cerebral neural dynamics during active visual exploration recorded by an electrocorticographic array covering the entire lateral surface of the marmoset cortex. We found that the dorsal stream was activated before the primary visual cortex with saccades and followed by the alteration of suppression and activation signals along the ventral stream. Similarly, the signal that propagated from the dorsal to ventral visual areas was accompanied by a travelling wave of low frequency oscillations. Such signal dynamics occurred at an average of 220 ms after saccades, which corresponded to the timing when whole-brain activation returned to background levels. We also demonstrated that saccades could occur at any point of signal flow, indicating the parallel computation of motor commands. Overall, this study reveals the neural dynamics of active vision, which are efficiently linked to the natural rhythms of visual exploration.

Human interictal epileptiform discharges are bidirectional traveling waves echoing ictal discharges

Interictal epileptiform discharges (IEDs), also known as interictal spikes, are large intermittent electrophysiological events observed between seizures in patients with epilepsy. Although they occur far more often than seizures, IEDs are less studied, and their relationship to seizures remains unclear. To better understand this relationship, we examined multi-day recordings of microelectrode arrays implanted in human epilepsy patients, allowing us to precisely observe the spatiotemporal propagation of IEDs, spontaneous seizures, and how they relate. These recordings showed that the majority of IEDs are traveling waves, traversing the same path as ictal discharges during seizures, and with a fixed direction relative to seizure propagation. Moreover, the majority of IEDs, like ictal discharges, were bidirectional, with one predominant and a second, less frequent antipodal direction. These results reveal a fundamental spatiotemporal similarity between IEDs and ictal discharges. These results also imply that most IEDs arise in brain tissue outside the site of seizure onset and propagate toward it, indicating that the propagation of IEDs provides useful information for localizing the seizure focus.

Human brain mapping with multithousand-channel PtNRGrids resolves spatiotemporal dynamics

Electrophysiological devices are critical for mapping eloquent and diseased brain regions and for therapeutic neuromodulation in clinical settings and are extensively used for research in brain-machine interfaces. However, the existing clinical and experimental devices are often limited in either spatial resolution or cortical coverage. Here, we developed scalable manufacturing processes with a dense electrical connection scheme to achieve reconfigurable thin-film, multithousand-channel neurophysiological recording grids using platinum nanorods (PtNRGrids). With PtNRGrids, we have achieved a multithousand-channel array of small (30 μm) contacts with low impedance, providing high spatial and temporal resolution over a large cortical area. We demonstrated that PtNRGrids can resolve submillimeter functional organization of the barrel cortex in anesthetized rats that captured the tissue structure. In the clinical setting, PtNRGrids resolved fine, complex temporal dynamics from the cortical surface in an awake human patient performing grasping tasks. In addition, the PtNRGrids identified the spatial spread and dynamics of epileptic discharges in a patient undergoing epilepsy surgery at 1-mm spatial resolution, including activity induced by direct electrical stimulation. Collectively, these findings demonstrated the power of the PtNRGrids to transform clinical mapping and research with brain-machine interfaces.

Local sleep: A new concept in brain plasticity

Traditionally, sleep and wakefulness have been considered as two global, mutually exclusive states. However, this view has been challenged by the discovery that sleep and wakefulness are actually locally regulated and that islands of these two states may often coexist in the same individual. Importantly, such a local regulation seems to be the key for many essential functions of sleep, including the maintenance of cognitive efficiency and the consolidation of new skills and memories. Indeed, local changes in sleep-related oscillations occur in brain areas that are used and involved in learning during wakefulness. In turn, these changes directly modulate experience-dependent brain adaptations and the consolidation of newly acquired memories. In line with these observations, alterations in the regional balance between wake- and sleep-like activity have been shown to accompany many pathologic conditions, including psychiatric and neurologic disorders. In the last decade, experimental research has started to shed light on the mechanisms involved in the local regulation of sleep and wakefulness. The results of this research have opened new avenues of investigation regarding the function of sleep and have revealed novel potential targets for the treatment of several pathologic conditions.

Cortical traveling waves reflect state-dependent hierarchical sequencing of local regions in the human connectome network

Recent human studies using electrocorticography have demonstrated that alpha and theta band oscillations form traveling waves on the cortical surface. According to neural synchronization theories, the cortical traveling waves may group local cortical regions and sequence them by phase synchronization; however these contributions have not yet been assessed. This study aimed to evaluate the functional contributions of traveling waves using connectome-based network modeling. In the simulation, we observed stable traveling waves on the entire cortical surface wherein the topographical pattern of these phases was substantially correlated with the empirically obtained resting-state networks, and local radial waves also appeared within the size of the empirical networks (< 50 mm). Importantly, individual regions in the entire network were instantaneously sequenced by their internal frequencies, and regions with higher intrinsic frequency were seen in the earlier phases of the traveling waves. Based on the communication-through-coherence theory, this phase configuration produced a hierarchical organization of each region by unidirectional communication between the arbitrarily paired regions. In conclusion, cortical traveling waves reflect the intrinsic frequency-dependent hierarchical sequencing of local regions, global traveling waves sequence the set of large-scale cortical networks, and local traveling waves sequence local regions within individual cortical networks.

Traveling waves in the prefrontal cortex during working memory

Neural oscillations are evident across cortex but their spatial structure is not well- explored. Are oscillations stationary or do they form "traveling waves", i.e., spatially organized patterns whose peaks and troughs move sequentially across cortex? Here, we show that oscillations in the prefrontal cortex (PFC) organized as traveling waves in the theta (4-8Hz), alpha (8-12Hz) and beta (12-30Hz) bands. Some traveling waves were planar but most rotated. The waves were modulated during performance of a working memory task. During baseline conditions, waves flowed bidirectionally along a specific axis of orientation. Waves in different frequency bands could travel in different directions. During task performance, there was an increase in waves in one direction over the other, especially in the beta band.

Whole-brain propagating patterns in human resting-state brain activities

Repetitive propagating activities in resting-state brain activities have been widely observed in various species and regions. Because they resemble the preceding brain activities during tasks, they are assumed to reflect past experiences embedded in neuronal circuits. "Whole-brain" propagating activities may also reflect a process that integrates information distributed over the entire brain, such as visual and motor information. Here we reveal whole-brain propagating activities from human resting-state magnetoencephalography (MEG) and electroencephalography (EEG) data. We simultaneously recorded the MEGs and EEGs and estimated the source currents from both measurements. Then using our recently proposed algorithm, we extracted repetitive spatiotemporal patterns from the source currents. The estimated patterns consisted of multiple frequency components, each of which transiently exhibited the frequency-specific resting-state networks (RSNs) of functional MRIs (fMRIs), such as the default mode and sensorimotor networks. A simulation test suggested that the spatiotemporal patterns reflected the phase alignment of the multiple frequency oscillators induced by the propagating activities along the anatomical connectivity. These results argue that whole-brain propagating activities transiently exhibited multiple RSNs in their multiple frequency components, suggesting that they reflected a process to integrate the information distributed over the frequencies and networks.

Decoding of cortex-wide brain activity from local recordings of neural potentials

Objective. Electrical recordings of neural activity from brain surface have been widely employed in basic neuroscience research and clinical practice for investigations of neural circuit functions, brain-computer interfaces, and treatments for neurological disorders. Traditionally, these surface potentials have been believed to mainly reflect local neural activity. It is not known how informative the locally recorded surface potentials are for the neural activities across multiple cortical regions.Approach. To investigate that, we perform simultaneous local electrical recording and wide-field calcium imaging in awake head-fixed mice. Using a recurrent neural network model, we try to decode the calcium fluorescence activity of multiple cortical regions from local electrical recordings.Main results. The mean activity of different cortical regions could be decoded from locally recorded surface potentials. Also, each frequency band of surface potentials differentially encodes activities from multiple cortical regions so that including all the frequency bands in the decoding model gives the highest decoding performance. Despite the close spacing between recording channels, surface potentials from different channels provide complementary information about the large-scale cortical activity and the decoding performance continues to improve as more channels are included. Finally, we demonstrate the successful decoding of whole dorsal cortex activity at pixel-level using locally recorded surface potentials.Significance. These results show that the locally recorded surface potentials indeed contain rich information of the large-scale neural activities, which could be further demixed to recover the neural activity across individual cortical regions. In the future, our cross-modality inference approach could be adapted to virtually reconstruct cortex-wide brain activity, greatly expanding the spatial reach of surface electrical recordings without increasing invasiveness. Furthermore, it could be used to facilitate imaging neural activity across the whole cortex in freely moving animals, without requirement of head-fixed microscopy configurations.

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.

Spontaneous traveling waves naturally emerge from horizontal fiber time delays and travel through locally asynchronous-irregular states

Studies of sensory-evoked neuronal responses often focus on mean spike rates, with fluctuations treated as internally-generated noise. However, fluctuations of spontaneous activity, often organized as traveling waves, shape stimulus-evoked responses and perceptual sensitivity. The mechanisms underlying these waves are unknown. Further, it is unclear whether waves are consistent with the low rate and weakly correlated “asynchronous-irregular” dynamics observed in cortical recordings. Here, we describe a large-scale computational model with topographically-organized connectivity and conduction delays relevant to biological scales. We find that spontaneous traveling waves are a general property of these networks. The traveling waves that occur in the model are sparse, with only a small fraction of neurons participating in any individual wave. Consequently, they do not induce measurable spike correlations and remain consistent with locally asynchronous irregular states. Further, by modulating local network state, they can shape responses to incoming inputs as observed in vivo.

Spontaneous traveling cortical waves shape neural responses. Using a large-scale computational model, the authors show that transmission delays shape locally asynchronous spiking dynamics into traveling waves without inducing correlations and boost responses to external input, as observed in vivo.

Replay in Deep Learning: Current Approaches and Missing Biological Elements

Replay is the reactivation of one or more neural patterns that are similar to the activation patterns experienced during past waking experiences. Replay was first observed in biological neural networks during sleep, and it is now thought to play a critical role in memory formation, retrieval, and consolidation. Replay-like mechanisms have been incorporated in deep artificial neural networks that learn over time to avoid catastrophic forgetting of previous knowledge. Replay algorithms have been successfully used in a wide range of deep learning methods within supervised, unsupervised, and reinforcement learning paradigms. In this letter, we provide the first comprehensive comparison between replay in the mammalian brain and replay in artificial neural networks. We identify multiple aspects of biological replay that are missing in deep learning systems and hypothesize how they could be used to improve artificial neural networks.

Brain Activity Fluctuations Propagate as Waves Traversing the Cortical Hierarchy

The brain exhibits highly organized patterns of spontaneous activity as measured by resting-state functional magnetic resonance imaging (fMRI) fluctuations that are being widely used to assess the brain's functional connectivity. Some evidence suggests that spatiotemporally coherent waves are a core feature of spontaneous activity that shapes functional connectivity, although this has been difficult to establish using fMRI given the temporal constraints of the hemodynamic signal. Here, we investigated the structure of spontaneous waves in human fMRI and monkey electrocorticography. In both species, we found clear, repeatable, and directionally constrained activity waves coursed along a spatial axis approximately representing cortical hierarchical organization. These cortical propagations were closely associated with activity changes in distinct subcortical structures, particularly those related to arousal regulation, and modulated across different states of vigilance. The findings demonstrate a neural origin of spatiotemporal fMRI wave propagation at rest and link it to the principal gradient of resting-state fMRI connectivity.

The impact of a closed-loop thalamocortical model on the spatiotemporal dynamics of cortical and thalamic traveling waves

Propagation of activity in spatially structured neuronal networks has been observed in awake, anesthetized, and sleeping brains. How these wave patterns emerge and organize across brain structures, and how network connectivity affects spatiotemporal neural activity remains unclear. Here, we develop a computational model of a two-dimensional thalamocortical network, which gives rise to emergent traveling waves similar to those observed experimentally. We illustrate how spontaneous and evoked oscillatory activity in space and time emerge using a closed-loop thalamocortical architecture, sustaining smooth waves in the cortex and staggered waves in the thalamus. We further show that intracortical and thalamocortical network connectivity, cortical excitation/inhibition balance, and thalamocortical or corticothalamic delay can independently or jointly change the spatiotemporal patterns (radial, planar and rotating waves) and characteristics (speed, direction, and frequency) of cortical and thalamic traveling waves. Computer simulations predict that increased thalamic inhibition induces slower cortical frequencies and that enhanced cortical excitation increases traveling wave speed and frequency. Overall, our results provide insight into the genesis and sustainability of thalamocortical spatiotemporal patterns, showing how simple synaptic alterations cause varied spontaneous and evoked wave patterns. Our model and simulations highlight the need for spatially spread neural recordings to uncover critical circuit mechanisms for brain functions.

Cross-participant prediction of vigilance stages through the combined use of wPLI and wSMI EEG functional connectivity metrics

Functional connectivity (FC) metrics describe brain inter-regional interactions and may complement information provided by common power-based analyses. Here, we investigated whether the FC-metrics weighted Phase Lag Index (wPLI) and weighted Symbolic Mutual Information (wSMI) may unveil functional differences across four stages of vigilance-wakefulness (W), NREM-N2, NREM-N3, and REM sleep-with respect to each other and to power-based features. Moreover, we explored their possible contribution in identifying differences between stages characterized by distinct levels of consciousness (REM+W vs. N2+N3) or sensory disconnection (REM vs. W). Overnight sleep and resting-state wakefulness recordings from 24 healthy participants (27 ± 6 years, 13F) were analyzed to extract power and FC-based features in six classical frequency bands. Cross-validated linear discriminant analyses (LDA) were applied to investigate the ability of extracted features to discriminate (1) the four vigilance stages, (2) W+REM vs. N2+N3, and (3) W vs. REM. For the four-way vigilance stages classification, combining features based on power and both connectivity metrics significantly increased accuracy relative to considering only power, wPLI, or wSMI features. Delta-power and connectivity (0.5-4 Hz) represented the most relevant features for all the tested classifications, in line with a possible involvement of slow waves in consciousness and sensory disconnection. Sigma-FC, but not sigma-power (12-16 Hz), was found to strongly contribute to the differentiation between states characterized by higher (W+REM) and lower (N2+N3) probabilities of conscious experiences. Finally, alpha-FC resulted as the most relevant FC-feature for distinguishing among wakefulness and REM sleep and may thus reflect the level of disconnection from the external environment.

Mean-Field Models for EEG/MEG: From Oscillations to Waves

Neural mass models have been used since the 1970s to model the coarse-grained activity of large populations of neurons. They have proven especially fruitful for understanding brain rhythms. However, although motivated by neurobiological considerations they are phenomenological in nature, and cannot hope to recreate some of the rich repertoire of responses seen in real neuronal tissue. Here we consider a simple spiking neuron network model that has recently been shown to admit an exact mean-field description for both synaptic and gap-junction interactions. The mean-field model takes a similar form to a standard neural mass model, with an additional dynamical equation to describe the evolution of within-population synchrony. As well as reviewing the origins of this next generation mass model we discuss its extension to describe an idealised spatially extended planar cortex. To emphasise the usefulness of this model for EEG/MEG modelling we show how it can be used to uncover the role of local gap-junction coupling in shaping large scale synaptic waves.