Multidimensional recording (MDR) and data sharing: an ecological open research and educational platform for neuroscience

The Thermodynamics of Mind

To not only survive, but also thrive, the brain must efficiently orchestrate distributed computation across space and time. This requires hierarchical organisation facilitating fast information transfer and processing at the lowest possible metabolic cost. Quantifying brain hierarchy is difficult but can be estimated from the asymmetry of information flow. Thermodynamics has successfully characterised hierarchy in many other complex systems. Here, we propose the 'Thermodynamics of Mind' framework as a natural way to quantify hierarchical brain orchestration and its underlying mechanisms. This has already provided novel insights into the orchestration of hierarchy in brain states including movie watching, where the hierarchy of the brain is flatter than during rest. Overall, this framework holds great promise for revealing the orchestration of cognition.

Identifying oscillatory brain networks with hidden Gaussian graphical spectral models of MEEG

Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflection of functional brain network activity with a temporal resolution that no other in vivo neuroimage may provide. A second step estimating functional connectivity from such activity pseudodata unveil the oscillatory brain networks that strongly correlate with all cognition and behavior. Simulations of such MEG or EEG inverse problem also reveal estimation errors of the functional connectivity determined by any of the state-of-the-art inverse solutions. We disclose a significant cause of estimation errors originating from misspecification of the functional network model incorporated into either inverse solution steps. We introduce the Bayesian identification of a Hidden Gaussian Graphical Spectral (HIGGS) model specifying such oscillatory brain networks model. In human EEG alpha rhythm simulations, the estimation errors measured as ROC performance do not surpass 2% in our HIGGS inverse solution and reach 20% in state-of-the-art methods. Macaque simultaneous EEG/ECoG recordings provide experimental confirmation for our results with 1/3 times larger congruence according to Riemannian distances than state-of-the-art methods.

Hierarchical fluctuation shapes a dynamic flow linked to states of consciousness

Consciousness arises from the spatiotemporal neural dynamics, however, its relationship with neural flexibility and regional specialization remains elusive. We identified a consciousness-related signature marked by shifting spontaneous fluctuations along a unimodal-transmodal cortical axis. This simple signature is sensitive to altered states of consciousness in single individuals, exhibiting abnormal elevation under psychedelics and in psychosis. The hierarchical dynamic reflects brain state changes in global integration and connectome diversity under task-free conditions. Quasi-periodic pattern detection revealed that hierarchical heterogeneity as spatiotemporally propagating waves linking to arousal. A similar pattern can be observed in macaque electrocorticography. Furthermore, the spatial distribution of principal cortical gradient preferentially recapitulated the genetic transcription levels of the histaminergic system and that of the functional connectome mapping of the tuberomammillary nucleus, which promotes wakefulness. Combining behavioral, neuroimaging, electrophysiological, and transcriptomic evidence, we propose that global consciousness is supported by efficient hierarchical processing constrained along a low-dimensional macroscale gradient.

Time-resolved correlation of distributed brain activity tracks E-I balance and accounts for diverse scale-free phenomena

Much of systems neuroscience posits the functional importance of brain activity patterns that lack natural scales of sizes, durations, or frequencies. The field has developed prominent, and sometimes competing, explanations for the nature of this scale-free activity. Here, we reconcile these explanations across species and modalities. First, we link estimates of excitation-inhibition (E-I) balance with time-resolved correlation of distributed brain activity. Second, we develop an unbiased method for sampling time series constrained by this time-resolved correlation. Third, we use this method to show that estimates of E-I balance account for diverse scale-free phenomena without need to attribute additional function or importance to these phenomena. Collectively, our results simplify existing explanations of scale-free brain activity and provide stringent tests on future theories that seek to transcend these explanations.

Structured tracking of alcohol reinforcement (STAR) for basic and translational alcohol research

There is inherent tension between methodologies developed to address basic research questions in model species and those intended for preclinical to clinical translation: basic investigations require flexibility of experimental design as hypotheses are rapidly tested and revised, whereas preclinical models emphasize standardized protocols and specific outcome measures. This dichotomy is particularly relevant in alcohol research, which spans a diverse range of basic sciences in addition to intensive efforts towards understanding the pathophysiology of alcohol use disorder (AUD). To advance these goals there is a great need for approaches that facilitate synergy across basic and translational areas of nonhuman alcohol research. In male and female mice, we establish a modular alcohol reinforcement paradigm: Structured Tracking of Alcohol Reinforcement (STAR). STAR provides a robust platform for quantitative assessment of AUD-relevant behavioral domains within a flexible framework that allows direct crosstalk between translational and mechanistically oriented studies. To achieve cross-study integration, despite disparate task parameters, a straightforward multivariate phenotyping analysis is used to classify subjects based on propensity for heightened alcohol consumption and insensitivity to punishment. Combining STAR with extant preclinical alcohol models, we delineate longitudinal phenotype dynamics and reveal putative neuro-biomarkers of heightened alcohol use vulnerability via neurochemical profiling of cortical and brainstem tissues. Together, STAR allows quantification of time-resolved biobehavioral processes essential for basic research questions simultaneous with longitudinal phenotyping of clinically relevant outcomes, thereby providing a framework to facilitate cohesion and translation in alcohol research.

Nine decades of electrocorticography: A comparison between epidural and subdural recordings

In recent years, electrocorticography (ECoG) has arisen as a neural signal recording tool in the development of clinically viable neural interfaces. ECoG electrodes are generally placed below the dura mater (subdural) but can also be placed on top of the dura (epidural). In deciding which of these modalities best suits long-term implants, complications and signal quality are important considerations. Conceptually, epidural placement may present a lower risk of complications as the dura is left intact but also a lower signal quality due to the dura acting as a signal attenuator. The extent to which complications and signal quality are affected by the dura, however, has been a matter of debate. To improve our understanding of the effects of the dura on complications and signal quality, we conducted a literature review. We inventorized the effect of the dura on signal quality, decodability and longevity of acute and chronic ECoG recordings in humans and non-human primates. Also, we compared the incidence and nature of serious complications in studies that employed epidural and subdural ECoG. Overall, we found that, even though epidural recordings exhibit attenuated signal amplitude over subdural recordings, particularly for high-density grids, the decodability of epidural recorded signals does not seem to be markedly affected. Additionally, we found that the nature of serious complications was comparable between epidural and subdural recordings. These results indicate that both epidural and subdural ECoG may be suited for long-term neural signal recordings, at least for current generations of clinical and high-density ECoG grids.

Spectral density-based clustering algorithms for complex networks

Introduction

Clustering is usually the first exploratory analysis step in empirical data. When the data set comprises graphs, the most common approaches focus on clustering its vertices. In this work, we are interested in grouping networks with similar connectivity structures together instead of grouping vertices of the graph. We could apply this approach to functional brain networks (FBNs) for identifying subgroups of people presenting similar functional connectivity, such as studying a mental disorder. The main problem is that real-world networks present natural fluctuations, which we should consider.

Methods

In this context, spectral density is an exciting feature because graphs generated by different models present distinct spectral densities, thus presenting different connectivity structures. We introduce two clustering methods: k-means for graphs of the same size and gCEM, a model-based approach for graphs of different sizes. We evaluated their performance in toy models. Finally, we applied them to FBNs of monkeys under anesthesia and a dataset of chemical compounds.

Results

We show that our methods work well in both toy models and real-world data. They present good results for clustering graphs presenting different connectivity structures even when they present the same number of edges, vertices, and degree of centrality.

Discussion

