Directional patterns of cross frequency phase and amplitude coupling within the resting state mimic patterns of fMRI functional connectivity

Temporospatial Nestedness in Consciousness: An Updated Perspective on the Temporospatial Theory of Consciousness

Time and space are fundamental elements that permeate the fabric of nature, and their significance in relation to neural activity and consciousness remains a compelling yet unexplored area of research. The Temporospatial Theory of Consciousness (TTC) provides a framework that links time, space, neural activity, and consciousness, shedding light on the intricate relationships among these dimensions. In this review, I revisit the fundamental concepts and mechanisms proposed by the TTC, with a particular focus on the central concept of temporospatial nestedness. I propose an extension of temporospatial nestedness by incorporating the nested relationship between the temporal circuit and functional geometry of the brain. To further unravel the complexities of temporospatial nestedness, future research directions should emphasize the characterization of functional geometry and the temporal circuit across multiple spatial and temporal scales. Investigating the links between these scales will yield a more comprehensive understanding of how spatial organization and temporal dynamics contribute to conscious states. This integrative approach holds the potential to uncover novel insights into the neural basis of consciousness and reshape our understanding of the world-brain dynamic.

Systematic evaluation of head motion on resting-state functional connectivity MRI in the neonate

Reliability and robustness of resting state functional connectivity MRI (rs-fcMRI) relies, in part, on minimizing the influence of head motion on measured brain signals. The confounding effects of head motion on functional connectivity have been extensively studied in adults, but its impact on newborn brain connectivity remains unexplored. Here, using a large newborn data set consisting of 159 rs-fcMRI scans acquired in the Developing Brain Institute at Children's National Hospital and 416 scans from The Developing Human Connectome Project (dHCP), we systematically investigated associations between head motion and rs-fcMRI. Head motion during the scan significantly affected connectivity at sensory-related networks and default mode networks, and at the whole brain scale; the direction of motion effects varied across the whole brain. Comparing high- versus low-head motion groups suggested that head motion can impact connectivity estimates across the whole brain. Censoring of high-motion volumes using frame-wise displacement significantly reduced the confounding effects of head motion on neonatal rs-fcMRI. Lastly, in the dHCP data set, we demonstrated similar persistent associations between head motion and network connectivity despite implementing a standard denoising strategy. Collectively, our results highlight the importance of using rigorous head motion correction in preprocessing neonatal rs-fcMRI to yield reliable estimates of brain activity.

Early adolescent psychological distress and cognition, correlates of resting-state EEG, interregional phase-amplitude coupling

Delineating neurobiological markers of youth mental health is crucial for early identification and treatment. One promising marker is phase-amplitude coupling (PAC), cross-frequency coupling between the phase of slower oscillatory activity and the amplitude of faster oscillatory activity in the brain. Prior research has demonstrated that PAC is associated with both cognition and mental health and can be modulated using neurostimulation. However, to date research investigating PAC has focused primarily on adults, and only within-region theta-gamma coupling in the context of mental health. We investigated associations between interregional resting-state PAC (posterior-anterior cortex), and cognition and psychological distress in N = 77 (M_age = 12.58 years, SD = 0.31; 51 % female) 12-year-olds. Firstly, while left theta-beta PAC showed a moderate positive correlation (r = 0.529, p < .01), right theta-gamma PAC showed a weak positive correlation, with psychological distress (r = 0.283, p < .05). In terms of cognition, moderate correlations were observed between: (i) increased left theta-beta PAC and increased psychomotor speed (r = -0.367, p < .05); (ii) increased left alpha-beta PAC and decreased attention (r = 0.355, p ≤0.01); and (iii) increased left alpha-beta PAC and decreased verbal learning and memory (r = -0.352, p < .01). Whereas weak associations were observed for: (i) increased left alpha-beta PAC and decreased executive functioning scores (r = 0.284, p < .05); and (ii) increased left alpha-gamma PAC and increased attention (r = -0.272, p < .05). The overall findings of this exploratory study are encouraging, although all the correlations were in the weak-to-moderate range and require replication. Further research may confirm interregional resting-state PAC as a biomarker that can help us better understand the link between mental health and cognition in adolescents and improve treatment of cognitive related deficits in mental illness.

Infraslow closed-loop brain training for anxiety and depression (ISAD): a protocol for a randomized, double-blind, sham-controlled pilot trial in adult females with internalizing disorders

BACKGROUND

The core intrinsic connectivity networks (core-ICNs), encompassing the default-mode network (DMN), salience network (SN) and central executive network (CEN), have been shown to be dysfunctional in individuals with internalizing disorders (IDs, e.g. major depressive disorder, MDD; generalized anxiety disorder, GAD; social anxiety disorder, SOC). As such, source-localized, closed-loop brain training of electrophysiological signals, also known as standardized low-resolution electromagnetic tomography (sLORETA) neurofeedback (NFB), targeting key cortical nodes within these networks has the potential to reduce symptoms associated with IDs and restore normal core ICN function. We intend to conduct a randomized, double-blind (participant and assessor), sham-controlled, parallel-group (3-arm) trial of sLORETA infraslow (<0.1 Hz) fluctuation neurofeedback (sLORETA ISF-NFB) 3 times per week over 4 weeks in participants (n=60) with IDs. Our primary objectives will be to examine patient-reported outcomes (PROs) and neurophysiological measures to (1) compare the potential effects of sham ISF-NFB to either genuine 1-region ISF-NFB or genuine 2-region ISF-NFB, and (2) assess for potential associations between changes in PRO scores and modifications of electroencephalographic (EEG) activity/connectivity within/between the trained regions of interest (ROIs). As part of an exploratory analysis, we will investigate the effects of additional training sessions and the potential for the potentiation of the effects over time.

