Axial variation of deoxyhemoglobin density as a source of the low-frequency time lag structure in blood oxygenation level-dependent signals

Brain-wide functional connectivity artifactually inflates throughout functional magnetic resonance imaging scans

Functional magnetic resonance imaging (fMRI) is a central tool for investigating human brain function, organization and disease. Here, we show that fMRI-based estimates of functional brain connectivity artifactually inflate at spatially heterogeneous rates during resting-state and task-based scans. This produces false positive connection strength changes and spatial distortion of brain connectivity maps. We demonstrate that this artefact is driven by temporal inflation of the non-neuronal, systemic low-frequency oscillation (sLFO) blood flow signal during fMRI scanning and is not addressed by standard denoising procedures. We provide evidence that sLFO inflation reflects perturbations in cerebral blood flow by respiration and heart rate changes that accompany diminishing arousal during scanning, although the mechanisms of this pathway are uncertain. Finally, we show that adding a specialized sLFO denoising procedure to fMRI processing pipelines mitigates the artifactual inflation of functional connectivity, enhancing the validity and within-scan reliability of fMRI findings.

Stronger influence of systemic than local hemodynamic-vascular factors on resting-state BOLD functional connectivity

Correlated fluctuations in the blood oxygenation level dependent (BOLD) signal of resting-state functional MRI (i.e., BOLD-functional connectivity, BOLD-FC) reflect a spectrum of neuronal and non-neuronal processes. In particular, there are multiple hemodynamic-vascular influences on BOLD-FC on both systemic (e.g., perfusion delay) and local levels (e.g., neurovascular coupling). While the influence of individual factors has been studied extensively, combined and comparative studies of systemic and local hemodynamic-vascular factors on BOLD-FC are scarce, notably in humans. We employed a multi-modal MRI approach to investigate and compare distinct hemodynamic-vascular processes and their impact on homotopic BOLD-FC in healthy controls and patients with unilateral asymptomatic internal carotid artery stenosis (ICAS). Asymptomatic ICAS is a cerebrovascular disorder, in which neuronal functioning is largely preserved but hemodynamic-vascular processes are impaired, mostly on the side of stenosis. Investigated indicators for local hemodynamic-vascular processes comprise capillary transit time heterogeneity (CTH) and cerebral blood volume (CBV) from dynamic susceptibility contrast (DSC) MRI, and cerebral blood flow (CBF) from pseudo-continuous arterial spin labeling (pCASL). Indicators for systemic processes are time-to-peak (TTP) from DSC MRI and BOLD lags from functional MRI. For each of these parameters, their influence on BOLD-FC was estimated by a comprehensive linear mixed model. Equally across groups, we found that individual mean BOLD-FC, local (CTH, CBV, and CBF) and systemic (TTP and BOLD lag) hemodynamic-vascular factors together explain 40.7% of BOLD-FC variance, with 20% of BOLD-FC variance explained by hemodynamic-vascular factors, with an about two-times larger contribution of systemic versus local factors. We conclude that regional differences in blood supply, i.e., systemic perfusion delays, exert a stronger influence on BOLD-FC than impairments in local neurovascular coupling.

Brain-wide functional connectivity artifactually inflates throughout fMRI scans: a problem and solution

The fMRI blood oxygen level-dependent (BOLD) signal is a mainstay of neuroimaging assessment of neuronal activity and functional connectivity in vivo. Thus, a chief priority is maximizing this signal's reliability and validity. To this end, the fMRI community has invested considerable effort into optimizing both experimental designs and physiological denoising procedures to improve the accuracy, across-scan reproducibility, and subject discriminability of BOLD-derived metrics like functional connectivity. Despite these advances, we discover that a substantial and ubiquitous defect remains in fMRI datasets: functional connectivity throughout the brain artifactually inflates during the course of fMRI scans - by an average of more than 70% in 15 minutes of scan time - at spatially heterogeneous rates, producing both spatial and temporal distortion of brain connectivity maps. We provide evidence that this inflation is driven by a previously unrecognized time-dependent increase of non-neuronal, systemic low-frequency oscillation (sLFO) blood flow signal during fMRI scanning. This signal is not removed by standard denoising procedures such as independent component analysis (ICA). However, we demonstrate that a specialized sLFO denoising procedure - Regressor Interpolation at Progressive Time Delays (RIPTiDe) - can be added to standard denoising pipelines to significantly attenuate functional connectivity inflation. We confirm the presence of sLFO-driven functional connectivity inflation in multiple independent fMRI datasets - including the Human Connectome Project - as well as across resting-state, task, and sleep-state conditions, and demonstrate its potential to produce false positive findings. Collectively, we present evidence for a previously unknown physiological phenomenon that spatiotemporally distorts estimates of brain connectivity in human fMRI datasets, and present a solution for mitigating this artifact.

