Origin, structure, and role of background EEG activity. Part 2. Analytic phase

Information-Theoretical Analysis of the Cycle of Creation of Knowledge and Meaning in Brains under Multiple Cognitive Modalities

It is of great interest to develop advanced sensory technologies allowing non-invasive monitoring of neural correlates of cognitive processing in people performing everyday tasks. A lot of progress has been reported in recent years in this research area using scalp EEG arrays, but the high level of noise in the electrode signals poses a lot of challenges. This study presents results of detailed statistical analysis of experimental data on the cycle of creation of knowledge and meaning in human brains under multiple cognitive modalities. We measure brain dynamics using a HydroCel Geodesic Sensor Net, 128-electrode dense-array electroencephalography (EEG). We compute a pragmatic information (PI) index derived from analytic amplitude and phase, by Hilbert transforming the EEG signals of 20 participants in six modalities, which combine various audiovisual stimuli, leading to different mental states, including relaxed and cognitively engaged conditions. We derive several relevant measures to classify different brain states based on the PI indices. We demonstrate significant differences between engaged brain states that require sensory information processing to create meaning and knowledge for intentional action, and relaxed-meditative brain states with less demand on psychophysiological resources. We also point out that different kinds of meanings may lead to different brain dynamics and behavioral responses.

The Electrome of a Parasitic Plant in a Putative State of Attention Increases the Energy of Low Band Frequency Waves: A Comparative Study with Neural Systems

Selective attention is an important cognitive phenomenon that allows organisms to flexibly engage with certain environmental cues or activities while ignoring others, permitting optimal behaviour. It has been proposed that selective attention can be present in many different animal species and, more recently, in plants. The phenomenon of attention in plants would be reflected in its electrophysiological activity, possibly being observable through electrophytographic (EPG) techniques. Former EPG time series obtained from the parasitic plant Cuscuta racemosa in a putative state of attention towards two different potential hosts, the suitable bean (Phaseolus vulgaris) and the unsuitable wheat (Triticum aestivum), were revisited. Here, we investigated the potential existence of different band frequencies (including low, delta, theta, mu, alpha, beta, and gamma waves) using a protocol adapted from neuroscientific research. Average band power (ABP) was used to analyse the energy distribution of each band frequency in the EPG signals, and time dispersion analysis of features (TDAF) was used to explore the variations in the energy of each band. Our findings indicated that most band waves were centred in the lower frequencies. We also observed that C. racemosa invested more energy in these low-frequency waves when suitable hosts were present. However, we also noted peaks of energy investment in all the band frequencies, which may be linked to extremely low oscillatory electrical signals in the entire tissue. Overall, the presence of suitable hosts induced a higher energy power, which supports the hypothesis of attention in plants. We further discuss and compare our results with generic neural systems.

Nonlinear dynamics of EEG responses to unmanned vehicle visual detection with different levels of task difficulty

The main objective of this study was to examine the presence of chaos in the EEG recordings of brain activity under simulated unmanned ground vehicle visual detection scenarios with different levels of task difficulty. One hundred and fifty people participated in the experiment and completed four visual detection task scenarios: (1) change detection, (2) a threat detection task, (3) a dual-task with different change detection task rates, and (4) a dual-task with different threat detection task rates. We used the largest Lyapunov exponent and correlation dimension of the EEG data and performed 0-1 tests on the EEG data. The results revealed a change in the level of nonlinearity in the EEG data corresponding to different levels of cognitive task difficulty. The differences in EEG nonlinearity measures among the studied levels of task difficulty, as well as between a single task scenario and a dual-task scenario, have also been assessed. The results increase our understanding of the nature of unmanned systems' operational requirements.

A Pilot Investigation of Visual Pathways in Patients with Mild Traumatic Brain Injury

In this study, we examined visual processing within primary visual areas (V1) in normal and visually impaired individuals who exhibit significant visual symptomology due to sports-related mild traumatic brain injury (mTBI). Five spatial frequency stimuli were applied to the right, left and both eyes in order to assess the visual processing of patients with sports-related mild traumatic brain injuries who exhibited visual abnormalities, i.e., photophobia, blurriness, etc., and controls. The measurement of the left/right eye and binocular integration was accomplished via the quantification of the spectral power and visual event-related potentials. The principal results have shown that the power spectral density (PSD) measurements display a distinct loss in the alpha band-width range, which corresponded to more instances of medium-sized receptive field loss. Medium-size receptive field loss may correspond to parvocellular (p-cell) processing deprecation. Our major conclusion provides a new measurement, using PSD analysis to assess mTBI conditions from primary V1 areas. The statistical analysis demonstrated significant differences between the mTBI and control cohort in the Visual Evoked Potentials (VEP) amplitude responses and PSD measurements. Additionally, the PSD measurements were able to assess the improvement in the mTBI primary visual areas over time through rehabilitation.

Classification of Motor Imagery EEG Signals Based on Data Augmentation and Convolutional Neural Networks

In brain-computer interface (BCI) systems, motor imagery electroencephalography (MI-EEG) signals are commonly used to detect participant intent. Many factors, including low signal-to-noise ratios and few high-quality samples, make MI classification difficult. In order for BCI systems to function, MI-EEG signals must be studied. In pattern recognition and other fields, deep learning approaches have recently been successfully applied. In contrast, few effective deep learning algorithms have been applied to BCI systems, especially MI-based systems. In this paper, we address these problems from two aspects based on the characteristics of EEG signals: first, we proposed a combined time-frequency domain data enhancement method. This method guarantees that the size of the training data is effectively increased while maintaining the intrinsic composition of the data. Second, our design consists of a parallel CNN that takes both raw EEG images and images transformed through continuous wavelet transform (CWT) as inputs. We conducted classification experiments on a public data set to verify the effectiveness of the algorithm. According to experimental results based on the BCI Competition IV Dataset2a, the average classification accuracy is 97.61%. A comparison of the proposed algorithm with other algorithms shows that it performs better in classification. The algorithm can be used to improve the classification performance of MI-based BCIs and BCI systems created for people with disabilities.