METHODS

We will randomly assign participants who meet the criteria for MDD, GAD, and/or SOC per the MINI (Mini International Neuropsychiatric Interview for DSM-5) to one of three groups: (1) 12 sessions of posterior cingulate cortex (PCC) ISF-NFB up-training (n=15), (2) 12 sessions of concurrent PCC ISF up-training and dorsal anterior cingulate cortex (dACC) ISF-NFB down-training (n=15), or (3) 6 sessions of yoked-sham training followed by 6 sessions genuine ISF-NFB (n=30). Transdiagnostic PROs (Hospital Anxiety and Depression Scale, HADS; Inventory of Depression and Anxiety Symptoms - Second Version, IDAS-II; Multidimensional Emotional Disorder Inventory, MEDI; Intolerance of Uncertainty Scale - Short Form, IUS-12; Repetitive Thinking Questionnaire, RTQ-10) as well as resting-state neurophysiological measures (full-band EEG and ECG) will be collected from all subjects during two baseline sessions (approximately 1 week apart) then at post 6 sessions, post 12 sessions, and follow-up (1 month later). We will employ Bayesian methods in R and advanced source-localisation software (i.e. exact low-resolution brain electromagnetic tomography; eLORETA) in our analysis.

DISCUSSION

This protocol will outline the rationale and research methodology for a clinical pilot trial of sLORETA ISF-NFB targeting key nodes within the core-ICNs in a female ID population with the primary aims being to assess its potential efficacy via transdiagnostic PROs and relevant neurophysiological measures.

TRIAL REGISTRATION

Our study was prospectively registered with the Australia New Zealand Clinical Trials Registry (ANZCTR; Trial ID: ACTRN12619001428156). Registered on October 15, 2019.

Personalizing repetitive transcranial magnetic stimulation for precision depression treatment based on functional brain network controllability and optimal control analysis

Brain neuromodulation effectively treats neurological diseases and psychiatric disorders such as Depression. However, due to patient heterogeneity, neuromodulation treatment outcomes are often highly variable, requiring patient-specific stimulation protocols throughout the recovery stages to optimize treatment outcomes. Therefore, it is critical to personalize neuromodulation protocol to optimize the patient-specific stimulation targets and parameters by accommodating inherent interpatient variability and intersession alteration during treatments. The study aims to develop a personalized repetitive transcranial magnetic stimulation (rTMS) protocol and evaluate its feasibility in optimizing the treatment efficiency using an existing dataset from an antidepressant experimental imaging study in depression. The personalization of the rTMS treatment protocol was achieved by personalizing both stimulation targets and parameters via a novel approach integrating the functional brain network controllability analysis and optimal control analysis. First, the functional brain network controllability analysis was performed to identify the optimal rTMS stimulation target from the effective connectivity network constructed from patient-specific resting-state functional magnetic resonance imaging data. The optimal control algorithm was then applied to optimize the rTMS stimulation parameters based on the optimized target. The performance of the proposed personalized rTMS technique was evaluated using datasets collected from a longitudinal antidepressant experimental imaging study in depression (n = 20). Simulation models demonstrated that the proposed personalized rTMS protocol outperformed the standard rTMS treatment by efficiently steering a depressive resting brain state to a healthy resting brain state, indicated by the significantly less control energy needed and higher model fitting accuracy achieved. The node with the maximum average controllability of each patient was designated as the optimal target region for the personalized rTMS protocol. Our results also demonstrated the theoretical feasibility of achieving comparable neuromodulation efficacy by stimulating a single node compared to stimulating multiple driver nodes. The findings support the feasibility of developing personalized neuromodulation protocols to more efficiently treat depression and other neurological diseases.

Disrupted Spatiotemporal Complexity of Resting-State Electroencephalogram Dynamics Is Associated With Adaptive and Maladaptive Rumination in Major Depressive Disorder

Patients with major depressive disorder (MDD) exhibit abnormal rumination, including both adaptive and maladaptive forms. However, the neural substrates of rumination in depression remain poorly understood. We hypothesize that divergent spatiotemporal complexity of brain oscillations would be associated with the levels of rumination in MDD. We employed the multi-scale entropy (MSE), power and phase-amplitude coupling (PAC) to estimate the complexity of rhythmic dynamics from the eye-closed high-density electroencephalographic (EEG) data in treatment-naive patients with MDD (n = 24) and healthy controls (n = 22). The depressive, brooding, and reflective subscales of the Ruminative Response Scale were assessed. MDD patients showed higher MSE in timescales finer than 5 (cluster P = 0.038) and gamma power (cluster P = 0.034), as well as lower PAC values between alpha/low beta and gamma bands (cluster P = 0.002- 0.021). Higher reflective rumination in MDD was region-specifically associated with the more localized EEG dynamics, including the greater MSE in scales finer than 8 (cluster P = 0.008), power in gamma (cluster P = 0.018) and PAC in low beta-gamma (cluster P = 0.042), as well as weaker alpha-gamma PAC (cluster P = 0.016- 0.029). Besides, the depressive and brooding rumination in MDD showed the lack of correlations with global long-range EEG variables. Our findings support the disturbed neural communications and point to the spatial reorganization of brain networks in a timescale-dependent migration toward local during adaptive and maladaptive rumination in MDD. These findings may provide potential implications on probing and modulating dynamic neuronal fluctuations during the rumination in depression.