We recommend using k-means-based clustering for graphs when graphs present the same number of vertices and the gCEM method when graphs present a different number of vertices.

Minimizing the distortions in electrophysiological source imaging of cortical oscillatory activity via Spectral Structured Sparse Bayesian Learning

Oscillatory processes at all spatial scales and on all frequencies underpin brain function. Electrophysiological Source Imaging (ESI) is the data-driven brain imaging modality that provides the inverse solutions to the source processes of the EEG, MEG, or ECoG data. This study aimed to carry out an ESI of the source cross-spectrum while controlling common distortions of the estimates. As with all ESI-related problems under realistic settings, the main obstacle we faced is a severely ill-conditioned and high-dimensional inverse problem. Therefore, we opted for Bayesian inverse solutions that posited a priori probabilities on the source process. Indeed, rigorously specifying both the likelihoods and a priori probabilities of the problem leads to the proper Bayesian inverse problem of cross-spectral matrices. These inverse solutions are our formal definition for cross-spectral ESI (cESI), which requires a priori of the source cross-spectrum to counter the severe ill-condition and high-dimensionality of matrices. However, inverse solutions for this problem were NP-hard to tackle or approximated within iterations with bad-conditioned matrices in the standard ESI setup. We introduce cESI with a joint a priori probability upon the source cross-spectrum to avoid these problems. cESI inverse solutions are low-dimensional ones for the set of random vector instances and not random matrices. We achieved cESI inverse solutions through the variational approximations via our Spectral Structured Sparse Bayesian Learning (ssSBL) algorithm https://github.com/CCC-members/Spectral-Structured-Sparse-Bayesian-Learning. We compared low-density EEG (10-20 system) ssSBL inverse solutions with reference cESIs for two experiments: (a) high-density MEG that were used to simulate EEG and (b) high-density macaque ECoG that were recorded simultaneously with EEG. The ssSBL resulted in two orders of magnitude with less distortion than the state-of-the-art ESI methods. Our cESI toolbox, including the ssSBL method, is available at https://github.com/CCC-members/BC-VARETA_Toolbox.

The INSIDEOUT framework provides precise signatures of the balance of intrinsic and extrinsic dynamics in brain states

Finding precise signatures of different brain states is a central, unsolved question in neuroscience. We reformulated the problem to quantify the 'inside out' balance of intrinsic and extrinsic brain dynamics in brain states. The difference in brain state can be described as differences in the detailed causal interactions found in the underlying intrinsic brain dynamics. We used a thermodynamics framework to quantify the breaking of the detailed balance captured by the level of asymmetry in temporal processing, i.e. the arrow of time. Specifically, the temporal asymmetry was computed by the time-shifted correlation matrices for the forward and reversed time series, reflecting the level of non-reversibility/non-equilibrium. We found precise, distinguishing signatures in terms of the reversibility and hierarchy of large-scale dynamics in three radically different brain states (awake, deep sleep and anaesthesia) in electrocorticography data from non-human primates. Significantly lower levels of reversibility were found in deep sleep and anaesthesia compared to wakefulness. Non-wakeful states also showed a flatter hierarchy, reflecting the diversity of the reversibility across the brain. Overall, this provides signatures of the breaking of detailed balance in different brain states, perhaps reflecting levels of conscious awareness.

Network Analysis of Time Series: Novel Approaches to Network Neuroscience

In the last two decades, there has been an explosion of interest in modeling the brain as a network, where nodes correspond variously to brain regions or neurons, and edges correspond to structural or statistical dependencies between them. This kind of network construction, which preserves spatial, or structural, information while collapsing across time, has become broadly known as "network neuroscience." In this work, we provide an alternative application of network science to neural data: network-based analysis of non-linear time series and review applications of these methods to neural data. Instead of preserving spatial information and collapsing across time, network analysis of time series does the reverse: it collapses spatial information, instead preserving temporally extended dynamics, typically corresponding to evolution through some kind of phase/state-space. This allows researchers to infer a, possibly low-dimensional, "intrinsic manifold" from empirical brain data. We will discuss three methods of constructing networks from nonlinear time series, and how to interpret them in the context of neural data: recurrence networks, visibility networks, and ordinal partition networks. By capturing typically continuous, non-linear dynamics in the form of discrete networks, we show how techniques from network science, non-linear dynamics, and information theory can extract meaningful information distinct from what is normally accessible in standard network neuroscience approaches.

Broadband Dynamics Rather than Frequency-Specific Rhythms Underlie Prediction Error in the Primate Auditory Cortex

Detection of statistical irregularities, measured as a prediction error response, is fundamental to the perceptual monitoring of the environment. We studied whether prediction error response is associated with neural oscillations or asynchronous broadband activity. Electrocorticography was conducted in three male monkeys, who passively listened to the auditory roving oddball stimuli. Local field potentials (LFPs) recorded over the auditory cortex underwent spectral principal component analysis, which decoupled broadband and rhythmic components of the LFP signal. We found that the broadband component captured the prediction error response, whereas none of the rhythmic components were associated with statistical irregularities of sounds. The broadband component displayed more stochastic, asymmetrical multifractal properties than the rhythmic components, which revealed more self-similar dynamics. We thus conclude that the prediction error response is captured by neuronal populations generating asynchronous broadband activity, defined by irregular dynamic states, which, unlike oscillatory rhythms, appear to enable the neural representation of auditory prediction error response.SIGNIFICANCE STATEMENT This study aimed to examine the contribution of oscillatory and asynchronous components of auditory local field potentials in the generation of prediction error responses to sensory irregularities, as this has not been directly addressed in the previous studies. Here, we show that mismatch negativity-an auditory prediction error response-is driven by the asynchronous broadband component of potentials recorded in the auditory cortex. This finding highlights the importance of nonoscillatory neural processes in the predictive monitoring of the environment. At a more general level, the study demonstrates that stochastic neural processes, which are often disregarded as neural noise, do have a functional role in the processing of sensory information.

Roving oddball paradigm elicits sensory gating, frequency sensitivity, and long-latency response in common marmosets

Mismatch negativity (MMN) is a candidate biomarker for neuropsychiatric disease. Understanding the extent to which it reflects cognitive deviance-detection or purely sensory processes will assist practitioners in making informed clinical interpretations. This study compares the utility of deviance-detection and sensory-processing theories for describing MMN-like auditory responses of a common marmoset monkey during roving oddball stimulation. The following exploratory analyses were performed on an existing dataset: responses during the transition and repetition sequence of the roving oddball paradigm (standard -> deviant/S1 -> S2 -> S3) were compared; long-latency potentials evoked by deviant stimuli were examined using a double-epoch waveform subtraction; effects of increasing stimulus repetitions on standard and deviant responses were analyzed; and transitions between standard and deviant stimuli were divided into ascending and descending frequency changes to explore contributions of frequency-sensitivity. An enlarged auditory response to deviant stimuli was observed. This decreased exponentially with stimulus repetition, characteristic of sensory gating. A slow positive deflection was viewed over approximately 300–800 ms after the deviant stimulus, which is more difficult to ascribe to afferent sensory mechanisms. When split into ascending and descending frequency transitions, the resulting difference waveforms were disproportionally influenced by descending frequency deviant stimuli. This asymmetry is inconsistent with the general deviance-detection theory of MMN. These findings tentatively suggest that MMN-like responses from common marmosets are predominantly influenced by rapid sensory adaptation and frequency preference of the auditory cortex, while deviance-detection may play a role in long-latency activity.

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.

Global waves synchronize the brain's functional systems with fluctuating arousal

