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

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.

Electrophysiological indicators of sleep-associated memory consolidation in 5- to 6-year-old children

In adults, the synchronized interplay of sleep spindles (SP) and slow oscillations (SO) supports memory consolidation. Given tremendous developmental changes in SP and SO morphology, it remains elusive whether across childhood the same mechanisms as identified in adults are functional. Based on topography and frequency, we characterize slow and fast SPs and their temporal coupling to SOs in 24 pre-school children. Further, we ask whether slow and fast SPs and their modulation during SOs are associated with behavioral indicators of declarative memory consolidation as suggested by the literature on adults. Employing an individually tailored approach, we reliably identify an inherent, development-specific fast centro-parietal SP type, nested in the adult-like slow SP frequency range, along with a dominant slow frontal SP type. Further, we provide evidence that the modulation of fast centro-parietal SPs during SOs is already present in pre-school children. However, the temporal coordination between fast centro-parietal SPs and SOs is weaker and less precise than expected from research on adults. While we do not find evidence for a critical contribution of SP-SO coupling for memory consolidation, crucially, slow frontal and fast centro-parietal SPs are each differentially related to sleep-associated consolidation of items of varying quality. Whereas a higher number of slow frontal SPs is associated with stronger maintenance of medium-quality memories, a higher number of fast centro-parietal SPs is linked to a greater gain of low-quality items. Our results demonstrate two functionally relevant inherent SP types in pre-school children although SP-SO coupling is not yet fully mature.

Travelling spindles create necessary conditions for spike-timing-dependent plasticity in humans

Sleep spindles facilitate memory consolidation in the cortex during mammalian non-rapid eye movement sleep. In rodents, phase-locked firing during spindles may facilitate spike-timing-dependent plasticity by grouping pre-then-post-synaptic cell firing within ~25 ms. Currently, microphysiological evidence in humans for conditions conducive for spike-timing-dependent plasticity during spindles is absent. Here, we analyze field potentials and unit firing from middle/upper layers during spindles from 10 × 10 microelectrode arrays at 400 μm pitch in humans. We report strong tonic and phase-locked increases in firing and co-firing within 25 ms during spindles, especially those co-occurring with down-to-upstate transitions. Co-firing, spindle co-occurrence, and spindle coherence are greatest within ~2 mm, and high co-firing of units on different contacts depends on high spindle coherence between those contacts. Spindles propagate at ~0.28 m/s in distinct patterns, with correlated cell co-firing sequences. Spindles hence organize spatiotemporal patterns of neuronal co-firing in ways that may provide pre-conditions for plasticity during non-rapid eye movement sleep.

Sleep spindles during non-rapid eye movement are important for memory consolidation and require specific neuronal firing conditions in non-human mammals. Here, the authors show these conditions are present in humans, potentially facilitating spike-timing-dependent plasticity.

Traveling waves in non-local pulse-coupled networks

Traveling phase waves are commonly observed in recordings of the cerebral cortex and are believed to organize behavior across different areas of the brain. We use this as motivation to analyze a one-dimensional network of phase oscillators that are nonlocally coupled via the phase response curve (PRC) and the Dirac delta function. Existence of waves is proven and the dispersion relation is computed. Using the theory of distributions enables us to write and solve an associated stability problem. First and second order perturbation theory is applied to get analytic insight and we show that long waves are stable while short waves are unstable. We apply the results to PRCs that come from mitral neurons. We extend the results to smooth pulse-like coupling by reducing the nonlocal equation to a local one and solving the associated boundary value problem.

Cross-Frequency Power-Power Coupling Analysis: A Useful Cross-Frequency Measure to Classify ICA-Decomposed EEG

Magneto-/Electro-encephalography (M/EEG) commonly uses (fast) Fourier transformation to compute power spectral density (PSD). However, the resulting PSD plot lacks temporal information, making interpretation sometimes equivocal. For example, consider two different PSDs: a central parietal EEG PSD with twin peaks at 10 Hz and 20 Hz and a central parietal PSD with twin peaks at 10 Hz and 50 Hz. We can assume the first PSD shows a mu rhythm and the second harmonic; however, the latter PSD likely shows an alpha peak and an independent line noise. Without prior knowledge, however, the PSD alone cannot distinguish between the two cases. To address this limitation of PSD, we propose using cross-frequency power-power coupling (PPC) as a post-processing of independent component (IC) analysis (ICA) to distinguish brain components from muscle and environmental artifact sources. We conclude that post-ICA PPC analysis could serve as a new data-driven EEG classifier in M/EEG studies. For the reader's convenience, we offer a brief literature overview on the disparate use of PPC. The proposed cross-frequency power-power coupling analysis toolbox (PowPowCAT) is a free, open-source toolbox, which works as an EEGLAB extension.

The unreasonable effectiveness of deep learning in artificial intelligence

Deep learning networks have been trained to recognize speech, caption photographs, and translate text between languages at high levels of performance. Although applications of deep learning networks to real-world problems have become ubiquitous, our understanding of why they are so effective is lacking. These empirical results should not be possible according to sample complexity in statistics and nonconvex optimization theory. However, paradoxes in the training and effectiveness of deep learning networks are being investigated and insights are being found in the geometry of high-dimensional spaces. A mathematical theory of deep learning would illuminate how they function, allow us to assess the strengths and weaknesses of different network architectures, and lead to major improvements. Deep learning has provided natural ways for humans to communicate with digital devices and is foundational for building artificial general intelligence. Deep learning was inspired by the architecture of the cerebral cortex and insights into autonomy and general intelligence may be found in other brain regions that are essential for planning and survival, but major breakthroughs will be needed to achieve these goals.

