Lia Maria Hocke

Low-frequency oscillations measured in the periphery with near-infrared spectroscopy are strongly correlated with blood oxygen level-dependent functional magnetic resonance imaging signals

Abstract. Low-frequency oscillations (LFOs) in the range of 0.01–0.15 Hz are commonly observed in functional imaging studies, such as blood oxygen level-dependent functional magnetic resonance imaging (BOLD fMRI) and functional near-infrared spectroscopy (fNIRS). Some of these LFOs are nonneuronal and are closely related to autonomic physiological processes. In the current study, we conducted a concurrent resting-state fMRI and NIRS experiment with healthy volunteers. LFO data was collected simultaneously at peripheral sites (middle fingertip and big toes) by NIRS, and centrally in the brain by BOLD fMRI. The cross-correlations of the LFOs collected from the finger, toes, and brain were calculated. Our data show that the LFOs measured in the periphery (NIRS signals) and in the brain (BOLD fMRI) were strongly correlated with varying time delays. This demonstrates that some portion of the LFOs actually reflect systemic physiological circulatory effects. Furthermore, we demonstrated that NIRS is effective for measuring the peripheral LFOs, and that these LFOs and the temporal shifts between them are consistent in healthy participants and may serve as useful biomarkers for detecting and monitoring circulatory dysfunction.

Evaluating the effects of systemic low frequency oscillations measured in the periphery on the independent component analysis results of resting state networks

Independent component analysis (ICA) is widely used in resting state functional connectivity studies. ICA is a data-driven method, which uses no a priori anatomical or functional assumptions. However, as a result, it still relies on the user to distinguish the independent components (ICs) corresponding to neuronal activation, peripherally originating signals (without directly attributable neuronal origin, such as respiration, cardiac pulsation and Mayer wave), and acquisition artifacts. In this concurrent near infrared spectroscopy (NIRS)/functional MRI (fMRI) resting state study, we developed a method to systematically and quantitatively identify the ICs that show strong contributions from signals originating in the periphery. We applied group ICA (MELODIC from FSL) to the resting state data of 10 healthy participants. The systemic low frequency oscillation (LFO) detected simultaneously at each participants fingertip by NIRS was used as a regressor to correlate with every subject-specific IC time course. The ICs that had high correlation with the systemic LFO were those closely associated with previously described sensorimotor, visual, and auditory networks. The ICs associated with the default mode and frontoparietal networks were less affected by the peripheral signals. The consistency and reproducibility of the results were evaluated using bootstrapping. This result demonstrates that systemic, low frequency oscillations in hemodynamic properties overlay the time courses of many spatial patterns identified in ICA analyses, which complicates the detection and interpretation of connectivity in these regions of the brain.

Can apparent resting state connectivity arise from systemic fluctuations

It is widely accepted that the fluctuations in resting state blood oxygenation level dependent (BOLD) functional MRI (fMRI) reflect baseline neuronal activation through neurovascular coupling; this data is used to infer functional connectivity in the human brain during rest. Consistent activation patterns, i.e., resting state networks (RSN) are seen across groups, conditions, and even species. In this study, we show that some of these patterns can also be generated from the dynamic, systemic, non-neuronal physiological low frequency oscillations (sLFOs) in the BOLD signal alone. We have previously used multimodal imaging to demonstrate the wide presence of the same sLFOs in the brain (BOLD) and periphery with different time delays. This study shows that these sLFOs from BOLD signals alone can give rise to stable spatial patterns, which can be detected during resting state analyses. We generated synthetic resting state data for 11 subjects based only on subject-specific, dynamic sLFO information obtained from resting state data using concurrent peripheral optical imaging or a novel recursive procedure. We compared the results obtained by performing a group independent component analysis (ICA) on this synthetic data (i.e., the result from simulation) to the results obtained from analysis of the real data. ICA detected most of the eight well-known RSNs, including visual, motor, and default mode networks (DMNs), in both the real and the synthetic data sets. These findings suggest that RSNs may reflect, to some extent, vascular anatomy associated with systemic fluctuations, rather than neuronal connectivity.

Short repetition time multiband echo-planar imaging with simultaneous pulse recording allows dynamic imaging of the cardiac pulsation signal

Recently developed simultaneous multislice echo‐planar imaging (EPI) sequences permit imaging of the whole brain at short repetition time (TR), allowing the cardiac fluctuations to be fully sampled in blood‐oxygen‐level dependent functional MRI (BOLD fMRI). A novel low computational analytical method was developed to dynamically map the passage of the pulsation signal through the brain and visualize the whole cerebral vasculature affected by the pulse signal. This algorithm is based on a simple combination of fast BOLD fMRI and the scanners own built‐in pulse oximeter.

Perfusion information extracted from resting state functional magnetic resonance imaging

It is widely known that blood oxygenation level dependent (BOLD) contrast in functional magnetic resonance imaging (fMRI) is an indirect measure for neuronal activations through neurovascular coupling. The BOLD signal is also influenced by many non-neuronal physiological fluctuations. In previous resting state (RS) fMRI studies, we have identified a moving systemic low frequency oscillation (sLFO) in BOLD signal and were able to track its passage through the brain. We hypothesized that this seemingly intrinsic signal moves with the blood, and therefore, its dynamic patterns represent cerebral blood flow. In this study, we tested this hypothesis by performing Dynamic Susceptibility Contrast (DSC) MRI scans (i.e. bolus tracking) following the RS scans on eight healthy subjects. The dynamic patterns of sLFO derived from RS data were compared with the bolus flow visually and quantitatively. We found that the flow of sLFO derived from RS fMRI does to a large extent represent the blood flow measured with DSC. The small differences, we hypothesize, are largely due to the difference between the methods in their sensitivity to different vessel types. We conclude that the flow of sLFO in RS visualized by our time delay method represents the blood flow in the capillaries and veins in the brain.