We propose and empirically support a parsimonious account of intrinsic, brain-wide spatiotemporal organization arising from traveling waves linked to arousal. We hypothesize that these waves are the predominant physiological process reflected in spontaneous functional magnetic resonance imaging (fMRI) signal fluctuations. The correlation structure ("functional connectivity") of these fluctuations recapitulates the large-scale functional organization of the brain. However, a unifying physiological account of this structure has so far been lacking. Here, using fMRI in humans, we show that ongoing arousal fluctuations are associated with global waves of activity that slowly propagate in parallel throughout the neocortex, thalamus, striatum, and cerebellum. We show that these waves can parsimoniously account for many features of spontaneous fMRI signal fluctuations, including topographically organized functional connectivity. Last, we demonstrate similar, cortex-wide propagation of neural activity measured with electrocorticography in macaques. These findings suggest that traveling waves spatiotemporally pattern brain-wide excitability in relation to arousal.

Topological analysis of differential effects of ketamine and propofol anaesthesia on brain dynamics

Research has found that the vividness of conscious experience is related to brain dynamics. Despite both being anaesthetics, propofol and ketamine produce different subjective states: we explore the different effects of these two anaesthetics on the structure of dynamic attractors reconstructed from electrophysiological activity recorded from cerebral cortex of two macaques. We used two methods: the first embeds the recordings in a continuous high-dimensional manifold on which we use topological data analysis to infer the presence of higher-order dynamics. The second reconstruction, an ordinal partition network embedding, allows us to create a discrete state-transition network, which is amenable to information-theoretic analysis and contains rich information about state-transition dynamics. We find that the awake condition generally had the 'richest' structure, visiting the most states, the presence of pronounced higher-order structures, and the least deterministic dynamics. By contrast, the propofol condition had the most dissimilar dynamics, transitioning to a more impoverished, constrained, low-structure regime. The ketamine condition, interestingly, seemed to combine aspects of both: while it was generally less complex than the awake condition, it remained well above propofol in almost all measures. These results provide deeper and more comprehensive insights than what is typically gained by using point-measures of complexity.

Robust Autoregression with Exogenous Input Model for System Identification and Predicting

Autoregression with exogenous input (ARX) is a widely used model to estimate the dynamic relationships between neurophysiological signals and other physiological parameters. Nevertheless, biological signals, such as electroencephalogram (EEG), arterial blood pressure (ABP), and intracranial pressure (ICP), are inevitably contaminated by unexpected artifacts, which may distort the parameter estimation due to the use of the L2 norm structure. In this paper, we defined the ARX in the Lp (p ≤ 1) norm space with the aim of resisting outlier influence and designed a feasible iteration procedure to estimate model parameters. A quantitative evaluation with various outlier conditions demonstrated that the proposed method could estimate ARX parameters more robustly than conventional methods. Testing with the resting-state EEG with ocular artifacts demonstrated that the proposed method could predict missing data with less influence from the artifacts. In addition, the results on ICP and ABP data further verified its efficiency for model fitting and system identification. The proposed Lp-ARX may help capture system parameters reliably with various input and output signals that are contaminated with artifacts.

MacaquePose: A Novel "In the Wild" Macaque Monkey Pose Dataset for Markerless Motion Capture

Video-based markerless motion capture permits quantification of an animal's pose and motion, with a high spatiotemporal resolution in a naturalistic context, and is a powerful tool for analyzing the relationship between the animal's behaviors and its brain functions. Macaque monkeys are excellent non-human primate models, especially for studying neuroscience. Due to the lack of a dataset allowing training of a deep neural network for the macaque's markerless motion capture in the naturalistic context, it has been challenging to apply this technology for macaques-based studies. In this study, we created MacaquePose, a novel open dataset with manually labeled body part positions (keypoints) for macaques in naturalistic scenes, consisting of >13,000 images. We also validated the application of the dataset by training and evaluating an artificial neural network with the dataset. The results indicated that the keypoint estimation performance of the trained network was close to that of a human-level. The dataset will be instrumental to train/test the neural networks for markerless motion capture of the macaques and developments of the algorithms for the networks, contributing establishment of an innovative platform for behavior analysis for non-human primates for neuroscience and medicine, as well as other fields using macaques as a model organism.

Differential effects of propofol and ketamine on critical brain dynamics

Whether the brain operates at a critical "tipping" point is a long standing scientific question, with evidence from both cellular and systems-scale studies suggesting that the brain does sit in, or near, a critical regime. Neuroimaging studies of humans in altered states of consciousness have prompted the suggestion that maintenance of critical dynamics is necessary for the emergence of consciousness and complex cognition, and that reduced or disorganized consciousness may be associated with deviations from criticality. Unfortunately, many of the cellular-level studies reporting signs of criticality were performed in non-conscious systems (in vitro neuronal cultures) or unconscious animals (e.g. anaesthetized rats). Here we attempted to address this knowledge gap by exploring critical brain dynamics in invasive ECoG recordings from multiple sessions with a single macaque as the animal transitioned from consciousness to unconsciousness under different anaesthetics (ketamine and propofol). We use a previously-validated test of criticality: avalanche dynamics to assess the differences in brain dynamics between normal consciousness and both drug-states. Propofol and ketamine were selected due to their differential effects on consciousness (ketamine, but not propofol, is known to induce an unusual state known as "dissociative anaesthesia"). Our analyses indicate that propofol dramatically restricted the size and duration of avalanches, while ketamine allowed for more awake-like dynamics to persist. In addition, propofol, but not ketamine, triggered a large reduction in the complexity of brain dynamics. All states, however, showed some signs of persistent criticality when testing for exponent relations and universal shape-collapse. Further, maintenance of critical brain dynamics may be important for regulation and control of conscious awareness.

Neuronal timescales are functionally dynamic and shaped by cortical microarchitecture

Complex cognitive functions such as working memory and decision-making require information maintenance over seconds to years, from transient sensory stimuli to long-term contextual cues. While theoretical accounts predict the emergence of a corresponding hierarchy of neuronal timescales, direct electrophysiological evidence across the human cortex is lacking. Here, we infer neuronal timescales from invasive intracranial recordings. Timescales increase along the principal sensorimotor-to-association axis across the entire human cortex, and scale with single-unit timescales within macaques. Cortex-wide transcriptomic analysis shows direct alignment between timescales and expression of excitation- and inhibition-related genes, as well as genes specific to voltage-gated transmembrane ion transporters. Finally, neuronal timescales are functionally dynamic: prefrontal cortex timescales expand during working memory maintenance and predict individual performance, while cortex-wide timescales compress with aging. Thus, neuronal timescales follow cytoarchitectonic gradients across the human cortex and are relevant for cognition in both short and long terms, bridging microcircuit physiology with macroscale dynamics and behavior.

Efficient search for informational cores in complex systems: Application to brain networks

An important step in understanding the nature of the brain is to identify "cores" in the brain network, where brain areas strongly interact with each other. Cores can be considered as essential sub-networks for brain functions. In the last few decades, an information-theoretic approach to identifying cores has been developed. In this approach, interactions between parts are measured by an information loss function, which quantifies how much information would be lost if interactions between parts were removed. Then, a core called a "complex" is defined as a subsystem wherein the amount of information loss is locally maximal. Although identifying complexes can be a novel and useful approach, its application is practically impossible because computation time grows exponentially with system size. Here we propose a fast and exact algorithm for finding complexes, called Hierarchical Partitioning for Complex search (HPC). HPC hierarchically partitions systems to narrow down candidates for complexes. The computation time of HPC is polynomial, enabling us to find complexes in large systems (up to several hundred) in a practical amount of time. We prove that HPC is exact when an information loss function satisfies a mathematical property, monotonicity. We show that mutual information is one such information loss function. We also show that a broad class of submodular functions can be considered as such information loss functions, indicating the expandability of our framework to the class. We applied HPC to electrocorticogram recordings from a monkey and demonstrated that HPC revealed temporally stable and characteristic complexes.

