This article is a commentary on:
Evaluation of Phase-Amplitude Coupling in Resting State Magnetoencephalographic Signals: Effect of Surrogates and Evaluation Approach
Esther Florin
Sylvain Baillet- 1Medical Faculty, Institute of Clinical Neuroscience and Medical Psychology, Heinrich-Heine University, Düsseldorf, Germany
- 2McConnell Brain Imaging Center, Montreal Neurological Institute, McGill University, Montreal, QC, Canada
A commentary on
Evaluation of Phase-Amplitude Coupling in Resting State Magnetoencephalographic Signals: Effect of Surrogates and Evaluation Approach
by Gohel, B., Lim, S., Kim,M.-Y., An, K., Kim, J.-E., Kwon, H., et al. (2016). Front. Comput. Neurosci. 10:120. doi: 10.3389/fncom.2016.00120
1. Introduction
Using Human Connectome Project (HCP) data, Gohel et al. (2016) (GEA) recently reported that “resting-state MEG signals failed to exhibit ubiquitous phase-amplitude coupling (PAC) phenomenon, contrary to what has been suggested” by Florin and Baillet (2015) (FB). GEA argued that the original PAC findings by FB were driven by false positives resulting from the use of inappropriate methods. In this commentary, we first correct GEA's mischaracterization of the approach actually used by FB. We then investigated the PAC computations in Gohel et al. (2016) (GEA) and demonstrate that with FB's approach, it is actually possible to detect PAC in the Human Connectome Project (HCP) resting-state data. Finally, when making the data processing as similar as possible to GEA we still found significant PAC across large portions of the brain.
2. Mischaracterization of Florin and Baillet (2015)
GEA claimed that FB's significance threshold was “generated after combining all surrogate PAC scores from all frequency pairs” and argued that such pooling “of PAC scores across the frequency pairs is not appropriate,” because it would bias the finding of significant couplings toward low frequencies. This claim is incorrect: The PAC detection threshold in FB was actually computed separately for each frequency pair (page 28 in FB: “the value at the 95th quantile of PAC scores for every low/high frequency pair was determined as the significance threshold at p < 0.05 for actual recordings”).
3. Technical Differences
3.1. Artifact Removal
From GEA's statement that at least 10 independent components were removed from the data, it seems that they used their own preprocessing pipeline, not the distributed preprocessed version of the HCP data, because this latter had fewer independent components removed. After direct communication with GEA, and in an attempt to reproduce their published method, we used the preprocessed HCP data set.
3.2. Source Modeling
GEA used a single-sphere head model, while FB used multiple-sphere forward modeling. Even though the latter is best-practice in MEG analysis (Gross et al., 2013), we tried to reproduce for the purpose of this commentary the approach of GEA by performing analyses with a single-sphere head model. However, their arbitrary selection of 58 sources where they computed PAC measures was not described with sufficient details to enable the identical reproduction of their analyses. In addition, such a small sample of sources can be expected to have biased GEA's results.
3.3. Frequency Selection
In the FB original analysis we used all frequency pairs (fφ, fa), with fφ ∈ [2, 48] Hz and fa ∈ [80, 150] Hz for the calculation of PAC. The binning width for fφ was 0.5 Hz from 2 to 12 Hz and 2 Hz from 12 to 48 Hz. The bins for fa were logarithmically spaced. Following Canolty et al. (2006), we used chirplets to extract the phase and amplitude with a chirp-factor of 0, hence yielding wavelet-like decompositions.
GEA chose three low frequencies (3, 6, 10 Hz) and as high frequencies the sum over two frequency ranges, 53:10:93 and 83:10:143 Hz. To extract phase and amplitude they used Morlet wavelets with a mother wavelet of full width at half maximum (FWHM) of 1s and a center frequency of 1 Hz. This choice leads to 46 Hz wide FWHM filter at 53 Hz. This results in considerable smearing of the high-frequency binning.
3.4. Statistical Detection of PAC
GEA used different statistical methods to determine statistical significance of PAC values. For the purpose of the present rebuttal, we only used one approach from GEA, which was also used in FB: We generated pink noise time series of the same lengths as the original data and calculated PAC on these surrogate data. We generated 500 pink noise time series as in GEA. For the determination of the significance threshold, GEA used a more conservative 99% quantile, while FB used 95%.
4. Method Comparison With HCP Data
We used the preprocessed data from four HCP subjects (100307, 108323, 113922, 877168). All subjects were included by GEA and we only used the one run described in the supplementary material of GEA. For preprocessing, as many ICA components were removed as described in the icaclass folder of the preprocessed data. The further analyses were conducted with Brainstorm (Tadel et al., 2011). For source reconstruction we used Fieldtrip's (Oostenveld et al., 2011) implementation to derive a single sphere head model and estimate source activations with a weighted minimum norm estimator. This is presumably the same approach as GEA used. We then computed PAC according to the method proposed by Özkurt et al. (2011) and used in FB and GEA at each of the 8,004 individual source locations.
As could be expected, using a higher significance threshold lead to fewer significant findings than in FB (Table 1). Comparing Chirplets to Morlets, fewer vertices with significant PAC were found with Morlet wavelets. Importantly the binning of low frequencies and summing of the gamma amplitudes as in GEA led to almost no PAC detection. Overall sparse sampling of the low frequencies had a greater influence on the detection of PAC than the sparse sampling in the gamma range.
TABLE 1
Table 1. Average percentage of significant PAC across the brain for HCP subjects 100307, 108323, 113922, and 877168 with different frequency binning options.
5. Discussion
We investigated the claim from GEA that the results of FB were only due to false positives. We did so by using four of the HCP datasets used by GEA, staying as close as possible to their approach.
Using the sparse low-frequency binning of GEA led to pronounced reduction in PAC detection. Thus, it is likely that the actual low-frequency phase component was missed using GEA's methodology. We believe that not restricting PAC analysis to few frequencies pairs provides a more encompassing picture of potential oscillatory coupling in the brain.
We used pink (1/f) Gaussian noise for the determination of the significance threshold as in FB. GEA demonstrated that using pink Gaussian noise is less conservative than the non-colored noise approach. While this might be the case, pink noise is the preferred option, because neural signals exhibit a 1/f frequency spectrum and non-parametric statistical testing requires that one should target the original data properties of non inferential interest with surrogate data. Of note for future research is that the non-sinusoidal shape of brain oscillations might affect the detection of cross-frequency coupling depending on the method used to extract the phase and amplitude coupling data (Yeh et al., 2016; Cole and Voytek, 2017). In principle, chirplets and wavelets can be constructed to account for the expected non-sinusoidal shape of brain oscillations (Mallat, 2009).
In summary, we can only restate that ubiquitous PAC in resting state MEG data is present across the brain. This was also the case for data from the HCP repository if all potential low-frequency high-frequency pairs are tested exhaustively.
Ethics Statment
We used publicly available, open access data from the human connectome project.
Author Contributions
All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.
Funding
EF gratefully acknowledges support from the Volkswagen Foundation (89387). SB was supported by a Discovery Grant from the National Science and Engineering Research Council of Canada and the NIH (2R01EB009048-05).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Canolty, R. T., Edwards, E., Dalal, S. S., Soltani, M., Nagarajan, S. S., Kirsch, H. E., et al. (2006). High gamma power is phase-locked to theta oscillations in human neocortex. Science 313, 1626–1628. doi: 10.1126/science.1128115
Cole, S. R., and Voytek, B. (2017). Brain oscillations and the importance of waveform shape. Trends Cogn. Sci. 21, 137–149. doi: 10.1016/j.tics.2016.12.008
Florin, E., and Baillet, S. (2015). The brain's resting-state activity is shaped by synchronized cross-frequency coupling of neural oscillations. Neuroimage 111, 26–35. doi: 10.1016/j.neuroimage.2015.01.054
Gohel, B., Lim, S., Kim, M. Y., An, K., Kim, J. E., Kwon, H., et al. (2016). Evaluation of phase-amplitude coupling in resting state magnetoencephalographic signals: effect of surrogates and evaluation approach. Front. Comput. Neurosci. 10:120. doi: 10.3389/fncom.2016.00120
Gross, J., Baillet, S., Barnes, G. R., Henson, R. N., Hillebrand, A., Jensen, O., et al. (2013). Good practice for conducting and reporting MEG research. Neuroimage 65, 349–363. doi: 10.1016/j.neuroimage.2012.10.001
Oostenveld, R., Pascal, F., Maris, E., and Schoffelen, J. M. (2011). FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. Comput. Intell. Neurosci. 2011:156869. doi: 10.1155/2011/156869
Özkurt, T. E., and Schnitzler, A. (2011). A critical note on the definition of phase-amplitude cross-frequency coupling. J. Neurosci. Methods 201, 438–443. doi: 10.1016/j.jneumeth.2011.08.014
Tadel, F., Baillet, S., Mosher, J. C., Pantazis, D., and Leahy, R. M. (2011). Brainstorm: a user-friendly application for MEG/EEG analysis. Comput. Intell. Neurosci. 2011:879716. doi: 10.1155/2011/879716
Keywords: phase-amplitude coupling (PAC), neural oscillations, resting state, magnetoencephalography (MEG), surrogate data
Citation: Florin E and Baillet S (2018) Commentary: Evaluation of Phase-Amplitude Coupling in Resting State Magnetoencephalographic Signals: Effect of Surrogates and Evaluation Approach. Front. Comput. Neurosci. 12:26. doi: 10.3389/fncom.2018.00026
Received: 21 December 2017; Accepted: 29 March 2018;
Published: 16 April 2018.
Edited by:
Xiaoli Li, Beijing Normal University, ChinaCopyright © 2018 Florin and Baillet. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Esther Florin, ZXN0aGVyLmZsb3JpbkBoaHUuZGU=