Electrocorticography to Investigate Age-Related Brain Lateralization on Pediatric Motor Inhibition

Response inhibition refers to the ability to suppress inappropriate actions that interfere with goal-driven behavior. The inferior frontal gyrus (IFG) is known to be associated with inhibition of a motor response by assuming executive control over motor cortex outputs. This study aimed to evaluate the pediatric development of response inhibition through subdural electrocorticography (ECoG) recording. Subdural ECoG recorded neural activities simultaneously during a Go/No-Go task, which was optimized for children. Different frequency power [theta: 4–8 Hz; beta: 12–40 Hz; high-gamma (HG): 70–200 Hz] was estimated within the IFG and motor cortex. Age-related analysis was computed by each bandpass power ratio between Go and No-Go conditions, and phase-amplitude coupling (PAC) over IFG by using the modulating index metric in two conditions. For all the eight pediatric patients, HG power was more activated in No-Go trials than in Go trials, in either right- or left-side IFG when available. In the IFG region, the power over theta and HG in No-Go conditions was higher than those in Go conditions, with significance over the right side (p < 0.05). The age-related lateralization from both sides to the right side was observed from the ratio of HG power and PAC value between the No-Go and Go trials. In the pediatric population, the role of motor inhibition was observed in both IFG, with age-related lateralization to the right side, which was proved in the previous functional magnetic resonance imaging studies. In this study, the evidence correlation of age and response inhibition was observed directly by the evidence of cortical recordings.

Connectomics of human electrophysiology

We present both a scientific overview and conceptual positions concerning the challenges and assets of electrophysiological measurements in the search for the nature and functions of the human connectome. We discuss how the field has been inspired by findings and approaches from functional magnetic resonance imaging (fMRI) and informed by a small number of significant multimodal empirical studies, which show that the canonical networks that are commonplace in fMRI are in fact rooted in electrophysiological processes. This review is also an opportunity to produce a brief, up-to-date critical survey of current data modalities and analytical methods available for deriving both static and dynamic connectomes from electrophysiology. We review hurdles that challenge the significance and impact of current electrophysiology connectome research. We then encourage the field to take a leap of faith and embrace the wealth of electrophysiological signals, despite their apparent, disconcerting complexity. Our position is that electrophysiology connectomics is poised to inform testable mechanistic models of information integration in hierarchical brain networks, constructed from observable oscillatory and aperiodic signal components and their polyrhythmic interactions.

Cross-Frequency Coupling in Childhood Absence Epilepsy

Objective: Absence seizures are the prototypic primarily generalized seizures, but there is incomplete understanding regarding their generation and maintenance. A core network for absence seizures has been defined, including focal cortical and thalamic regions that have frequency-dependent interactions. The purpose of this study was to investigate within-frequency coupling and cross-frequency coupling (CFC) during human absence seizures, to identify key regions (hubs) within the absence network that contribute to propagation and maintenance. Methods: Thirteen children with new-onset and untreated childhood absence epilepsy had over 60 typical absence seizures during both electroencephalography-functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) recordings. The spatial map of the ictal network was defined using fMRI and used as prior information for MEG connectivity. A multilayer network approach was used to investigate within-frequency coupling and CFC for canonical frequency bands. A rigorous null-modeling approach was used to determine connections outside the noise floor. Results: Strong coupling between beta and gamma frequencies, within the left frontal cortex, and between the left frontal and right parietal regions was observed. There was also strong connectivity between left frontal and right parietal nodes within the gamma band. Multilayer versatility analysis identified a cluster of network hubs in the left frontal region. Interpretation: Cortical regions commonly identified as being critical for absence seizure generation (frontal cortex, precuneus) have strong CFC and within-frequency coupling between beta and gamma bands. As nonpharmacologic treatments, such as neuromodulation, become available for generalized epilepsies, detailed mechanistic understanding of how "diffuse" seizures are generated and maintained will be necessary to provide optimal outcomes.

A novel beamformer-based imaging of phase-amplitude coupling (BIPAC) unveiling the inter-regional connectivity of emotional prosody processing in women with primary dysmenorrhea