Pre-Stimulus Power but Not Phase Predicts Prefrontal Cortical Excitability in TMS-EEG

The cortical response to transcranial magnetic stimulation (TMS) has notable inter-trial variability. One source of this variability can be the influence of the phase and power of pre-stimulus neuronal oscillations on single-trial TMS responses. Here, we investigate the effect of brain oscillatory activity on TMS response in 49 distinct healthy participants (64 datasets) who had received single-pulse TMS over the left dorsolateral prefrontal cortex. Across all frequency bands of theta (4-7 Hz), alpha (8-13 Hz), and beta (14-30 Hz), there was no significant effect of pre-TMS phase on single-trial cortical evoked activity. After high-powered oscillations, whether followed by a TMS pulse or not, the subsequent activity was larger than after low-powered oscillations. We further defined a measure, corrected_effect, to enable us to investigate brain responses to the TMS pulse disentangled from the power of ongoing (spontaneous) oscillations. The corrected_effect was significantly different from zero (meaningful added effect of TMS) only in theta and beta bands. Our results suggest that brain state prior to stimulation might play some role in shaping the subsequent TMS-EEG response. Specifically, our findings indicate that the power of ongoing oscillatory activity, but not phase, can influence brain responses to TMS. Aligning the TMS pulse with specific power thresholds of an EEG signal might therefore reduce variability in neurophysiological measurements and also has the potential to facilitate more robust therapeutic effects of stimulation.

Dual EEG alignment between participants during shared intentionality experiments

Electroencephalograph (EEG) analysis from human subjects have demonstrated that beta oscillations carried perceptual information across the cortex featuring amplitude and phase modulation occurrences when subjects are engaged in task-oriented activities. A hypothesis was tested that synchronized patterns could be found in the scalp EEG of two human subjects engaged in similar intentional activity. Signals were recorded from scalp electrodes and band-pass filtered. The Hilbert transform decomposes the EEG signals into the analytic phase and amplitude. With these components of the EEG signal, a systematic search of the alpha, beta, delta, gamma, and theta spectrum is executed to locate temporal patterns. The amplitude and phase modulation were classified with respect to task intervals. Temporal patterns were found in the alpha-beta range (15-30 Hz). Our results suggest that the scalp EEG can yield information about the timing of episodically synchronized brain activity in higher cognitive function between two individuals engaged in similar task-oriented activities.

Interpretable functional specialization emerges in deep convolutional networks trained on brain signals

OBJECTIVE

Functional specialization is fundamental to neural information processing. Here, we study whether and how functional specialization emerges in artificial deep convolutional neural networks (CNNs) during a brain-computer interfacing (BCI) task.

APPROACH

We trained CNNs to predict hand movement speed from intracranial EEG (iEEG) and delineated how units across the different CNN hidden layers learned to represent the iEEG signal.

MAIN RESULTS

We show that distinct, functionally interpretable neural populations emerged as a result of the training process. While some units became sensitive to either iEEG amplitude or phase, others showed bimodal behavior with significant sensitivity to both features. Pruning of highly-sensitive units resulted in a steep drop of decoding accuracy not observed for pruning of less sensitive units, highlighting the functional relevance of the amplitude- and phase-specialized populations.

SIGNIFICANCE

We anticipate that emergent functional specialization as uncovered here will become a key concept in research towards interpretable deep learning for neuroscience and BCI applications.

Long-Range Respiratory and Theta Oscillation Networks Depend on Spatial Sensory Context

Neural oscillations can couple networks of brain regions, especially at lower frequencies. The nasal respiratory rhythm, which elicits robust olfactory bulb oscillations, has been linked to episodic memory, locomotion, and exploration, along with widespread oscillatory coherence. The piriform cortex is implicated in propagating the olfactory-bulb-driven respiratory rhythm, but this has not been tested explicitly in the context of both hippocampal theta and nasal respiratory rhythm during exploratory behaviors. We investigated systemwide interactions during foraging behavior, which engages respiratory and theta rhythms. Local field potentials from the olfactory bulb, piriform cortex, dentate gyrus, and CA1 of hippocampus, primary visual cortex, and nasal respiration were recorded simultaneously from male rats. We compared interactions among these areas while rats foraged using either visual or olfactory spatial cues. We found high coherence during foraging compared with home cage activity in two frequency bands that matched slow and fast respiratory rates. Piriform cortex and hippocampus maintained strong coupling at theta frequency during periods of slow respiration, whereas other pairs showed coupling only at the fast respiratory frequency. Directional analysis shows that the modality of spatial cues was matched to larger influences in the network by the respective primary sensory area. Respiratory and theta rhythms also coupled to faster oscillations in primary sensory and hippocampal areas. These data provide the first evidence of widespread interactions among nasal respiration, olfactory bulb, piriform cortex, and hippocampus in awake freely moving rats, and support the piriform cortex as an integrator of respiratory and theta activity.SIGNIFICANCE STATEMENT Recent studies have shown widespread interactions between the nasally driven respiratory rhythm and neural oscillations in hippocampus and neocortex. With this study, we address how the respiratory rhythm interacts with ongoing slow brain rhythms across olfactory, hippocampal, and visual systems in freely moving rats. Patterns of network connectivity change with behavioral state, with stronger interactions at fast and slow respiratory frequencies during foraging as compared with home cage activity. Routing of interactions between sensory cortices depends on the modality of spatial cues present during foraging. Functional connectivity and cross-frequency coupling analyses suggest strong bidirectional interactions between olfactory and hippocampal systems related to respiration and point to the piriform cortex as a key area for mediating respiratory and theta rhythms.

