Improved 7 Tesla resting-state fMRI connectivity measurements by cluster-based modeling of respiratory volume and heart rate effects

Comparing end-tidal CO2, respiration volume per time (RVT), and average gray matter signal for mapping cerebrovascular reactivity amplitude and delay with breath-hold task BOLD fMRI

Cerebrovascular reactivity (CVR), defined as the cerebral blood flow response to a vasoactive stimulus, is an imaging biomarker with demonstrated utility in a range of diseases and in typical development and aging processes. A robust and widely implemented method to map CVR involves using a breath-hold task during a BOLD fMRI scan. Recording end-tidal CO₂ (P_ETCO₂) changes during the breath-hold task is recommended to be used as a reference signal for modeling CVR amplitude in standard units (%BOLD/mmHg) and CVR delay in seconds. However, obtaining reliable P_ETCO₂ recordings requires equipment and task compliance that may not be achievable in all settings. To address this challenge, we investigated two alternative reference signals to map CVR amplitude and delay in a lagged general linear model (lagged-GLM) framework: respiration volume per time (RVT) and average gray matter BOLD response (GM-BOLD). In 8 healthy adults with multiple scan sessions, we compare spatial agreement of CVR maps from RVT and GM-BOLD to those generated with P_ETCO₂. We define a threshold to determine whether a P_ETCO₂ recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient P_ETCO₂ and 6 datasets with insufficient P_ETCO₂. When P_ETCO₂ quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from P_ETCO₂ (after accounting for differences in scale), with the caveat they are not in standard units to facilitate between-group comparisons. CVR delays are comparable to P_ETCO₂ with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when P_ETCO₂ quality is insufficient, RVT and GM-BOLD CVR recover reasonable CVR amplitude and delay maps, provided the participant attempted the breath-hold task. Therefore, our framework offers a solution for achieving high quality CVR maps in both retrospective and prospective studies where sufficient P_ETCO₂ recordings are not available and especially in populations where obtaining reliable measurements is a known challenge (e.g., children). Our results have the potential to improve the accessibility of CVR mapping and to increase the prevalence of this promising metric of vascular health.

Depression Classification Using Frequent Subgraph Mining Based on Pattern Growth of Frequent Edge in Functional Magnetic Resonance Imaging Uncertain Network

The brain network structure is highly uncertain due to the noise in imaging signals and evaluation methods. Recent works have shown that uncertain brain networks could capture uncertain information with regards to functional connections. Most of the existing research studies covering uncertain brain networks used graph mining methods for analysis; for example, the mining uncertain subgraph patterns (MUSE) method was used to mine frequent subgraphs and the discriminative feature selection for uncertain graph classification (DUG) method was used to select discriminant subgraphs. However, these methods led to a lack of effective discriminative information; this reduced the classification accuracy for brain diseases. Therefore, considering these problems, we propose an approximate frequent subgraph mining algorithm based on pattern growth of frequent edge (unFEPG) for uncertain brain networks and a novel discriminative feature selection method based on statistical index (dfsSI) to perform graph mining and selection. Results showed that compared with the conventional methods, the unFEPG and dfsSI methods achieved a higher classification accuracy. Furthermore, to demonstrate the efficacy of the proposed method, we used consistent discriminative subgraph patterns based on thresholding and weighting approaches to compare the classification performance of uncertain networks and certain networks in a bidirectional manner. Results showed that classification performance of the uncertain network was superior to that of the certain network within a defined sparsity range. This indicated that if a better classification performance is to be achieved, it is necessary to select a certain brain network with a higher threshold or an uncertain brain network model. Moreover, if the uncertain brain network model was selected, it is necessary to make full use of the uncertain information of its functional connection.

Advances in resting state fMRI acquisitions for functional connectomics

