Extracting electrophysiological correlates of functional magnetic resonance imaging data using the canonical polyadic decomposition. Issue 13 (14th May 2022)

• Record Type: Journal Article
• Authors: Mann‐Krzisnik, Dylan; Mitsis, Georgios D.
• Is Part Of: Human brain mapping. Volume 43:Issue 13(2022)
• Page Start: 4045
• Page End: 4073
• Publication Date: 2022-05-14
• Subjects: computational modelling -- electrophysiology -- magnetic resonance imaging -- motor‐imagery -- neurovascular coupling -- resting‐state -- tensor decomposition; Brain mapping -- Periodicals; 611.81
• DOI: 10.1002/hbm.25902
• Languages: English
• ISSNs: 1065-9471

Abstract

Abstract: The relation between electrophysiology and BOLD‐fMRI requires further elucidation. One approach for studying this relation is to find time‐frequency features from electrophysiology that explain the variance of BOLD time‐series. Convolution of these features with a canonical hemodynamic response function (HRF) is often required to model neurovascular coupling mechanisms and thus account for time shifts between electrophysiological and BOLD‐fMRI data. We propose a framework for extracting the spatial distribution of these time‐frequency features while also estimating more flexible, region‐specific HRFs. The core component of this method is the decomposition of a tensor containing impulse response functions using the Canonical Polyadic Decomposition. The outputs of this decomposition provide insight into the relation between electrophysiology and BOLD‐fMRI and can be used to construct estimates of BOLD time‐series. We demonstrated the performance of this method on simulated data while also examining the effects of simulated measurement noise and physiological confounds. Afterwards, we validated our method on publicly available task‐based and resting‐state EEG‐fMRI data. We adjusted our method to accommodate the multisubject nature of these datasets, enabling the investigation of inter‐subject variability with regards to EEG‐to‐BOLD neurovascular coupling mechanisms. We thus also demonstrate how EEG features for modelling the BOLD signal differ across subjects. Abstract : We show that EEG features of BOLD‐fMRI dynamics can be obtained using the canonical polyadic decomposition. Hemodynamic response functions are also estimated in the process, which allows one to model BOLD‐fMRI signals using EEG signals. We employ our method on simulated data, as well as publicly available task‐based and resting‐state EEG‐fMRI data.