Turning the Stimulus On and Off Changes the Direction of α Traveling Waves

Traveling waves have been studied to characterize the complex spatiotemporal dynamics of the brain. Several studies have suggested that the propagation direction of α traveling waves can be task dependent. For example, a recent electroencephalography (EEG) study from our group found that forward waves (i.e., occipital to frontal, FW waves) were observed during visual processing, whereas backward waves (i.e., frontal to occipital, BW waves) mostly occurred in the absence of sensory input. These EEG recordings, however, were obtained from different experimental sessions and different groups of subjects. To further examine how the waves’ direction changes between task conditions, 13 human participants were tested on a target detection task while EEG signals were recorded simultaneously. We alternated visual stimulation (5-s display of visual luminance sequences) and resting state (5 s of black screen) within each single trial, allowing us to monitor the moment-to-moment progression of traveling waves. As expected, the direction of α waves was closely linked with task conditions. First, FW waves from occipital to frontal regions, absent during rest, emerged as a result of visual processing, while BW waves in the opposite direction dominated in the absence of visual inputs, and were reduced (but not eliminated) by external visual inputs. Second, during visual stimulation (but not rest), both waves coexisted on average, but were negatively correlated. In summary, we conclude that the functional role of α traveling waves is closely related with their propagating direction, with stimulus-evoked FW waves supporting visual processing and spontaneous BW waves involved more in top-down control.

Brain Waves: Emergence of Localized, Persistent, Weakly Evanescent Cortical Loops

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of nonevanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full-brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the commonly agreed paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation because of the cortex wave modes in the highly folded areas will exhibit no apparent correlation with the fiber directions. Nonlinear coupling of those linear weakly evanescent wave modes then provides a universal mechanism for the emergence of synchronized brain wave field activity. The resonant and nonresonant terms of nonlinear coupling between multiple modes produce both synchronous spiking-like high-frequency wave activity as well as low-frequency wave rhythms. Numerical simulation of forced multiple-mode dynamics shows that, as forcing increases, there is a transition from damped to oscillatory regime that can then transition quickly to a nonoscillatory state when a critical excitation threshold is reached. The resonant nonlinear coupling results in the emergence of low-frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these weakly evanescent cortical wave modes have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and magnetoencephalography, and for the understanding of brain activity in general, including mechanisms of memory.

Language prediction mechanisms in human auditory cortex

Spoken language, both perception and production, is thought to be facilitated by an ensemble of predictive mechanisms. We obtain intracranial recordings in 37 patients using depth probes implanted along the anteroposterior extent of the supratemporal plane during rhythm listening, speech perception, and speech production. These reveal two predictive mechanisms in early auditory cortex with distinct anatomical and functional characteristics. The first, localized to bilateral Heschl's gyri and indexed by low-frequency phase, predicts the timing of acoustic events. The second, localized to planum temporale only in language-dominant cortex and indexed by high-gamma power, shows a transient response to acoustic stimuli that is uniquely suppressed during speech production. Chronometric stimulation of Heschl's gyrus selectively disrupts speech perception, while stimulation of planum temporale selectively disrupts speech production. This work illuminates the fundamental acoustic infrastructure-both architecture and function-for spoken language, grounding cognitive models of speech perception and production in human neurobiology.

Coupling of gamma band activity to sleep spindle oscillations - a combined EEG/MEG study

Sleep spindles are crucial to memory consolidation. Cortical gamma oscillations (30-100 Hz) are considered to reflect processing of memory in local cortical networks. The temporal and regulatory relationship between spindles and gamma activity might therefore provide clues into how sleep strengthens cortical memory representations. Here, combining EEG with MEG recordings during sleep in healthy humans (n = 12), we investigated the temporal relationships of cortical gamma band activity, always measured by MEG, during fast (12-16 Hz) and slow (8-12 Hz) sleep spindles detected in the EEG or MEG. Time-frequency distributions did not show a consistent coupling of gamma to the spindle oscillation, although activity in the low gamma (30-40 Hz) and neighboring beta range (<30 Hz) was generally increased during spindles. However, more fine-grained analyses of cross-frequency interactions revealed that both low and high gamma power (30-100 Hz) was coupled to the phase of slow and fast EEG spindles, importantly, with this coupling at a fixed phase only for the oscillations within an individual spindle, but with variable phase across spindles. We did not observe any coupling of gamma activity for spindles detected solely in the MEG and not in parallel EEG recordings, raising the possibility that these are more local spindles of different quality. Similar to fast spindle activity, low gamma band power followed a ~0.025 Hz infraslow rhythm during sleep whose frequency, however, was significantly faster than that of spindle activity. Our findings suggest a general function of fast and slow spindles that by spanning larger cortical networks might serve to synchronize gamma band activity occurring in more local but distributed networks. Thereby, spindles might help linking local memory processing between distributed networks.

Spontaneous travelling cortical waves gate perception in behaving primates