Resting state functional magnetic resonance imaging (rs-fMRI) is based on spontaneous fluctuations in the blood oxygen level dependent (BOLD) signal, which occur simultaneously in different brain regions, without the subject performing an explicit task. The low-frequency oscillations of the rs-fMRI signal demonstrate an intrinsic spatiotemporal organization in the brain (brain networks) that may relate to the underlying neural activity. In this review article, we briefly describe the current acquisition techniques for rs-fMRI data, from the most common approaches for resting state acquisition strategies, to more recent investigations with dedicated hardware and ultra-high fields. Specific sequences that allow very fast acquisitions, or multiple echoes, are discussed next. We then consider how acquisition methods weighted towards specific parts of the BOLD signal, like the Cerebral Blood Flow (CBF) or Volume (CBV), can provide more spatially specific network information. These approaches are being developed alongside the commonly used BOLD-weighted acquisitions. Finally, specific applications of rs-fMRI to challenging regions such as the laminae in the neocortex, and the networks within the large areas of subcortical white matter regions are discussed. We finish the review with recommendations for acquisition strategies for a range of typical applications of resting state fMRI.

Physiological and motion signatures in static and time-varying functional connectivity and their subject identifiability

Human brain connectivity yields significant potential as a noninvasive biomarker. Several studies have used fMRI-based connectivity fingerprinting to characterize individual patterns of brain activity. However, it is not clear whether these patterns mainly reflect neural activity or the effect of physiological and motion processes. To answer this question, we capitalize on a large data sample from the Human Connectome Project and rigorously investigate the contribution of the aforementioned processes on functional connectivity (FC) and time-varying FC, as well as their contribution to subject identifiability. We find that head motion, as well as heart rate and breathing fluctuations, induce artifactual connectivity within distinct resting-state networks and that they correlate with recurrent patterns in time-varying FC. Even though the spatiotemporal signatures of these processes yield above-chance levels in subject identifiability, removing their effects at the preprocessing stage improves identifiability, suggesting a neural component underpinning the inter-individual differences in connectivity.

Cerebrovascular Reactivity Mapping Without Gas Challenges: A Methodological Guide

Cerebrovascular reactivity (CVR) is defined as the ability of vessels to alter their caliber in response to vasoactive factors, by means of dilating or constricting, in order to increase or decrease regional cerebral blood flow (CBF). Importantly, CVR may provide a sensitive biomarker for pathologies where vasculature is compromised. Furthermore, the spatiotemporal dynamics of CVR observed in healthy subjects, reflecting regional differences in cerebral vascular tone and response, may also be important in functional MRI studies based on neurovascular coupling mechanisms. Assessment of CVR is usually based on the use of a vasoactive stimulus combined with a CBF measurement technique. Although transcranial Doppler ultrasound has been frequently used to obtain global flow velocity measurements, MRI techniques are being increasingly employed for obtaining CBF maps. For the vasoactive stimulus, vasodilatory hypercapnia is usually induced through the manipulation of respiratory gases, including the inhalation of increased concentrations of carbon dioxide. However, most of these methods require an additional apparatus and complex setups, which not only may not be well-tolerated by some populations but are also not widely available. For these reasons, strategies based on voluntary breathing fluctuations without the need for external gas challenges have been proposed. These include the task-based methodologies of breath holding and paced deep breathing, as well as a new generation of methods based on spontaneous breathing fluctuations during resting-state. Despite the multitude of alternatives to gas challenges, existing literature lacks definitive conclusions regarding the best practices for the vasoactive modulation and associated analysis protocols. In this work, we perform an extensive review of CVR mapping techniques based on MRI and CO2 variations without gas challenges, focusing on the methodological aspects of the breathing protocols and corresponding data analysis. Finally, we outline a set of practical guidelines based on generally accepted practices and available data, extending previous reports and encouraging the wider application of CVR mapping methodologies in both clinical and academic MRI settings.

Impact of physiological noise in characterizing the functional MRI default-mode network in Alzheimer's disease

The functional connectivity of the default-mode network (DMN) monitored by functional magnetic resonance imaging (fMRI) in Alzheimer's disease (AD) patients has been found weaker than that in healthy participants. Since breathing and heart beating can cause fluctuations in the fMRI signal, these physiological activities may affect the fMRI data differently between AD patients and healthy participants. We collected resting-state fMRI data from AD patients and age-matched healthy participants. With concurrent cardiac and respiratory recordings, we estimated both physiological responses phase-locked and non-phase-locked to heart beating and breathing. We found that the cardiac and respiratory physiological responses in AD patients were 3.00 ± 0.51 s and 3.96 ± 0.52 s later (both p < 0.0001) than those in healthy participants, respectively. After correcting the physiological noise in the resting-state fMRI data by population-specific physiological response functions, the DMN estimated by seed-correlation was more localized to the seed region. The DMN difference between AD patients and healthy controls became insignificant after suppressing physiological noise. Our results indicate the importance of controlling physiological noise in the resting-state fMRI analysis to obtain clinically related characterizations in AD.