Classification of Phonological Categories in Imagined Speech using Phase Synchronization Measure

Phonological categories in articulated speech are defined based on the place and manner of articulation. In this work, we investigate whether the phonological categories of the prompts imagined during speech imagery lead to differences in phase synchronization in various cortical regions that can be discriminated from the EEG captured during the imagination. Nasal and bilabial consonant are the two phonological categories considered due to their differences in both place and manner of articulation. Mean phase coherence (MPC) is used for measuring the phase synchronization and shallow neural network (NN) is used as the classifier. As a benchmark, we have also designed another NN based on statistical parameters extracted from imagined speech EEG. The NN trained on MPC values in the beta band gives classification results superior to NN trained on alpha band MPC values, gamma band MPC values and statistical parameters extracted from the EEG.Clinical relevance: Brain-computer interface (BCI) is a promising tool for aiding differently-abled people and for neurorehabilitation. One of the challenges in designing speech imagery based BCI is the identification of speech prompts that can lead to distinct neural activations. We have shown that nasal and blilabial consonants lead to dissimilar activations. Hence prompts orthogonal in these phonological categories are good choices as speech imagery prompts.

Can we identify the category of imagined phoneme from EEG?

Phonemes are classified into different categories based on the place and manner of articulation. We investigate the differences between the neural correlates of imagined nasal and bilabial consonants (distinct phonological categories). Mean phase coherence is used as a metric for measuring the phase synchronisation between pairs of electrodes in six cortical regions (auditory, motor, prefrontal, sensorimotor, so-matosensory and premotor) during the imagery of nasal and bilabial consonants. Statistically significant difference at 95% confidence interval is observed in beta and lower-gamma bands in various cortical regions. Our observations are inline with the directions into velocities of articulators and dual stream prediction models and support the hypothesis that phonological categories not only exist in articulated speech but can also be distinguished from the EEG of imagined speech.

Prefrontal Theta-Phase Synchronized Brain Stimulation With Real-Time EEG-Triggered TMS

Background

Theta-band neuronal oscillations in the prefrontal cortex are associated with several cognitive functions. Oscillatory phase is an important correlate of excitability and phase synchrony mediates information transfer between neuronal populations oscillating at that frequency. The ability to extract and exploit the prefrontal theta rhythm in real time in humans would facilitate insight into neurophysiological mechanisms of cognitive processes involving the prefrontal cortex, and development of brain-state-dependent stimulation for therapeutic applications.

Objectives

We investigate individual source-space beamforming-based estimation of the prefrontal theta oscillation as a method to target specific phases of the ongoing theta oscillations in the human dorsomedial prefrontal cortex (DMPFC) with real-time EEG-triggered transcranial magnetic stimulation (TMS). Different spatial filters for extracting the prefrontal theta oscillation from EEG signals are compared and additional signal quality criteria are assessed to take into account the dynamics of this cortical oscillation.

Methods

Twenty two healthy participants were recruited for anatomical MRI scans and EEG recordings with 18 composing the final analysis. We calculated individual spatial filters based on EEG beamforming in source space. The extracted EEG signal was then used to simulate real-time phase-detection and quantify the accuracy as compared to post-hoc phase estimates. Different spatial filters and triggering parameters were compared. Finally, we validated the feasibility of this approach by actual real-time triggering of TMS pulses at different phases of the prefrontal theta oscillation.

Results

Higher phase-detection accuracy was achieved using individualized source-based spatial filters, as compared to an average or standard Laplacian filter, and also by detecting and avoiding periods of low theta amplitude and periods containing a phase reset. Using optimized parameters, prefrontal theta-phase synchronized TMS of DMPFC was achieved with an accuracy of ±55°.

Conclusion

This study demonstrates the feasibility of triggering TMS pulses during different phases of the ongoing prefrontal theta oscillation in real time. This method is relevant for brain state-dependent stimulation in human studies of cognition. It will also enable new personalized therapeutic repetitive TMS protocols for more effective treatment of neuropsychiatric disorders.

Scale-Free Coupled Dynamics in Brain Networks Captured by Bivariate Focus-Based Multifractal Analysis

While most connectivity studies investigate functional connectivity (FC) in a scale-dependent manner, coupled neural processes may also exhibit broadband dynamics, manifesting as power-law scaling of their measures of interdependence. Here we introduce the bivariate focus-based multifractal (BFMF) analysis as a robust tool for capturing such scale-free relations and use resting-state electroencephalography (EEG) recordings of 12 subjects to demonstrate its performance in reconstructing physiological networks. BFMF was employed to characterize broadband FC between 62 cortical regions in a pairwise manner, with all investigated connections being tested for true bivariate multifractality. EEG channels were also grouped to represent the activity of six resting-state networks (RSNs) in the brain, thus allowing for the analysis of within- and between- RSNs connectivity, separately. Most connections featured true bivariate multifractality, which could be attributed to the genuine scale-free coupling of neural dynamics. Bivariate multifractality showed a characteristic topology over the cortex that was highly concordant among subjects. Long-term autocorrelation was higher in within-RSNs, while the degree of multifractality was generally found stronger in between-RSNs connections. These results offer statistical evidence of the bivariate multifractal nature of functional coupling in the brain and validate BFMF as a robust method to capture such scale-independent coupled dynamics.

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.

