No abstract

Despite being the focus of a thriving field of research, the biological mechanisms that underlie information integration in the brain are not yet fully understood. A theory that has gained a lot of traction in recent years suggests that multi-scale integration is regulated by a hierarchy of mutually interacting neural oscillations. In particular, there is accumulating evidence that phase-amplitude coupling (PAC), a specific form of cross-frequency interaction, plays a key role in numerous cognitive processes. Current research in the field is not only hampered by the absence of a gold standard for PAC analysis, but also by the computational costs of running exhaustive computations on large and high-dimensional electrophysiological brain signals. In addition, various signal properties and analyses parameters can lead to spurious PAC. Here, we present Tensorpac, an open-source Python toolbox dedicated to PAC analysis of neurophysiological data. The advantages of Tensorpac include (1) higher computational efficiency thanks to software design that combines tensor computations and parallel computing, (2) the implementation of all most widely used PAC methods in one package, (3) the statistical analysis of PAC measures, and (4) extended PAC visualization capabilities. Tensorpac is distributed under a BSD-3-Clause license and can be launched on any operating system (Linux, OSX and Windows). It can be installed directly via pip or downloaded from Github (https://github.com/EtienneCmb/tensorpac). By making Tensorpac available, we aim to enhance the reproducibility and quality of PAC research, and provide open tools that will accelerate future method development in neuroscience.

We describe a method to detect brain activation in cortically constrained maps of current density computed from MEG data using multivariate statistical inference. We apply time-frequency (wavelet) analysis to individual epochs to produce dynamic images of brain signal power on the cerebral cortex in multiple time-frequency bands. We form vector observations by concatenating the power in each frequency band, and fit them into separate multivariate linear models for each time band and cortical location with experimental conditions as predictor variables. The resulting Roy's maximum statistic maps are thresholded for significance using permutation tests and the maximum statistic approach. A source is considered significant if it exceeds a statistical threshold, which is chosen to control the familywise error rate, or the probability of at least one false positive, across the cortical surface. We compare and evaluate the multivariate approach with existing univariate approaches to time-frequency MEG signal analysis, both on simulated data and experimental data from an MEG visuomotor task study. Our results indicate that the multivariate method is more powerful than the univariate approach in detecting experimental effects when correlations exist between power across frequency bands. We further describe protected F-tests and linear discriminant analysis to identify individual frequencies that contribute significantly to experimental effects.

No abstract

Cross-frequency phase-amplitude coupling (PAC) within large neuronal populations is hypothesized to play a functional role in information processing in a range of cognitive tasks. The goal of our study was to examine the putative role of PAC in the brain networks that mediate continuous visuomotor control. We estimated the cortical activity that mediates visuomotor control via magnetoencephalography (MEG) recordings in 15 healthy volunteers. We extracted the cortical signal amplitudes and phases at the frequencies of interest by means of band-pass filtering followed by Hilbert transforms. To quantify task-related changes of PAC, we implemented a technique based on the Kullback-Leibler divergence. The choice of this technique among others was based on the results of comparisons performed on simulations of coupled sources in various noise conditions. The application of PAC to the MEG data revealed a significant task-related increase in coupling between the phase of delta (2-5 Hz) and the amplitude of high-gamma (60-90 Hz) oscillations in the occipital and parietal cortices as well as in the cerebellum. Remarkably, when comparing PAC in the early trials to PAC recorded towards the end of the experiment we found a significant increase in delta-high-gamma coupling over time in the superior parietal lobule, possibly reflecting visuomotor adaptation processes. Our results suggest that, in addition to power modulations, cross-frequency interactions play a key role in visuomotor behavior.