Hemodynamic timing in resting-state and breathing-task BOLD fMRI

The blood flow response to a vasoactive stimulus demonstrates regional heterogeneity across both the healthy brain and in cerebrovascular pathology. The timing of a regional hemodynamic response is emerging as an important biomarker of cerebrovascular dysfunction, as well as a confound within fMRI analyses. Previous research demonstrated that hemodynamic timing is more robustly characterized when a larger systemic vascular response is evoked by a breathing challenge, compared to when only spontaneous fluctuations in vascular physiology are present (i.e., in resting-state data). However, it is not clear whether hemodynamic delays in these two conditions are physiologically interchangeable, and how methodological signal-to-noise factors may limit their agreement. To address this, we generated whole-brain maps of hemodynamic delays in nine healthy adults. We assessed the agreement of voxel-wise gray matter (GM) hemodynamic delays between two conditions: resting-state and breath-holding. We found that delay values demonstrated poor agreement when considering all GM voxels, but increasingly greater agreement when limiting analyses to voxels showing strong correlation with the GM mean time-series. Voxels showing the strongest agreement with the GM mean time-series were primarily located near large venous vessels, however these voxels explain some, but not all, of the observed agreement in timing. Increasing the degree of spatial smoothing of the fMRI data enhanced the correlation between individual voxel time-series and the GM mean time-series. These results suggest that signal-to-noise factors may be limiting the accuracy of voxel-wise timing estimates and hence their agreement between the two data segments. In conclusion, caution must be taken when using voxel-wise delay estimates from resting-state and breathing-task data interchangeably, and additional work is needed to evaluate their relative sensitivity and specificity to aspects of vascular physiology and pathology.

Reproducibility of cerebral perfusion measurements using BOLD delay

BOLD delay is an emerging, noninvasive method for assessing cerebral perfusion that does not require the use of intravenous contrast agents and is thus particularly suited for longitudinal monitoring. In this study, we assess the reproducibility of BOLD delay using data from 136 subjects with normal cerebral perfusion scanned on two separate occasions with scanners, sequence parameters, and intervals between scans varying between subjects. The effects of various factors on the reproducibility of BOLD delay, defined here as the differences in BOLD delay values between the scanning sessions, were investigated using a linear mixed model. Reproducibility was additionally assessed using the intraclass correlation coefficient of BOLD delay between sessions. Reproducibility was highest in the posterior cerebral artery territory. The mean BOLD delay test-retest difference after accounting for the aforementioned factors was 1.2 s (95% CI = 1.0 to 1.4 s). Overall, BOLD delay shows good reproducibility, but care should be taken when interpreting longitudinal BOLD delay changes that are either very small or are located in certain brain regions.

Human CSF movement influenced by vascular low frequency oscillations and respiration

