Granger causal time-dependent source connectivity in the somatosensory network

Exploration of transient Granger causal interactions in neural sources of electrophysiological activities provides deeper insights into brain information processing mechanisms. However, the underlying neural patterns are confounded by time-dependent dynamics, non-stationarity and observational noise contamination. Here we investigate transient Granger causal interactions using source time-series of somatosensory evoked magnetoencephalographic (MEG) elicited by air puff stimulation of right index finger and recorded using 306-channel MEG from 21 healthy subjects. A new time-varying connectivity approach, combining renormalised partial directed coherence with state space modelling, is employed to estimate fast changing information flow among the sources. Source analysis confirmed that somatosensory evoked MEG was mainly generated from the contralateral primary somatosensory cortex (SI) and bilateral secondary somatosensory cortices (SII). Transient Granger causality shows a serial processing of somatosensory information, 1) from contralateral SI to contralateral SII, 2) from contralateral SI to ipsilateral SII, 3) from contralateral SII to contralateral SI, and 4) from contralateral SII to ipsilateral SII. These results are consistent with established anatomical connectivity between somatosensory regions and previous source modeling results, thereby providing empirical validation of the time-varying connectivity analysis. We argue that the suggested approach provides novel information regarding transient cortical dynamic connectivity, which previous approaches could not assess.

[1]  F. Varela,et al.  Neuromagnetic imaging of cortical oscillations accompanying tactile stimulation. , 2003, Brain research. Cognitive brain research.

[2]  Riitta Salmelin,et al.  Comparison of BOLD fMRI and MEG characteristics to vibrotactile stimulation , 2003, NeuroImage.

[3]  N. Forss,et al.  Functional characterization of human second somatosensory cortex by magnetoencephalography , 2002, Behavioural Brain Research.

[4]  Tongsheng Zhang,et al.  Facilitating Neuronal Connectivity Analysis of Evoked Responses by Exposing Local Activity with Principal Component Analysis Preprocessing: Simulation of Evoked MEG , 2013, Brain Topography.

[5]  A. Mouraux,et al.  Nociceptive laser-evoked brain potentials do not reflect nociceptive-specific neural activity. , 2009, Journal of neurophysiology.

[6]  R. Ilmoniemi,et al.  Signal-space projection method for separating MEG or EEG into components , 1997, Medical and Biological Engineering and Computing.

[7]  竹安 数博,et al.  Time series analysis and its applications , 2007 .

[8]  Ernst Nennig,et al.  Time-efficient localization of the human secondary somatosensory cortex by functional magnetic resonance imaging , 2005, Neuroscience Letters.

[9]  Gian Domenico Iannetti,et al.  A multisensory investigation of the functional significance of the “pain matrix” , 2011, NeuroImage.

[10]  D. Haar,et al.  Statistical Physics , 1971, Nature.

[11]  J. Martinerie,et al.  The brainweb: Phase synchronization and large-scale integration , 2001, Nature Reviews Neuroscience.

[12]  G. T. Coleman,et al.  Functional characteristics of the parallel SI- and SII-projecting neurons of the thalamic ventral posterior nucleus in the marmoset. , 2001, Journal of neurophysiology.

[13]  Ralph G. Andrzejak,et al.  Detecting event-related time-dependent directional couplings , 2006 .

[14]  Hiroshi Shibasaki,et al.  Attention modulates both primary and second somatosensory cortical activities in humans: a magnetoencephalographic study. , 1998, Journal of neurophysiology.

[15]  J. Downar,et al.  A multimodal cortical network for the detection of changes in the sensory environment , 2000, Nature Neuroscience.

[16]  Klaus Lehnertz,et al.  Inferring directional interactions from transient signals with symbolic transfer entropy. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  Paolo Maria Rossini,et al.  Topographic Organization of the Human Primary and Secondary Somatosensory Cortices: Comparison of fMRI and MEG Findings , 2002, NeuroImage.

