48 research outputs found
Removing motion and physiological artifacts from intrinsic BOLD fluctuations using short echo data
AbstractDiffering noise variance across study populations has been shown to cause artifactual group differences in functional connectivity measures. In this study, we investigate the use of short echo time functional MRI data to correct for these noise sources in blood oxygenation level dependent (BOLD)-weighted time series. A dualâecho sequence was used to simultaneously acquire data at both a short (TE=3.3ms) and a BOLD-weighted (TE=35ms) echo time. This approach is effectively âfree,â using dead-time in the pulse sequence to collect an additional echo without affecting overall scan time or temporal resolution. The proposed correction method uses voxelwise regression of the short TE data from the BOLD-weighted data to remove noise variance. In addition to a typical resting state scan, non-compliant behavior associated with patient groups was simulated via increased head motion or physiological fluctuations in 10 subjects. Short TE data showed significant correlations with the traditional motion-related and physiological noise regressors used in current connectivity analyses. Following traditional preprocessing, the extent of significant additional variance explained by the short TE data regressors was significantly correlated with the average head motion across the scan in the resting data (r2=0.93, p<0.0001). The reduction in data variance following the inclusion of short TE regressors was also correlated with scan head motion (r2=0.48, p=0.027). Task-related data were used to demonstrate the effects of the short TE correction on BOLD activation time series with known temporal structure; the size and strength of the activation were significantly decreased, but it is not clear whether this reflects BOLD contamination in the short TE data or correlated changes in blood volume. Finally, functional connectivity maps of the default mode network were constructed using a seed correlation approach. The effects of short TE correction and low-pass filtering on the resulting correlations maps were compared. Results suggest that short TE correction more accurately differentiates artifactual correlations from the correlations of interest in conditions of amplified noise
Potential pitfalls when denoising resting state fMRI data using nuisance regression
In resting state fMRI, it is necessary to remove signal variance associated with noise sources, leaving cleaned fMRI time-series that more accurately reflect the underlying intrinsic brain fluctuations of interest. This is commonly achieved through nuisance regression, in which the fit is calculated of a noise model of head motion and physiological processes to the fMRI data in a General Linear Model, and the âcleanedâ residuals of this fit are used in further analysis. We examine the statistical assumptions and requirements of the General Linear Model, and whether these are met during nuisance regression of resting state fMRI data. Using toy examples and real data we show how pre-whitening, temporal filtering and temporal shifting of regressors impact model fit. Based on our own observations, existing literature, and statistical theory, we make the following recommendations when employing nuisance regression: pre-whitening should be applied to achieve valid statistical inference of the noise model fit parameters; temporal filtering should be incorporated into the noise model to best account for changes in degrees of freedom; temporal shifting of regressors, although merited, should be achieved via optimisation and validation of a single temporal shift. We encourage all readers to make simple, practical changes to their fMRI denoising pipeline, and to regularly assess the appropriateness of the noise model used. By negotiating the potential pitfalls described in this paper, and by clearly reporting the details of nuisance regression in future manuscripts, we hope that the field will achieve more accurate and precise noise models for cleaning the resting state fMRI time-series
Is fMRI ânoiseâ really noise? Resting state nuisance regressors remove variance with network structure
Noise correction is a critical step towards accurate mapping of resting state BOLD fMRI connectivity. Noise sources related to head motion or physiology are typically modelled by nuisance regressors, and a generalised linear model is applied to regress out the associated signal variance. In this study, we use independent component analysis (ICA) to characterise the data variance typically discarded in this pre-processing stage in a cohort of 12 healthy volunteers. The signal variance removed by 24, 12, 6, or only 3 head motion parameters demonstrated network structure typically associated with functional connectivity, and certain networks were discernable in the variance extracted by as few as 2 physiologic regressors. Simulated nuisance regressors, unrelated to the true data noise, also removed variance with network structure, indicating that any group of regressors that randomly sample variance may remove highly structured âsignalâ as well as ânoise.â Furthermore, to support this we demonstrate that random sampling of the original data variance continues to exhibit robust network structure, even when as few as 10% of the original volumes are considered. Finally, we examine the diminishing returns of increasing the number of nuisance regressors used in pre-processing, showing that excessive use of motion regressors may do little better than chance in removing variance within a functional network. It remains an open challenge to understand the balance between the benefits and confounds of noise correction using nuisance regressors
Vascular physiology drives functional brain networks
We present the first evidence for vascular regulation driving fMRI signals in specific functional brain networks. Using concurrent neuronal and vascular stimuli, we collected 30 BOLD fMRI datasets in 10 healthy individuals: a working memory task, flashing checkerboard stimulus, and CO2 inhalation challenge were delivered in concurrent but orthogonal paradigms. The resulting imaging data were averaged together and decomposed using independent component analysis, and three âneuronal networksâ were identified as demonstrating maximum temporal correlation with the neuronal stimulus paradigms: Default Mode Network, Task Positive Network, and Visual Network. For each of these, we observed a second network component with high spatial overlap. Using dual regression in the original 30 datasets, we extracted the time-series associated with these network pairs and calculated the percent of variance explained by the neuronal or vascular stimuli using a normalized R2 parameter. In each pairing, one network was dominated by the appropriate neuronal stimulus, and the other was dominated by the vascular stimulus as represented by the end-tidal CO2 time-series recorded in each scan. We acquired a second dataset in 8 of the original participants, where no CO2 challenge was delivered and CO2 levels fluctuated naturally with breathing variations. Although splitting of functional networks was not robust in these data, performing dual regression with the network maps from the original analysis in this new dataset successfully replicated our observations. Thus, in addition to responding to localized metabolic changes, the brainâs vasculature may be regulated in a coordinated manner that mimics (and potentially supports) specific functional brain networks. Multi-modal imaging and advances in fMRI acquisition and analysis could facilitate further study of the dual nature of functional brain networks. It will be critical to understand network-specific vascular function, and the behavior of a coupled vascular-neural network, in future studies of brain pathology
Multiparametric measurement of cerebral physiology using calibrated fMRI
The ultimate goal of calibrated fMRI is the quantitative imaging of oxygen metabolism (CMRO2), and this has been the focus of numerous methods and approaches. However, one underappreciated aspect of this quest is that in the drive to measure CMRO2, many other physiological parameters of interest are often acquired along the way. This can significantly increase the value of the dataset, providing greater information that is clinically relevant, or detail that can disambiguate the cause of signal variations. This can also be somewhat of a double-edged sword: calibrated fMRI experiments combine multiple parameters into a physiological model that requires multiple steps, thereby providing more opportunity for error propagation and increasing the noise and error of the final derived values. As with all measurements, there is a trade-off between imaging time, spatial resolution, coverage, and accuracy. In this review, we provide a brief overview of the benefits and pitfalls of extracting multiparametric measurements of cerebral physiology through calibrated fMRI experiments
Arterial CO2 fluctuations modulate neuronal rhythmicity: Implications for MEG and fMRI studies of resting-state networks
A fast emerging technique for studying human resting state networks (RSNs) is based on spontaneous temporal fluctuations in neuronal oscillatory power, as measured by magnetoencephalography. However, it has been demonstrated recently that this power is sensitive to modulations in arterial CO2 concentration. Arterial CO2 can be modulated by natural fluctuations in breathing pattern, as might typically occur during the acquisition of an RSN experiment. Here, we demonstrate for the first time the fine-scale dependence of neuronal oscillatory power on arterial CO2 concentration, showing that reductions in alpha, beta, and gamma power are observed with even very mild levels of hypercapnia (increased arterial CO2). We use a graded hypercapnia paradigm and participant feedback to rule out a sensory cause, suggesting a predominantly physiological origin. Furthermore, we demonstrate that natural fluctuations in arterial CO2, without administration of inspired CO2, are of a sufficient level to influence neuronal oscillatory power significantly in the delta-, alpha-, beta-, and gamma-frequency bands. A more thorough understanding of the relationship between physiological factors and cortical rhythmicity is required. In light of these findings, existing results, paradigms, and analysis techniques for the study of resting-state brain data should be revisited
Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow
Published: 15 June 2022Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health,
is commonly studied with the Blood Oxygenation Level Dependent functional MRI
(BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence
suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude
and may influence BOLD-CVR estimates. We address how acquisition and modeling
choices affect the relationship between baseline cerebral blood flow (bCBF) and
BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task,
and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We
assessed between-subject correlations of average GM values and within-subject spatial
correlations across cortical regions. Our results suggest that a breathing task addition to
a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can
improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely
because these CVR estimates are more physiologically accurate. We report positive
correlations between bCBF and BOLD-CVR, both between- and within-subjects. The
physiological explanation of this positive correlation is unclear; research with larger
samples and tightly controlled vasoactive stimuli is needed. Insights into what drives
variability in BOLD-CVR measurements and related measurements of cerebrovascular
function are particularly relevant when interpreting results in populations with altered
vascular and/or metabolic baselines or impaired cerebrovascular reserve.This work was supported by the Center for Translational Imaging
at Northwestern University. The authors disclosed receipt of
the following financial support for the research, authorship,
and/or publication of this article: This research was supported by
the Eunice Kennedy Shriver National Institute of Child Health
and Human Development of the National Institutes of Health
[K12HD073945]. KZ was supported by an NIH-funded training
program [T32EB025766]. SM was supported by the European
Unionâs Horizon 2020 research and innovation program [Marie
SkĆodowska-Curie grant agreement No. 713673] and a fellowship
from La Caixa Foundation [ID 100010434, fellowship code
LCF/BQ/IN17/11620063]. CC-G was supported by the Spanish
Ministry of Economy and Competitiveness [Ramon y Cajal
Fellowship, RYC2017-21845], the Basque Government [BERC
2018-2021 and PIBA_2019_104], and the Spanish Ministry
of Science, Innovation and Universities [MICINN; PID2019-
105520GB-100]
Changes in arterial cerebral blood volume during lower body negative pressure measured with MRI
Cerebral Autoregulation (CA), defined as the ability of the cerebral vasculature to maintain stable levels of blood flow despite changes in systemic blood pressure, is a critical factor in neurophysiological health. Magnetic resonance imaging (MRI) is a powerful technique for investigating cerebrovascular function, offering high spatial resolution and wide fields of view (FOV), yet it is relatively underutilized as a tool for assessment of CA. The aim of this study was to demonstrate the potential of using MRI to measure changes in cerebrovascular resistance in response to lower body negative pressure (LBNP). A Pulsed Arterial Spin Labeling (PASL) approach with short inversion times (TI) was used to estimate cerebral arterial blood volume (CBVa) in eight healthy subjects at baseline and -40mmHg LBNP. We estimated group mean CBVa values of 3.13±1.00 and 2.70±0.38 for baseline and lbnp respectively, which were the result of a differential change in CBVa during -40 mmHg LBNP that was dependent on baseline CBVa. These data suggest that the PASL CBVa estimates are sensitive to the complex cerebrovascular response that occurs during the moderate orthostatic challenge delivered by LBNP, which we speculatively propose may involve differential changes in vascular tone within different segments of the arterial vasculature. These novel data provide invaluable insight into the mechanisms that regulate perfusion of the brain, and establishes the use of MRI as a tool for studying CA in more detail
A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function
Available online 24 June 2021.Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO 2 pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemo- dynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural ac- tivity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting- state data segments, and in data segments which added a 2â3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO 2 pressure: a breath-hold task to induce hypercapnia (CO 2 increase) and a cued deep breathing task to induce hypocapnia (CO 2 decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO 2 by systematically shifting the CO 2 regressor in time to optimize the model fit. This optimization inher- ently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO 2 and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemody- namic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.This research was supported by the Eunice Kennedy Shriver Na- tional Institute of Child Health and Human Development of the Na- tional Institutes of Health under award number K12HD073945. The pediatric dataset and cerebral palsy dataset were collected with sup- port of National Institutes of Health award R03 HD094615â01A1. The authors would like to acknowledge Marie Wasielewski and Carson Ingo for their support in acquiring these data. K.Z. was supported by an NIH-funded training program (T32EB025766). S.M. was supported by the European Unionâs Horizon 2020 research and innovation pro- gram (Marie Sk Ć odowska-Curie grant agreement No. 713673), a fel- lowship from La Caixa Foundation (ID 100010434, fellowship code LCF/BQ/IN17/11620063) and C.C.G was supported by the Spanish Ministry of Economy and Competitiveness (Ramon y Cajal Fellowship, RYC-2017- 21845), the Basque Government (BERC 2018â2021 and PIBA_2019_104) and the Spanish Ministry of Science, Innovation and Universities (MICINN; PID2019â105520GB-100)
ICA-based denoising strategies in breath-hold induced cerebrovascular reactivity mapping with multi echo BOLD fMRI
Available online 6 March 2021.Performing a BOLD functional MRI (fMRI) acquisition during breath-hold (BH) tasks is a non-invasive, robust method to estimate cerebrovascular reactivity (CVR). However, movement and breathing-related artefacts caused by the BH can substantially hinder CVR estimates due to their high temporal collinearity with the effect of interest, and attention has to be paid when choosing which analysis model should be applied to the data. In this study, we evaluate the performance of multiple analysis strategies based on lagged general linear models applied on multi- echo BOLD fMRI data, acquired in ten subjects performing a BH task during ten sessions, to obtain subject-specific CVR and haemodynamic lag estimates. The evaluated approaches range from conventional regression models, i.e. including drifts and motion timecourses as nuisance regressors, applied on single-echo or optimally-combined data, to more complex models including regressors obtained from multi-echo independent component analysis with different grades of orthogonalization in order to preserve the effect of interest, i.e. the CVR. We compare these models in terms of their ability to make signal intensity changes independent from motion, as well as the reliability as measured by voxelwise intraclass correlation coefficients of both CVR and lag maps over time. Our results reveal that a conservative independent component analysis model applied on the optimally-combined multi-echo fMRI signal offers the largest reduction of motion-related effects in the signal, while yielding reliable CVR amplitude and lag estimates, although a conventional regression model applied on the optimally-combined data results in similar estimates. This work demonstrates the usefulness of multi-echo based fMRI acquisitions and independent component analysis denoising for precision mapping of CVR in single subjects based on BH paradigms, fostering its potential as a clinically-viable neuroimaging tool for individual patients. It also proves that the way in which data-driven regressors should be incorporated in the analysis model is not straight-forward due to their complex interaction with the BH-induced BOLD response.This research was supported by the European Unionâs Horizon 2020 research and innovation program ( Marie Sk Ć odowska-Curie grant agreement No. 713673 ), a fellowship from La Caixa Foundation (ID 100010434 , fellowship code LCF/BQ/IN17/11620063 ), the Spanish Ministry of Economy and Competitiveness ( Ramon y Cajal Fellowship, RYC-2017- 21845 ), the Spanish State Research Agency (BCBL âSevero Ochoa âexcellence accreditation, SEV- 2015-490 ), the Basque Govern- ment ( BERC 2018-2021 and PIBA_2019_104 ), the Spanish Ministry of Science, Innovation and Universities (MICINN; PID2019-105520GB-100 and FJCI-2017-31814 ), and the Eunice Kennedy Shriver National Insti- tute of Child Health and Human Development of the National Institutes of Health under award number K12HD073945