A multi-measure approach for assessing the performance of fMRI preprocessing strategies in resting-state functional connectivity

Comparing data-driven physiological denoising approaches for resting-state fMRI: implications for the study of aging

Introduction

Physiological nuisance contributions by cardiac and respiratory signals have a significant impact on resting-state fMRI data quality. As these physiological signals are often not recorded, data-driven denoising methods are commonly used to estimate and remove physiological noise from fMRI data. To investigate the efficacy of these denoising methods, one of the first steps is to accurately capture the cardiac and respiratory signals, which requires acquiring fMRI data with high temporal resolution.

Methods

In this study, we used such high-temporal resolution fMRI data to evaluate the effectiveness of several data-driven denoising methods, including global-signal regression (GSR), white matter and cerebrospinal fluid regression (WM-CSF), anatomical (aCompCor) and temporal CompCor (tCompCor), ICA-AROMA. Our analysis focused on the consequence of changes in low-frequency, cardiac and respiratory signal power, as well as age-related differences in terms of functional connectivity (fcMRI).

Results

Our results confirm that the ICA-AROMA and GSR removed the most physiological noise but also more low-frequency signals. These methods are also associated with substantially lower age-related fcMRI differences. On the other hand, aCompCor and tCompCor appear to be better at removing high-frequency physiological signals but not low-frequency signal power. These methods are also associated with relatively higher age-related fcMRI differences, whether driven by neuronal signal or residual artifact. These results were reproduced in data downsampled to represent conventional fMRI sampling frequency. Lastly, methods differ in performance depending on the age group.

Discussion

While this study cautions direct comparisons of fcMRI results based on different denoising methods in the study of aging, it also enhances the understanding of different denoising methods in broader fcMRI applications.

Cortical timescales and the modular organization of structural and functional brain networks

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

Test-retest reliability and predictive utility of a macroscale principal functional connectivity gradient

Mapping individual differences in brain function has been hampered by poor reliability as well as limited interpretability. Leveraging patterns of brain-wide functional connectivity (FC) offers some promise in this endeavor. In particular, a macroscale principal FC gradient that recapitulates a hierarchical organization spanning molecular, cellular, and circuit level features along a sensory-to-association cortical axis has emerged as both a parsimonious and interpretable measure of individual differences in behavior. However, the measurement reliabilities of this FC gradient have not been fully evaluated. Here, we assess the reliabilities of both global and regional principal FC gradient measures using test-retest data from the young adult Human Connectome Project (HCP-YA) and the Dunedin Study. Analyses revealed that the reliabilities of principal FC gradient measures were (1) consistently higher than those for traditional edge-wise FC measures, (2) higher for FC measures derived from general FC (GFC) in comparison with resting-state FC, and (3) higher for longer scan lengths. We additionally examined the relative utility of these principal FC gradient measures in predicting cognition and aging in both datasets as well as the HCP-aging dataset. These analyses revealed that regional FC gradient measures and global gradient range were significantly associated with aging in all three datasets, and moderately associated with cognition in the HCP-YA and Dunedin Study datasets, reflecting contractions and expansions of the cortical hierarchy, respectively. Collectively, these results demonstrate that measures of the principal FC gradient, especially derived using GFC, effectively capture a reliable feature of the human brain subject to interpretable and biologically meaningful individual variation, offering some advantages over traditional edge-wise FC measures in the search for brain-behavior associations.

Brain Connectivity Correlates of Breathing and Cardiac Irregularities in SUDEP: A Resting-State fMRI Study