[18]  C Braun,et al.  Unilaterally applied stimuli in a frequency discrimination task are represented bilaterally in primary somatosensory cortex. , 2004, Neurology & clinical neurophysiology : NCN.

[19]  Masahito Yamada,et al.  Structural Time Series Models and the Kalman Filter , 1989 .

[20]  S. Taulu,et al.  The Signal Space Separation method , 2004, physics/0401166.

[21]  W. Backes,et al.  Somatosensory cortex responses to median nerve stimulation: fMRI effects of current amplitude and selective attention , 2000, Clinical Neurophysiology.

[22]  Celso Grebogi,et al.  Overarching framework for data-based modelling , 2014 .

[23]  G. Innocenti.  General Organization of Callosal Connections in the Cerebral Cortex , 1986 .

[24]  Yoshio Okada,et al.  Recursive artifact windowed–single tone extraction method (RAW–STEM) as periodic noise filter for electrophysiological signals with interfering transients , 2006, Journal of Neuroscience Methods.

[25]  C. C. Wood,et al.  Human cortical potentials evoked by stimulation of the median nerve. I. Cytoarchitectonic areas generating short-latency activity. , 1989, Journal of neurophysiology.

[26]  S. Taulu,et al.  Presentation of electromagnetic multichannel data: The signal space separation method , 2005 .

[27]  David C. Alsop,et al.  Mapping of secondary somatosensory cortex activation induced by vibrational stimulation: an fMRI study , 1999, Brain Research.

[28]  Maria Blatow,et al.  fMRI reflects functional connectivity of human somatosensory cortex , 2007, NeuroImage.

[29]  Andrew Harvey,et al.  Forecasting, Structural Time Series Models and the Kalman Filter. , 1991 .

[30]  L. Lauronen,et al.  Trigeminal somatosensory evoked magnetic fields to tactile stimulation , 2006, Clinical Neurophysiology.

[31]  Matthew T. Sutherland,et al.  Reliable detection of bilateral activation in human primary somatosensory cortex by unilateral median nerve stimulation , 2006, NeuroImage.

[32]  V. Jousmäki,et al.  Activation of a distributed somatosensory cortical network in the human brain. A dipole modelling study of magnetic fields evoked by median nerve stimulation. Part I: Location and activation timing of SEF sources. , 1997, Electroencephalography and clinical neurophysiology.

[33]  M Mishkin,et al.  Serial and parallel processing of tactual information in somatosensory cortex of rhesus monkeys. , 1992, Journal of neurophysiology.

[34]  J. Kaas,et al.  The relation of corpus callosum connections to architectonic fields and body surface maps in sensorimotor cortex of new and old world monkeys , 1983, The Journal of comparative neurology.

[35]  Riitta Hari,et al.  Transient Suppression of Ipsilateral Primary Somatosensory Cortex during Tactile Finger Stimulation , 2006, The Journal of Neuroscience.

[36]  Lars Arendt-Nielsen,et al.  Dynamic Brain Topography of Somatosensory Evoked Potentials and Equivalent Dipoles in Response to Graded Painful Skin and Muscle Stimulation , 2004, Brain Topography.

[37]  Gian Domenico Iannetti,et al.  Representation of different trigeminal divisions within the primary and secondary human somatosensory cortex , 2003, NeuroImage.

[38]  M A Nicolelis,et al.  Bilateral Integration of Whisker Information in the Primary Somatosensory Cortex of Rats , 2001, The Journal of Neuroscience.

[39]  J S George,et al.  MUSIC seeded multi-dipole MEG modeling using the Constrained Start Spatio-Temporal modeling procedure. , 2004, Neurology & clinical neurophysiology : NCN.

[40]  R. Hari,et al.  Functional Organization of the Human First and Second Somatosensory Cortices: a Neuromagnetic Study , 1993, The European journal of neuroscience.

[41]  M. Eichler,et al.  Assessing the strength of directed influences among neural signals using renormalized partial directed coherence , 2009, Journal of Neuroscience Methods.