Propofol Anesthesia Increases Long-range Frontoparietal Corticocortical Interaction in the Oculomotor Circuit in Macaque Monkeys

WHAT WE ALREADY KNOW ABOUT THIS TOPIC

A decrease in frontoparietal functional connectivity has been demonstrated with multiple anesthetic agents, and this decrease has been proposed as a final common functional pathway to produce anesthesia.Two alternative measures of long-range cortical interaction are coherence and phase-amplitude coupling. Although phase-amplitude coupling within frontal cortex changes with propofol administration, the effects of propofol on phase-amplitude coupling between different cortical areas have not previously been reported.

WHAT THIS ARTICLE TELLS US THAT IS NEW

Using a previously published monkey electrocorticography data set, it was found that propofol induced coherent slow oscillations in visual and oculomotor networks made up of cortical areas with strong anatomic projections.Frontal eye field within-area phase-amplitude coupling increased.Contrary to expectations from previous functional connectivity studies, interareal phase-amplitude coupling also increased with propofol.

BACKGROUND

Frontoparietal functional connectivity decreases with multiple anesthetics using electrophysiology and functional imaging. This decrease has been proposed as a final common functional pathway to produce anesthesia. Two alternative measures of long-range cortical interaction are coherence and phase-amplitude coupling. Although phase-amplitude coupling within frontal cortex changes with propofol administration, the effects of propofol on phase-amplitude coupling between different cortical areas have not previously been reported. Based on phase-amplitude coupling observed within frontal lobe during the anesthetized period, it was hypothesized that between-lead phase-amplitude coupling analysis should decrease between frontal and parietal leads during propofol anesthesia.

METHODS

A published monkey electrocorticography data set (N = 2 animals) was used to test for interactions in the cortical oculomotor circuit, which is robustly interconnected in primates, and in the visual system during propofol anesthesia using coherence and interarea phase-amplitude coupling.

RESULTS

Propofol induces coherent slow oscillations in visual and oculomotor networks made up of cortical areas with strong anatomic projections. Frontal eye field within-area phase-amplitude coupling increases with a time course consistent with a bolus response to intravenous propofol (modulation index increase of 12.6-fold). Contrary to the hypothesis, interareal phase-amplitude coupling also increases with propofol, with the largest increase in phase-amplitude coupling in frontal eye field low-frequency phase modulating lateral intraparietal area β-power (27-fold increase) and visual area 2 low-frequency phase altering visual area 1 β-power (19-fold increase).

CONCLUSIONS

Propofol anesthesia induces coherent oscillations and increases certain frontoparietal interactions in oculomotor cortices. Frontal eye field and lateral intraparietal area show increased coherence and phase-amplitude coupling. Visual areas 2 and 1, which have similar anatomic projection patterns, show similar increases in phase-amplitude coupling, suggesting higher order feedback increases in influence during propofol anesthesia relative to wakefulness. This suggests that functional connectivity between frontal and parietal areas is not uniformly decreased by anesthetics.

Complexity of brain activity and connectivity in functional neuroimaging

Understanding the complexity of human brain dynamics and brain connectivity across the repertoire of functional neuroimaging and various conditions, is of paramount importance. Novel measures should be designed tailored to the input focusing on multichannel activity and dynamic functional brain connectivity (DFBC). Here, we defined a novel complexity index (CI) from the field of symbolic dynamics that quantifies patterns of different words up to a length from a symbolic sequence. The CI characterizes the complexity of the brain activity. We analysed DFBC by adopting the sliding window approach using imaginary part of phase locking value (iPLV) for EEG/ECoG/MEG and wavelet coherence (WC) for fMRI. Both intra and cross-frequency couplings (CFC) namely phase-to-amplitude were estimated using iPLV/WC at every snapshot of the DFBC. Using proper surrogate analysis, we defined the dominant intrinsic coupling mode (DICM) per pair of regions-of-interest (ROI). The spatiotemporal probability distribution of DICM were reported to reveal the most prominent coupling modes per condition and modality. Finally, a novel flexibility index is defined that quantifies the transition of DICM per pair of ROIs between consecutive time windows. The whole methodology was demonstrated using four neuroimaging datasets (EEG/ECoG/MEG/fMRI). Finally, we succeeded to totally discriminate healthy controls from schizophrenic using FI and dynamic reconfiguration of DICM. Anaesthesia independently of the drug caused a global decreased of complexity in all frequency bands with the exception in δ and alters the dynamic reconfiguration of DICM. CI and DICM of MEG/fMRI resting-state recordings in two spatial scales were high reliable.

The generation and propagation of the human alpha rhythm

The alpha rhythm is the longest-studied brain oscillation and has been theorized to play a key role in cognition. Still, its physiology is poorly understood. In this study, we used microelectrodes and macroelectrodes in surgical epilepsy patients to measure the intracortical and thalamic generators of the alpha rhythm during quiet wakefulness. We first found that alpha in both visual and somatosensory cortex propagates from higher-order to lower-order areas. In posterior cortex, alpha propagates from higher-order anterosuperior areas toward the occipital pole, whereas alpha in somatosensory cortex propagates from associative regions toward primary cortex. Several analyses suggest that this cortical alpha leads pulvinar alpha, complicating prevailing theories of a thalamic pacemaker. Finally, alpha is dominated by currents and firing in supragranular cortical layers. Together, these results suggest that the alpha rhythm likely reflects short-range supragranular feedback, which propagates from higher- to lower-order cortex and cortex to thalamus. These physiological insights suggest how alpha could mediate feedback throughout the thalamocortical system.

EECoG-Comp: An Open Source Platform for Concurrent EEG/ECoG Comparisons-Applications to Connectivity Studies