Sudden unexpected death in epilepsy (SUDEP) is the leading cause of premature mortality among people with epilepsy. Evidence from witnessed and monitored SUDEP cases indicate seizure-induced cardiovascular and respiratory failures; yet, the underlying mechanisms remain obscure. SUDEP occurs often during the night and early morning hours, suggesting that sleep or circadian rhythm-induced changes in physiology contribute to the fatal event. Resting-state fMRI studies have found altered functional connectivity between brain structures involved in cardiorespiratory regulation in later SUDEP cases and in individuals at high-risk of SUDEP. However, those connectivity findings have not been related to changes in cardiovascular or respiratory patterns. Here, we compared fMRI patterns of brain connectivity associated with regular and irregular cardiorespiratory rhythms in SUDEP cases with those of living epilepsy patients of varying SUDEP risk, and healthy controls. We analysed resting-state fMRI data from 98 patients with epilepsy (9 who subsequently succumbed to SUDEP, 43 categorized as low SUDEP risk (no tonic-clonic seizures (TCS) in the year preceding the fMRI scan), and 46 as high SUDEP risk (>3 TCS in the year preceding the scan)) and 25 healthy controls. The global signal amplitude (GSA), defined as the moving standard deviation of the fMRI global signal, was used to identify periods with regular ('low state') and irregular ('high state') cardiorespiratory rhythms. Correlation maps were derived from seeds in twelve regions with a key role in autonomic or respiratory regulation, for the low and high states. Following principal component analysis, component weights were compared between the groups. We found widespread alterations in connectivity of precuneus/posterior cingulate cortex in epilepsy compared to controls, in the low state (regular cardiorespiratory activity). In the low state, and to a lesser degree in the high state, reduced anterior insula connectivity (mainly with anterior and posterior cingulate cortex) in epilepsy appeared, relative to healthy controls. For SUDEP cases, the insula connectivity differences were inversely related to the interval between the fMRI scan and death. The findings suggest that anterior insula connectivity measures may provide a biomarker of SUDEP risk. The neural correlates of autonomic brain structures associated with different cardiorespiratory rhythms may shed light on the mechanisms underlying terminal apnea observed in SUDEP.

Extracting electrophysiological correlates of functional magnetic resonance imaging data using the canonical polyadic decomposition

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.

Prediction of Obsessive-Compulsive Disorder: Importance of Neurobiology-Aided Feature Design and Cross-Diagnosis Transfer Learning

BACKGROUND

Machine learning applications using neuroimaging provide a multidimensional, data-driven approach that captures the level of complexity necessary for objectively aiding diagnosis and prognosis in psychiatry. However, models learned from small training samples often have limited generalizability, which continues to be a problem with automated diagnosis of mental illnesses such as obsessive-compulsive disorder (OCD). Earlier studies have shown that features incorporating prior neurobiological knowledge of brain function and combining brain parcellations from various sources can potentially improve the overall prediction. However, it is unknown whether such knowledge-driven methods can provide a performance that is comparable to state-of-the-art approaches based on neural networks.

METHODS

In this study, we apply a transparent and explainable multiparcellation ensemble learning framework EMPaSchiz (Ensemble algorithm with Multiple Parcellations for Schizophrenia prediction) to the task of predicting OCD, based on a resting-state functional magnetic resonance imaging dataset of 350 subjects. Furthermore, we apply transfer learning using the features found effective for schizophrenia to OCD to leverage the commonality in brain alterations across these psychiatric diagnoses.

RESULTS

We show that our knowledge-based approach leads to a prediction performance of 80.3% accuracy for OCD diagnosis that is better than domain-agnostic and automated feature design using neural networks. Furthermore, we show that a selection of reduced feature sets can be transferred from schizophrenia to the OCD prediction model without significant loss in prediction performance.

CONCLUSIONS

This study presents a machine learning framework for OCD prediction with neurobiology-aided feature design using resting-state functional magnetic resonance imaging that is generalizable and reasonably interpretable.

Future Directions for Chemosensory Connectomes: Best Practices and Specific Challenges

Ecological chemosensory stimuli almost always evoke responses in more than one sensory system. Moreover, any sensory processing takes place along a hierarchy of brain regions. So far, the field of chemosensory neuroimaging is dominated by studies that examine the role of brain regions in isolation. However, to completely understand neural processing of chemosensation, we must also examine interactions between regions. In general, the use of connectivity methods has increased in the neuroimaging field, providing important insights to physical sensory processing, such as vision, audition, and touch. A similar trend has been observed in chemosensory neuroimaging, however, these established techniques have largely not been rigorously applied to imaging studies on the chemical senses, leaving network insights overlooked. In this article, we first highlight some recent work in chemosensory connectomics and we summarize different connectomics techniques. Then, we outline specific challenges for chemosensory connectome neuroimaging studies. Finally, we review best practices from the general connectomics and neuroimaging fields. We recommend future studies to develop or use the following methods we perceive as key to improve chemosensory connectomics: (1) optimized study designs, (2) reporting guidelines, (3) consensus on brain parcellations, (4) consortium research, and (5) data sharing.