The shaky ground truth of real-time phase estimation

Instantaneous phase of brain oscillations in electroencephalography (EEG) is a measure of brain state that is relevant to neuronal processing and modulates evoked responses. However, determining phase at the time of a stimulus with standard signal processing methods is not possible due to the stimulus artifact masking the future part of the signal. Here, we quantify the degree to which signal-to-noise ratio and instantaneous amplitude of the signal affect the variance of phase estimation error and the precision with which "ground truth" phase is even defined, using both the variance of equivalent estimators and realistic simulated EEG data with known synthetic phase. Necessary experimental conditions are specified in which pre-stimulus phase estimation is meaningfully possible based on instantaneous amplitude and signal-to-noise ratio of the oscillation of interest. An open source toolbox is made available for causal (using pre-stimulus signal only) phase estimation along with a EEG dataset consisting of recordings from 140 participants and a best practices workflow for algorithm optimization and benchmarking. As an illustration, post-hoc sorting of open-loop transcranial magnetic stimulation (TMS) trials according to pre-stimulus sensorimotor μ-rhythm phase is performed to demonstrate modulation of corticospinal excitability, as indexed by the amplitude of motor evoked potentials.

The rhythm of memory: how breathing shapes memory function

The mammalian olfactory bulb displays a prominent respiratory rhythm, which is linked to the sniff cycle and is driven by sensory input from olfactory receptors in the nasal sensory epithelium. In rats and mice, respiratory frequencies occupy the same band as the hippocampal θ-rhythm, which has been shown to be a key player in memory processes. Hippocampal and olfactory bulb rhythms were previously found to be uncorrelated except in specific odor-contingency learning circumstances. However, many recent electrophysiological studies in both rodents and humans reveal a surprising cycle-by-cycle influence of nasal respiration on neuronal activity throughout much of the cerebral cortex beyond the olfactory system, including the prefrontal cortex, hippocampus, and subcortical structures. In addition, respiratory phase has been shown to influence higher-frequency oscillations associated with cognitive functions, including attention and memory, such as the power of γ-rhythms and the timing of hippocampal sharp wave ripples. These new findings support respiration's role in cognitive function, which is supported by studies in human subjects, in which nasal respiration has been linked to memory processes. Here, we review recent reports from human and rodent experiments that link respiration to the modulation of memory function and the neurophysiological processes involved in memory in rodents and humans. We argue that respiratory influence on the neuronal activity of two key memory structures, the hippocampus and prefrontal cortex, provides a potential neuronal mechanism behind respiratory modulation of memory.

Cortico-Striatal Cross-Frequency Coupling and Gamma Genesis Disruptions in Huntington's Disease Mouse and Computational Models

Abnormal gamma band power across cortex and striatum is an important phenotype of Huntington's disease (HD) in both patients and animal models, but neither the origin nor the functional relevance of this phenotype is well understood. Here, we analyzed local field potential (LFP) activity in freely behaving, symptomatic R6/2 and Q175 mouse models and corresponding wild-type (WT) controls. We focused on periods of quiet rest, which show strong γ activity in HD mice. Simultaneous recording from motor cortex and its target area in dorsal striatum in the R6/2 model revealed exaggerated functional coupling over that observed in WT between the phase of delta frequencies (1-4 Hz) in cortex and striatum and striatal amplitude modulation of low γ frequencies (25-55 Hz; i.e., phase-amplitude coupling, PAC), but no evidence that abnormal cortical activity alone can account for the increase in striatal γ power. Both HD mouse models had stronger coupling of γ amplitude to δ phase and more unimodal phase distributions than their WT counterparts. To assess the possible role of striatal fast-spiking interneurons (FSIs) in these phenomena, we developed a computational model based on additional striatal recordings from Q175 mice. Changes in peak γ frequency and power ratio were readily reproduced by our computational model, accounting for several experimental findings reported in the literature. Our results suggest that HD is characterized by both a reorganization of cortico-striatal drive and specific population changes related to intrastriatal synaptic coupling.

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.

1/f electrophysiological spectra in resting and drug-induced states can be explained by the dynamics of multiple oscillatory relaxation processes

Neurophysiological recordings are dominated by arhythmical activity whose spectra can be characterised by power-law functions, and on this basis are often referred to as reflecting scale-free brain dynamics (1/f^β). Relatively little is known regarding the neural generators and temporal dynamics of this arhythmical behaviour compared to rhythmical behaviour. Here we used Irregularly Resampled AutoSpectral Analysis (IRASA) to quantify β, in both the high (5-100 Hz, β_hf) and low frequency bands (0.1-2.5 Hz, β_lf) in MEG/EEG/ECoG recordings and to separate arhythmical from rhythmical modes of activity, such as, alpha rhythms. In MEG/EEG/ECoG data, we demonstrate that oscillatory alpha power dynamically correlates over time with β_hf and similarly, participants with higher rhythmical alpha power have higher β_hf. In a series of pharmaco-MEG investigations using the GABA reuptake inhibitor tiagabine, the glutamatergic AMPA receptor antagonist perampanel, the NMDA receptor antagonist ketamine and the mixed partial serotonergic agonist LSD, a variety of effects on both β_hf and β_lf were observed. Additionally, strong modulations of β_hf were seen in monkey ECoG data during general anaesthesia using propofol and ketamine. We develop and test a unifying model which can explain, the 1/f nature of electrophysiological spectra, their dynamic interaction with oscillatory rhythms as well as the sensitivity of 1/f activity to drug interventions by considering electrophysiological spectra as being generated by a collection of stochastically perturbed damped oscillators having a distribution of relaxation rates.