Electrophysiological Source Imaging (ESI) is hampered by lack of "gold standards" for model validation. Concurrent electroencephalography (EEG) and electrocorticography (ECoG) experiments (EECoG) are useful for this purpose, especially primate models due to their flexibility and translational value for human research. Unfortunately, there is only one EECoG experiments in the public domain that we know of: the Multidimensional Recording (MDR) is based on a single monkey ( www.neurotycho.org ). The mining of this type of data is hindered by lack of specialized procedures to deal with: (1) Severe EECoG artifacts due to the experimental produces; (2) Sophisticated forward models that account for surgery induced skull defects and implanted ECoG electrode strips; (3) Reliable statistical procedures to estimate and compare source connectivity (partial correlation). We provide solutions to the processing issues just mentioned with EECoG-Comp: an open source platform ( https://github.com/Vincent-wq/EECoG-Comp ). EECoG lead fields calculated with FEM (Simbio) for MDR data are also provided and were used in other papers of this special issue. As a use case with the MDR, we show: (1) For real MDR data, 4 popular ESI methods (MNE, LCMV, eLORETA and SSBL) showed significant but moderate concordance with a usual standard, the ECoG Laplacian (standard partial [Formula: see text]); (2) In both monkey and human simulations, all ESI methods as well as Laplacian had a significant but poor correspondence with the true source connectivity. These preliminary results may stimulate the development of improved ESI connectivity estimators but require the availability of more EECoG data sets to obtain neurobiologically valid inferences.

Information integration in large brain networks

An outstanding problem in neuroscience is to understand how information is integrated across the many modules of the brain. While classic information-theoretic measures have transformed our understanding of feedforward information processing in the brain's sensory periphery, comparable measures for information flow in the massively recurrent networks of the rest of the brain have been lacking. To address this, recent work in information theory has produced a sound measure of network-wide "integrated information", which can be estimated from time-series data. But, a computational hurdle has stymied attempts to measure large-scale information integration in real brains. Specifically, the measurement of integrated information involves a combinatorial search for the informational "weakest link" of a network, a process whose computation time explodes super-exponentially with network size. Here, we show that spectral clustering, applied on the correlation matrix of time-series data, provides an approximate but robust solution to the search for the informational weakest link of large networks. This reduces the computation time for integrated information in large systems from longer than the lifespan of the universe to just minutes. We evaluate this solution in brain-like systems of coupled oscillators as well as in high-density electrocortigraphy data from two macaque monkeys, and show that the informational "weakest link" of the monkey cortex splits posterior sensory areas from anterior association areas. Finally, we use our solution to provide evidence in support of the long-standing hypothesis that information integration is maximized by networks with a high global efficiency, and that modular network structures promote the segregation of information.

Adaptive, locally linear models of complex dynamics

The dynamics of complex systems generally include high-dimensional, nonstationary, and nonlinear behavior, all of which pose fundamental challenges to quantitative understanding. To address these difficulties, we detail an approach based on local linear models within windows determined adaptively from data. While the dynamics within each window are simple, consisting of exponential decay, growth, and oscillations, the collection of local parameters across all windows provides a principled characterization of the full time series. To explore the resulting model space, we develop a likelihood-based hierarchical clustering, and we examine the eigenvalues of the linear dynamics. We demonstrate our analysis with the Lorenz system undergoing stable spiral dynamics and in the standard chaotic regime. Applied to the posture dynamics of the nematode Caenorhabditis elegans, our approach identifies fine-grained behavioral states and model dynamics which fluctuate about an instability boundary, and we detail a bifurcation in a transition from forward to backward crawling. We analyze whole-brain imaging in C. elegans and show that global brain dynamics is damped away from the instability boundary by a decrease in oxygen concentration. We provide additional evidence for such near-critical dynamics from the analysis of electrocorticography in monkey and the imaging of a neural population from mouse visual cortex at single-cell resolution.

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.

Efficient Algorithms for Searching the Minimum Information Partition in Integrated Information Theory

The ability to integrate information in the brain is considered to be an essential property for cognition and consciousness. Integrated Information Theory (IIT) hypothesizes that the amount of integrated information (Φ) in the brain is related to the level of consciousness. IIT proposes that, to quantify information integration in a system as a whole, integrated information should be measured across the partition of the system at which information loss caused by partitioning is minimized, called the Minimum Information Partition (MIP). The computational cost for exhaustively searching for the MIP grows exponentially with system size, making it difficult to apply IIT to real neural data. It has been previously shown that, if a measure of Φ satisfies a mathematical property, submodularity, the MIP can be found in a polynomial order by an optimization algorithm. However, although the first version of Φ is submodular, the later versions are not. In this study, we empirically explore to what extent the algorithm can be applied to the non-submodular measures of Φ by evaluating the accuracy of the algorithm in simulated data and real neural data. We find that the algorithm identifies the MIP in a nearly perfect manner even for the non-submodular measures. Our results show that the algorithm allows us to measure Φ in large systems within a practical amount of time.

A piecewise probabilistic regression model to decode hand movement trajectories from epidural and subdural ECoG signals

OBJECTIVE

The primary concern of this study is to develop a probabilistic regression method that would improve the decoding of the hand movement trajectories from epidural ECoG as well as from subdural ECoG signals.

APPROACH

The model is characterized by the conditional expectation of the hand position given the ECoG signals. The conditional expectation of the hand position is then modeled by a linear combination of the conditional probability density functions defined for each segment of the movement. Moreover, a spatial linear filter is proposed for reducing the dimension of the feature space. The spatial linear filter is applied to each frequency band of the ECoG signals and extract the features with highest decoding performance.

MAIN RESULTS

For evaluating the proposed method, a dataset including 28 ECoG recordings from four adult Japanese macaques is used. The results show that the proposed decoding method outperforms the results with respect to the state of the art methods using this dataset. The relative kinematic information of each frequency band is also investigated using mutual information and decoding performance. The decoding performance shows that the best performance was obtained for high gamma bands from 50 to 200 Hz as well as high frequency ECoG band from 200 to 400 Hz for subdural recordings. However, the decoding performance was decreased for these frequency bands using epidural recordings. The mutual information shows that, on average, the high gamma band from 50 to 200 Hz and high frequency ECoG band from 200 to 400 Hz contain significantly more information than the average of the rest of the frequency bands [Formula: see text] for both subdural and epidural recordings. The results of high resolution time-frequency analysis show that ERD/ERS patterns in all frequency bands could reveal the dynamics of the ECoG responses during the movement. The onset and offset of the movement can be clearly identified by the ERD/ERS patterns.

SIGNIFICANCE

Reliable decoding the kinematic information from the brain signals paves the way for robust control of external devices.

Efficient communication dynamics on macro-connectome, and the propagation speed

Global communication dynamics in the brain can be captured using fMRI, MEG, or electrocorticography (ECoG), and the global slow dynamics often represent anatomical constraints. Complementary single-/multi-unit recordings have described local fast temporal dynamics. However, global fast temporal dynamics remain incompletely understood with considering of anatomical constraints. Therefore, we compared temporal aspects of cross-area propagations of single-unit recordings and ECoG, and investigated their anatomical bases. First, we demonstrated how both evoked and spontaneous ECoGs can accurately predict latencies of single-unit recordings. Next, we estimated the propagation velocity (1.0–1.5 m/s) from brain-wide data and found that it was fairly stable among different conscious levels. We also found that the shortest paths in anatomical topology strongly predicted the latencies. Finally, we demonstrated that Communicability, a novel graph-theoretic measure, is able to quantify that more than 90% of paths should use shortest paths and the remaining are non-shortest walks. These results revealed that macro-connectome is efficiently wired for detailed communication dynamics in the brain.

Mapping Brain Activity with Electrocorticography: Resolution Properties and Robustness of Inverse Solutions

Electrocorticography (ECoG) is an electrophysiological technique that records brain activity directly from the cortical surface with high temporal (ms) and spatial (mm) resolution. Its major limitations are in the high invasiveness and in the restricted field-of-view of the electrode grid, which partially covers the cortex. To infer brain activity at locations different from just below the electrodes, it is necessary to solve the electromagnetic inverse problem. Limitations in the performance of source reconstruction algorithms from ECoG have been, to date, only partially addressed in the literature, and a systematic evaluation is still lacking. The main goal of this study is to provide a quantitative evaluation of resolution properties of widely used inverse methods (eLORETA and MNE) for various ECoG grid sizes, in terms of localization error, spatial dispersion, and overall amplitude. Additionally, this study aims at evaluating how the use of simultaneous electroencephalography (EEG) affects the above properties. For these purposes, we take advantage of a unique dataset in which a monkey underwent a simultaneous recording with a 128 channel ECoG grid and an 18 channel EEG grid. Our results show that, in general conditions, the reconstruction of cortical activity located more than 1 cm away from the ECoG grid is not accurate, since the localization error increases linearly with the distance from the electrodes. This problem can be partially overcome by recording simultaneously ECoG and EEG. However, this analysis enlightens the necessity to design inverse algorithms specifically targeted at taking into account the limited field-of-view of the ECoG grid.

Recursive Exponentially Weighted N-way Partial Least Squares Regression with Recursive-Validation of Hyper-Parameters in Brain-Computer Interface Applications

A tensor-input/tensor-output Recursive Exponentially Weighted N-Way Partial Least Squares (REW-NPLS) regression algorithm is proposed for high dimension multi-way (tensor) data treatment and adaptive modeling of complex processes in real-time. The method unites fast and efficient calculation schemes of the Recursive Exponentially Weighted PLS with the robustness of tensor-based approaches. Moreover, contrary to other multi-way recursive algorithms, no loss of information occurs in the REW-NPLS. In addition, the Recursive-Validation method for recursive estimation of the hyper-parameters is proposed instead of conventional cross-validation procedure. The approach was then compared to state-of-the-art methods. The efficiency of the methods was tested in electrocorticography (ECoG) and magnetoencephalography (MEG) datasets. The algorithms are implemented in software suitable for real-time operation. Although the Brain-Computer Interface applications are used to demonstrate the methods, the proposed approaches could be efficiently used in a wide range of tasks beyond neuroscience uniting complex multi-modal data structures, adaptive modeling, and real-time computational requirements.

Loss of consciousness is related to hyper-correlated gamma-band activity in anesthetized macaques and sleeping humans

Loss of consciousness can result from a wide range of causes, including natural sleep and pharmacologically induced anesthesia. Important insights might thus come from identifying neuronal mechanisms of loss and re-emergence of consciousness independent of a specific manipulation. Therefore, to seek neuronal signatures of loss of consciousness common to sleep and anesthesia we analyzed spontaneous electrophysiological activity recorded in two experiments. First, electrocorticography (ECoG) acquired from 4 macaque monkeys anesthetized with different anesthetic agents (ketamine, medetomidine, propofol) and, second, stereo-electroencephalography (sEEG) from 10 epilepsy patients in different wake-sleep stages (wakefulness, NREM, REM). Specifically, we investigated co-activation patterns among brain areas, defined as correlations between local amplitudes of gamma-band activity. We found that resting wakefulness was associated with intermediate levels of gamma-band coupling, indicating neither complete dependence, nor full independence among brain regions. In contrast, loss of consciousness during NREM sleep and propofol anesthesia was associated with excessively correlated brain activity, as indicated by a robust increase of number and strength of positive correlations. However, such excessively correlated brain signals were not observed during REM sleep, and were present only to a limited extent during ketamine anesthesia. This might be related to the fact that, despite suppression of behavioral responsiveness, REM sleep and ketamine anesthesia often involve presence of dream-like conscious experiences. We conclude that hyper-correlated gamma-band activity might be a signature of loss of consciousness common across various manipulations and independent of behavioral responsiveness.

A Chronically Implantable Bidirectional Neural Interface for Non-human Primates

Optogenetics has potential applications in the study of epilepsy and neuroprostheses, and for studies on neural circuit dynamics. However, to achieve translation to clinical usage, optogenetic interfaces that are capable of chronic stimulation and monitoring with minimal brain trauma are required. We aimed to develop a chronically implantable device for photostimulation of the brain of non-human primates. We used a micro-light-emitting diode (LED) array with a flexible polyimide film. The array was combined with a whole-cortex electrocorticographic (ECoG) electrode array for simultaneous photostimulation and recording. Channelrhodopsin-2 (ChR2) was virally transduced into the cerebral cortex of common marmosets, and then the device was epidurally implanted into their brains. We recorded the neural activity during photostimulation of the awake monkeys for 4 months. The neural responses gradually increased after the virus injection for ~8 weeks and remained constant for another 8 weeks. The micro-LED and ECoG arrays allowed semi-invasive simultaneous stimulation and recording during long-term implantation in the brains of non-human primates. The development of this device represents substantial progress in the field of optogenetic applications.

Broadband Electrophysiological Dynamics Contribute to Global Resting-State fMRI Signal

Spontaneous activity observed with resting-state fMRI is used widely to uncover the brain's intrinsic functional networks in health and disease. Although many networks appear modular and specific, global and nonspecific fMRI fluctuations also exist and both pose a challenge and present an opportunity for characterizing and understanding brain networks. Here, we used a multimodal approach to investigate the neural correlates to the global fMRI signal in the resting state. Like fMRI, resting-state power fluctuations of broadband and arrhythmic, or scale-free, macaque electrocorticography and human magnetoencephalography activity were correlated globally. The power fluctuations of scale-free human electroencephalography (EEG) were coupled with the global component of simultaneously acquired resting-state fMRI, with the global hemodynamic change lagging the broadband spectral change of EEG by ∼5 s. The levels of global and nonspecific fluctuation and synchronization in scale-free population activity also varied across and depended on arousal states. Together, these results suggest that the neural origin of global resting-state fMRI activity is the broadband power fluctuation in scale-free population activity observable with macroscopic electrical or magnetic recordings. Moreover, the global fluctuation in neurophysiological and hemodynamic activity is likely modulated through diffuse neuromodulation pathways that govern arousal states and vigilance levels.

SIGNIFICANCE STATEMENT

This study provides new insights into the neural origin of resting-state fMRI. Results demonstrate that the broadband power fluctuation of scale-free electrophysiology is globally synchronized and directly coupled with the global component of spontaneous fMRI signals, in contrast to modularly synchronized fluctuations in oscillatory neural activity. These findings lead to a new hypothesis that scale-free and oscillatory neural processes account for global and modular patterns of functional connectivity observed with resting-state fMRI, respectively.

Structure Shapes Dynamics and Directionality in Diverse Brain Networks: Mathematical Principles and Empirical Confirmation in Three Species

Identifying how spatially distributed information becomes integrated in the brain is essential to understanding higher cognitive functions. Previous computational and empirical studies suggest a significant influence of brain network structure on brain network function. However, there have been few analytical approaches to explain the role of network structure in shaping regional activities and directionality patterns. In this study, analytical methods are applied to a coupled oscillator model implemented in inhomogeneous networks. We first derive a mathematical principle that explains the emergence of directionality from the underlying brain network structure. We then apply the analytical methods to the anatomical brain networks of human, macaque, and mouse, successfully predicting simulation and empirical electroencephalographic data. The results demonstrate that the global directionality patterns in resting state brain networks can be predicted solely by their unique network structures. This study forms a foundation for a more comprehensive understanding of how neural information is directed and integrated in complex brain networks.

Using pupil size and heart rate to infer affective states during behavioral neurophysiology and neuropsychology experiments

BACKGROUND

Nonhuman primates (NHPs) are a valuable research model because of their behavioral, physiological and neuroanatomical similarities to humans. In the absence of language, autonomic activity can provide crucial information about cognitive and affective states during single-unit recording, inactivation and lesion studies. Methods standardized for use in humans are not easily adapted to NHPs and detailed guidance has been lacking.

NEW METHOD

We provide guidance for monitoring heart rate and pupil size in the behavioral neurophysiology setting by addressing the methodological issues, pitfalls and solutions for NHP studies. The methods are based on comparative physiology to establish a rationale for each solution. We include examples from both electrophysiological and lesion studies.

RESULTS

Single-unit recording, pupil responses and heart rate changes represent a range of decreasing temporal resolution, a characteristic that impacts experimental design and analysis. We demonstrate the unexpected result that autonomic measures acquired before and after amygdala lesions are comparable despite disruption of normal autonomic function.

COMPARISON WITH EXISTING METHODS

Species and study design differences can render standard techniques used in human studies inappropriate for NHP studies. We show how to manage data from small groups typical of NHP studies, data from the short behavioral trials typical of neurophysiological studies, issues associated with longitudinal studies, and differences in anatomy and physiology.

CONCLUSIONS

Autonomic measurement to infer cognitive and affective states in NHP is neither off-the-shelf nor onerous. Familiarity with the issues and solutions will broaden the use of autonomic signals in NHP single unit and lesion studies.

A new method for quantifying the performance of EEG blind source separation algorithms by referencing a simultaneously recorded ECoG signal

Blind source separation (BSS) algorithms extract neural signals from electroencephalography (EEG) data. However, it is difficult to quantify source separation performance because there is no criterion to dissociate neural signals and noise in EEG signals. This study develops a method for evaluating BSS performance. The idea is neural signals in EEG can be estimated by comparison with simultaneously measured electrocorticography (ECoG). Because the ECoG electrodes cover the majority of the lateral cortical surface and should capture most of the original neural sources in the EEG signals. We measured real EEG and ECoG data and developed an algorithm for evaluating BSS performance. First, EEG signals are separated into EEG components using the BSS algorithm. Second, the EEG components are ranked using the correlation coefficients of the ECoG regression and the components are grouped into subsets based on their ranks. Third, canonical correlation analysis estimates how much information is shared between the subsets of the EEG components and the ECoG signals. We used our algorithm to compare the performance of BSS algorithms (PCA, AMUSE, SOBI, JADE, fastICA) via the EEG and ECoG data of anesthetized nonhuman primates. The results (Best case >JADE = fastICA >AMUSE = SOBI ≥ PCA >random separation) were common to the two subjects. To encourage the further development of better BSS algorithms, our EEG and ECoG data are available on our Web site (http://neurotycho.org/) as a common testing platform.

Penalized Multi-Way Partial Least Squares for Smooth Trajectory Decoding from Electrocorticographic (ECoG) Recording

In the current paper the decoding algorithms for motor-related BCI systems for continuous upper limb trajectory prediction are considered. Two methods for the smooth prediction, namely Sobolev and Polynomial Penalized Multi-Way Partial Least Squares (PLS) regressions, are proposed. The methods are compared to the Multi-Way Partial Least Squares and Kalman Filter approaches. The comparison demonstrated that the proposed methods combined the prediction accuracy of the algorithms of the PLS family and trajectory smoothness of the Kalman Filter. In addition, the prediction delay is significantly lower for the proposed algorithms than for the Kalman Filter approach. The proposed methods could be applied in a wide range of applications beyond neuroscience.

Social Suppressive Behavior Is Organized by the Spatiotemporal Integration of Multiple Cortical Regions in the Japanese Macaque

Under social conflict, monkeys develop hierarchical positions through social interactions. Once the hierarchy is established, the dominant monkey dominates the space around itself and the submissive monkey tries not to violate this space. Previous studies have shown the contributions of the frontal and parietal cortices in social suppression, but the contributions of other cortical areas to suppressive functions remain elusive. We recorded neural activity in large cortical areas using electrocorticographic (ECoG) arrays while monkeys performed a social food-grab task in which a target monkey was paired with either a dominant or a submissive monkey. If the paired monkey was dominant, the target monkey avoided taking food in the shared conflict space, but not in other areas. By contrast, when the paired monkey was submissive, the target monkey took the food freely without hesitation. We applied decoding analysis to the ECoG data to see when and which cortical areas contribute to social behavioral suppression. Neural information discriminating the social condition was more evident when the conflict space was set in the area contralateral to the recording hemisphere. We found that the information increased as the social pressure increased during the task. Before food presentation, when the pressure was relatively low, the parietal and somatosensory-motor cortices showed sustained discrimination of the social condition. After food presentation, when the monkey faced greater pressure to make a decision as to whether it should take the food, the prefrontal and visual cortices started to develop buildup responses. The social representation was found in a sustained form in the parietal and somatosensory-motor regions, followed by additional buildup form in the visual and prefrontal cortices. The representation was less influenced by reward expectation. These findings suggest that social adaptation is achieved by a higher-order self-regulation process (incorporating motor preparation/execution processes) in accordance with the embodied social contexts.

Facilitative effect of repetitive presentation of one stimulus on cortical responses to other stimuli in macaque monkeys--a possible neural mechanism for mismatch negativity

The event-related potential 'mismatch negativity' (MMN) is an indicator of a perceiver's ability to detect deviations in sensory signal streams. MMN and its homologue in animals, mismatch activity (MMA), are differential neural responses to a repeatedly presented stimulus and a subsequent deviant stimulus (oddball). Because neural mechanisms underlying MMN and MMA remain unclear, there is a controversy as to whether MMN and MMA arise solely from stimulus-specific adaptation (SSA), in which the response to a stimulus cumulatively attenuates with its repetitive presentation. To address this issue, we used electrocorticography and the auditory roving-oddball paradigm in two awake macaque monkeys. We examined the effect of stimulus repetition number on MMA and on responses to repeated stimuli and oddballs across the cerebral cortex in the time-frequency domain. As the repetition number increased, MMA spread across the temporal, frontal and parietal cortices, and each electrode yielded a larger MMA. Surprisingly, this increment in MMA largely depended on response augmentation to the oddball rather than on SSA to the repeated stimulus. Following sufficient repetition, the oddball evoked a spectral power increment in some electrodes on the frontal cortex that had shown no power increase to the stimuli with less or no preceding repetition. We thereby revealed that repetitive presentation of one stimulus not only leads to SSA but also facilitates the cortical response to oddballs involving a wide range of cortical regions. This facilitative effect might underlie the generation of MMN-like scalp potentials in macaques that potentially shares similar neural mechanisms with MMN in humans.

Untangling Brain-Wide Dynamics in Consciousness by Cross-Embedding

Brain-wide interactions generating complex neural dynamics are considered crucial for emergent cognitive functions. However, the irreducible nature of nonlinear and high-dimensional dynamical interactions challenges conventional reductionist approaches. We introduce a model-free method, based on embedding theorems in nonlinear state-space reconstruction, that permits a simultaneous characterization of complexity in local dynamics, directed interactions between brain areas, and how the complexity is produced by the interactions. We demonstrate this method in large-scale electrophysiological recordings from awake and anesthetized monkeys. The cross-embedding method captures structured interaction underlying cortex-wide dynamics that may be missed by conventional correlation-based analysis, demonstrating a critical role of time-series analysis in characterizing brain state. The method reveals a consciousness-related hierarchy of cortical areas, where dynamical complexity increases along with cross-area information flow. These findings demonstrate the advantages of the cross-embedding method in deciphering large-scale and heterogeneous neuronal systems, suggesting a crucial contribution by sensory-frontoparietal interactions to the emergence of complex brain dynamics during consciousness.

Mismatch negativity in common marmosets: Whole-cortical recordings with multi-channel electrocorticograms

Mismatch negativity (MMN) is a component of event-related potentials (ERPs) evoked by violations of regularity in sensory stimulus-series in humans. Recently, the MMN has received attention as a clinical and translatable biomarker of psychiatric disorders such as schizophrenia, and for the development animal models of these psychiatric disorders. In this study, we investigated the generation of MMN in common marmosets, which are an important non-human primate model with genetic manipulability. We recorded the electrocorticograms (ECoGs) from two common marmosets with epidurally implanted electrodes covering a wide range of cortical regions. ECoG recordings were conducted in a passive listening condition with a roving oddball paradigm. We compared the ERPs evoked by repeatedly presented standard stimuli and those evoked by the deviant stimuli. Significant differences in the ERPs were observed in several cortical areas. In particular, deviant stimuli elicited larger negative activity than standard stimuli in the temporal area. In addition, the latency and polarity of the activity were comparable to human MMNs. This is thus the first report of MMN-like activity in common marmosets. Our findings have the potential to advance future gene-manipulation studies that aim to establish non-human primate models of schizophrenia.

Cortical network architecture for context processing in primate brain

Context is information linked to a situation that can guide behavior. In the brain, context is encoded by sensory processing and can later be retrieved from memory. How context is communicated within the cortical network in sensory and mnemonic forms is unknown due to the lack of methods for high-resolution, brain-wide neuronal recording and analysis. Here, we report the comprehensive architecture of a cortical network for context processing. Using hemisphere-wide, high-density electrocorticography, we measured large-scale neuronal activity from monkeys observing videos of agents interacting in situations with different contexts. We extracted five context-related network structures including a bottom-up network during encoding and, seconds later, cue-dependent retrieval of the same network with the opposite top-down connectivity. These findings show that context is represented in the cortical network as distributed communication structures with dynamic information flows. This study provides a general methodology for recording and analyzing cortical network neuronal communication during cognition.

DOI:http://dx.doi.org/10.7554/eLife.06121.001

If we see someone looking frightened, the way we respond to the situation is influenced by other information, referred to as the ‘context’. For example, if the person is frightened because another individual is shouting at them, we might try to intervene. However, if the person is watching a horror video we may decide that they don't need our help and leave them to it. Nevertheless, it is not clear how the brain processes the context of a situation to inform our response.

Here, Chao et al. developed a new method to study electrical activity across the whole of the brain and used it to study how monkeys process context in response to several different social situations. In the experiments, monkeys were shown video clips in which one monkey—known as the ‘video monkey’—was threatened by a human or another monkey, or in which the video monkey is facing an empty wall (i.e., in three different contexts). Afterwards, the video monkey either displays a frightened expression or a neutral one. Chao et al. found that if the video monkey looked frightened by the context, the monkeys watching the video clip shifted their gaze to observe the apparent threat. How these monkeys shifted their gaze depended on the context, but this behavior was absent when the video monkey gave a neutral expression.

The experiments used an array of electrodes that covered a wide area of the monkeys' brains to record electrical activity of nerve cells as the monkeys watched the videos. Chao et al. investigated how brain regions communicated with each other in response to different contexts, and found that the information of contexts was presented in the interactions between distant brain regions. The monkeys' brains sent information from a region called the temporal cortex (which is involved in processing sensory and social information), to another region called the prefrontal cortex (which is involved in functions such as reasoning, attention, and memory). Seconds later, the flow of information was reversed as the monkeys utilized information about the context to guide their behavior.

Chao et al.'s findings reveal how information about the context of a situation is transmitted around the brain to inform a response. The next challenge is to experimentally manipulate the identified brain circuits to investigate if problems in context processing could lead to the inappropriate responses that contribute to schizophrenia, post-traumatic stress disorder and other psychiatric disorders.

DOI:http://dx.doi.org/10.7554/eLife.06121.002

Frequency-dependent spatiotemporal profiles of visual responses recorded with subdural ECoG electrodes in awake monkeys: Differences between high- and low-frequency activity

Electrocorticography (ECoG) constitutes a powerful and promising neural recording modality in humans and animals. ECoG signals are often decomposed into several frequency bands, among which the so-called high-gamma band (80-250Hz) has been proposed to reflect local cortical functions near the cortical surface below the ECoG electrodes. It is typically assumed that the lower the frequency bands, the lower the spatial resolution of the signals; thus, there is not much to gain by analyzing the event-related changes of the ECoG signals in the lower-frequency bands. However, differences across frequency bands have not been systematically investigated. To address this issue, we recorded ECoG activity from two awake monkeys performing a retinotopic mapping task. We characterized the spatiotemporal profiles of the visual responses in the time-frequency domain. We defined the preferred spatial position, receptive field (RF), and response latencies of band-limited power (BLP) (i.e., alpha [3.9-11.7Hz], beta [15.6-23.4Hz], low [30-80Hz] and high [80-250Hz] gamma) for each electrode and compared them across bands and time-domain visual evoked potentials (VEPs). At the population level, we found that the spatial preferences were comparable across bands and VEPs. The high-gamma power showed a smaller RF than the other bands and VEPs. The response latencies for the alpha band were always longer than the latencies for the other bands and fastest in VEPs. Comparing the response profiles in both space and time for each cortical region (V1, V4+, and TEO/TE) revealed regional idiosyncrasies. Although the latencies of visual responses in the beta, low-, and high-gamma bands were almost identical in V1 and V4+, beta and low-gamma BLP occurred about 17ms earlier than high-gamma power in TEO/TE. Furthermore, TEO/TE exhibited a unique pattern in the spatial response profile: the alpha and high-gamma responses tended to prefer the foveal regions, whereas the beta and low-gamma responses preferred the peripheral visual fields with larger RFs. This suggests that neurons in TEO/TE first receive less selective spatial information via beta and low-gamma BLP but later receive more fine-tuned spatial foveal information via high-gamma power. This result is consistent with a hypothesis previously proposed by Nakamura et al. (1993) that states that visual processing in TEO/TE starts with coarse-grained information, which primes subsequent fine-grained information. Collectively, our results demonstrate that ECoG can be a potent tool for investigating the nature of the neural computations in each cortical region that cannot be fully understood by measuring only the spiking activity, through the incorporation of the knowledge of the spatiotemporal characteristics across all frequency bands.

Estimating Directed Connectivity from Cortical Recordings and Reconstructed Sources

In cognitive neuroscience, electrical brain activity is most commonly recorded at the scalp. In order to infer the contributions and connectivity of underlying neuronal sources within the brain, it is necessary to reconstruct sensor data at the source level. Several approaches to this reconstruction have been developed, thereby solving the so-called implicit inverse problem Michel et al. (Clin Neurophysiol 115:2195–2222, 2004). However, a unifying premise against which to validate these source reconstructions is seldom available. The dataset provided in this work, in which brain activity is simultaneously recorded on the scalp (non-invasively) by electroencephalography (EEG) and on the cortex (invasively) by electrocorticography (ECoG), can be of a great help in this direction. These multimodal recordings were obtained from a macaque monkey under wakefulness and sedation. Our primary goal was to establish the connectivity architecture between two sources of interest (frontal and parietal), and to assess how their coupling changes over the conditions. We chose these sources because previous studies have shown that the connections between them are modified by anaesthesia Boly et al. (J Neurosci 32:7082–7090, 2012). Our secondary goal was to evaluate the consistency of the connectivity results when analyzing sources recorded from invasive data (128 implanted ECoG sources) and source activity reconstructed from scalp recordings (19 EEG sensors) at the same locations as the ECoG sources. We conclude that the directed connectivity in the frequency domain between cortical sources reconstructed from scalp EEG is qualitatively similar to the connectivity inferred directly from cortical recordings, using both data-driven (directed transfer function) and biologically grounded (dynamic causal modelling) methods. Furthermore, the connectivity changes identified were consistent with previous findings Boly et al. (J Neurosci 32:7082–7090, 2012). Our findings suggest that inferences about directed connectivity based upon non-invasive electrophysiological data have construct validity in relation to invasive recordings.