Controlling for lesions, kinematics and physiological noise: impact on fMRI results of spastic post-stroke patients

Functional magnetic resonance imaging (fMRI) is a widely used technique for assessing brain function in both healthy and pathological populations. Some factors, such as motion, physiological noise and lesion presence, can contribute to signal change and confound the fMRI data, but fMRI data processing techniques have been developed to correct for these confounding effects. Fifteen spastic subacute stroke patients underwent fMRI while performing a highly controlled task (i.e. passive extension of their affected and unaffected wrists). We investigated the impact on activation maps of lesion masking during preprocessing and first- and second-level analyses, and of adding wrist extension amplitudes and physiological data as regressors using the Statistical Parametric Mapping toolbox (SPM12).

We observed a significant decrease in sensorimotor region activation after the addition of lesion masks and movement/physiological regressors during the processing of stroke patients’ fMRI data. Our results demonstrate that:•

The unified segmentation routine results in good normalization accuracy when dealing with stroke lesions regardless of their size;

•

Adding a group lesion mask during the second-level analysis seems to be a suitable option when none of the patients have lesions in target regions. Otherwise, no masking is acceptable;

•

Movement amplitude is a significant contributor to the sensorimotor activation observed during passive wrist extension in spastic stroke patients;

•

Movement features and physiological noise are relevant factors when interpreting for sensorimotor activation in studies of the motor system in patients with brain lesions. They can be added as nuisance covariates during large patient groups’ analyses.

Resting-state "physiological networks"

Slow changes in systemic brain physiology can elicit large fluctuations in fMRI time series, which manifest as structured spatial patterns of temporal correlations between distant brain regions. Here, we investigated whether such "physiological networks"-sets of segregated brain regions that exhibit similar responses following slow changes in systemic physiology-resemble patterns associated with large-scale networks typically attributed to remotely synchronized neuronal activity. By analyzing a large group of subjects from the 3T Human Connectome Project (HCP) database, we demonstrate brain-wide and noticeably heterogenous dynamics tightly coupled to either respiratory variation or heart rate changes. We show, using synthesized data generated from physiological recordings across subjects, that these physiologically-coupled fluctuations alone can produce networks that strongly resemble previously reported resting-state networks, suggesting that, in some cases, the "physiological networks" seem to mimic the neuronal networks. Further, we show that such physiologically-relevant connectivity estimates appear to dominate the overall connectivity observations in multiple HCP subjects, and that this apparent "physiological connectivity" cannot be removed by the use of a single nuisance regressor for the entire brain (such as global signal regression) due to the clear regional heterogeneity of the physiologically-coupled responses. Our results challenge previous notions that physiological confounds are either localized to large veins or globally coherent across the cortex, therefore emphasizing the necessity to consider potential physiological contributions in fMRI-based functional connectivity studies. The rich spatiotemporal patterns carried by such "physiological" dynamics also suggest great potential for clinical biomarkers that are complementary to large-scale neuronal networks.

EEG-Informed fMRI: A Review of Data Analysis Methods

The simultaneous acquisition of electroencephalography (EEG) with functional magnetic resonance imaging (fMRI) is a very promising non-invasive technique for the study of human brain function. Despite continuous improvements, it remains a challenging technique, and a standard methodology for data analysis is yet to be established. Here we review the methodologies that are currently available to address the challenges at each step of the data analysis pipeline. We start by surveying methods for pre-processing both EEG and fMRI data. On the EEG side, we focus on the correction for several MR-induced artifacts, particularly the gradient and pulse artifacts, as well as other sources of EEG artifacts. On the fMRI side, we consider image artifacts induced by the presence of EEG hardware inside the MR scanner, and the contamination of the fMRI signal by physiological noise of non-neuronal origin, including a review of several approaches to model and remove it. We then provide an overview of the approaches specifically employed for the integration of EEG and fMRI when using EEG to predict the blood oxygenation level dependent (BOLD) fMRI signal, the so-called EEG-informed fMRI integration strategy, the most commonly used strategy in EEG-fMRI research. Finally, we systematically review methods used for the extraction of EEG features reflecting neuronal phenomena of interest.