Perceptual sensitivity varies from moment to moment. One potential source of variability is spontaneous fluctuations in cortical activity that can travel as a wave¹. Spontaneous traveling waves have been reported during anesthesia^2–7, but it is not known whether spontaneous traveling waves play a role during waking perception. Using newly developed analytic techniques to characterize the moment-to-moment dynamics of noisy multielectrode data, we find spontaneous waves of activity in extrastriate visual cortex of awake, behaving marmosets (Callithrix jacchus). In monkeys trained to detect faint visual targets, the timing and position of spontaneous traveling waves, prior to target onset, predict the magnitude of target-evoked activity and the likelihood of target detection. In contrast, spatially disorganized fluctuations of neural activity are much less predictive. These results reveal an important role for spontaneous traveling waves in sensory processing through modulating neural and perceptual sensitivity.

EEG microstate periodicity explained by rotating phase patterns of resting-state alpha oscillations

Spatio-temporal patterns in electroencephalography (EEG) can be described by microstate analysis, a discrete approximation of the continuous electric field patterns produced by the cerebral cortex. Resting-state EEG microstates are largely determined by alpha frequencies (8-12 Hz) and we recently demonstrated that microstates occur periodically with twice the alpha frequency. To understand the origin of microstate periodicity, we analyzed the analytic amplitude and the analytic phase of resting-state alpha oscillations independently. In continuous EEG data we found rotating phase patterns organized around a small number of phase singularities which varied in number and location. The spatial rotation of phase patterns occurred with the underlying alpha frequency. Phase rotors coincided with periodic microstate motifs involving the four canonical microstate maps. The analytic amplitude showed no oscillatory behaviour and was almost static across time intervals of 1-2 alpha cycles, resulting in the global pattern of a standing wave. In n=23 healthy adults, time-lagged mutual information analysis of microstate sequences derived from amplitude and phase signals of awake eyes-closed EEG records showed that only the phase component contributed to the periodicity of microstate sequences. Phase sequences showed mutual information peaks at multiples of 50 ms and the group average had a main peak at 100 ms (10 Hz), whereas amplitude sequences had a slow and monotonous information decay. This result was confirmed by an independent approach combining temporal principal component analysis (tPCA) and autocorrelation analysis. We reproduced our observations in a generic model of EEG oscillations composed of coupled non-linear oscillators (Stuart-Landau model). Phase-amplitude dynamics similar to experimental EEG occurred when the oscillators underwent a supercritical Hopf bifurcation, a common feature of many computational models of the alpha rhythm. These findings explain our previous description of periodic microstate recurrence and its relation to the time scale of alpha oscillations. Moreover, our results corroborate the predictions of computational models and connect experimentally observed EEG patterns to properties of critical oscillator networks.

Local Sleep and Alzheimer's Disease Pathophysiology

Even prior to the onset of the prodromal stages of Alzheimer's disease (AD), a constellation of sleep disturbances are apparent. A series of epidemiological studies indicate that multiple forms of these sleep disturbances are associated with increased risk for developing mild cognitive impairment (MCI) and AD, even triggering disease onset at an earlier age. Through the combination of causal manipulation studies in humans and rodents, as well as targeted examination of sleep disturbance with respect to AD biomarkers, mechanisms linking sleep disturbance to AD are beginning to emerge. In this review, we explore recent evidence linking local deficits in brain oscillatory function during sleep with local AD pathological burden and circuit-level dysfunction and degeneration. In short, three deficits in the local expression of sleep oscillations have been identified in relation to AD pathophysiology: (1) frequency-specific frontal deficits in slow wave expression during non-rapid eye movement (NREM) sleep, (2) deficits in parietal sleep spindle expression, and (3) deficits in the quality of electroencephalographic (EEG) desynchrony characteristic of REM sleep. These deficits are noteworthy since they differ from that seen in normal aging, indicating the potential presence of an abnormal aging process. How each of these are associated with β-amyloid (Aβ) and tau pathology, as well as neurodegeneration of circuits sensitive to AD pathophysiology, are examined in the present review, with a focus on the role of dysfunction within fronto-hippocampal and subcortical sleep-wake circuits. It is hypothesized that each of these local sleep deficits arise from distinct network-specific dysfunctions driven by regionally-specific accumulation of AD pathologies, as well as their associated neurodegeneration. Overall, the evolution of these local sleep deficits offer unique windows into the circuit-specific progression of distinct AD pathophysiological processes prior to AD onset, as well as their impact on brain function. This includes the potential erosion of sleep-dependent memory mechanisms, which may contribute to memory decline in AD. This review closes with a discussion of the remaining critical knowledge gaps and implications of this work for future mechanistic studies and studies implementing sleep-based treatment interventions.

A methodological framework for inverse-modeling of propagating cortical activity using MEG/EEG

