Localizing bicoherence from EEG and MEG

We propose a new method for the localization of nonlinear cross-frequency coupling in EEG and MEG data analysis, based on the estimation of bicoherences at the source level. While for the analysis of rhythmic brain activity, source directions are commonly chosen to maximize power, we suggest to maximize bicoherence instead. The resulting nonlinear cost function can be minimized effectively using a gradient approach. We argue, that bicoherence is also a generally useful tool to analyze phase-amplitude coupling (PAC), by deriving formal relations between PAC and bispectra. This is illustrated in simulated and empirical LFP data. The localization method is applied to EEG resting state data, where the most prominent bicoherence signatures originate from the occipital alpha rhythm and the mu rhythm. While the latter is hardly visible using power analysis, we observe clear bicoherence peaks in the high alpha range of sensorymotor areas. We additionally apply our method to resting-state data of subjects with schizophrenia and healthy controls and observe significant bicoherence differences in motor areas which could not be found from analyzing power differences.