Analysis of functional connectivity and oscillatory power using DICS: from raw MEG data to group-level statistics in Python

Communication between brain regions is thought to be facilitated by the synchronization of oscillatory activity. Hence, large-scale functional networks within the brain may be estimated by measuring synchronicity between regions. Neurophysiological recordings, such as magnetoencephalography (MEG) and electroencephalography (EEG), provide a direct measure of oscillatory neural activity with millisecond temporal resolution. In this paper, we describe a full data analysis pipeline for functional connectivity analysis based on dynamic imaging of coherent sources (DICS) of MEG data. DICS is a beamforming technique in the frequency-domain that enables the study of the cortical sources of oscillatory activity and synchronization between brain regions. All the analysis steps, starting from the raw MEG data up to publication-ready group-level statistics and visualization, are discussed in depth, including methodological considerations, rules of thumb and tradeoffs. We start by computing cross-spectral density (CSD) matrices using a wavelet approach in several frequency bands (alpha, theta, beta, gamma). We then provide a way to create comparable source spaces across subjects and discuss the cortical mapping of spectral power. For connectivity analysis, we present a canonical computation of coherence that facilitates a stable estimation of all-to-all connectivity. Finally, we use group-level statistics to limit the network to cortical regions for which significant differences between experimental conditions are detected and produce vertex-and parcel-level visualizations of the different brain networks. Code examples using the MNE-Python package are provided at each step, guiding the reader through a complete analysis of the freely available openfMRI ds000117 “familiar vs. unfamiliar vs. scrambled faces” dataset. The goal is to educate both novice and experienced data analysts with the “tricks of the trade” necessary to successfully perform this type of analysis on their own data.

[1]  Martin Luessi,et al.  MEG and EEG data analysis with MNE-Python , 2013, Front. Neuroinform..

[2]  R. Ilmoniemi,et al.  Magnetoencephalography-theory, instrumentation, and applications to noninvasive studies of the working human brain , 1993 .

[3]  Joachim Gross,et al.  Investigating causality between interacting brain areas with multivariate autoregressive models of MEG sensor data , 2009, NeuroImage.

[4]  Mia Liljeström,et al.  Large-scale functional networks connect differently for processing words and symbol strings , 2018, PloS one.

[5]  A. Dale,et al.  High‐resolution intersubject averaging and a coordinate system for the cortical surface , 1999, Human brain mapping.

[6]  J. Schoffelen,et al.  Source connectivity analysis with MEG and EEG , 2009, Human brain mapping.

[7]  W. El-Deredy,et al.  Reconstructing Coherent Networks from Electroencephalography and Magnetoencephalography with Reduced Contamination from Volume Conduction or Magnetic Field Spread , 2013, PloS one.

[8]  A. Dale,et al.  Distributed current estimates using cortical orientation constraints , 2006, Human brain mapping.

[9]  J Gross,et al.  REPRINTS , 1962, The Lancet.

[10]  Richard N Henson,et al.  A multi-subject, multi-modal human neuroimaging dataset , 2015, Scientific Data.

[11]  P. Fries.  A mechanism for cognitive dynamics: neuronal communication through neuronal coherence , 2005, Trends in Cognitive Sciences.

[12]  Mia Liljeström,et al.  Task- and stimulus-related cortical networks in language production: Exploring similarity of MEG- and fMRI-derived functional connectivity , 2015, NeuroImage.

[13]  J. Palva,et al.  Neuronal synchrony reveals working memory networks and predicts individual memory capacity , 2010, Proceedings of the National Academy of Sciences.

[14]  Hannu Laaksonen,et al.  Cortical rhythms as markers of neural processing , 2012 .

[15]  Jan Kujala,et al.  A method for spatiotemporal mapping of event-related modulation of cortical rhythmic activity , 2008, NeuroImage.

[16]  Jan Kujala,et al.  Localization of correlated network activity at the cortical level with MEG , 2008, NeuroImage.

[17]  Riitta Salmelin,et al.  Neural representation of language: activation versus long-range connectivity , 2006, Trends in Cognitive Sciences.