The prevailing view on the dynamics of large-scale electrical activity in the human cortex is that it constitutes a functional network of discrete and localized circuits. Within this view, a natural way to analyse magnetoencephalographic (MEG) and electroencephalographic (EEG) data is by adopting methods from network theory. Invasive recordings, however, demonstrate that cortical activity is spatially continuous, rather than discrete, and exhibits propagation behavior. Furthermore, human cortical activity is known to propagate under a variety of conditions such as non-REM sleep, general anesthesia, and coma. Although several MEG/EEG studies have investigated propagating cortical activity, not much is known about the conditions under which such activity can be successfully reconstructed from MEG/EEG sensor-data. This study provides a methodological framework for inverse-modeling of propagating cortical activity. Within this framework, cortical activity is represented in the spatial frequency domain, which is more natural than the dipole domain when dealing with spatially continuous activity. We define angular power spectra, which show how the power of cortical activity is distributed across spatial frequencies, angular gain/phase spectra, which characterize the spatial filtering properties of linear inverse operators, and angular resolution matrices, which summarize how linear inverse operators leak signal within and across spatial frequencies. We adopt the framework to provide insight into the performance of several linear inverse operators in reconstructing propagating cortical activity from MEG/EEG sensor-data. We also describe how prior spatial frequency information can be incorporated into the inverse-modeling to obtain better reconstructions.

Global enhancement of cortical excitability following coactivation of large neuronal populations

Correlated activation of cortical neurons often occurs in the brain and repetitive correlated neuronal firing could cause long-term modifications of synaptic efficacy and intrinsic excitability. We found that repetitive optogenetic activation of neuronal populations in the mouse cortex caused enhancement of optogenetically evoked firing of local coactivated neurons as well as distant cortical neurons in both ipsilateral and contralateral hemispheres. This global enhancement of evoked responses required coactivation of a sufficiently large population of neurons either within one cortical area or distributed in several areas. Enhancement of neuronal firing was saturable after repeated episodes of coactivation, diminished by inhibition of N-methyl-d-aspartic acid receptors, and accompanied by elevated excitatory postsynaptic potentials, all consistent with activity-induced synaptic potentiation. Chemogenetic inhibition of neuronal activity of the thalamus decreased the enhancement effect, suggesting thalamic involvement. Thus, correlated excitation of large neuronal populations leads to global enhancement of neuronal excitability.

Universal theory of brain waves: from linear loops to nonlinear synchronized spiking and collective brain rhythms

An inhomogeneous anisotropic physical model of the brain cortex is presented that predicts the emergence of non-evanescent (weakly damped) wave-like modes propagating in the thin cortex layers transverse to both the mean neural fiber direction and to the cortex spatial gradient. Although the amplitude of these modes stays below the typically observed axon spiking potential, the lifetime of these modes may significantly exceed the spiking potential inverse decay constant. Full brain numerical simulations based on parameters extracted from diffusion and structural MRI confirm the existence and extended duration of these wave modes. Contrary to the standard paradigm that the neural fibers determine the pathways for signal propagation in the brain, the signal propagation due to the cortex wave modes in highly folded areas will exhibit no apparent correlation with the fiber directions. The results are consistent with numerous recent experimental animal and human brain studies demonstrating the existence of electrostatic field activity in the form of traveling waves (including studies where neuronal connections were severed) and with wave loop induced peaks observed in EEG spectra. In addition, we demonstrate that the resonant and non-resonant terms of the nonlinear coupling between multiple modes produce both synchronous spiking-like high frequency wave activity as well as low frequency wave rhythms as a result of their unique dispersion properties. Numerical simulation of forced multiple mode dynamics shows that as forcing increases there is a transition from damped to oscillatory regime that subsequently decays away as over-excitation is reached. The resonant nonlinear coupling results in the emergence of low frequency rhythms with frequencies that are several orders of magnitude below the linear frequencies of modes taking part in the coupling. The localization and persistence of these cortical wave modes, and this new mechanism for understanding the nature of spiking behavior, have significant implications in particular for neuroimaging methods that detect electromagnetic physiological activity, such as EEG and MEG, and in general for the understanding of brain activity, including mechanisms of memory.

Signal Propagation via Open-Loop Intrathalamic Architectures: A Computational Model

Propagation of signals across the cerebral cortex is a core component of many cognitive processes and is generally thought to be mediated by direct intracortical connectivity. The thalamus, by contrast, is considered to be devoid of internal connections and organized as a collection of parallel inputs to the cortex. Here, we provide evidence that "open-loop" intrathalamic pathways involving the thalamic reticular nucleus (TRN) can support propagation of oscillatory activity across the cortex. Recent studies support the existence of open-loop thalamo-reticulo-thalamic (TC-TRN-TC) synaptic motifs in addition to traditional closed-loop architectures. We hypothesized that open-loop structural modules, when connected in series, might underlie thalamic and, therefore cortical, signal propagation. Using a supercomputing platform to simulate thousands of permutations of a thalamocortical network based on physiological data collected in mice, rats, ferrets, and cats and in which select synapses were allowed to vary both by class and individually, we evaluated the relative capacities of closed-loop and open-loop TC-TRN-TC synaptic configurations to support both propagation and oscillation. We observed that (1) signal propagation was best supported in networks possessing strong open-loop TC-TRN-TC connectivity; (2) intrareticular synapses were neither primary substrates of propagation nor oscillation; and (3) heterogeneous synaptic networks supported more robust propagation of oscillation than their homogeneous counterparts. These findings suggest that open-loop, heterogeneous intrathalamic architectures might complement direct intracortical connectivity to facilitate cortical signal propagation.

One-Shot Tagging During Wake and Cueing During Sleep With Spatiotemporal Patterns of Transcranial Electrical Stimulation Can Boost Long-Term Metamemory of Individual Episodes in Humans