Objective. Neural communication or the interactions of brain regions play a key role in the formation of functional neural networks. A type of neural communication can be measured in the form of phase-amplitude coupling (PAC), which is the coupling between the phase of low-frequency oscillations and the amplitude of high-frequency oscillations. This paper presents a beamformer-based imaging method, beamformer-based imaging of PAC (BIPAC), to quantify the strength of PAC between a seed region and other brain regions.Approach. A dipole is used to model the ensemble of neural activity within a group of nearby neurons and represents a mixture of multiple source components of cortical activity. From ensemble activity at each brain location, the source component with the strongest coupling to the seed activity is extracted, while unrelated components are suppressed to enhance the sensitivity of coupled-source estimation.Main results. In evaluations using simulation data sets, BIPAC proved advantageous with regard to estimation accuracy in source localization, orientation, and coupling strength. BIPAC was also applied to the analysis of magnetoencephalographic signals recorded from women with primary dysmenorrhea in an implicit emotional prosody experiment. In response to negative emotional prosody, auditory areas revealed strong PAC with the ventral auditory stream and occipitoparietal areas in the theta-gamma and alpha-gamma bands, which may respectively indicate the recruitment of auditory sensory memory and attention reorientation. Moreover, patients with more severe pain experience appeared to have stronger coupling between auditory areas and temporoparietal regions.Significance. Our findings indicate that the implicit processing of emotional prosody is altered by menstrual pain experience. The proposed BIPAC is feasible and applicable to imaging inter-regional connectivity based on cross-frequency coupling estimates. The experimental results also demonstrate that BIPAC is capable of revealing autonomous brain processing and neurodynamics, which are more subtle than active and attended task-driven processing.

Concurrent Deep Brain Stimulation Reduces the Direct Cortical Stimulation Necessary for Motor Output

BACKGROUND

Converging literatures suggest that deep brain stimulation (DBS) in Parkinson's disease affects multiple circuit mechanisms. One proposed mechanism is the normalization of primary motor cortex (M1) pathophysiology via effects on the hyperdirect pathway.

OBJECTIVES

We hypothesized that DBS would reduce the current intensity necessary to modulate motor-evoked potentials from focally applied direct cortical stimulation (DCS).

METHODS

Intraoperative subthalamic DBS, DCS, and preoperative diffusion tensor imaging data were acquired in 8 patients with Parkinson's disease.

RESULTS

In 7 of 8 patients, DBS significantly reduced the M1 DCS current intensity required to elicit motor-evoked potentials. This neuromodulation was specific to select DBS bipolar configurations. In addition, the volume of activated tissue models of these configurations were significantly associated with overlap of the hyperdirect pathway.

CONCLUSIONS

DBS reduces the current necessary to elicit a motor-evoked potential using DCS. This supports a circuit mechanism of DBS effectiveness, potentially involving the hyperdirect pathway that speculatively may underlie reductions in hypokinetic abnormalities in Parkinson's disease. © 2020 International Parkinson and Movement Disorder Society.

Dissociation as a disorder of integration - On the footsteps of Pierre Janet

At the end of the 19th century Pierre Janet described dissociation as an altered state of consciousness manifested in disrupted integration of psychological functions. Clinically, such disruption comprises compartmentalization symptoms like amnesia, detachment symptoms like depersonalization/derealization, and structural dissociation of personality with changes in the sense of self. The exact neuronal mechanisms leading to these different symptoms remain unclear. We here suggest to put Janet's original account of dissociation as disrupted integration of psychological functions into a novel context, that is, a neuronal context as related to current brain imaging. This requires a combined theoretical and empirical approach on data supporting such neuronal reframing of Janet. For that, we here review (i) past and (ii) recent psychological and neuronal views on dissociation together with neuroscientific theories of integration, which (iii) are supported and complemented by preliminary fMRI data. We propose three neuronal mechanisms of dynamic integration operating at different levels of the brain's spontaneous activity - temporo-spatial binding on the regional level, temporo-spatial synchronization on the network level, and temporo-spatial globalization on the global level. These neuronal mechanisms, in turn, may be related to different symptomatic manifestation of dissociation operating at different levels, e.g., compartmentalization, detachment, and structural, which, as we suggest, can all be traced to disrupted integration of neuronal and psychological functions as originally envisioned by Janet.

Differentiation of epileptic regions from voluntary high-gamma activation via interictal cross-frequency windowed power-power correlation

OBJECTIVE

Electrocorticography is an indispensable tool in identifying the epileptogenic zone in the presurgical evaluation of many epilepsy patients. Traditional electrocorticographic features (spikes, ictal onset changes, and recently high-frequency oscillations [HFOs]) rely on the presence of transient features that occur within or near epileptogenic cortex. Here the authors report on a novel corticography feature of epileptogenic cortex-covariation of high-gamma and beta frequency band power profiles. Band-limited power was measured from each recording site based on native physiological signal differences without relying on clinical ictal or interictal epileptogenic features. In this preliminary analysis, frequency windowed power correlation appears to be a specific marker of the epileptogenic zone. The authors' overall aim was to validate this observation with the location of the eventual resection and outcome.

METHODS

The authors conducted a retrospective analysis of 13 adult patients who had undergone electrocorticography for surgical planning at their center. They quantified the correlation of high-gamma (70-200 Hz) and beta (12-18 Hz) band frequency power per electrode site during a cognitive task. They used a sliding window method to correlate the power of smoothed, Hilbert-transformed high-gamma and beta bands. They then compared positive and negative correlations between power in the high-gamma and beta bands in the setting of a hand versus a tongue motor task as well as within the resting state. Significant positive correlations were compared to surgically resected areas and outcomes based on reviewed records.

RESULTS

Positive high-gamma and beta correlations appeared to predict the area of eventual resection and, preliminarily, surgical outcome independent of spike detection. In general, patients with the best outcomes had well-localized positive correlations (high-gamma and beta activities) to areas of eventual resection, while those with poorer outcomes displayed more diffuse patterns.