Time delay processing of hypercapnic fMRI allows quantitative parameterization of cerebrovascular reactivity and blood flow delays

Blood oxygenation level-dependent fMRI contrast depends on the volume and oxygenation of blood flowing through the circulatory system. The effects on image intensity depend temporally on the arrival of blood within a voxel, and signal can be monitored during the time course of such blood flow. It has been previously shown that the passage of global endogenous variations in blood volume and oxygenation can be tracked as blood passes through the brain by determining the strength and peak time lag of their cross-correlation with blood oxygenation level-dependent data. By manipulating blood composition using transient hypercarbia and hyperoxia, we can induce much larger oxygenation and volume changes in the blood oxygenation level-dependent signal than result from natural endogenous fluctuations. This technique was used to examine cerebrovascular parameters in healthy subjects (n = 8) and subjects with intracranial stenosis (n = 22), with a subgroup of intracranial stenosis subjects scanned before and after surgical revascularization (n = 6). The halfwidth of cross-correlation lag times in the brain was larger in IC stenosis subjects (21.21 ± 14.22 s) than in healthy control subjects (8.03 ± 3.67), p < 0.001, and was subsequently reduced in regions that co-localized with surgical revascularization. These data show that blood circulatory timing can be measured robustly and longitudinally throughout the brain using simple respiratory challenges.

Isolating the sources of widespread physiological fluctuations in functional near-infrared spectroscopy signals

Physiological fluctuations at low frequency (<0.1 Hz) are prominent in functional near-infrared spectroscopy (fNIRS) measurements in both resting state and functional task studies. In this study, we used the high spatial resolution and full brain coverage of functional magnetic resonance imaging (fMRI) to understand the origins and commonalities of these fluctuations. Specifically, we applied a newly developed method, regressor interpolation at progressive time delays, to analyze concurrently recorded fNIRS and fMRI data acquired both in a resting state study and in a finger tapping study. The method calculates the voxelwise correlations between blood oxygen level dependent (BOLD) fMRI and fNIRS signals with different time shifts and localizes the areas in the brain that highly correlate with the fNIRS signal recorded at the surface of the head. The results show the wide spatial distribution of this physiological fluctuation in BOLD data, both in task and resting states. The brain areas that are highly correlated with global physiological fluctuations observed by fNIRS have a pattern that resembles the venous system of the brain, indicating the blood fluctuation from veins on the brain surface might strongly contribute to the overall fNIRS signal.

Correcting for Blood Arrival Time in Global Mean Regression Enhances Functional Connectivity Analysis of Resting State fMRI-BOLD Signals

Resting state functional connectivity analysis is a widely used method for mapping intrinsic functional organization of the brain. Global signal regression (GSR) is commonly employed for removing systemic global variance from resting state BOLD-fMRI data; however, recent studies have demonstrated that GSR may introduce spurious negative correlations within and between functional networks, calling into question the meaning of anticorrelations reported between some networks. In the present study, we propose that global signal from resting state fMRI is composed primarily of systemic low frequency oscillations (sLFOs) that propagate with cerebral blood circulation throughout the brain. We introduce a novel systemic noise removal strategy for resting state fMRI data, “dynamic global signal regression” (dGSR), which applies a voxel-specific optimal time delay to the global signal prior to regression from voxel-wise time series. We test our hypothesis on two functional systems that are suggested to be intrinsically organized into anticorrelated networks: the default mode network (DMN) and task positive network (TPN). We evaluate the efficacy of dGSR and compare its performance with the conventional “static” global regression (sGSR) method in terms of (i) explaining systemic variance in the data and (ii) enhancing specificity and sensitivity of functional connectivity measures. dGSR increases the amount of BOLD signal variance being modeled and removed relative to sGSR while reducing spurious negative correlations introduced in reference regions by sGSR, and attenuating inflated positive connectivity measures. We conclude that incorporating time delay information for sLFOs into global noise removal strategies is of crucial importance for optimal noise removal from resting state functional connectivity maps.

Comparison of peripheral near‐infrared spectroscopy low‐frequency oscillations to other denoising methods in resting state functional MRI with ultrahigh temporal resolution

Functional MRI (fMRI) blood–oxygen level–dependent (BOLD) signals result not only from neuronal activation, but also from nonneuronal physiological processes. These changes, especially in the low‐frequency domain (0.01–0.2 Hz), can significantly confound inferences about neuronal processes. It is crucial to effectively identify these nuisance low‐frequency oscillations (LFOs).

Systemic Low-Frequency Oscillations in BOLD Signal Vary with Tissue Type.

Blood-oxygen-level dependent (BOLD) signals are widely used in functional magnetic resonance imaging (fMRI) as a proxy measure of brain activation. However, because these signals are blood-related, they are also influenced by other physiological processes. This is especially true in resting state fMRI, during which no experimental stimulation occurs. Previous studies have found that the amplitude of resting state BOLD is closely related to regional vascular density. In this study, we investigated how some of the temporal fluctuations of the BOLD signal also possibly relate to regional vascular density. We began by identifying the blood-bound systemic low-frequency oscillation (sLFO). We then assessed the distribution of all voxels based on their correlations with this sLFO. We found that sLFO signals are widely present in resting state BOLD signals and that the proportion of these sLFOs in each voxel correlates with different tissue types, which vary significantly in underlying vascular density. These results deepen our understanding of the BOLD signal and suggest new imaging biomarkers based on fMRI data, such as amplitude of low-frequency fluctuation (ALFF) and sLFO, a combination of both, for assessing vascular density.