Targeted memory reactivation (TMR) during slow-wave oscillations (SWOs) in sleep has been demonstrated with sensory cues to achieve about 5-12% improvement in post-nap memory performance on simple laboratory tasks. But prior work has not yet addressed the one-shot aspect of episodic memory acquisition, or dealt with the presence of interference from ambient environmental cues in real-world settings. Further, TMR with sensory cues may not be scalable to the multitude of experiences over one's lifetime. We designed a novel non-invasive non-sensory paradigm that tags one-shot experiences of minute-long naturalistic episodes in immersive virtual reality (VR) with unique spatiotemporal amplitude-modulated patterns (STAMPs) of transcranial electrical stimulation (tES). In particular, we demonstrated that these STAMPs can be re-applied as brief pulses during SWOs in sleep to achieve about 10-20% improvement in the metamemory of targeted episodes compared to the control episodes at 48 hours after initial viewing. We found that STAMPs can not only facilitate but also impair metamemory for the targeted episodes based on an interaction between pre-sleep metamemory and the number of STAMP applications during sleep. Overnight metamemory improvements were mediated by spectral power increases following the offset of STAMPs in the slow-spindle band (8-12 Hz) for left temporal areas in the scalp electroencephalography (EEG) during sleep. These results prescribe an optimal strategy to leverage STAMPs for boosting metamemory and suggest that real-world episodic memories can be modulated in a targeted manner even with coarser, non-invasive spatiotemporal stimulation.

Sleep Spindles: Mechanisms and Functions

Sleep spindles are burstlike signals in the electroencephalogram (EEG) of the sleeping mammalian brain and electrical surface correlates of neuronal oscillations in thalamus. As one of the most inheritable sleep EEG signatures, sleep spindles probably reflect the strength and malleability of thalamocortical circuits that underlie individual cognitive profiles. We review the characteristics, organization, regulation, and origins of sleep spindles and their implication in non-rapid-eye-movement sleep (NREMS) and its functions, focusing on human and rodent. Spatially, sleep spindle-related neuronal activity appears on scales ranging from small thalamic circuits to functional cortical areas, and generates a cortical state favoring intracortical plasticity while limiting cortical output. Temporally, sleep spindles are discrete events, part of a continuous power band, and elements grouped on an infraslow time scale over which NREMS alternates between continuity and fragility. We synthesize diverse and seemingly unlinked functions of sleep spindles for sleep architecture, sensory processing, synaptic plasticity, memory formation, and cognitive abilities into a unifying sleep spindle concept, according to which sleep spindles 1) generate neural conditions of large-scale functional connectivity and plasticity that outlast their appearance as discrete EEG events, 2) appear preferentially in thalamic circuits engaged in learning and attention-based experience during wakefulness, and 3) enable a selective reactivation and routing of wake-instated neuronal traces between brain areas such as hippocampus and cortex. Their fine spatiotemporal organization reflects NREMS as a physiological state coordinated over brain and body and may indicate, if not anticipate and ultimately differentiate, pathologies in sleep and neurodevelopmental, -degenerative, and -psychiatric conditions.

Interictal epileptiform discharges shape large-scale intercortical communication

Dynamic interactions between remote but functionally specialized brain regions enable complex information processing. This intercortical communication is disrupted in the neural networks of patients with focal epilepsy, and epileptic activity can exert widespread effects within the brain. Using large-scale human intracranial electroencephalography recordings, we show that interictal epileptiform discharges (IEDs) are significantly coupled with spindles in discrete, individualized brain regions outside of the epileptic network. We found that a substantial proportion of these localized spindles travel across the cortical surface. Brain regions that participate in this IED-driven oscillatory coupling express spindles that have a broader spatial extent and higher tendency to propagate than spindles occurring in uncoupled regions. These altered spatiotemporal oscillatory properties identify areas that are shaped by epileptic activity independent of IED or seizure detection. Our findings suggest that IED-spindle coupling may be an important mechanism of interictal global network dysfunction that could be targeted to prevent disruption of normal neural activity.

Atypicalities in sleep and semantic consolidation in autism

Sleep is known to support the neocortical consolidation of declarative memory, including the acquisition of new language. Autism spectrum disorder (ASD) is often characterized by both sleep and language learning difficulties, but few studies have explored a potential connection between the two. Here, 54 children with and without ASD (matched on age, nonverbal ability and vocabulary) were taught nine rare animal names (e.g., pipa). Memory was assessed via definitions, naming and speeded semantic decision tasks immediately after learning (pre‐sleep), the next day (post‐sleep, with a night of polysomnography between pre‐ and post‐sleep tests) and roughly 1 month later (follow‐up). Both groups showed comparable performance at pre‐test and similar levels of overnight change on all tasks; but at follow‐up children with ASD showed significantly greater forgetting of the unique features of the new animals (e.g., pipa is a flat frog). Children with ASD had significantly lower central non‐rapid eye movement (NREM) sigma power. Associations between spindle properties and overnight changes in speeded semantic decisions differed by group. For the TD group, spindle duration predicted overnight changes in responses to novel animals but not familiar animals, reinforcing a role for sleep in the stabilization of new semantic knowledge. For the ASD group, sigma power and spindle duration were associated with improvements in responses to novel and particularly familiar animals, perhaps reflecting more general sleep‐associated improvements in task performance. Plausibly, microstructural sleep atypicalities in children with ASD and differences in how information is prioritized for consolidation may lead to cumulative consolidation difficulties, compromising the quality of newly formed semantic representations in long‐term memory.