Mesoscopic neuron population modeling of normal/epileptic brain dynamics

Simulations of EEG data provide the understanding of how the limbic system exhibits normal and abnormal states of the electrical activity of the brain. While brain activity exhibits a type of homeostasis of excitatory and inhibitory mesoscopic neuron behavior, abnormal neural firings found in the seizure state exhibits brain instability due to runaway oscillatory entrained neural behavior. We utilize a model of mesoscopic brain activity, the KIV model, where each network represents the areas of the limbic system, i.e., hippocampus, sensory cortex, and the amygdala. Our model initially demonstrates oscillatory entrained neural behavior as the epileptogenesis, and then by increasing the external weights that join the three networks that represent the areas of the limbic system, seizure activity entrains the entire system. By introducing an external signal into the model, simulating external electrical titration therapy, the modeled seizure behavior can be 'rebalanced' back to its normal state.

Power Spectral Density Analysis of Electrocorticogram Recordings during Cerebral Hypothermia in Neonatal Seizures

BACKGROUND

Neonatal seizures (NS) are the most common form of neurological dysfunction observed in newborns.

PURPOSE

The purpose of this study in newborn piglets was to determine the effect of cerebral hypothermia (CH) on neural activity during pharmacologically induced NS. We hypothesized that the neuroprotective effects of CH would preserve higher frequencies observed in electrocorticogram (ECoG) recordings.

METHODS

Power spectral density was employed to determine the levels of brain activity in ECoGs to quantitatively assess the power of each frequency observed in neurological brain states of delta, theta, alpha, and beta-gamma frequencies.

RESULT

The most significant reduction of power occurs in the lower frequency band of delta-theta-alpha of CH cohorts, while t score probabilities imply that high-frequency brain activity in the beta-gamma range is preserved in the CH population.

CONCLUSION

While the overall power density decreases over time in both groups, the decrease is to a lesser degree in the CH population.

Cinematic Operation of the Cerebral Cortex Interpreted via Critical Transitions in Self-Organized Dynamic Systems

Measurements of local field potentials over the cortical surface and the scalp of animals and human subjects reveal intermittent bursts of beta and gamma oscillations. During the bursts, narrow-band metastable amplitude modulation (AM) patters emerge for a fraction of a second and ultimately dissolve to the broad-band random background activity. The burst process depends on previously learnt conditioned stimuli (CS), thus different AM patterns may emerge in response to different CS. This observation leads to our cinematic theory of cognition when perception happens in discrete steps manifested in the sequence of AM patterns. Our article summarizes findings in the past decades on experimental evidence of cinematic theory of cognition and relevant mathematical models. We treat cortices as dissipative systems that self-organize themselves near a critical level of activity that is a non-equilibrium metastable state. Criticality is arguably a key aspect of brains in their rapid adaptation, reconfiguration, high storage capacity, and sensitive response to external stimuli. Self-organized criticality (SOC) became an important concept to describe neural systems. We argue that transitions from one AM pattern to the other require the concept of phase transitions, extending beyond the dynamics described by SOC. We employ random graph theory (RGT) and percolation dynamics as fundamental mathematical approaches to model fluctuations in the cortical tissue. Our results indicate that perceptions are formed through a phase transition from a disorganized (high entropy) to a well-organized (low entropy) state, which explains the swiftness of the emergence of the perceptual experience in response to learned stimuli.

Breathing as a Fundamental Rhythm of Brain Function

Ongoing fluctuations of neuronal activity have long been considered intrinsic noise that introduces unavoidable and unwanted variability into neuronal processing, which the brain eliminates by averaging across population activity (Georgopoulos et al., 1986; Lee et al., 1988; Shadlen and Newsome, 1994; Maynard et al., 1999). It is now understood, that the seemingly random fluctuations of cortical activity form highly structured patterns, including oscillations at various frequencies, that modulate evoked neuronal responses (Arieli et al., 1996; Poulet and Petersen, 2008; He, 2013) and affect sensory perception (Linkenkaer-Hansen et al., 2004; Boly et al., 2007; Sadaghiani et al., 2009; Vinnik et al., 2012; Palva et al., 2013). Ongoing cortical activity is driven by proprioceptive and interoceptive inputs. In addition, it is partially intrinsically generated in which case it may be related to mental processes (Fox and Raichle, 2007; Deco et al., 2011). Here we argue that respiration, via multiple sensory pathways, contributes a rhythmic component to the ongoing cortical activity. We suggest that this rhythmic activity modulates the temporal organization of cortical neurodynamics, thereby linking higher cortical functions to the process of breathing.

Reflections on a giant of brain science: How lucky we are having Walter J. Freeman as our beacon in cognitive neurodynamics research

Walter J. Freeman was a giant of the field of neuroscience whose visionary work contributed various experimental and theoretical breakthroughs to brain research in the past 60 years. He has pioneered a number of Electroencephalogram and Electrocorticogram tools and approaches that shaped the field, while "Freeman Neurodynamics" is a theoretical concept that is widely known, used, and respected among neuroscientists all over the world. His recent death is a profound loss to neuroscience and biomedical engineering. Many of his revolutionary ideas on brain dynamics have been ahead of their time by decades. We summarize his following groundbreaking achievements: (1) Mass Action in the Nervous System, from microscopic (single cell) recordings, through mesoscopic populations, to large-scale collective brain patterns underlying cognition; (2) Freeman-Kachalsky model of multi-scale, modular brain dynamics; (3) cinematic theory of cognitive dynamics; (4) phase transitions in cortical dynamics modeled with random graphs and quantum field theory; (5) philosophical aspects of intentionality, consciousness, and the unity of brain-mind-body. His work has been admired by many of his neuroscientist colleagues and followers. At the same time, his multidisciplinary approach combining advanced concepts of control theory and the mathematics of nonlinear systems and chaos, poses significant challenges to those who wish to thoroughly understand his message. The goal of this commemorative paper is to review key aspects of Freeman's neurodynamics and to provide some handles to gain better understanding about Freeman's extraordinary intellectual achievement.

