1
|
Hindriks R, Broeders TAA, Schoonheim MM, Douw L, Santos F, van Wieringen W, Tewarie PKB. Higher-order functional connectivity analysis of resting-state functional magnetic resonance imaging data using multivariate cumulants. Hum Brain Mapp 2024; 45:e26663. [PMID: 38520377 PMCID: PMC10960559 DOI: 10.1002/hbm.26663] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Revised: 02/12/2024] [Accepted: 03/08/2024] [Indexed: 03/25/2024] Open
Abstract
Blood-level oxygenation-dependent (BOLD) functional magnetic resonance imaging (fMRI) is the most common modality to study functional connectivity in the human brain. Most research to date has focused on connectivity between pairs of brain regions. However, attention has recently turned towards connectivity involving more than two regions, that is, higher-order connectivity. It is not yet clear how higher-order connectivity can best be quantified. The measures that are currently in use cannot distinguish between pairwise (i.e., second-order) and higher-order connectivity. We show that genuine higher-order connectivity can be quantified by using multivariate cumulants. We explore the use of multivariate cumulants for quantifying higher-order connectivity and the performance of block bootstrapping for statistical inference. In particular, we formulate a generative model for fMRI signals exhibiting higher-order connectivity and use it to assess bias, standard errors, and detection probabilities. Application to resting-state fMRI data from the Human Connectome Project demonstrates that spontaneous fMRI signals are organized into higher-order networks that are distinct from second-order resting-state networks. Application to a clinical cohort of patients with multiple sclerosis further demonstrates that cumulants can be used to classify disease groups and explain behavioral variability. Hence, we present a novel framework to reliably estimate genuine higher-order connectivity in fMRI data which can be used for constructing hyperedges, and finally, which can readily be applied to fMRI data from populations with neuropsychiatric disease or cognitive neuroscientific experiments.
Collapse
Affiliation(s)
- Rikkert Hindriks
- Department of Mathematics, Faculty of ScienceVrije Universiteit AmsterdamAmsterdamThe Netherlands
| | - Tommy A. A. Broeders
- Department of Anatomy and Neurosciences, Amsterdam NeuroscienceAmsterdam UMC, Vrije Universiteit AmsterdamAmsterdamThe Netherlands
| | - Menno M. Schoonheim
- Department of Anatomy and Neurosciences, Amsterdam NeuroscienceAmsterdam UMC, Vrije Universiteit AmsterdamAmsterdamThe Netherlands
| | - Linda Douw
- Department of Anatomy and Neurosciences, Amsterdam NeuroscienceAmsterdam UMC, Vrije Universiteit AmsterdamAmsterdamThe Netherlands
| | - Fernando Santos
- Dutch Institute for Emergent Phenomena (DIEP)Institute for Advanced Studies, University of AmsterdamAmsterdamThe Netherlands
- Korteweg de Vries Institute for MathematicsUniversity of AmsterdamAmsterdamthe Netherlands
| | - Wessel van Wieringen
- Department of Epidemiology and BiostatisticsAmsterdam UMC, Vrije Universiteit AmsterdamAmsterdamThe Netherlands
| | - Prejaas K. B. Tewarie
- Sir Peter Mansfield Imaging CenterSchool of Physics, University of NottinghamNottinghamUnited Kingdom
- Clinical Neurophysiology GroupUniversity of TwenteEnschedeThe Netherlands
| |
Collapse
|
2
|
Tewarie PKB, Hindriks R, Lai YM, Sotiropoulos SN, Kringelbach M, Deco G. Non-reversibility outperforms functional connectivity in characterisation of brain states in MEG data. Neuroimage 2023; 276:120186. [PMID: 37268096 DOI: 10.1016/j.neuroimage.2023.120186] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2023] [Revised: 04/27/2023] [Accepted: 05/22/2023] [Indexed: 06/04/2023] Open
Abstract
Characterising brain states during tasks is common practice for many neuroscientific experiments using electrophysiological modalities such as electroencephalography (EEG) and magnetoencephalography (MEG). Brain states are often described in terms of oscillatory power and correlated brain activity, i.e. functional connectivity. It is, however, not unusual to observe weak task induced functional connectivity alterations in the presence of strong task induced power modulations using classical time-frequency representation of the data. Here, we propose that non-reversibility, or the temporal asymmetry in functional interactions, may be more sensitive to characterise task induced brain states than functional connectivity. As a second step, we explore causal mechanisms of non-reversibility in MEG data using whole brain computational models. We include working memory, motor, language tasks and resting-state data from participants of the Human Connectome Project (HCP). Non-reversibility is derived from the lagged amplitude envelope correlation (LAEC), and is based on asymmetry of the forward and reversed cross-correlations of the amplitude envelopes. Using random forests, we find that non-reversibility outperforms functional connectivity in the identification of task induced brain states. Non-reversibility shows especially better sensitivity to capture bottom-up gamma induced brain states across all tasks, but also alpha band associated brain states. Using whole brain computational models we find that asymmetry in the effective connectivity and axonal conduction delays play a major role in shaping non-reversibility across the brain. Our work paves the way for better sensitivity in characterising brain states during both bottom-up as well as top-down modulation in future neuroscientific experiments.
Collapse
Affiliation(s)
- Prejaas K B Tewarie
- Center for Brain and Cognition, Computational Neuroscience Group, Department of Information and Communication Technologies, Universitat Pompeu Fabra, Spain; Clinical Neurophysiology Group, University of Twente, Enschede, The Netherlands; Department of Neurology, Amsterdam UMC, Amsterdam, the Netherlands; Sir Peter Mansfield Imaging Centre, School of Physics, University of Nottingham, Nottingham, United Kingdom.
| | - Rikkert Hindriks
- Department of Mathematics, Faculty of Science, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands
| | - Yi Ming Lai
- Sir Peter Mansfield Imaging Centre, School of Medicine, University of Nottingham, United Kingdom
| | - Stamatios N Sotiropoulos
- Sir Peter Mansfield Imaging Centre, School of Medicine, University of Nottingham, United Kingdom; NIHR Biomedical Research Centre, University of Nottingham, Nottingham University Hospitals NHS Trust, Nottingham, UK
| | - Morten Kringelbach
- Centre for Eudaimonia and Human Flourishing, Linacre College, University of Oxford, Oxford, UK; Center for Music in the Brain, Department of Clinical Medicine, Aarhus University, Aarhus, Denmark; Department of Psychiatry, University of Oxford, Oxford, UK
| | - Gustavo Deco
- Center for Brain and Cognition, Computational Neuroscience Group, Department of Information and Communication Technologies, Universitat Pompeu Fabra, Spain; Institució Catalana de la Recerca i Estudis Avançats (ICREA), Barcelona, Spain; Department of Neuropsychology, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany
| |
Collapse
|
3
|
Tewarie PKB, Tjepkema-Cloostermans MC, Abeysuriya RG, Hofmeijer J, van Putten MJAM. Preservation of thalamocortical circuitry is essential for good recovery after cardiac arrest. PNAS Nexus 2023; 2:pgad119. [PMID: 37143862 PMCID: PMC10153639 DOI: 10.1093/pnasnexus/pgad119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Revised: 03/10/2023] [Accepted: 04/04/2023] [Indexed: 05/06/2023]
Abstract
Continuous electroencephalographam (EEG) monitoring contributes to prediction of neurological outcome in comatose cardiac arrest survivors. While the phenomenology of EEG abnormalities in postanoxic encephalopathy is well known, the pathophysiology, especially the presumed role of selective synaptic failure, is less understood. To further this understanding, we estimate biophysical model parameters from the EEG power spectra from individual patients with a good or poor recovery from a postanoxic encephalopathy. This biophysical model includes intracortical, intrathalamic, and corticothalamic synaptic strengths, as well as synaptic time constants and axonal conduction delays. We used continuous EEG measurements from hundred comatose patients recorded during the first 48 h postcardiac arrest, 50 with a poor neurological outcome [cerebral performance category ( CPC = 5 ) ] and 50 with a good neurological outcome ( CPC = 1 ). We only included patients that developed (dis-)continuous EEG activity within 48 h postcardiac arrest. For patients with a good outcome, we observed an initial relative excitation in the corticothalamic loop and corticothalamic propagation that subsequently evolved towards values observed in healthy controls. For patients with a poor outcome, we observed an initial increase in the cortical excitation-inhibition ratio, increased relative inhibition in the corticothalamic loop, delayed corticothalamic propagation of neuronal activity, and severely prolonged synaptic time constants that did not return to physiological values. We conclude that the abnormal EEG evolution in patients with a poor neurological recovery after cardiac arrest may result from persistent and selective synaptic failure that includes corticothalamic circuitry and also delayed corticothalamic propagation.
Collapse
Affiliation(s)
| | - Marleen C Tjepkema-Cloostermans
- Clinical Neurophysiology Group, University of Twente, 7522 NH Enschede, Netherlands
- Department of Neurology and Clinical Neurophysiology, Medisch Spectrum Twente, 7512 KZ Enschede, Netherlands
| | - Romesh G Abeysuriya
- Computational Epidemic Modelling, Burnet Institute, 3004 Melbourne, Australia
| | - Jeannette Hofmeijer
- Clinical Neurophysiology Group, University of Twente, 7522 NH Enschede, Netherlands
- Department of Neurology, Rijnstate Hospital, 6815 AD Arnhem, Netherlands
| | | |
Collapse
|
4
|
Hindriks R, Tewarie PKB. Dissociation between phase and power correlation networks in the human brain is driven by co-occurrent bursts. Commun Biol 2023; 6:286. [PMID: 36934153 PMCID: PMC10024695 DOI: 10.1038/s42003-023-04648-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2022] [Accepted: 03/02/2023] [Indexed: 03/20/2023] Open
Abstract
Well-known haemodynamic resting-state networks are better mirrored in power correlation networks than phase coupling networks in electrophysiological data. However, what do these power correlation networks reflect? We address this long-outstanding question in neuroscience using rigorous mathematical analysis, biophysical simulations with ground truth and application of these mathematical concepts to empirical magnetoencephalography (MEG) data. Our mathematical derivations show that for two non-Gaussian electrophysiological signals, their power correlation depends on their coherence, cokurtosis and conjugate-coherence. Only coherence and cokurtosis contribute to power correlation networks in MEG data, but cokurtosis is less affected by artefactual signal leakage and better mirrors haemodynamic resting-state networks. Simulations and MEG data show that cokurtosis may reflect co-occurrent bursting events. Our findings shed light on the origin of the complementary nature of power correlation networks to phase coupling networks and suggests that the origin of resting-state networks is partly reflected in co-occurent bursts in neuronal activity.
Collapse
Affiliation(s)
- Rikkert Hindriks
- Department of Mathematics, Faculty of Science, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.
| | - Prejaas K B Tewarie
- Clinical Neurophysiology Group, University of Twente, Enschede, The Netherlands
- Sir Peter Mansfield Imaging Center, School of Physics, University of Nottingham, Nottingham, UK
| |
Collapse
|
5
|
Tewarie PKB, Beernink TMJ, Eertman-Meyer CJ, Cornet AD, Beishuizen A, van Putten MJAM, Tjepkema-Cloostermans MC. Early EEG monitoring predicts clinical outcome in patients with moderate to severe traumatic brain injury. Neuroimage Clin 2023; 37:103350. [PMID: 36801601 PMCID: PMC9984683 DOI: 10.1016/j.nicl.2023.103350] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Revised: 01/23/2023] [Accepted: 02/11/2023] [Indexed: 02/16/2023]
Abstract
There is a need for reliable predictors in patients with moderate to severe traumatic brain injury to assist clinical decision making. We assess the ability of early continuous EEG monitoring at the intensive care unit (ICU) in patients with traumatic brain injury (TBI) to predict long term clinical outcome and evaluate its complementary value to current clinical standards. We performed continuous EEG measurements in patients with moderate to severe TBI during the first week of ICU admission. We assessed the Extended Glasgow Outcome Scale (GOSE) at 12 months, dichotomized into poor (GOSE 1-3) and good (GOSE 4-8) outcome. We extracted EEG spectral features, brain symmetry index, coherence, aperiodic exponent of the power spectrum, long range temporal correlations, and broken detailed balance. A random forest classifier using feature selection was trained to predict poor clinical outcome based on EEG features at 12, 24, 48, 72 and 96 h after trauma. We compared our predictor with the IMPACT score, the best available predictor, based on clinical, radiological and laboratory findings. In addition we created a combined model using EEG as well as the clinical, radiological and laboratory findings. We included hundred-seven patients. The best prediction model using EEG parameters was found at 72 h after trauma with an AUC of 0.82 (0.69-0.92), specificity of 0.83 (0.67-0.99) and sensitivity of 0.74 (0.63-0.93). The IMPACT score predicted poor outcome with an AUC of 0.81 (0.62-0.93), sensitivity of 0.86 (0.74-0.96) and specificity of 0.70 (0.43-0.83). A model using EEG and clinical, radiological and laboratory parameters resulted in a better prediction of poor outcome (p < 0.001) with an AUC of 0.89 (0.72-0.99), sensitivity of 0.83 (0.62-0.93) and specificity of 0.85 (0.75-1.00). EEG features have potential use for predicting clinical outcome and decision making in patients with moderate to severe TBI and provide complementary information to current clinical standards.
Collapse
Affiliation(s)
- Prejaas K B Tewarie
- Clinical Neurophysiology Group, University of Twente, Enschede, the Netherlands; Department of Neurology and Clinical Neurophysiology, Medisch Spectrum Twente, Enschede, the Netherlands; Department of Neurology, Amsterdam UMC/VUmc, Amsterdam, the Netherlands.
| | - Tim M J Beernink
- Department of Neurology and Clinical Neurophysiology, Medisch Spectrum Twente, Enschede, the Netherlands
| | - Carin J Eertman-Meyer
- Department of Neurology and Clinical Neurophysiology, Medisch Spectrum Twente, Enschede, the Netherlands
| | - Alexander D Cornet
- Intensive Care Center, Medisch Spectrum Twente, Enschede, the Netherlands
| | | | - Michel J A M van Putten
- Clinical Neurophysiology Group, University of Twente, Enschede, the Netherlands; Department of Neurology and Clinical Neurophysiology, Medisch Spectrum Twente, Enschede, the Netherlands
| | - Marleen C Tjepkema-Cloostermans
- Clinical Neurophysiology Group, University of Twente, Enschede, the Netherlands; Department of Neurology and Clinical Neurophysiology, Medisch Spectrum Twente, Enschede, the Netherlands
| |
Collapse
|