[42]  Nina Forss,et al.  Characteristics of the human contra- versus ipsilateral SII cortex , 2000, Clinical Neurophysiology.

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

[44]  N. M. Maurits,et al.  Original article/Article originalMultichannel recording of median nerve somatosensory evoked potentialsEnregistrement multicanal des PES du nerf médian , 2008 .

[45]  Y. Iwamura.  Hierarchical somatosensory processing , 1998, Current Opinion in Neurobiology.

[46]  N. Maurits,et al.  Multichannel recording of median nerve somatosensory evoked potentials , 2008, Neurophysiologie Clinique/Clinical Neurophysiology.

[47]  C. Granger.  Investigating causal relations by econometric models and cross-spectral methods , 1969 .

[48]  Li Hu,et al.  A time-varying source connectivity approach to reveal human somatosensory information processing , 2012, NeuroImage.

[49]  N. M. Maurits,et al.  Multichannel recording of tibial-nerve somatosensory evoked potentials , 2008, Neurophysiologie Clinique/Clinical Neurophysiology.

[50]  C C Wood,et al.  Human cortical potentials evoked by stimulation of the median nerve. II. Cytoarchitectonic areas generating long-latency activity. , 1989, Journal of neurophysiology.

[51]  Felix Blankenburg,et al.  Stimulus-dependent EEG activity reflects internal updating of tactile working memory in humans , 2011, Proceedings of the National Academy of Sciences.

[52]  R. Kakigi,et al.  Serial processing in the human somatosensory system. , 2004, Cerebral cortex.

[53]  C. Grebogi,et al.  Inference of Granger causal time-dependent influences in noisy multivariate time series , 2012, Journal of Neuroscience Methods.

[54]  Karin Schwab,et al.  Comparison of linear signal processing techniques to infer directed interactions in multivariate neural systems , 2005, Signal Process..

[55]  Luiz A. Baccalá,et al.  Partial directed coherence: a new concept in neural structure determination , 2001, Biological Cybernetics.

[56]  Rodrigo Quian Quiroga,et al.  Nonlinear multivariate analysis of neurophysiological signals , 2005, Progress in Neurobiology.

[57]  Jonathan D. Cryer,et al.  Time Series Analysis , 1986 .

[58]  R. Ilmoniemi,et al.  Interpreting magnetic fields of the brain: minimum norm estimates , 2006, Medical and Biological Engineering and Computing.

[59]  T. Elbert,et al.  The separation of overlapping neuromagnetic sources in first and second somatosensory cortices , 2005, Brain Topography.

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

[61]  A. Schnitzler,et al.  Normal and pathological oscillatory communication in the brain , 2005, Nature Reviews Neuroscience.

[62]  N Forss,et al.  Functional overlap of finger representations in human SI and SII cortices. , 2001, Journal of neurophysiology.

[63]  G. Buzsáki,et al.  Neuronal Oscillations in Cortical Networks , 2004, Science.

[64]  Takeo Ishigaki,et al.  Contralateral and ipsilateral responses in primary somatosensory cortex following electrical median nerve stimulation—an fMRI study , 2005, Clinical Neurophysiology.

[65]  M. Alexander,et al.  Principles of Neural Science , 1981 .

[66]  Shubhodeep Chakrabarti,et al.  Differential origin of projections from SI barrel cortex to the whisker representations in SII and MI , 2006, The Journal of comparative neurology.

[67]  Andrew Harvey,et al.  Forecasting, Structural Time Series Models and the Kalman Filter , 1990 .

[68]  Koji Inui,et al.  Somatosensory off-response in humans: an ERP study , 2008, Experimental Brain Research.

[69]  Juha Virtanen,et al.  Activation of multiple cortical areas in response to somatosensory stimulation: Combined magnetoencephalographic and functional magnetic resonance imaging , 1999, Human brain mapping.

[70]  Otto W. Witte,et al.  Modelling and analysis of time-variant directed interrelations between brain regions based on BOLD-signals , 2009, NeuroImage.