Cerebrospinal fluid (CSF) movement through the pathways within the central nervous system is of high significance for maintaining normal brain health and function. Low frequency hemodynamics and respiration have been shown to drive CSF in humans independently. Here, we hypothesize that CSF movement may be driven simultaneously (and in synchrony) by both mechanisms and study their independent and coupled effects on CSF movement using novel neck fMRI scans. Caudad CSF movement at the fourth ventricle and hemodynamics of the major neck blood vessels (internal carotid arteries and internal jugular veins) was measured from 11 young, healthy volunteers using novel neck fMRI scans with simultaneous measurement of respiration. Two distinct models of CSF movement (1. Low-frequency hemodynamics and 2. Respiration) and possible coupling between them were investigated. We show that the dynamics of brain fluids can be assessed from the neck by studying the interrelationships between major neck blood vessels and the CSF movement in the fourth ventricle. We also demonstrate that there exists a cross-frequency coupling between these two separable mechanisms. The human CSF system can respond to multiple coupled physiological forces at the same time. This information may help inform the pathological mechanisms behind CSF movement-related disorders.

Quantifying cerebral blood arrival times using hypoxia-mediated arterial BOLD contrast

Cerebral blood arrival and tissue transit times are sensitive measures of the efficiency of tissue perfusion and can provide clinically meaningful information on collateral blood flow status. We exploit the arterial blood oxygen level dependent (BOLD) signal contrast established by precisely decreasing, and then increasing, arterial hemoglobin saturation using respiratory re-oxygenation challenges to quantify arterial blood arrival times throughout the brain. We term this approach the Step Hemoglobin re-Oxygenation Contrast Stimulus (SHOCS). Carpet plot analysis yielded measures of signal onset (blood arrival), global transit time (gTT) and calculations of relative total blood volume. Onset times averaged across 12 healthy subjects were 1.1 ± 0.4 and 1.9 ± 0.6 for cortical gray and deep white matter, respectively. The average whole brain gTT was 4.5 ± 0.9 seconds. The SHOCS response was 1.7 fold higher in grey versus white matter; in line with known differences in tissue-specific blood volume fraction. SHOCS was also applied in a patient with unilateral carotid artery occlusion revealing ipsilateral prolonged signal onset with normal perfusion in the unaffected hemisphere. We anticipate that SHOCS will further inform on the extent of collateral blood flow in patients with upstream steno-occlusive vascular disease, including those already known to manifest reductions in vasodilatory reserve capacity or vascular steal.

The relationship between attention and avoidance coping in anorexia nervosa: functional magnetic resonance imaging study

BACKGROUND

Numerous studies have demonstrated attentional control difficulties and high avoidance coping in patients with anorexia nervosa. Attention is a critical coping resource because it enables individuals to demonstrate self-control and complete goal-directed behaviours.

AIMS

We aimed to examine whether attentional control difficulty is related to high avoidance coping, and investigate the neural underpinnings of attentional control difficulties in individuals with anorexia nervosa.

METHOD

Twenty-three patients with anorexia nervosa and 17 healthy controls completed questionnaires that assessed attention and coping, and underwent functional magnetic resonance imaging while performing a go/no-go task.

RESULTS

Patients with anorexia nervosa showed weaker attentional control, higher omission error rates and higher avoidance coping compared with healthy controls. Attentional control difficulty was associated with higher avoidance coping in both groups. Functional magnetic resonance imaging analysis showed less deactivation in regions representing internal mental processing, such as the praecuneus, cuneus and left lingual gyrus, during the no-go condition. Moreover, weakened deactivation of the left lingual gyrus was associated with higher commission error rate in the anorexia nervosa group.

CONCLUSIONS

Our results suggest that patients with anorexia nervosa may have difficulty in maintaining attention to external ongoing events because of disturbance from internal self-related thought, and support the notion that attentional control difficulties underlie the frequent use of avoidance coping in anorexia nervosa.

The Role of Cerebrovascular-Reactivity Mapping in Functional MRI: Calibrated fMRI and Resting-State fMRI

