Phase-based measures of cross-frequency coupling in brain electrical dynamics under general anesthesia

The state of general anesthesia (GA) is associated with an increase in spectral power in scalp electroencephalogram (EEG) at frequencies below 40 Hz, including spectral peaks in the slow oscillation (SO, 0.1-1 Hz) and α (8-14 Hz) bands. Because conventional power spectral analyses are insensitive to possible cross-frequency coupling, the relationships among the oscillations at different frequencies remain largely unexplored. Quantifying such coupling is essential for improving clinical monitoring of anesthesia and understanding the neuroscience of this brain state. We tested the usefulness of two measures of cross-frequency coupling: the bispectrum-derived SynchFastSlow, which is sensitive to phase-phase coupling in different frequency bands, and modulogram analysis of coupling between SO phase and α rhythm amplitude. SynchFastSlow, a metric that is used in clinical depth-of-anesthesia monitors, showed a robust correlation with the loss of consciousness at the induction of propofol GA, but this could be largely explained by power spectral changes without considering cross-frequency coupling. Modulogram analysis revealed two distinct modes of cross-frequency coupling under GA. The waking and two distinct states under GA could be discriminated by projecting in a two-dimensional phase space defined by the SynchFastSlow and the preferred SO phase of α activity. Our results show that a stereotyped pattern of phase-amplitude coupling accompanies multiple stages of anesthetic-induced unconsciousness. These findings suggest that modulogram analysis can improve EEG based monitoring of brain state under GA.National Institutes of Health (U.S.) (Grant DP2-OD006454)National Institutes of Health (U.S.) (Grant K25-NS057580)National Institutes of Health (U.S.) (Grant DP1-OD003646)National Institutes of Health (U.S.) (Grant R01-EB006385)National Institutes of Health (U.S.) (Grant R01-MH071847

Robust time-varying multivariate coherence estimation: Application to electroencephalogram recordings during general anesthesia

Coherence analysis characterizes frequency-dependent covariance between signals, and is useful for multivariate oscillatory data often encountered in neuroscience. The global coherence provides a summary of coherent behavior in high-dimensional multivariate data by quantifying the concentration of variance in the first mode of an eigenvalue decomposition of the cross-spectral matrix. Practical application of this useful method is sensitive to noise, and can confound coherent activity in disparate neural populations or spatial locations that have a similar frequency structure. In this paper we describe two methodological enhancements to the global coherence procedure that increase robustness of the technique to noise, and that allow characterization of how power within specific coherent modes change through time.National Institutes of Health (U.S.) (Grant DP2-OD006454)National Institutes of Health (U.S.) (Grant K25-NS057580)National Institutes of Health (U.S.) (Grant DP1-OD003646)National Institutes of Health (U.S.) (Grant R01-EB006385)National Institutes of Health (U.S.) (Grant R01-MH071847

State space methods for phase amplitude coupling analysis

Phase amplitude coupling (PAC) is thought to play a fundamental role in the dynamic coordination of brain circuits and systems. There are however growing concerns that existing methods for PAC analysis are prone to error and misinterpretation. Improper frequency band selection can render true PAC undetectable, while non-linearities or abrupt changes in the signal can produce spurious PAC. Current methods require large amounts of data and lack formal statistical inference tools. We describe here a novel approach for PAC analysis that substantially addresses these problems. We use a state space model to estimate the component oscillations, avoiding problems with frequency band selection, nonlinearities, and sharp signal transitions. We represent cross-frequency coupling in parametric and time-varying forms to further improve statistical efficiency and estimate the posterior distribution of the coupling parameters to derive their credible intervals. We demonstrate the method using simulated data, rat local field potentials (LFP) data, and human EEG data.P01GM118269 - NIH HHS; R01AG056015 - NIH HHS; R01AG054081 - NIH HHS; R21DA048323 - NIH HHSPublished versio