Trajectory of frequency stability in typical development

This work explores a feature of brain dynamics, metastability, by which transients are observed in functional brain data. Metastability is a balance between static (stable) and dynamic (unstable) tendencies in electrophysiological brain activity. Furthermore, metastability is a theoretical mechanism underlying the rapid synchronization of cell assemblies that serve as neural substrates for cognitive states, and it has been associated with cognitive flexibility. While much previous research has sought to characterize metastability in the adult human brain, few studies have examined metastability in early development, in part because of the challenges of acquiring adequate, noise free continuous data in young children. To accomplish this endeavor, we studied a new method for characterizing the stability of EEG frequency in early childhood, as inspired by prior approaches for describing cortical phase resets in the scalp EEG of healthy adults. Specifically, we quantified the variance of the rate of change of the signal phase (i.e., frequency) as a proxy for phase resets (signal instability), given that phase resets occur almost simultaneously across large portions of the scalp. We tested our method in a cohort of 39 preschool age children (age =53 ± 13.6 months). We found that our outcome variable of interest, frequency variance, was a promising marker of signal stability, as it increased with the number of phase resets in surrogate (artificial) signals. In our cohort of children, frequency variance decreased cross-sectionally with age (r = -0.47, p = 0.0028). EEG signal stability, as quantified by frequency variance, increases with age in preschool age children. Future studies will relate this biomarker with the development of executive function and cognitive flexibility in children, with the overarching goal of understanding metastability in atypical development.

Transcranial focused ultrasound modulates intrinsic and evoked EEG dynamics

BACKGROUND

The integration of EEG recordings and transcranial neuromodulation has provided a useful construct for noninvasively investigating the modification of human brain circuit activity. Recent evidence has demonstrated that focused ultrasound can be targeted through the human skull to affect the amplitude of somatosensory evoked potentials and its associated spectral content.

OBJECTIVE/HYPOTHESIS

The present study tests whether focused ultrasound transmitted through the human skull and targeted to somatosensory cortex can affect the phase and phase rate of cortical oscillatory dynamics.

METHODS

A computational model was developed to gain insight regarding the insertion behavior of ultrasound induced pressure waves in the human head. The instantaneous phase and phase rate of EEG recordings before, during, and after transmission of transcranial focused ultrasound (tFUS) to human somatosensory cortex were examined to explore its effects on phase dynamics.

RESULTS

Computational modeling results show the skull effectively reinforces the focusing of tFUS due to curvature of material interfaces. Neurophysiological recordings show that tFUS alters the phase distribution of intrinsic brain activity for beta frequencies, but not gamma. This modulation was accompanied by a change in phase rate of both beta and gamma frequencies. Additionally, tFUS modulated phase distributions in the beta band of early sensory-evoked activity but did not affect late sensory-evoked activity, lending support to the spatial specificity of tFUS for neuromodulation. This spatial specificity was confirmed through an additional experiment where the ultrasound transducer was moved 1 cm laterally from the original cortical target.

CONCLUSIONS

Focused ultrasonic energy can alter EEG oscillatory dynamics through local mechanical perturbation of discrete cortical circuits.

The use of many-body physics and thermodynamics to describe the dynamics of rhythmic generators in sensory cortices engaged in memory and learning

The problem of the transition from the molecular and cellular level to the macroscopic level of observed assemblies of myriads of neurons is the subject addressed in this report. The great amount of detailed information available at molecular and cellular level seems not sufficient to account for the high effectiveness and reliability observed in the brain macroscopic functioning. It is suggested that the dissipative many-body model and thermodynamics might offer the dynamical frame underlying the rich phenomenology observed at microscopic and macroscopic level and help in the understanding on how to fill the gap between the bio-molecular and cellular level and the one of brain macroscopic functioning.

Review of the BCI Competition IV

The BCI competition IV stands in the tradition of prior BCI competitions that aim to provide high quality neuroscientific data for open access to the scientific community. As experienced already in prior competitions not only scientists from the narrow field of BCI compete, but scholars with a broad variety of backgrounds and nationalities. They include high specialists as well as students. The goals of all BCI competitions have always been to challenge with respect to novel paradigms and complex data. We report on the following challenges: (1) asynchronous data, (2) synthetic, (3) multi-class continuous data, (4) session-to-session transfer, (5) directionally modulated MEG, (6) finger movements recorded by ECoG. As after past competitions, our hope is that winning entries may enhance the analysis methods of future BCIs.

The origin of extracellular fields and currents--EEG, ECoG, LFP and spikes

Neuronal activity in the brain gives rise to transmembrane currents that can be measured in the extracellular medium. Although the major contributor of the extracellular signal is the synaptic transmembrane current, other sources--including Na(+) and Ca(2+) spikes, ionic fluxes through voltage- and ligand-gated channels, and intrinsic membrane oscillations--can substantially shape the extracellular field. High-density recordings of field activity in animals and subdural grid recordings in humans, combined with recently developed data processing tools and computational modelling, can provide insight into the cooperative behaviour of neurons, their average synaptic input and their spiking output, and can increase our understanding of how these processes contribute to the extracellular signal.