New perspectives on dimensionality and variability from large-scale cortical dynamics

The neocortex is a multi-scale network, with intricate local circuitry interwoven into a global mesh of long-range connections. Neural activity propagates within this network on a wide range of temporal and spatial scales. At the micro scale, neurophysiological recordings reveal coordinated dynamics in local neural populations, which support behaviorally relevant computations. At the macro scale, neuroimaging modalities measure global activity fluctuations organized into spatiotemporal patterns across the entire brain. Here we review recent advances linking the local and global scales of cortical dynamics and their relationship to behavior. We argue that diverse experimental observations on the dimensionality and variability of neural activity can be reconciled by considering how activity propagates in space and time on multiple spatial scales.

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

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

Engineering Adenylate Cyclase Activated by Near-Infrared Window Light for Mammalian Optogenetic Applications

Light in the near-infrared optical window (NIRW) penetrates deep through mammalian tissues, including the skull and brain tissue. Here we engineered an adenylate cyclase (AC) activated by NIRW light (NIRW-AC) and suitable for mammalian applications. To accomplish this goal, we constructed fusions of several bacteriophytochrome photosensory and bacterial AC modules using guidelines for designing chimeric homodimeric bacteriophytochromes. One engineered NIRW-AC, designated IlaM5, has significantly higher activity at 37 °C, is better expressed in mammalian cells, and can mediate cAMP-dependent photoactivation of gene expression in mammalian cells, in favorable contrast to the NIRW-ACs engineered earlier. The ilaM5 gene expressed from an AAV vector was delivered into the ventral basal thalamus region of the mouse brain, resulting in the light-controlled suppression of the cAMP-dependent wave pattern of the sleeping brain known as spindle oscillations. Reversible spindle oscillation suppression was observed in sleeping mice exposed to light from an external light source. This study confirms the robustness of principles of homodimeric bacteriophytochrome engineering, describes a NIRW-AC suitable for mammalian optogenetic applications, and demonstrates the feasibility of controlling brain activity via NIRW-ACs using transcranial irradiation.

Electric field dynamics in the brain during multi-electrode transcranial electric stimulation

Neural oscillations play a crucial role in communication between remote brain areas. Transcranial electric stimulation with alternating currents (TACS) can manipulate these brain oscillations in a non-invasive manner. Recently, TACS using multiple electrodes with phase shifted stimulation currents were developed to alter long-range connectivity. Typically, an increase in coordination between two areas is assumed when they experience an in-phase stimulation and a disorganization through an anti-phase stimulation. However, the underlying biophysics of multi-electrode TACS has not been studied in detail. Here, we leverage direct invasive recordings from two non-human primates during multi-electrode TACS to characterize electric field magnitude and phase as a function of the phase of stimulation currents. Further, we report a novel "traveling wave" stimulation where the location of the electric field maximum changes over the stimulation cycle. Our results provide a mechanistic understanding of the biophysics of multi-electrode TACS and enable future developments of novel stimulation protocols.

Bird-like propagating brain activity in anesthetized Nile crocodiles

Study Objectives

The changes in electroencephalogram (EEG) activity that characterize sleep and its sub-states-slow-wave sleep (SWS) and rapid eye movement (REM) sleep-are similar in mammals and birds. SWS is characterized by EEG slow waves resulting from the synchronous alternation of neuronal membrane potentials between hyperpolarized down-states with neuronal quiescence and depolarized up-states associated with action potentials. By contrast, studies of non-avian reptiles report the presence of high-voltage sharp waves (HShW) during sleep. How HShW relate to EEG phenomena occurring during mammalian and avian sleep is unclear. We investigated the spatiotemporal patterns of electrophysiological phenomena in Nile crocodiles (Crocodylus niloticus) anesthetized with isoflurane to determine whether they share similar spatiotemporal patterns to mammalian and avian slow waves.

Methods

Recordings of anesthetized crocodiles were made using 64-channel penetrating arrays with electrodes arranged in an 8 × 8 equally spaced grid. The arrays were placed in the dorsal ventricular ridge (DVR), a region implicated in the genesis of HShW. Various aspects of the spatiotemporal distribution of recorded signals were investigated.

Results

Recorded signals revealed the presence of HShW resembling those reported in earlier studies of naturally sleeping reptiles. HShW propagated in complex and variable patterns across the DVR.

Conclusions

We demonstrate that HShW within the DVR propagate in complex patterns similar to those observed for avian slow waves recorded from homologous brain regions. Consequently, sleep with HShW may represent an ancestral form of SWS, characterized by up-states occurring less often and for a shorter duration than in mammals and birds.

Synchronization dependent on spatial structures of a mesoscopic whole-brain network