[18]  Viviana Betti,et al.  Natural Scenes Viewing Alters the Dynamics of Functional Connectivity in the Human Brain , 2013, Neuron.

[19]  R. Oostenveld,et al.  Nonparametric statistical testing of EEG- and MEG-data , 2007, Journal of Neuroscience Methods.

[20]  Matti Hämäläinen,et al.  MEG/EEG group study with MNE: recommendations, quality assessments and best practices , 2017, bioRxiv.

[21]  C. Rosazza,et al.  Resting-state brain networks: literature review and clinical applications , 2011, Neurological Sciences.

[22]  M. Siegel,et al.  A framework for local cortical oscillation patterns , 2011, Trends in Cognitive Sciences.

[23]  Peter A. Bandettini,et al.  Task-based dynamic functional connectivity: Recent findings and open questions , 2017, NeuroImage.

[24]  Riitta Salmelin,et al.  Task‐sensitive reconfiguration of corticocortical 6–20 Hz oscillatory coherence in naturalistic human performance , 2015, Human brain mapping.

[25]  Robert Oostenveld,et al.  An improved index of phase-synchronization for electrophysiological data in the presence of volume-conduction, noise and sample-size bias , 2011, NeuroImage.

[26]  Anders M. Dale,et al.  Cortical Surface-Based Analysis I. Segmentation and Surface Reconstruction , 1999, NeuroImage.

[27]  Jan Kujala,et al.  Neural interactions at the core of phonological and semantic priming of written words. , 2012, Cerebral cortex.

[28]  A. Dale,et al.  Improved Localizadon of Cortical Activity by Combining EEG and MEG with MRI Cortical Surface Reconstruction: A Linear Approach , 1993, Journal of Cognitive Neuroscience.

[29]  R. Salmelin,et al.  Dynamic reconfiguration of the language network preceding onset of speech in picture naming , 2014, Human brain mapping.

[30]  M. Corbetta,et al.  Large-scale cortical correlation structure of spontaneous oscillatory activity , 2012, Nature Neuroscience.

[31]  Olaf Sporns,et al.  Mechanisms of Zero-Lag Synchronization in Cortical Motifs , 2013, PLoS Comput. Biol..

[32]  J. Kelso,et al.  Cortical coordination dynamics and cognition , 2001, Trends in Cognitive Sciences.

[33]  Peter Hagoort,et al.  Frequency-specific directed interactions in the human brain network for language , 2017 .

[34]  Nikos Makris,et al.  Automatically parcellating the human cerebral cortex. , 2004, Cerebral cortex.

[35]  W. Drongelen,et al.  Localization of brain electrical activity via linearly constrained minimum variance spatial filtering , 1997, IEEE Transactions on Biomedical Engineering.

[36]  O. Bertrand,et al.  Oscillatory gamma-band (30-70 Hz) activity induced by a visual search task in humans. , 1997, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[37]  J. Pernier,et al.  Oscillatory γ-Band (30–70 Hz) Activity Induced by a Visual Search Task in Humans , 1997, The Journal of Neuroscience.

[38]  R. Desimone,et al.  Laminar differences in gamma and alpha coherence in the ventral stream , 2011, Proceedings of the National Academy of Sciences.

[39]  M. Hallett,et al.  Identifying true brain interaction from EEG data using the imaginary part of coherency , 2004, Clinical Neurophysiology.

[40]  Mia Liljeström,et al.  MEG evoked responses and rhythmic activity provide spatiotemporally complementary measures of neural activity in language production , 2012, NeuroImage.

[41]  C. Ciulla,et al.  MEG Characterization of Spontaneous Alpha Rhythm in the Human Brain , 2004, Brain Topography.

[42]  Jan Kujala,et al.  The right hemisphere is highlighted in connected natural speech production and perception , 2017, NeuroImage.

[43]  Mia Liljeström,et al.  Multivariate analysis of correlation between electrophysiological and hemodynamic responses during cognitive processing , 2014, NeuroImage.