Scale-free modulation of resting-state neuronal oscillations reflects prolonged brain maturation in humans

Human neuronal circuits undergo life-long functional reorganization with profound effects on cognition and behavior. Well documented prolonged development of anatomical brain structures includes white and gray matter changes that continue into the third decade of life. We investigated resting-state EEG oscillations in 1433 subjects from 5 to 71 years. Neuronal oscillations exhibit scale-free amplitude modulation as reflected in power-law decay of autocorrelations--also known as long-range temporal correlations (LRTC)--which was assessed by detrended fluctuation analysis. We observed pronounced increases in LRTC from childhood to adolescence, during adolescence, and even into early adulthood (∼25 years of age) after which the temporal structure stabilized. A principal component analysis of the spatial distribution of LRTC revealed increasingly uniform scores across the scalp. Together, these findings indicate that the scale-free modulation of resting-state oscillations reflects brain maturation, and suggests that scaling analysis may prove useful as a biomarker of pathophysiology in neurodevelopmental disorders such as attention deficit hyperactivity disorder and schizophrenia.

CORTICAL PHASE TRANSITIONS: PROPERTIES DEMONSTRATED IN CONTINUUM SIMULATIONS AT MESOSCOPIC AND MACROSCOPIC SCALES

Continuum simulations of cortical dynamics permit consistent simulations to be performed at different spatial scales, using scale-adjusted parameter values. Properties of the simulations described here accord with Freeman's experimental and theoretical findings on gamma synchrony, phase transition, phase cones, and null spikes. State equations include effects of retrograde action potential propagation into dendritic trees, and kinetics of AMPA, GABA, and NMDA receptors. Realistic field potentials and pulse rates, gamma resonance and oscillation, and 1/f2 background activity are obtained. Zero-lag synchrony and traveling waves occur as complementary aspects of cortical transmission, and lead/lag relations between excitatory and inhibitory cell populations vary systematically around transition to autonomous gamma oscillation. Autonomous gamma is initiated by focal excitation of excitatory cells and suppressed by laterally spreading trans-cortical excitation. By implication, patches of cortex excited to gamma oscillation can mutually synchronize into larger fields, self-organized into sequences by mutual negative feedback relations, while the sequence of synchronous fields is regulated both by cortical/subcortical interactions and by traveling waves in the cortex — the latter observable as phase cones. At a critical level of cortical excitation, just before transition to autonomous gamma, patches of cortex exhibit selective sensitivity to action potential pulse trains modulated in the gamma band, while autonomous gamma releases pulse trains modulated in the same band, implying coupling of input and output modes. Transition between input and output modes may be heralded by phase slips and null spikes. Synaptic segregation by retrograde action potential propagation implies state-specific synaptic information storage.

COHERENT STATES, FRACTALS AND BRAIN WAVES

I show that a functional representation of self-similarity (as the one occurring in fractals) is provided by squeezed coherent states. In this way, the dissipative model of brain is shown to account for the self-similarity in brain background activity suggested by power-law distributions of power spectral densities of electrocorticograms. I also briefly discuss the action-perception cycle in the dissipative model with reference to intentionality in terms of trajectories in the memory state space.

Estimation of electrode location in a rat motor cortex by laminar analysis of electrophysiology and intracortical electrical stimulation

While the development of microelectrode arrays has enabled access to disparate regions of a cortex for neurorehabilitation, neuroprosthetic and basic neuroscience research, accurate interpretation of the signals and manipulation of the cortical neurons depend upon the anatomical placement of the electrode arrays in a layered cortex. Toward this end, this report compares two in vivo methods for identifying the placement of electrodes in a linear array spaced 100 µm apart based on in situ laminar analysis of (1) ketamine-xylazine-induced field potential oscillations in a rat motor cortex and (2) an intracortical electrical stimulation-induced movement threshold. The first method is based on finding the polarity reversal in laminar oscillations which is reported to appear at the transition between layers IV and V in laminar 'high voltage spindles' of the rat cortical column. Analysis of histological images in our dataset indicates that polarity reversal is detected 150.1 ± 104.2 µm below the start of layer V. The second method compares the intracortical microstimulation currents that elicit a physical movement for anodic versus cathodic stimulation. It is based on the hypothesis that neural elements perpendicular to the electrode surface are preferentially excited by anodic stimulation while cathodic stimulation excites those with a direction component parallel to its surface. With this method, we expect to see a change in the stimulation currents that elicits a movement at the beginning of layer V when comparing anodic versus cathodic stimulation as the upper cortical layers contain neuronal structures that are primarily parallel to the cortical surface and lower layers contain structures that are primarily perpendicular. Using this method, there was a 78.7 ± 68 µm offset in the estimate of the depth of the start of layer V. The polarity reversal method estimates the beginning of layer V within ±90 µm with 95% confidence and the intracortical stimulation method estimates it within ±69.3 µm. We propose that these methods can be used to estimate the in situ location of laminar electrodes implanted in the rat motor cortex.

Critical issues in state-of-the-art brain-computer interface signal processing

This paper reviews several critical issues facing signal processing for brain-computer interfaces (BCIs) and suggests several recent approaches that should be further examined. The topics were selected based on discussions held during the 4th International BCI Meeting at a workshop organized to review and evaluate the current state of, and issues relevant to, feature extraction and translation of field potentials for BCIs. The topics presented in this paper include the relationship between electroencephalography and electrocorticography, novel features for performance prediction, time-embedded signal representations, phase information, signal non-stationarity, and unsupervised adaptation.