Complex structural connectivity of the mammalian brain is believed to underlie the versatility of neural computations. Many previous studies have investigated properties of small subsystems or coarse connectivity among large brain regions that are often binarized and lack spatial information. Yet little is known about spatial embedding of the detailed whole-brain connectivity and its functional implications. We focus on closing this gap by analyzing how spatially-constrained neural connectivity shapes synchronization of the brain dynamics based on a system of coupled phase oscillators on a mammalian whole-brain network at the mesoscopic level. This was made possible by the recent development of the Allen Mouse Brain Connectivity Atlas constructed from viral tracing experiments together with a new mapping algorithm. We investigated whether the network can be compactly represented based on the spatial dependence of the network topology. We found that the connectivity has a significant spatial dependence, with spatially close brain regions strongly connected and distal regions weakly connected, following a power law. However, there are a number of residuals above the power-law fit, indicating connections between brain regions that are stronger than predicted by the power-law relationship. By measuring the sensitivity of the network order parameter, we show how these strong connections dispersed across multiple spatial scales of the network promote rapid transitions between partial synchronization and more global synchronization as the global coupling coefficient changes. We further demonstrate the significance of the locations of the residual connections, suggesting a possible link between the network complexity and the brain’s exceptional ability to swiftly switch computational states depending on stimulus and behavioral context.

In a previous study, a data-driven large-scale model of mouse brain connectivity was constructed. This mouse brain connectivity model is estimated by a simplified model which only takes in account anatomy and distance dependence of connection strength which is best fit by a power law. The distance dependence model captures the connection strengths of the mouse whole-brain network well. But can it capture the dynamics? In this study, we show that a small number of connections which are missed by the simple spatial model lead to significant differences in dynamics. The presence of a small number of strong connections over longer distances increases sensitivity of synchronization to perturbations. Unlike the data-driven network, the network without these long-range connections, as well as the network in which these long range connections are shuffled, lose global synchronization while maintaining localized synchrony, underlining the significance of the exact topology of these distal connections in the data-driven brain network.

Simulating human sleep spindle MEG and EEG from ion channel and circuit level dynamics

BACKGROUND

Although they form a unitary phenomenon, the relationship between extracranial M/EEG and transmembrane ion flows is understood only as a general principle rather than as a well-articulated and quantified causal chain.

METHOD

We present an integrated multiscale model, consisting of a neural simulation of thalamus and cortex during stage N2 sleep and a biophysical model projecting cortical current densities to M/EEG fields. Sleep spindles were generated through the interactions of local and distant network connections and intrinsic currents within thalamocortical circuits. 32,652 cortical neurons were mapped onto the cortical surface reconstructed from subjects' MRI, interconnected based on geodesic distances, and scaled-up to current dipole densities based on laminar recordings in humans. MRIs were used to generate a quasi-static electromagnetic model enabling simulated cortical activity to be projected to the M/EEG sensors.

RESULTS

The simulated M/EEG spindles were similar in amplitude and topography to empirical examples in the same subjects. Simulated spindles with more core-dominant activity were more MEG weighted.

COMPARISON WITH EXISTING METHODS

Previous models lacked either spindle-generating thalamic neural dynamics or whole head biophysical modeling; the framework presented here is the first to simultaneously capture these disparate scales.

CONCLUSIONS

This multiscale model provides a platform for the principled quantitative integration of existing information relevant to the generation of sleep spindles, and allows the implications of future findings to be explored. It provides a proof of principle for a methodological framework allowing large-scale integrative brain oscillations to be understood in terms of their underlying channels and synapses.

Metastable brain waves

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

Precise Slow Oscillation-Spindle Coupling Promotes Memory Consolidation in Younger and Older Adults

Memory consolidation during sleep relies on the precisely timed interaction of rhythmic neural events. Here, we investigate differences in slow oscillations (SO; 0.5-1 Hz), sleep spindles (SP), and their coupling across the adult human lifespan and ask whether observed alterations relate to the ability to retain associative memories across sleep. We demonstrate that older adults do not show the fine-tuned coupling of fast SPs (12.5-16 Hz) to the SO peak present in younger adults but, instead, are characterized most by a slow SP power increase (9-12.5 Hz) at the end of the SO up-state. This slow SP power increase, typical for older adults, coincides with worse memory consolidation in young age already, whereas the tight precision of SO-fast SP coupling promotes memory consolidation across younger and older adults. Crucially, brain integrity in source regions of SO and SP generation, including the medial prefrontal cortex, thalamus, hippocampus and entorhinal cortex, reinforces this beneficial SO-SP coupling in old age. Our results reveal that cognitive functioning is not only determined by maintaining structural brain integrity across the adult lifespan, but also by the preservation of precisely timed neural interactions during sleep that enable the consolidation of declarative memories.

The Hidden Spatial Dimension of Alpha: 10-Hz Perceptual Echoes Propagate as Periodic Traveling Waves in the Human Brain

EEG reverse-correlation techniques have revealed that visual information processing entails a ∼10-Hz (alpha) occipital response that reverberates sensory inputs up to 1 s. However, the spatial distribution of these perceptual echoes remains unknown: are they synchronized across the brain, or do they propagate like a traveling wave? Here, in two experiments with varying stimulus locations, we demonstrate the systematic phase propagation of perceptual echoes. A single stimulation in the upper visual field produced an "echo traveling wave" propagating from posterior to frontal sensors. The simultaneous presentation of two independent stimuli in separate visual hemifields produced two superimposed traveling waves propagating in opposite directions. Strikingly, in each sensor, the phase of the two echoes differed, with a phase advance for the contralateral stimulus. Thus, alpha traveling waves sweep across the human brain, encoding stimulus position in the phase domain, in line with the 70-year-old "cortical scanning" hypothesis (Pitts and McCulloch, 1947).