CONCLUSIONS

Data in this study suggest that positive high-gamma and beta correlations independent of any behavioral metric may have clinical applicability in surgical decision-making. Further studies are needed to evaluate the clinical potential of this methodology. Additional work is also needed to relate these results to other methods, such as HFO detection or connectivity with other cortical areas.

Intrinsic connectome organization across temporal scales: New insights from cross-modal approaches

The discovery of a stable, whole-brain functional connectivity organization that is largely independent of external events has drastically extended our view of human brain function. However, this discovery has been primarily based on functional magnetic resonance imaging (fMRI). The role of this whole-brain organization in fast oscillation-based connectivity as measured, for example, by electroencephalography (EEG) and magnetoencephalography (MEG) is only beginning to emerge. Here, we review studies of intrinsic connectivity and its whole-brain organization in EEG, MEG, and intracranial electrophysiology with a particular focus on direct comparisons to connectome studies in fMRI. Synthesizing this literature, we conclude that irrespective of temporal scale over four orders of magnitude, intrinsic neurophysiological connectivity shows spatial similarity to the connectivity organization commonly observed in fMRI. A shared structural connectivity basis and cross-frequency coupling are possible mechanisms contributing to this similarity. Acknowledging that a stable whole-brain organization governs long-range coupling across all timescales of neural processing motivates researchers to take "baseline" intrinsic connectivity into account when investigating brain-behavior associations, and further encourages more widespread exploration of functional connectomics approaches beyond fMRI by using EEG and MEG modalities.

Intact Auditory Cortical Cross-Frequency Coupling in Early and Chronic Schizophrenia

BACKGROUND

Previous work has identified a hierarchical organization of neural oscillations that supports performance of complex cognitive and perceptual tasks, and can be indexed with phase-amplitude coupling (PAC) between low- and high-frequency oscillations. Our aim was to employ enhanced source localization afforded by magnetoencephalography (MEG) to expand on earlier reports of intact auditory cortical PAC in schizophrenia and to investigate how PAC may evolve over the early and chronic phases of the illness.

METHODS

Individuals with early schizophrenia (n=12) (≤5 years of illness duration), chronic schizophrenia (n=16) (>5 years of illness duration) and healthy comparators (n = 17) performed the auditory steady state response (ASSR) to 40, 30, and 20 Hz stimuli during MEG recordings. We estimated amplitude and PAC on the MEG ASSR source localized to the auditory cortices.

RESULTS

Gamma amplitude during 40-Hz ASSR exhibited a significant group by hemisphere interaction, with both patient groups showing reduced right hemisphere amplitude and no overall lateralization in contrast to the right hemisphere lateralization demonstrated in controls. We found significant PAC in the right auditory cortex during the 40-Hz entrainment condition relative to baseline, however, PAC did not differ significantly between groups.

CONCLUSIONS

In the current study, we demonstrated an apparent sparing of ASSR-related PAC across phases of the illness, in contrast with impaired cortical gamma oscillation amplitudes. The distinction between our PAC and evoked ASSR findings supports the notion of separate but interacting circuits for the generation and maintenance of sensory gamma oscillations. The apparent sparing of PAC in both early and chronic schizophrenia patients could imply that the neuropathology of schizophrenia differentially affects these mechanisms across different stages of the disease. Future studies should investigate the distinction between PAC during passive tasks and more cognitively demanding task such as working memory so that we can begin to understand the influence of schizophrenia neuropathology on the larger framework for modulating neurocomputational capacity.

Spontaneous Variation in Electrocorticographic Resting-State Connectivity