Single pulse TMS-induced modulations of resting brain neurodynamics encoded in EEG phase

Integration of electroencephalographic (EEG) recordings and transcranial magnetic stimulation (TMS) provides a useful framework for quantifying stimulation-induced modulations of neural dynamics. Amplitude and frequency modulations by different TMS protocols have been previously investigated, but the study of stimulation-induced effects on EEG phase has been more limited. We examined changes in resting brain dynamics following single TMS pulses, focusing on measures in the phase domain, to assess their sensitivity to stimulation effects. We observed a significant, approximately global increase in EEG relative phase following prolonged (>20 min) single-pulse TMS. In addition, we estimated higher rates of phase fluctuation from the slope of estimated phase curves, and higher numbers of phase resetting intervals following TMS over motor cortex, particularly in frontal and centro-parietal/parietal channels. Phase changes were only significantly different from their pre-TMS values at the end of the stimulation session, which suggests that prolonged single-pulse TMS may result in cumulative changes in neural activity reflected in the phase of the EEG. This is a novel result, as prior studies have reported only transient stimulation-related effects in the amplitude and frequency domains following single-pulse TMS.

Behavioral states may be associated with distinct spatial patterns in electrocorticogram

To determine if behavioral states are associated with unique spatial electrocorticographic (ECoG) patterns, we obtained recordings with a microgrid electrode array applied to the cortical surface of a human subject. The array was constructed with the intent of extracting maximal spatial information by optimizing interelectrode distances. A 34-year-old patient with intractable epilepsy underwent intracranial ECoG monitoring after standard methods failed to reveal localization of seizures. During the 8-day period of invasive recording, in addition to standard clinical electrodes a square 1 × 1 cm microgrid array with 64 electrodes (1.25 mm separation) was placed on the right inferior temporal gyrus. Careful review of video recordings identified four extended naturalistic behaviors: reading, conversing on the telephone, looking at photographs, and face-to-face interactions. ECoG activity recorded with the microgrid that corresponded to these behaviors was collected and ECoG spatial patterns were analyzed. During periods of ECoG selected for analysis, no electrographic seizures or epileptiform patterns were present. Moments of maximal spatial variance are shown to cluster by behavior. Comparisons between conditions using a permutation test reveal significantly different spatial patterns for each behavior. We conclude that ECoG recordings obtained on the cortical surface with optimal high spatial frequency resolution reveal distinct local spatial patterns that reflect different behavioral states, and we predict that similar patterns will be found in many if not most cortical areas on which a microgrid is placed.

Utility of independent component analysis for interpretation of intracranial EEG

Electrode arrays are sometimes implanted in the brains of patients with intractable epilepsy to better localize seizure foci before epilepsy surgery. Analysis of intracranial EEG (iEEG) recordings is typically performed in the electrode channel domain without explicit separation of the sources that generate the signals. However, intracranial EEG signals, like scalp EEG signals, could be linear mixtures of local activity and volume-conducted activity arising in multiple source areas. Independent component analysis (ICA) has recently been applied to scalp EEG data, and shown to separate the signal mixtures into independently generated brain and non-brain source signals. Here, we applied ICA to unmix source signals from intracranial EEG recordings from four epilepsy patients during a visually cued finger movement task in the presence of background pathological brain activity. This ICA decomposition demonstrated that the iEEG recordings were not maximally independent, but rather are linear mixtures of activity from multiple sources. Many of the independent component (IC) projections to the iEEG recording grid were consistent with sources from single brain regions, including components exhibiting classic movement-related dynamics. Notably, the largest IC projection to each channel accounted for no more than 20–80% of the channel signal variance, implying that in general intracranial recordings cannot be accurately interpreted as recordings of independent brain sources. These results suggest that ICA can be used to identify and monitor major field sources of local and distributed functional networks generating iEEG data. ICA decomposition methods are useful for improving the fidelity of source signals of interest, likely including distinguishing the sources of pathological brain activity.

Attractor dynamics and thermodynamic analogies in the cerebral cortex: synchronous oscillation, the background EEG, and the regulation of attention

Ongoing changes in attention and cognition depend upon cortical/subcortical interactions, which select sequences of different spatial patterns of activation in the cortex. It is proposed that each pattern of cortical activation permits evolution of electrocortical wave activity toward statistically stationary states, analogous to thermodynamic equilibrium. In each steady-state, neurons fire with an intrinsic Poisson spike probability and also with a bursting pattern related to network oscillations. Excitatory cell dendrites act as a regenerative reservoir in which pulse generation is balanced against dissipations. Equilibria exhibit contrasting limits. One limit, at high cortical activation, generates widespread zero-lag synchrony among excitatory cells, with partial suppression of noise. Excitatory and inhibitory cells approach zero-lag local correlation, with 1/4 cycle lag-correlation at greater distances of separation. The high-activation limit defines a correlated system of attractor basins, capable of co-ordinating synaptic modifications and intracortical signal generation. Suppression of noise would enhance convergence about attractor basins in the manner of simulated annealing, while, conversely, the persistence of some noise prevents network paralysis by phase locking. At the opposite limit-that of low activation-spikes and waves have low cross- and auto-correlation, but have wide-spectrum sensitivity to inputs. It is hypothesised that cortical regions, transiently at equilibrium near these extremes, engage in interaction with each other and with subcortical systems, to generate ongoing sequences of attention and cognition. This account is compatible with classical and recently observed experimental phenomena. The principle features inferred from a simplified linear mathematical account are reproduced in a more physiologically realistic and non-linear numerical simulation.