Individual spindle detection and analysis in high-density recordings across the night and in thalamic stroke

Sleep spindles are thalamocortical oscillations associated with several behavioural and clinical phenomena. In clinical populations, spindle activity has been shown to be reduced in schizophrenia, as well as after thalamic stroke. Automatic spindle detection algorithms present the only feasible way to systematically examine individual spindle characteristics. We took an established algorithm for spindle detection, and adapted it to high-density EEG sleep recordings. To illustrate the detection and analysis procedure, we examined how spindle characteristics changed across the night and introduced a linear mixed model approach applied to individual spindles in adults (n = 9). Next we examined spindle characteristics between a group of paramedian thalamic stroke patients (n = 9) and matched controls. We found a high spindle incidence rate and that, from early to late in the night, individual spindle power increased with the duration and globality of spindles; despite decreases in spindle incidence and peak-to-peak amplitude. In stroke patients, we found that only left-sided damage reduced individual spindle power. Furthermore, reduction was specific to posterior/fast spindles. Altogether, we demonstrate how state-of-the-art spindle detection techniques, applied to high-density recordings, and analysed using advanced statistical approaches can yield novel insights into how both normal and pathological circumstances affect sleep.

Detection and analysis of spatiotemporal patterns in brain activity

There is growing evidence that population-level brain activity is often organized into propagating waves that are structured in both space and time. Such spatiotemporal patterns have been linked to brain function and observed across multiple recording methodologies and scales. The ability to detect and analyze these patterns is thus essential for understanding the working mechanisms of neural circuits. Here we present a mathematical and computational framework for the identification and analysis of multiple classes of wave patterns in neural population-level recordings. By drawing a conceptual link between spatiotemporal patterns found in the brain and coherent structures such as vortices found in turbulent flows, we introduce velocity vector fields to characterize neural population activity. These vector fields are calculated for both phase and amplitude of oscillatory neural signals by adapting optical flow estimation methods from the field of computer vision. Based on these velocity vector fields, we then introduce order parameters and critical point analysis to detect and characterize a diverse range of propagating wave patterns, including planar waves, sources, sinks, spiral waves, and saddle patterns. We also introduce a novel vector field decomposition method that extracts the dominant spatiotemporal structures in a recording. This enables neural data to be represented by the activity of a small number of independent spatiotemporal modes, providing an alternative to existing dimensionality reduction techniques which separate space and time components. We demonstrate the capabilities of the framework and toolbox with simulated data, local field potentials from marmoset visual cortex and optical voltage recordings from whole mouse cortex, and we show that pattern dynamics are non-random and are modulated by the presence of visual stimuli. These methods are implemented in a MATLAB toolbox, which is freely available under an open-source licensing agreement.

Structured activity such as propagating wave patterns at the level of neural circuits can arise from highly variable firing activity of individual neurons. This property makes the brain, a quintessential example of a complex system, analogous to other complex physical systems such as turbulent fluids, in which structured patterns like vortices similarly emerge from molecules that behave irregularly. In this study, by uniquely adapting techniques for the identification of coherent structures in fluid turbulence, we develop new analytical and computational methods for the reliable detection of a diverse range of propagating wave patterns in large-scale neural recordings, for comprehensive analysis and visualization of these patterns, and for analysis of their dominant spatiotemporal modes. We demonstrate that these methods can be used to uncover the essential spatiotemporal properties of neural population activity recorded by different modalities, thus offering new insights into understanding the working mechanisms of neural systems.

High-Density Porous Graphene Arrays Enable Detection and Analysis of Propagating Cortical Waves and Spirals

Cortical propagating waves have recently attracted significant attention by the neuroscience community. These travelling waves have been suggested to coordinate different brain areas and play roles in assisting neural plasticity and learning. However, it is extremely challenging to record them with very fine spatial scales over large areas to investigate their effect on neural dynamics or network connectivity changes. In this work, we employ high-density porous graphene microelectrode arrays fabricated using laser pyrolysis on flexible substrates to study the functional network connectivity during cortical propagating waves. The low-impedance porous graphene arrays are used to record cortical potentials during theta oscillations and drug-induced seizures in vivo. Spatiotemporal analysis on the neural recordings reveal that theta oscillations and epileptiform activities have distinct characteristics in terms of both synchronization and resulting propagating wave patterns. To investigate the network connectivity during the propagating waves, we perform network analysis. The results show that the propagating waves are consistent with the functional connectivity changes in the neural circuits, suggesting that the underlying network states are reflected by the cortical potential propagation patterns.

Investigation of Propagating Cortical Waves and Spirals Recorded by High Density Porous Graphene Arrays

Propagating waves along the cortical surface have recently attracted significant attention by the neuroscience community. However, whether these propagating waves imply network connectivity changes for the neural circuits is not known. In this work, we employ a high density porous graphene microelectrode array and perform in vivo experiments with rodents to investigate network connectivity during cortical propagating waves. The spatial-temporal analysis of the cortical recordings reveals various types of propagating waves across the recording area. Network analysis results show that these propagating waves are consistent with the functional connectivity changes in the neural circuits, suggesting that the underlying network states are reflected by the cortical potential propagation patterns.