Task and resting-state functional MRI (fMRI) is primarily based on the same blood-oxygenation level-dependent (BOLD) phenomenon that MRI-based cerebrovascular reactivity (CVR) mapping has most commonly relied upon. This technique is finding an ever-increasing role in neuroscience and clinical research as well as treatment planning. The estimation of CVR has unique applications in and associations with fMRI. In particular, CVR estimation is part of a family of techniques called calibrated BOLD fMRI, the purpose of which is to allow the mapping of cerebral oxidative metabolism (CMRO2) using a combination of BOLD and cerebral-blood flow (CBF) measurements. Moreover, CVR has recently been shown to be a major source of vascular bias in computing resting-state functional connectivity, in much the same way that it is used to neutralize the vascular contribution in calibrated fMRI. Furthermore, due to the obvious challenges in estimating CVR using gas challenges, a rapidly growing field of study is the estimation of CVR without any form of challenge, including the use of resting-state fMRI for that purpose. This review addresses all of these aspects in which CVR interacts with fMRI and the role of CVR in calibrated fMRI, provides an overview of the physiological biases and assumptions underlying hypercapnia-based CVR and calibrated fMRI, and provides a view into the future of non-invasive CVR measurement.

An MRI protocol for anatomical and functional evaluation of the California sea lion brain

BACKGROUND

Domoic acid (DOM) is a neurotoxin produced by some harmful algae blooms in coastal waters. California sea lions (Zalophus californianus) exposed to DOM often strand on beaches where they exhibit a variety of symptoms, including seizures. These animals typically show hippocampal atrophy on MRI scans.

NEW METHOD

We describe an MRI protocol for comprehensive evaluation of DOM toxicosis in the sea lion brain. We intend to study brain development in pups exposed in utero. The protocol depicts the hippocampal formation as the primary region of interest. We include scans for quantitative morphometry, functional and structural connectivity, and a cerebral blood flow map.

RESULTS

High-resolution 3D anatomical scans facilitate post hoc slicing in arbitrary planes and accurate morphometry. We demonstrate the first cerebral blood flow map using MRI, and the first structural tractography from a live sea lion brain.

COMPARISON WITH EXISTING METHODS

Scans were compared to prior anatomical and functional studies in live sea lions, and structural connectivity in post mortem specimens. Hippocampal volumes were broadly in line with prior studies, with differences likely attributable to the 3D approach used here. Functional connectivity of the dorsal left hippocampus matched that found in a prior study conducted at a lower magnetic field, while structural connectivity in the live brain agreed with findings observed in post mortem studies.

CONCLUSIONS

Our protocol provides a comprehensive, longitudinal view of the functional and anatomical changes expected to result from DOM toxicosis. It can also screen for other common neurological pathologies and is suitable for any pinniped that can fit inside an MRI scanner.

A venous mechanism of ventriculomegaly shared between traumatic brain injury and normal ageing

Recently, age-related timing dissociation between the superficial and deep venous systems has been observed; this was particularly pronounced in patients with normal pressure hydrocephalus, suggesting a common mechanism of ventriculomegaly. Establishing the relationship between venous drainage and ventricular enlargement would be clinically relevant and could provide insight into the mechanisms underlying brain ageing. To investigate a possible link between venous drainage and ventriculomegaly in both normal ageing and pathological conditions, we compared 225 healthy subjects (137 males and 88 females) and 71 traumatic brain injury patients of varying ages (53 males and 18 females) using MRI-based volumetry and a novel perfusion-timing analysis. Volumetry, focusing on the CSF space, revealed that the sulcal space and ventricular size presented different lifespan profiles with age; the latter presented a quadratic, rather than linear, pattern of increase. The venous timing shift slightly preceded this change, supporting a role for venous drainage in ventriculomegaly. In traumatic brain injury, a small but significant disease effect, similar to idiopathic normal pressure hydrocephalus, was found in venous timing, but it tended to decrease with age at injury, suggesting an overlapping mechanism with normal ageing. Structural bias due to, or a direct causative role of ventriculomegaly was unlikely to play a dominant role, because of the low correlation between venous timing and ventricular size after adjustment for age in both patients and controls. Since post-traumatic hydrocephalus can be asymptomatic and occasionally overlooked, the observation suggested a link between venous drainage and CSF accumulation. Thus, hydrocephalus, involving venous insufficiency, may be a part of normal ageing, can be detected non-invasively, and is potentially treatable. Further investigation into the clinical application of this new marker of venous function is therefore warranted.

Rude mechanicals in brain haemodynamics: non-neural actors that influence blood flow