Prior studies using functional magnetic resonance imaging, electroencephalography, and magnetoencephalography have observed both structured patterns in resting-state functional connectivity and spontaneous longitudinal variation in connectivity patterns independent of a task. In this first study using electrocorticography (ECoG), we characterized spontaneous, intersession variation in resting-state functional connectivity not linked to a task. We evaluated pairwise connectivity between electrodes using three measures (phase locking value [PLV], amplitude correlation, and coherence) for six canonical frequency bands, capturing different characteristics of time-evolving signals. We grouped electrodes into 10 functional regions and used intraclass correlation (ICC) to estimate pairwise longitudinal stability. We found that stronger PLV (PLV ≥0.4) in theta through gamma bands and strong correlation in all bands (R²'s ≥0.6) are linked to substantial stability (ICC ≥0.6), but that stability does not imply strong phase locking or amplitude correlation. There was no notable link between strong coherence and high ICC. All within-region PLVs are markedly stable across frequencies. In addition, we highlight interaction patterns across several regions: parahippocampal/entorhinal cortex is characterized by stable, weak functional connectivity except self-connections. Dorsolateral prefrontal cortex connectivity is weak and unstable, except self-connections. Inferior parietal lobule has little stability despite narrow connectivity bounds. We confirm prior studies linking functional connectivity strength and intersession variability, extending into higher frequencies than other modalities, with greater spatial specificity than scalp electrophysiology. We suggest further studies quantitatively compare ECoG to other modalities and/or use these findings as a baseline to capture functional connectivity and dynamics linked to perturbations with a task or disease state.

Consistent Gradient of Performance and Decoding of Stimulus Type and Valence From Local and Network Activity

The individual differences approach focuses on the variation of behavioral and neural signatures across subjects. In this context, we searched for intracranial neural markers of performance in three individuals with distinct behavioral patterns (efficient, borderline, and inefficient) in a dual-valence task assessing facial and lexical emotion recognition. First, we performed a preliminary study to replicate well-established evoked responses in relevant brain regions. Then, we examined time series data and network connectivity, combined with multivariate pattern analyses and machine learning, to explore electrophysiological differences in resting-state versus task-related activity across subjects. Next, using the same methodological approach, we assessed the neural decoding of performance for different dimensions of the task. The classification of time series data mirrored the behavioral gradient across subjects for stimulus type but not for valence. However, network-based measures reflected the subjects' hierarchical profiles for both stimulus types and valence. Therefore, this measure serves as a sensitive marker for capturing distributed processes such as emotional valence discrimination, which relies on an extended set of regions. Network measures combined with classification methods may offer useful insights to study single subjects and understand inter-individual performance variability. Promisingly, this approach could eventually be extrapolated to other neuroscientific techniques.

A Blind Module Identification Approach for Predicting Effective Connectivity Within Brain Dynamical Subnetworks

Model-based network discovery measures, such as the brain effective connectivity, require fitting of generative process models to measurements obtained from key areas across the network. For distributed dynamic phenomena, such as generalized seizures and slow-wave sleep, studying effective connectivity from real-time recordings is significantly complicated since (i) outputs from only a subnetwork can be practically measured, and (ii) exogenous subnetwork inputs are unobservable. Model fitting, therefore, constitutes a challenging blind module identification or model inversion problem for finding both the parameters and the many unknown inputs of the subnetwork. We herein propose a novel estimation framework for identifying nonlinear dynamic subnetworks in the case of slowly-varying, otherwise unknown local inputs. Starting with approximate predictions obtained using Cubature Kalman filtering, residuals of local output predictions are utilized to improve upon local input estimates. The algorithm performance is tested on both simulated and clinical EEG of induced seizures under electroconvulsive therapy (ECT). For the simulated network, the algorithm significantly boosted the estimation accuracy for inputs and connections from noisy EEG. For the clinical data, the algorithm predicted increased subnetwork inputs during the pre-stimulus anesthesia condition. Importantly, it predicted an increased frontocentral connectivity during the generalized seizure that is commensurate with electrode placement and that corroborates the clinical hypothesis of increased frontal focality of therapeutic ECT seizures. The proposed framework can be extended to account for several input configurations and can in principle be applied to study effective connectivity within brain subnetworks defined at the microscale (cortical lamina interaction) or at the macroscale (sensory integration).

An interspecies comparative study of invasive electrophysiological functional connectivity

INTRODUCTION

Resting-state connectivity patterns have been observed in humans and other mammal species, and can be recorded using a variety of different technologies. Functional connectivity has been previously compared between species using resting-state fMRI, but not in electrophysiological studies.

METHODS

We compared connectivity with implanted electrodes in humans (electrocorticography) to macaques and sheep (microelectrocorticography), which are capable of recording neural data at high frequencies with spatial precision. We specifically examined synchrony, implicated in functional integration between regions.

RESULTS

We found that connectivity strength was overwhelmingly similar in humans and monkeys for pairs of two different brain regions (prefrontal, motor, premotor, parietal), but differed more often within single brain regions. The two connectivity measures, correlation and phase locking value, were similar in most comparisons. Connectivity strength agreed more often between the species at higher frequencies. Where the species differed, monkey synchrony was stronger than human in all but one case. In contrast, human and sheep connectivity within somatosensory cortex diverged in almost all frequencies, with human connectivity stronger than sheep.

DISCUSSION

Our findings imply greater heterogeneity within regions in humans than in monkeys, but comparable functional interactions between regions in the two species. This suggests that monkeys may be effectively used to probe resting-state connectivity in humans, and that such findings can then be validated in humans. Although the discrepancy between humans and sheep is larger, we suggest that findings from sheep in highly invasive studies may be used to provide guidance for studies in other species.

Intrinsic frequency specific brain networks for identification of MCI individuals using resting-state fMRI

Numerous brain oscillations are well organized into several brain rhythms to support complex brain activities within distinct frequency bands. These rhythms temporally coexist in the same or different brain areas and may interact with each other with specific properties and physiological functions. However, the identification and evaluation of these various brain rhythms derived from BOLD-fMRI signals are obscure. To address this issue, we introduced a data-driven method named Complementary Ensemble Empirical Mode Decomposition (CEEMD) to automatically decompose the BOLD oscillations into several brain rhythms within distinct frequency bands. Thereafter, in order to evaluate the performance of CEEMD in the detection of subtle BOLD signals, a novel CEEMD-based high-dimensional pattern classification framework was proposed to accurately identify mild cognitive impairment individuals from the healthy controls. Our results showed CEEMD is a stable frequency decomposition method. Furthermore, CEEMD-based frequency specific topological profiles provided a classification accuracy of 93.33%, which was saliently higher than that of the conventional frequency separation based scheme. Importantly, our findings demonstrated that CEEMD could provide an effective means for brain oscillation separation, by which a more meaningful frequency bins could be used to detect the subtle changes embedded in the BOLD signals.

BCI Use and Its Relation to Adaptation in Cortical Networks

Brain-computer interfaces (BCIs) carry great potential in the treatment of motor impairments. As a new motor output, BCIs interface with the native motor system, but acquisition of BCI proficiency requires a degree of learning to integrate this new function. In this review, we discuss how BCI designs often take advantage of the brain's motor system infrastructure as sources of command signals. We highlight a growing body of literature examining how this approach leads to changes in activity across cortex, including beyond motor regions, as a result of learning the new skill of BCI control. We discuss the previous research identifying patterns of neural activity associated with BCI skill acquisition and use that closely resembles those associated with learning traditional native motor tasks. We then discuss recent work in animals probing changes in connectivity of the BCI control site, which were linked to BCI skill acquisition, and use this as a foundation for our original work in humans. We present our novel work showing changes in resting state connectivity across cortex following the BCI learning process. We find substantial, heterogeneous changes in connectivity across regions and frequencies, including interactions that do not involve the BCI control site. We conclude from our review and original work that BCI skill acquisition may potentially lead to significant changes in evoked and resting state connectivity across multiple cortical regions. We recommend that future studies of BCIs look beyond motor regions to fully describe the cortical networks involved and long-term adaptations resulting from BCI skill acquisition.

The Detection of Phase Amplitude Coupling during Sensory Processing

There is increasing interest in understanding how the phase and amplitude of distinct neural oscillations might interact to support dynamic communication within the brain. In particular, previous work has demonstrated a coupling between the phase of low frequency oscillations and the amplitude (or power) of high frequency oscillations during certain tasks, termed phase amplitude coupling (PAC). For instance, during visual processing in humans, PAC has been reliably observed between ongoing alpha (8–13 Hz) and gamma-band (>40 Hz) activity. However, the application of PAC metrics to electrophysiological data can be challenging due to numerous methodological issues and lack of coherent approaches within the field. Therefore, in this article we outline the various analysis steps involved in detecting PAC, using an openly available MEG dataset from 16 participants performing an interactive visual task. Firstly, we localized gamma and alpha-band power using the Fieldtrip toolbox, and extracted time courses from area V1, defined using a multimodal parcelation scheme. These V1 responses were analyzed for changes in alpha-gamma PAC, using four common algorithms. Results showed an increase in alpha (7–13 Hz)–gamma (40–100 Hz) PAC in response to the visual grating stimulus, though specific patterns of coupling were somewhat dependent upon the algorithm employed. Additionally, post-hoc analyses showed that these results were not driven by the presence of non-sinusoidal oscillations, and that trial length was sufficient to obtain reliable PAC estimates. Finally, throughout the article, methodological issues and practical guidelines for ongoing PAC research will be discussed.

Nasal Respiration Entrains Human Limbic Oscillations and Modulates Cognitive Function

The need to breathe links the mammalian olfactory system inextricably to the respiratory rhythms that draw air through the nose. In rodents and other small animals, slow oscillations of local field potential activity are driven at the rate of breathing (∼2-12 Hz) in olfactory bulb and cortex, and faster oscillatory bursts are coupled to specific phases of the respiratory cycle. These dynamic rhythms are thought to regulate cortical excitability and coordinate network interactions, helping to shape olfactory coding, memory, and behavior. However, while respiratory oscillations are a ubiquitous hallmark of olfactory system function in animals, direct evidence for such patterns is lacking in humans. In this study, we acquired intracranial EEG data from rare patients (Ps) with medically refractory epilepsy, enabling us to test the hypothesis that cortical oscillatory activity would be entrained to the human respiratory cycle, albeit at the much slower rhythm of ∼0.16-0.33 Hz. Our results reveal that natural breathing synchronizes electrical activity in human piriform (olfactory) cortex, as well as in limbic-related brain areas, including amygdala and hippocampus. Notably, oscillatory power peaked during inspiration and dissipated when breathing was diverted from nose to mouth. Parallel behavioral experiments showed that breathing phase enhances fear discrimination and memory retrieval. Our findings provide a unique framework for understanding the pivotal role of nasal breathing in coordinating neuronal oscillations to support stimulus processing and behavior.

SIGNIFICANCE STATEMENT

Animal studies have long shown that olfactory oscillatory activity emerges in line with the natural rhythm of breathing, even in the absence of an odor stimulus. Whether the breathing cycle induces cortical oscillations in the human brain is poorly understood. In this study, we collected intracranial EEG data from rare patients with medically intractable epilepsy, and found evidence for respiratory entrainment of local field potential activity in human piriform cortex, amygdala, and hippocampus. These effects diminished when breathing was diverted to the mouth, highlighting the importance of nasal airflow for generating respiratory oscillations. Finally, behavioral data in healthy subjects suggest that breathing phase systematically influences cognitive tasks related to amygdala and hippocampal functions.

Phase-amplitude coupling at the organism level: The amplitude of spontaneous alpha rhythm fluctuations varies with the phase of the infra-slow gastric basal rhythm

A fundamental feature of the temporal organization of neural activity is phase-amplitude coupling between brain rhythms at different frequencies, where the amplitude of a higher frequency varies according to the phase of a lower frequency. Here, we show that this rule extends to brain-organ interactions. We measured both the infra-slow (~0.05Hz) rhythm intrinsically generated by the stomach - the gastric basal rhythm - using electrogastrography, and spontaneous brain dynamics with magnetoencephalography during resting-state with eyes open. We found significant phase-amplitude coupling between the infra-slow gastric phase and the amplitude of the cortical alpha rhythm (10-11Hz), with gastric phase accounting for 8% of the variance of alpha rhythm amplitude fluctuations. Gastric-alpha coupling was localized to the right anterior insula, and bilaterally to occipito-parietal regions. Transfer entropy, a measure of directionality of information transfer, indicates that gastric-alpha coupling is due to an ascending influence from the stomach to both the right anterior insula and occipito-parietal regions. Our results show that phase-amplitude coupling so far only observed within the brain extends to brain-viscera interactions. They further reveal that the temporal structure of spontaneous brain activity depends not only on neuron and network properties endogenous to the brain, but also on the slow electrical rhythm generated by the stomach.

From intentions to actions: Neural oscillations encode motor processes through phase, amplitude and phase-amplitude coupling

Goal-directed motor behavior is associated with changes in patterns of rhythmic neuronal activity across widely distributed brain areas. In particular, movement initiation and execution are mediated by patterns of synchronization and desynchronization that occur concurrently across distinct frequency bands and across multiple motor cortical areas. To date, motor-related local oscillatory modulations have been predominantly examined by quantifying increases or suppressions in spectral power. However, beyond signal power, spectral properties such as phase and phase-amplitude coupling (PAC) have also been shown to carry information with regards to the oscillatory dynamics underlying motor processes. Yet, the distinct functional roles of phase, amplitude and PAC across the planning and execution of goal-directed motor behavior remain largely elusive. Here, we address this question with unprecedented resolution thanks to multi-site intracerebral EEG recordings in human subjects while they performed a delayed motor task. To compare the roles of phase, amplitude and PAC, we monitored intracranial brain signals from 748 sites across six medically intractable epilepsy patients at movement execution, and during the delay period where motor intention is present but execution is withheld. In particular, we used a machine-learning framework to identify the key contributions of various neuronal responses. We found a high degree of overlap between brain network patterns observed during planning and those present during execution. Prominent amplitude increases in the delta (2-4Hz) and high gamma (60-200Hz) bands were observed during both planning and execution. In contrast, motor alpha (8-13Hz) and beta (13-30Hz) power were suppressed during execution, but enhanced during the delay period. Interestingly, single-trial classification revealed that low-frequency phase information, rather than spectral power change, was the most discriminant feature in dissociating action from intention. Additionally, despite providing weaker decoding, PAC features led to statistically significant classification of motor states, particularly in anterior cingulate cortex and premotor brain areas. These results advance our understanding of the distinct and partly overlapping involvement of phase, amplitude and the coupling between them, in the neuronal mechanisms underlying motor intentions and executions.

Evaluation of Phase-Amplitude Coupling in Resting State Magnetoencephalographic Signals: Effect of Surrogates and Evaluation Approach

Phase-amplitude coupling (PAC) plays an important role in neural communication and computation. Interestingly, recent studies have indicated the presence of ubiquitous PAC phenomenon even during the resting state. Despite the importance of PAC phenomenon, estimation of significant physiological PAC is challenging because of the lack of appropriate surrogate measures to control false positives caused by non-physiological PAC. Therefore, in the present study, we evaluated PAC phenomenon during resting-state magnetoencephalography (MEG) signal and considered various surrogate measures and computational approaches widely used in the literature in addition to proposing new ones. We evaluated PAC phenomenon over the entire length of the MEG signal and for multiple shorter time segments. The results indicate that the extent of PAC phenomenon mainly depends on the surrogate measures and PAC computational methods used, as well as the evaluation approach. After a careful and critical evaluation, we found that resting-state MEG signals failed to exhibit ubiquitous PAC phenomenon, contrary to what has been suggested previously.

Regional Patterns of Cortical Phase Synchrony in the Resting State

Synchronized phase estimates between oscillating neuronal signals at the macroscale level reflect coordinated activities between neuronal assemblies. Recent electrophysiological evidence suggests the presence of significant spontaneous phase synchrony within the resting state. The purpose of this study was to investigate phase synchrony, including directional interactions, in resting state subdural electrocorticographic recordings to better characterize patterns of regional phase interactions across the lateral cortical surface during the resting state. We estimated spontaneous phase locking value (PLV) as a measure of functional connectivity, and phase slope index (PSI) as a measure of pseudo-causal phase interactions, across a broad range of canonical frequency bands and the modulation of the amplitude envelope of high gamma (amHG), a band that is believed to best reflect the physiological processes giving rise to the functional magnetic resonance imaging BOLD signal. Long-distance interactions had higher PLVs in slower frequencies (≤theta) than in higher ones (≥beta) with amHG behaving more like slow frequencies, and a general trend of increasing frequency band of significant PLVs when moving across the lateral surface along an anterior-posterior axis. Moreover, there was a strong trend of frontal-to-parietal directional phase synchronization, measured by PSI across multiple frequencies. These findings, which are likely indicative of coordinated and structured spontaneous cortical interactions, are important in the study of time scales and directional nature of resting state functional connectivity, and may ultimately contribute to a better understanding of how spontaneous synchrony is linked to variation in regional architecture across the lateral cortical surface.