Fluctuations in blood oxygenation and flow are widely used to infer brain activity during resting-state functional magnetic resonance imaging (fMRI). However, there are strong systemic and vascular contributions to resting-state signals that are unrelated to ongoing neural activity. Importantly, these non-neural contributions to haemodynamic signals (or 'rude mechanicals') can be as large as or larger than the neurally evoked components. Here, we review the two broad classes of drivers of these signals. One is systemic and is tied to fluctuations in external drivers such as heart rate and breathing, and the robust autoregulatory mechanisms that try to maintain a constant milieu in the brain. The other class comprises local, active fluctuations that appear to be intrinsic to vascular tissue and are likely similar to active local fluctuations seen in vasculature all over the body. In this review, we describe these non-neural fluctuations and some of the tools developed to correct for them when interpreting fMRI recordings. However, we also emphasize the links between these vascular fluctuations and brain physiology and point to ways in which fMRI measurements can be used to exploit such links to gain valuable information about neurovascular health and about internal brain states. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Timescales of spontaneous fMRI fluctuations relate to structural connectivity in the brain

Intrinsic timescales of activity fluctuations vary hierarchically across the brain. This variation reflects a broad gradient of functional specialization in information storage and processing, with integrative association areas displaying slower timescales that are thought to reflect longer temporal processing windows. The organization of timescales is associated with cognitive function, distinctive between individuals, and disrupted in disease, but we do not yet understand how the temporal properties of activity dynamics are shaped by the brain's underlying structural connectivity network. Using resting-state fMRI and diffusion MRI data from 100 healthy individuals from the Human Connectome Project, here we show that the timescale of resting-state fMRI dynamics increases with structural connectivity strength, matching recent results in the mouse brain. Our results hold at the level of individuals, are robust to parcellation schemes, and are conserved across a range of different timescale- related statistics. We establish a comprehensive BOLD dynamical signature of structural connectivity strength by comparing over 6,000 time series features, highlighting a range of new temporal features for characterizing BOLD dynamics, including measures of stationarity and symbolic motif frequencies. Our findings indicate a conserved property of mouse and human brain organization in which a brain region's spontaneous activity fluctuations are closely related to their surrounding structural scaffold.

The Effect of Scan Length on the Assessment of BOLD Delay in Ischemic Stroke

Objectives: To evaluate the impact of resting-state functional MRI scan length on the diagnostic accuracy, image quality and lesion volume estimation of BOLD delay maps used for brain perfusion assessment in acute ischemic stroke.

Methods: Sixty-three acute ischemic stroke patients received a 340 s resting-state functional MRI within 24 h of stroke symptom onset. BOLD delay maps were calculated from the full scan and four shortened versions (68 s, 136 s, 204 s, 272 s). The BOLD delay lesions on these maps were compared in terms of spatial overlap and volumetric agreement with the lesions derived from the full scans and with time-to-maximum (Tmax) lesions derived from DSC-MRI in a subset of patients (n = 10). In addition, the interpretability and quality of these maps were compared across different scan lengths using mixed models.

Results: Shortened BOLD delay scans showed a small volumetric bias (ranging from 0.05 to 5.3 mL; between a 0.13% volumetric underestimation and a 7.7% overestimation relative to the mean of the volumes, depending on scan length) compared to the full scan. Decreased scan length was associated with decreased spatial overlap with both the BOLD delay lesions derived from the full scans and with Tmax lesions. Only the two shortest scan lengths (68 and 136 s) were associated with substantially decreased interpretability, decreased structure clarity, and increased noisiness of BOLD delay maps.

Conclusions: BOLD delay maps derived from resting-state fMRI scans lasting 272 and 204 s provide sufficient diagnostic quality and adequate assessment of perfusion lesion volumes. Such shortened scans may be helpful in situations where quick clinical decisions need to be made.

Correction: Axial variation of deoxyhemoglobin density as a source of the low-frequency time lag structure in blood oxygenation level-dependent signals

[This corrects the article DOI: 10.1371/journal.pone.0222787.].