Spontaneous fluctuations in cerebral oxygen supply. An introduction

The impact of vasomotion on analysis of rodent fMRI data

Introduction

Small animal fMRI is an essential part of translational research in the cognitive neurosciences. Due to small dimensions and animal physiology preclinical fMRI is prone to artifacts that may lead to misinterpretation of the data. To reach unbiased translational conclusions, it is, therefore, crucial to identify potential sources of experimental noise and to develop correction methods for contributions that cannot be avoided such as physiological noise. Aim of this study was to assess origin and prevalence of hemodynamic oscillations (HDO) in preclinical fMRI in rat, as well as their impact on data analysis.

Methods

Following the development of algorithms for HDO detection and suppression, HDO prevalence in fMRI measurements was investigated for different anesthetic regimens, comprising isoflurane and medetomidine, and for both gradient echo and spin echo fMRI sequences. In addition to assessing the effect of vasodilation on HDO, it was studied if HDO have a direct neuronal correlate using local field potential (LFP) recordings. Finally, the impact of HDO on analysis of fMRI data was assessed, studying both the impact on calculation of activation maps as well as the impact on brain network analysis. Overall, 303 fMRI measurements and 32 LFP recordings were performed in 71 rats.

Results

In total, 62% of the fMRI measurements showed HDO with a frequency of (0.20 ± 0.02) Hz. This frequent occurrence indicated that HDO cannot be generally neglected in fMRI experiments. Using the developed algorithms, HDO were detected with a specificity of 95%, and removed efficiently from the signal time courses. HDO occurred brain-wide under vasoconstrictive conditions in both small and large blood vessels. Vasodilation immediately interrupted HDO, which, however, returned within 1 h under vasoconstrictive conditions. No direct neuronal correlate of HDO was observed in LFP recordings. HDO significantly impacted analysis of fMRI data, leading to altered cluster sizes and F-values for activated voxels, as well as altered brain networks, when comparing data with and without HDO.

Discussion

We therefore conclude that HDO are caused by vasomotion under certain anesthetic conditions and should be corrected during fMRI data analysis to avoid bias.

Function and development of interneurons involved in brain tissue oxygen regulation

The regulation of oxygen in brain tissue is one of the most important fundamental questions in neuroscience and medicine. The brain is a metabolically demanding organ, and its health directly depends on maintaining oxygen concentrations within a relatively narrow range that is both sufficiently high to prevent hypoxia, and low enough to restrict the overproduction of oxygen species. Neurovascular interactions, which are responsible for oxygen delivery, consist of neuronal and glial components. GABAergic interneurons play a particularly important role in neurovascular interactions. The involvement of interneurons extends beyond the perspective of inhibition, which prevents excessive neuronal activity and oxygen consumption, and includes direct modulation of the microvasculature depending upon their sub-type. Namely, nitric oxide synthase-expressing (NOS), vasoactive intestinal peptide-expressing (VIP), and somatostatin-expressing (SST) interneurons have shown modulatory effects on microvessels. VIP interneurons are known to elicit vasodilation, SST interneurons typically cause vasoconstriction, and NOS interneurons have to propensity to induce both effects. Given the importance and heterogeneity of interneurons in regulating local brain tissue oxygen concentrations, we review their differing functions and developmental trajectories. Importantly, VIP and SST interneurons display key developmental milestones in adolescence, while NOS interneurons mature much earlier. The implications of these findings point to different periods of critical development of the interneuron-mediated oxygen regulatory systems. Such that interference with normal maturation processes early in development may effect NOS interneuron neurovascular interactions to a greater degree, while insults later in development may be more targeted toward VIP- and SST-mediated mechanisms of oxygen regulation.

Brain tissue oxygen dynamics while mimicking the functional deficiency of interneurons

The dynamic interaction between excitatory and inhibitory activity in the brain is known as excitatory-inhibitory balance (EIB). A significant shift in EIB toward excitation has been observed in numerous pathological states and diseases, such as autism or epilepsy, where interneurons may be dysfunctional. The consequences of this on neurovascular interactions remains to be elucidated. Specifically, it is not known if there is an elevated metabolic consumption of oxygen due to increased excitatory activity. To investigate this, we administered microinjections of picrotoxin, a gamma aminobutyric acid (GABA) antagonist, to the rabbit cortex in the awake state to mimic the functional deficiency of GABAergic interneurons. This caused an observable shift in EIB toward excitation without the induction of seizures. We used chronically implanted electrodes to measure both neuronal activity and brain tissue oxygen concentrations (PO₂) simultaneously and in the same location. Using a high-frequency recording rate for PO₂, we were able to detect two important phenomena, (1) the shift in EIB led to a change in the power spectra of PO₂ fluctuations, such that higher frequencies (8-15 cycles per minute) were suppressed and (2) there were brief periods (dips with a duration of less than 100 ms associated with neuronal bursts) when PO₂ dropped below 10 mmHg, which we defined as the threshold for hypoxia. The dips were followed by an overshoot, which indicates either a rapid vascular response or decrease in oxygen consumption. Our results point to the essential role of interneurons in brain tissue oxygen regulation in the resting state.

Origins of 1/f-like tissue oxygenation fluctuations in the murine cortex

The concentration of oxygen in the brain spontaneously fluctuates, and the distribution of power in these fluctuations has a 1/f-like spectra, where the power present at low frequencies of the power spectrum is orders of magnitude higher than at higher frequencies. Though these oscillations have been interpreted as being driven by neural activity, the origin of these 1/f-like oscillations is not well understood. Here, to gain insight of the origin of the 1/f-like oxygen fluctuations, we investigated the dynamics of tissue oxygenation and neural activity in awake behaving mice. We found that oxygen signal recorded from the cortex of mice had 1/f-like spectra. However, band-limited power in the local field potential did not show corresponding 1/f-like fluctuations. When local neural activity was suppressed, the 1/f-like fluctuations in oxygen concentration persisted. Two-photon measurements of erythrocyte spacing fluctuations and mathematical modeling show that stochastic fluctuations in erythrocyte flow could underlie 1/f-like dynamics in oxygenation. These results suggest that the discrete nature of erythrocytes and their irregular flow, rather than fluctuations in neural activity, could drive 1/f-like fluctuations in tissue oxygenation.

Behavior and Regional Cortical BOLD Signal Fluctuations Are Altered in Adult Rabbits After Neonatal Volatile Anesthetic Exposure

Neonatal and infant exposure to volatile anesthetics has been associated with long-term learning, memory, and behavioral deficits. Although early anesthesia exposure has been linked to a number of underlying structural abnormalities, functional changes associated with these impairments remain poorly understood. To investigate the relationship between functional alteration in neuronal circuits and learning deficiency, resting state functional MRI (rsfMRI) connectivity was examined in adolescent rabbits exposed to general anesthesia as neonates (1 MAC isoflurane for 2 h on postnatal days P8, P11, and P14) and unanesthetized controls before and after training with a trace eyeblink classical conditioning (ECC) paradigm. Long-range connectivity was measured between several key regions of interest (ROIs), including primary and secondary somatosensory cortices, thalamus, hippocampus, and cingulate. In addition, metrics of regional BOLD fluctuation amplitudes and coherence, amplitude of low-frequency fluctuation (ALFF), fractional ALFF (fALFF), and regional homogeneity (ReHo) were calculated. Our results showed that the trace ECC learning rate was significantly lower in the anesthesia-exposed group. No anesthesia-related changes in long-range connectivity, fALFF, or ReHo were found between any ROIs. However, ALFF was significantly higher in anesthesia-exposed rabbits in the primary and secondary somatosensory cortices, and ALFF in those areas was a significant predictor of the learning performance for trace ECC. The absence of anesthesia-related changes in long-range thalamocortical connectivity indicates that functional thalamocortical input is not affected. Higher ALFF in the somatosensory cortex may indicate the developmental disruption of cortical neuronal circuits after neonatal anesthesia exposure, including excessive neuronal synchronization that may underlie the observed cognitive deficits.

The glymphatic system and its role in cerebral homeostasis

The brain's high bioenergetic state is paralleled by high metabolic waste production. Authentic lymphatic vasculature is lacking in brain parenchyma. Cerebrospinal fluid (CSF) flow has long been thought to facilitate central nervous system detoxification in place of lymphatics, but the exact processes involved in toxic waste clearance from the brain remain incompletely understood. Over the past 8 yr, novel data in animals and humans have begun to shed new light on these processes in the form of the "glymphatic system," a brain-wide perivascular transit passageway dedicated to CSF transport and interstitial fluid exchange that facilitates metabolic waste drainage from the brain. Here we will discuss glymphatic system anatomy and methods to visualize and quantify glymphatic system (GS) transport in the brain and also discuss physiological drivers of its function in normal brain and in neurodegeneration.

Oscillations and concentration dynamics of brain tissue oxygen in neonates and adults

The brain is a metabolically demanding organ and its health directly depends on brain oxygen dynamics to prevent hypoxia and ischemia. Localized brain tissue oxygen is characterized by a baseline level combined with spontaneous oscillations. These oscillations are attributed to spontaneous changes of vascular tone at the level of arterioles and their frequencies depend on age. Specifically, lower frequencies are more typical for neonates than for adults. We have built a mathematical model which analyses the diffusion abilities of oxygen based on the frequency of source brain oxygen oscillations and neuronal demand. We have found that a lower frequency of spontaneous oscillations of localized brain tissue oxygen can support higher amplitudes of oxygen concentration at areas distant from a source relative to oscillations at higher frequencies. Since hypoxia and ischemia are very common events during early development and the neurovascular unit is underdeveloped in neonates, our results indicate that lower frequency oxygen oscillations can represent an effective passive method of neonatal brain protection against hypoxia. These results can have a potential impact on future studies aiming to find new treatment strategies for brain ischemia.

Origins of the Resting-State Functional MRI Signal: Potential Limitations of the "Neurocentric" Model

Resting-state functional connectivity (rsFC) is emerging as a research tool for systems and clinical neuroscience. The mechanism underlying resting-state functional MRI (rsfMRI) signal, however, remains incompletely understood. A widely held assumption is that the spontaneous fluctuations in blood oxygenation level-dependent (BOLD) signal reflect ongoing neuronal processes (herein called “neurocentric” model). In support of this model, evidence from human and animal studies collectively reveals that the spatial synchrony of spontaneously occurring electrophysiological signal recapitulates BOLD rsFC networks. Two recent experiments from independent labs designed to specifically examine neuronal origins of rsFC, however, suggest that spontaneously occurring neuronal events, as assessed by multiunit activity or local field potential (LFP), although statistically significant, explain only a small portion (∼10%) of variance in resting-state BOLD fluctuations. These two studies, although each with its own limitations, suggest that the spontaneous fluctuations in rsfMRI, may have complex cellular origins, and the “neurocentric” model may not apply to all brain regions.

Detailed view on slow sinusoidal, hemodynamic oscillations on the human brain cortex by Fourier transforming oxy/deoxy hyperspectral images

Slow sinusoidal, hemodynamic oscillations (SSHOs) around 0.1 Hz are frequently seen in mammalian and human brains. In four patients undergoing epilepsy surgery, subtle but robust fluctuations in oxy‐ and deoxyhemoglobin were detected using hyperspectral imaging of the cortex. These SSHOs were stationary during the entire 4 to 10 min acquisition time. By Fourier filtering the oxy‐ and deoxyhemoglobin time signals with a small bandwidth, SSHOs became visible within localized regions of the brain, with distinctive frequencies and a continuous phase variation within that region. SSHOs of deoxyhemoglobin appeared to have an opposite phase and 11% smaller amplitude with respect to the oxyhemoglobin SSHOs. Although the origin of SSHOs remains unclear, we find indications that the observed SSHOs may embody a local propagating hemodynamic wave with velocities in line with capillary blood velocities, and conceivably related to vasomotion and maintenance of adequate tissue perfusion. Hyperspectral imaging of the human cortex during surgery allow in‐depth characterization of SSHOs, and may give further insight in the nature and potential (clinical) use of SSHOs.

Contribution of systemic vascular effects to fMRI activity in white matter

To investigate a potential contribution of systemic physiology to recently reported BOLD fMRI signals in white matter, we compared photo-plethysmography (PPG) and whole-brain fMRI signals recorded simultaneously during long resting-state scans from an overnight sleep study. We found that intermittent drops in the amplitude of the PPG signal exhibited strong and widespread correlations with the fMRI signal, both in white matter (WM) and in gray matter (GM). The WM signal pattern resembled that seen in previous resting-state fMRI studies and closely tracked the location of medullary veins. Its temporal cross-correlation with the PPG amplitude was bipolar, with an early negative value. In GM, the correlation was consistently positive. Consistent with previous studies comparing physiological signals with fMRI, these findings point to a systemic vascular contribution to WM fMRI signals. The PPG drops are interpreted as systemic vasoconstrictive events, possibly related to intermittent increases in sympathetic tone related to fluctuations in arousal state. The counter-intuitive polarity of the WM signal is explained by long blood transit times in the medullary vasculature of WM, which cause blood oxygenation loss and a substantial timing mismatch between blood volume and blood oxygenation effects. A similar mechanism may explain previous findings of negative WM signals around large draining veins during both task- and resting-state fMRI.

Brain tissue oxygen regulation in awake and anesthetized neonates

Inhaled general anesthetics are used commonly in adults and children, and a growing body of literature from animals and humans suggests that exposure to anesthesia at an early age can impact brain development. While the origin of these effects is not well understood, it is known that anesthesia can disrupt oxygen regulation in the brain, which is critically important for maintaining healthy brain function. Here we investigated how anesthesia affected brain tissue oxygen regulation in neonatal rabbits by comparing brain tissue oxygen and single unit activity in the awake and anesthetized states. We tested two common general anesthetics, isoflurane and sevoflurane, delivered in both air and 80% oxygen. Our findings show that general anesthetics can greatly increase brain tissue PO₂ in neonates, especially when combined with supplemental oxygen. Although isoflurane and sevoflurane belong to the same class of anesthetics, notable differences were observed in their effects upon neuronal activity and spontaneous respiration. Our findings point to the need to consider the potential effects of hyperoxia when supplemental oxygen is utilized, particularly in children and neonates.

Co-activation patterns in resting-state fMRI signals

The brain is a complex system that integrates and processes information across multiple time scales by dynamically coordinating activities over brain regions and circuits. Correlations in resting-state functional magnetic resonance imaging (rsfMRI) signals have been widely used to infer functional connectivity of the brain, providing a metric of functional associations that reflects a temporal average over an entire scan (typically several minutes or longer). Not until recently was the study of dynamic brain interactions at much shorter time scales (seconds to minutes) considered for inference of functional connectivity. One method proposed for this objective seeks to identify and extract recurring co-activation patterns (CAPs) that represent instantaneous brain configurations at single time points. Here, we review the development and recent advancement of CAP methodology and other closely related approaches, as well as their applications and associated findings. We also discuss the potential neural origins and behavioral relevance of CAPs, along with methodological issues and future research directions in the analysis of fMRI co-activation patterns.

Effects of Glutamate Receptor Activation on Local Oxygen Changes

Glutamate is ubiquitous throughout the brain and serves as the primary excitatory neurotransmitter. Neurons require energy to fire, and energetic substrates (i.e., O₂, glucose) are renewed via cerebral blood flow (CBF) to maintain metabolic homeostasis. Magnetic resonance brain functionality studies rely on the assumption that CBF and neuronal activity are coupled consistently throughout the brain; however, the origin of neuronal activity does not always coincide with signals indicative of energy consumption (e.g., O₂ decreases) at high spatial resolutions. Therefore, relationships between excitatory neurotransmission and energy use must be evaluated at higher resolutions. In this study, we showed that both endogenously released and exogenously ejected glutamate decrease local tissue O₂ concentrations, but whether hyperemic O₂ restoration followed depended on the stimulus method. Electrically stimulating the glutamatergic corticostriatal pathway evoked biphasic O₂ responses at striatal terminals: first O₂ decreased, then concentrations increased above baseline. Using iontophoresis to locally eject ionotropic glutamate receptor antagonists revealed that these receptors only influenced the O₂ decrease. We compared electrical stimulation to iontophoretic glutamate stimulation, and measured concurrent single-unit activity and O₂ to limit both stimulation and recordings to <50 μm radius from our sensor. Similarly, iontophoretic glutamate delivery elicited monophasic O₂ decreases without subsequent increases.

Retinal oxygen: from animals to humans

This article discusses retinal oxygenation and retinal metabolism by focusing on measurements made with two of the principal methods used to study O₂ in the retina: measurements of PO₂ with oxygen-sensitive microelectrodes in vivo in animals with a retinal circulation similar to that of humans, and oximetry, which can be used non-invasively in both animals and humans to measure O₂ concentration in retinal vessels. Microelectrodes uniquely have high spatial resolution, allowing the mapping of PO₂ in detail, and when combined with mathematical models of diffusion and consumption, they provide information about retinal metabolism. Mathematical models, grounded in experiments, can also be used to simulate situations that are not amenable to experimental study. New methods of oximetry, particularly photoacoustic ophthalmoscopy and visible light optical coherence tomography, provide depth-resolved methods that can separate signals from blood vessels and surrounding tissues, and can be combined with blood flow measures to determine metabolic rate. We discuss the effects on retinal oxygenation of illumination, hypoxia and hyperoxia, and describe retinal oxygenation in diabetes, retinal detachment, arterial occlusion, and macular degeneration. We explain how the metabolic measurements obtained from microelectrodes and imaging are different, and how they need to be brought together in the future. Finally, we argue for revisiting the clinical use of hyperoxia in ophthalmology, particularly in retinal arterial occlusions and retinal detachment, based on animal research and diffusion theory.

Noise and non-neuronal contributions to the BOLD signal: applications to and insights from animal studies

The BOLD signal reflects hemodynamic events within the brain, which in turn are driven by metabolic changes and neural activity. However, the link between BOLD changes and neural activity is indirect and can be influenced by a number of non-neuronal processes. Motion and physiological cycles have long been known to affect the BOLD signal and are present in both humans and animal models. Differences in physiological baseline can also contribute to intra- and inter-subject variability. The use of anesthesia, common in animal studies, alters neural activity, vascular tone, and neurovascular coupling. Most intriguing, perhaps, are the contributions from other processes that do not appear to be neural in origin but which may provide information about other aspects of neurophysiology. This review discusses different types of noise and non-neuronal contributors to the BOLD signal, sources of variability for animal studies, and insights to be gained from animal models.

Spontaneous Fluctuations of PO2 in the Rabbit Somatosensory Cortex

In many tissues, PO2 fluctuates spontaneously with amplitudes of a few mmHg. Here we further characterized these oscillations. PO2 recordings were made from the whisker barrel cortex of six rabbits with acutely or chronically placed polarographic electrodes. Measurements were made while rabbits were awake and while anesthetized with isoflurane, during air breathing, and during 100% oxygen inspiration. In awake rabbits, 90% of the power was between 0 and 20 cycles per minute (cpm), not uniformly distributed over this range, but with a peak frequently near 10 cpm. This was much slower than heart or respiratory rhythms and is similar to the frequency content observed in other tissues. During hyperoxia, total power was higher than during air-breathing, and the dominant frequencies tended to shift toward lower values (0-10 cpm). These observations suggest that at least the lower frequency fluctuations represent efforts by the circulation to regulate local PO2. There were no consistent changes in total power during 0.5 or 1.5% isoflurane anesthesia, but the power shifted to lower frequencies. Thus, both hyperoxia and anesthesia cause characteristic, but distinct, changes in spontaneous fluctuations. These PO2 fluctuations may be caused by vasomotion, but other factors cannot be ruled out.

Time-frequency dynamics of the sum of intra- and extracerebral hemodynamic functional connectivity during resting-state and respiratory challenges assessed by multimodal functional near-infrared spectroscopy

Monitoring respiratory processes is important for evaluating neuroimaging data, given their influence on time-frequency dynamics of intra- and extracerebral hemodynamics. Here we investigated the time-frequency dynamics of the sum of intra- and extracerebral hemodynamic functional connectivity states during hypo- and hypercapnia by using three different respiratory challenge tasks (i.e., hyperventilation, breath-holding, and rebreathing) compared to resting-state. The sum of intra- and extracerebral hemodynamic responses were assessed using functional near-infrared spectroscopy (fNIRS) within two regions of interest (i.e., the dorsolateral and the medial prefrontal cortex). Time-frequency fNIRS analysis was performed based on wavelet transform coherence to quantify functional connectivity in terms of positive and negative phase-coupling within each region of interest. Physiological measures were assessed in the form of partial end-tidal carbon dioxide, heart rate, arterial tissue oxygen saturation, and respiration rate. We found that the three respiration challenges modulated time-frequency dynamics differently with respect to resting-state: 1) Hyperventilation and breath-holding exhibited inverse patterns of positive and negative phase-coupling. 2) In contrast, rebreathing had no significant effect. 3) Low-frequency oscillations contributed to a greater extent to time-frequency dynamics compared to high-frequency oscillations. The results highlight that there exist distinct differences in time-frequency dynamics of the sum of intra- and extracerebral functional connectivity not only between hypo- (hyperventilation) and hypercapnia but also between different states of hypercapnia (breath-holding versus rebreathing). This suggests that a multimodal assessment of intra-/extracerebral and systemic physiological changes during respiratory challenges compared to resting-state may have potential use in the differentiation between physiological and pathological respiratory behavior accompanied by the psycho-physiological state of a human.

Frequency-specific network topologies in the resting human brain

A community is a set of nodes with dense inter-connections, while there are sparse connections between different communities. A hub is a highly connected node with high centrality. It has been shown that both “communities” and “hubs” exist simultaneously in the brain's functional connectivity network (FCN), as estimated by correlations among low-frequency spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signal changes (0.01–0.10 Hz). This indicates that the brain has a spatial organization that promotes both segregation and integration of information. Here, we demonstrate that frequency-specific network topologies that characterize segregation and integration also exist within this frequency range. In investigating the coherence spectrum among 87 brain regions, we found that two frequency bands, 0.01–0.03 Hz (very low frequency [VLF] band) and 0.07–0.09 Hz (low frequency [LF] band), mainly contributed to functional connectivity. Comparing graph theoretical indices for the VLF and LF bands revealed that the network in the former had a higher capacity for information segregation between identified communities than the latter. Hubs in the VLF band were mainly located within the anterior cingulate cortices, whereas those in the LF band were located in the posterior cingulate cortices and thalamus. Thus, depending on the timescale of brain activity, at least two distinct network topologies contributed to information segregation and integration. This suggests that the brain intrinsically has timescale-dependent functional organizations.

Early evaluation of nerve regeneration after nerve injury and repair using functional connectivity MRI

Resting state functional connectivity magnetic resonance imaging studies in rat brain show brain reorganization caused by nerve injury and repair. In this study, distinguishable differences were found in healthy, nerve transection without repair (R-) and nerve transection with repair (R+) groups in the subacute stage (2 weeks after initial injury). Only forepaw on the healthy side was used to determine seed voxel regions in this study. Disturbance of neuronal network in the primary sensory region of cortex occurs within two hours after initial injury, and the network pattern was restored in R+ group in subacute stage, while the disturbed pattern remained in R- group. These are the central findings of the study. This technique provides a novel way of detecting and monitoring the effectiveness of peripheral nerve injury treatment in the early stage and potentially offers a tool for clinicians to avoid poor clinical outcomes.

Detecting resting-state brain activity by spontaneous cerebral blood volume fluctuations using whole brain vascular space occupancy imaging

Resting-state brain activity has been investigated extensively using BOLD contrast. However, BOLD signal represents the combined effects of multiple physiological processes and its spatial localization is less accurate than that of cerebral blood flow and volume (CBF and CBF, respectively). In this study, we demonstrate that resting-state brain activity can be reliably detected by spontaneous fluctuations of CBV-weighted signal using whole-brain gradient and spin echo (GRASE) based vascular space occupancy (VASO) imaging. Specifically, using independent component analysis, intrinsic brain networks, including default mode, salience, executive control, visual, auditory, and sensorimotor networks were revealed robustly by the VASO technique. We further demonstrate that task-evoked VASO signal aligned well with expected gray matter areas, while blood-oxygenation level dependent (BOLD) signal extended outside of these areas probably due to their different spatial specificity. The improved spatial localization of VASO is consistent with previous studies using animal models. Moreover, we showed that the 3D-GRASE VASO images had reduced susceptibility-induced signal voiding, compared to the BOLD technique. This is attributed to the fact that VASO does not require T2* weighting, thus the acquisition can use a shorter TE and can employ spin-echo scheme. Consequently VASO-based functional connectivity signals were well preserved in brain regions that tend to suffer from signal loss and geometric distortion in BOLD, such as orbital prefrontal cortex. Our study suggests that 3D-GRASE VASO imaging, with its improved spatial specificity and less sensitivity to susceptibility artifacts, may have advantages in resting-state fMRI studies.

Resting-state fMRI confounds and cleanup

The goal of resting-state functional magnetic resonance imaging (fMRI) is to investigate the brain's functional connections by using the temporal similarity between blood oxygenation level dependent (BOLD) signals in different regions of the brain "at rest" as an indicator of synchronous neural activity. Since this measure relies on the temporal correlation of fMRI signal changes between different parts of the brain, any non-neural activity-related process that affects the signals will influence the measure of functional connectivity, yielding spurious results. To understand the sources of these resting-state fMRI confounds, this article describes the origins of the BOLD signal in terms of MR physics and cerebral physiology. Potential confounds arising from motion, cardiac and respiratory cycles, arterial CO₂ concentration, blood pressure/cerebral autoregulation, and vasomotion are discussed. Two classes of techniques to remove confounds from resting-state BOLD time series are reviewed: 1) those utilising external recordings of physiology and 2) data-based cleanup methods that only use the resting-state fMRI data itself. Further methods that remove noise from functional connectivity measures at a group level are also discussed. For successful interpretation of resting-state fMRI comparisons and results, noise cleanup is an often over-looked but essential step in the analysis pipeline.

The effects of respiratory CO2 fluctuations in the resting-state BOLD signal differ between eyes open and eyes closed

Resting fluctuations in arterial CO2 (a cerebral vasodilator) are believed to be an important source of low-frequency blood oxygenation level dependent (BOLD) signal fluctuations. In this study we focus on the two commonly used resting-states in functional magnetic resonance imaging experiments, eyes open and eyes closed, and quantify the degree to which measured spontaneous fluctuations in the partial pressure of end-tidal CO2 (Petco2) relate to BOLD signal time series. A significantly longer latency of BOLD signal changes following Petco2 fluctuations was found in the eyes closed condition compared to with eyes open, which may reveal different intrinsic vascular response delays in CO2 reactivity or an alteration in the net BOLD signal arising from Petco2 fluctuations and altered neural activity with eyes closed. By allowing a spatially varying time delay for the compensation of this temporal difference, a more spatially consistent CO2 correlation map can be obtained. Finally, Granger-causality analysis demonstrated a "causal" relationship between Petco2 and BOLD. The identified dominant Petco2→BOLD directional coupling supports the notion that Petco2 fluctuations are indeed a cause of resting BOLD variance in the majority of subjects.

The role of physiological noise in resting-state functional connectivity

Functional connectivity between different brain regions can be estimated from MRI data by computing the temporal correlation of low frequency (<0.1Hz) fluctuations in the MRI signal. These correlated fluctuations occur even when the subject is "at rest" (not asked to perform any particular task) and result from spontaneous neuronal activity synchronized within multiple distinct networks of brain regions. This estimate of connectivity, however, can be influenced by physiological noise, such as cardiac and respiratory fluctuations. This brief review looks at the effect of physiological noise on estimates of resting-state functional connectivity, discusses ways to remove physiological noise, and provides a personal recollection of the early developments in these approaches. This review also discusses the importance of physiological noise correction and provides a summary of evidence demonstrating that functional connectivity does have a neuronal underpinning and cannot purely be the result of physiological noise.

Rapid assessment of cerebral autoregulation by near-infrared spectroscopy and a single dose of phenylephrine

Rapid bedside determination of cerebral blood pressure autoregulation (AR) may improve clinical utility. We tested the hypothesis that cerebral Hb oxygenation (HbDiff) and cerebral Hb volume (HbTotal) measured by near-infrared spectroscopy (NIRS) would correlate with cerebral blood flow (CBF) after single dose phenylephrine (PE). Critically ill patients requiring artificial ventilation and arterial lines were eligible. During rapid blood pressure rise induced by i.v. PE bolus, ΔHbDiff and ΔHbTotal were calculated by subtracting values at baseline (normotension) from values at peak blood pressure elevation (hypertension). With the aid of NIRS and bolus injection of indocyanine green, relative measures of CBF, called blood flow index (BFI), were determined during normotension and during hypertension. BFI during hypertension was expressed as percentage from BFI during normotension (BFI%). Autoregulation indices (ARIs) were calculated by dividing BFI%, ΔHbDiff, and ΔHbTotal by the concomitant change in blood pressure. In 24 patients (11 newborns and 13 children), significant correlations between BFI% and ΔHbDiff (or ΔHbTotal) were found. In addition, the associations between Hb-based ARI and BFI%-based ARI were significant with correlation coefficients of 0.73 (or 0.72). Rapid determination of dynamic AR with the aid of cerebral Hb signals and PE bolus seems to be reliable.

Hypersensitivity to thromboxane receptor mediated cerebral vasomotion and CBF oscillations during acute NO-deficiency in rats

BACKGROUND

Low frequency (4-12 cpm) spontaneous fluctuations of the cerebrovascular tone (vasomotion) and oscillations of the cerebral blood flow (CBF) have been reported in diseases associated with endothelial dysfunction. Since endothelium-derived nitric oxide (NO) suppresses constitutively the release and vascular effects of thromboxane A(2) (TXA(2)), NO-deficiency is often associated with activation of thromboxane receptors (TP). In the present study we hypothesized that in the absence of NO, overactivation of the TP-receptor mediated cerebrovascular signaling pathway contributes to the development of vasomotion and CBF oscillations.

METHODOLOGY/PRINCIPAL FINDINGS

Effects of pharmacological modulation of TP-receptor activation and its downstream signaling pathway have been investigated on CBF oscillations (measured by laser-Doppler flowmetry in anesthetized rats) and vasomotion (measured by isometric tension recording in isolated rat middle cerebral arteries, MCAs) both under physiological conditions and after acute inhibition of NO synthesis. Administration of the TP-receptor agonist U-46619 (1 µg/kg i.v.) to control animals failed to induce any changes of the systemic or cerebral circulatory parameters. Inhibition of the NO synthesis by nitro-L-arginine methyl ester (L-NAME, 100 mg/kg i.v.) resulted in increased mean arterial blood pressure and a decreased CBF accompanied by appearance of CBF-oscillations with a dominant frequency of 148±2 mHz. U-46619 significantly augmented the CBF-oscillations induced by L-NAME while inhibition of endogenous TXA(2) synthesis by ozagrel (10 mg/kg i.v.) attenuated it. In isolated MCAs U-46619 in a concentration of 100 nM, which induced weak and stable contraction under physiological conditions, evoked sustained vasomotion in the absence of NO, which effect could be completely reversed by inhibition of Rho-kinase by 10 µM Y-27632.

CONCLUSION/SIGNIFICANCE

These results suggest that hypersensitivity of the TP-receptor-Rho-kinase signaling pathway contributes to the development of low frequency cerebral vasomotion which may propagate to vasospasm in pathophysiological states associated with NO-deficiency.

Age-Related Differences in Functional Nodes of the Brain Cortex - A High Model Order Group ICA Study

Functional MRI measured with blood oxygen dependent (BOLD) contrast in the absence of intermittent tasks reflects spontaneous activity of so-called resting state networks (RSN) of the brain. Group level independent component analysis (ICA) of BOLD data can separate the human brain cortex into 42 independent RSNs. In this study we evaluated age-related effects from primary motor and sensory, and, higher level control RSNs. One hundred sixty-eight healthy subjects were scanned and divided into three groups: 55 adolescents (ADO, 13.2 ± 2.4 years), 59 young adults (YA, 22.2 ± 0.6 years), and 54 older adults (OA, 42.7 ± 0.5 years), all with normal IQ. High model order group probabilistic ICA components (70) were calculated and dual-regression analysis was used to compare 21 RSN's spatial differences between groups. The power spectra were derived from individual ICA mixing matrix time series of the group analyses for frequency domain analysis. We show that primary sensory and motor networks tend to alter more in younger age groups, whereas associative and higher level cognitive networks consolidate and re-arrange until older adulthood. The change has a common trend: both spatial extent and the low frequency power of the RSN's reduce with increasing age. We interpret these result as a sign of normal pruning via focusing of activity to less distributed local hubs.

Global cerebral oxidative metabolism during hypercapnia and hypocapnia in humans: implications for BOLD fMRI

The effect of carbon dioxide (CO(2)) on cerebral metabolism is of tremendous interest to functional imaging. In particular, mild-to-moderate hypercapnia is routinely used in calibrated blood oxygen-level dependent (BOLD)-functional magnetic resonance imaging (fMRI)-based quantification of cerebral oxidative metabolism changes (DeltaCMRO(2)), and relies on the assumption of a stable CMRO(2) during CO(2) challenges. However, this assumption has been challenged by certain animal studies, necessitating its verification in humans and under conditions customary to fMRI. We report, for the first time, on global DeltaCMRO(2) measurements made noninvasively in humans during graded hypercapnia and hypocapnia. We used computerized end-tidal CO(2) modulation to minimize undesired concurrent changes in oxygen pressure, and our findings suggest that no significant change in global CMRO(2) is expected at the levels of end-tidal CO(2) changes customary to calibrated BOLD.

Mechanism of asynchronous Ca(2+) waves underlying agonist-induced contraction in the rat basilar artery

BACKGROUND AND PURPOSE

Uridine 5'-triphosphate (UTP) is a potent vasoconstrictor of cerebral arteries and induces Ca(2+) waves in vascular smooth muscle cells (VSMCs). This study aimed to determine the mechanisms underlying UTP-induced Ca(2+) waves in VSMCs of the rat basilar artery.

EXPERIMENTAL APPROACH

Isometric force and intracellular Ca(2+) ([Ca(2+)](i)) were measured in endothelium-denuded rat basilar artery using wire myography and confocal microscopy respectively.

KEY RESULTS

Uridine 5'-triphosphate (0.1-1000 micromol.L(-1)) concentration-dependently induced tonic contraction (pEC(50) = 4.34 +/- 0.13), associated with sustained repetitive oscillations in [Ca(2+)](i) propagating along the length of the VSMCs as asynchronized Ca(2+) waves. Inhibition of Ca(2+) reuptake in sarcoplasmic reticulum (SR) by cyclopiazonic acid abolished the Ca(2+) waves and resulted in a dramatic drop in tonic contraction. Nifedipine reduced the frequency of Ca(2+) waves by 40% and tonic contraction by 52%, and the nifedipine-insensitive component was abolished by SKF-96365, an inhibitor of receptor- and store-operated channels, and KB-R7943, an inhibitor of reverse-mode Na(+)/Ca(2+) exchange. Ongoing Ca(2+) waves and tonic contraction were also abolished after blockade of inositol-1,4,5-triphosphate-sensitive receptors by 2-aminoethoxydiphenylborate, but not by high concentrations of ryanodine or tetracaine. However, depletion of ryanodine-sensitive SR Ca(2+) stores prior to UTP stimulation prevented Ca(2+) waves.

CONCLUSIONS AND IMPLICATIONS

Uridine 5'-triphosphate-induced Ca(2+) waves may underlie tonic contraction and appear to be produced by repetitive cycles of regenerative Ca(2+) release from the SR through inositol-1,4,5-triphosphate-sensitive receptors. Maintenance of Ca(2+) waves requires SR Ca(2+) reuptake from Ca(2+) entry across the plasma membrane via L-type Ca(2+) channels, receptor- and store-operated channels, and reverse-mode Na(+)/Ca(2+) exchange.

Resting-state functional connectivity in animal models: modulations by exsanguination

We studied the spatiotemporal characteristics of the resting state low frequency fluctuations in functional MRI (fMRI), blood oxygenation level dependent (BOLD) signal in isoflurane-anesthetized rats. fMRI-BOLD measurements at 9.4 Telsa were made during normal and exsanguinated condition previously known to alter cerebral blood flow (CBF) fluctuations in anesthetized rats. fMRI signal time series were low-pass filtered and studied by spectral analysis. During normal conditions, baseline mean arterial pressure (MAP) was 110 +/- 10 mm Hg and low-frequency fluctuations in BOLD signal were observed in the frequency range of 0.01 - 0.125 Hz. Following blood withdrawal (exsanguination), MAP decreased to 68 +/- 7 mm Hg, resulting in an increase in the amplitude of the low-frequency fluctuations in BOLD signal time series and an increase in power at several frequencies between 0.01 and 0.125 Hz. Spatially, the BOLD fluctuations were confined to the cortex and thalamus spanning both hemispheres with sparse presence in the caudate putamen and hippocampus during both normal and exsanguinated states. Spatial distribution of the low frequency fluctuations in BOLD signal, from cross correlation analysis, indicates substantial inter-hemispheric synchrony similar to that observed in the conscious human brain. The behavior of the resting state BOLD signal fluctuations similar to CBF fluctuations during exsanguination indicates a myogenic dependence. Also, a high inter-hemispheric synchrony combined with different phase characteristics of the low frequency BOLD fluctuations particularly in the hippocampus relative to the cortex emphasizes distinct functional networks.

Physiological noise in MR images: an indicator of the tissue response to ischemia?

PURPOSE

To determine whether measuring signal intensity (SI) fluctuations in MRI time series data from acute stroke patients would identify ischemic tissue.

MATERIALS AND METHODS

Prebolus perfusion-weighted MRI data from 32 acute ischemic stroke patients (N = 32) was analyzed as a time series. Ischemic and normal tissue regions were outlined and compared.

RESULTS

The magnitude of the measured SI fluctuations was significantly lower in ischemic regions relative to normal tissue. Spatial differences in these fluctuations occurred in a manner that was different than other perfusion-based metrics.

CONCLUSION

Prior studies have shown that SI fluctuations in MRI time series data correspond to the presence of physiological "noise," which includes vasomotion, an autoregulatory phenomenon that affects the tissue response to ischemia. In this study, SI fluctuations were found to decrease in ischemia, consistent with the notion that small vessels will remain open (fluctuations in vessel diameter will decrease) when there is a challenge to flow. Spatial variation in SI fluctuations appeared to be different from spatial variation seen on other perfusion-based metrics, suggesting that a separate contrast mechanism is responsible, one that might be of diagnostic and prognostic value in acute stroke in which the ability of tissue to withstand ischemia is currently not well visualized.

Source of low-frequency fluctuations in functional MRI signal

PURPOSE

To investigate the source of native low-frequency fluctuations (LFF) in functional MRI (fMRI) signal.

MATERIALS AND METHODS

Phase analysis was performed on tissue-segmented fMRI data acquired at systematically varying sampling rates.

RESULTS

LFF in fMRI signal were both native and aliased in origin. Scanner instability did not contribute to native or aliased LFF. Aliased LFF arose from cardiorespiratory processes and head motion. Native LFF did not arise from cardiorespiratory processes, but did so, at least in part, from head motion. Motion correction reduced native LFF, but did not eliminate them. The residual native LFF in motion-corrected fMRI data showed a systematic phase difference among different tissue structures. The native LFF in fMRI signals of cerebral blood vessels and CSF were synchronous, and preceded those of gray and white matter, indicating that the vascular fluctuations lead the metabolic fluctuations.

CONCLUSION

The primary physiologic source of native LFF in fMRI signal is vasomotion.

Endogenous brain fluctuations and diagnostic imaging

Much of the rising health care costs in aging populations can be attributed to congenital disease and psychiatric and neurologic disorders. Early detection of changes related to these diseases can promote the development of new therapeutic strategies and effective treatments. Changes in tissue, such as damage resulting from continued functional abnormality, often exhibit a time-delay before detection is possible. Methods for detecting functional alterations in endogenous brain fluctuations allow for an early diagnosis before tissue damage occurs, enabling early treatment and a more likely positive outcome. A literature review and comprehensive overview of the current state of knowledge about endogenous brain fluctuations is presented here. Recent findings of the association between various pathological conditions and endogenous fluctuations are discussed. A particular emphasis is placed on research showing the relationship between clinical measures and pathological findings to the dynamics of endogenous fluctuations of the brain. Recent discoveries of methods for detecting abnormal functional connectivity are discussed and future research directions explored.

Spatio-temporal characteristics of low-frequency BOLD signal fluctuations in isoflurane-anesthetized rat brain

We studied the spatio-temporal characteristics of the resting state low-frequency fluctuations in fMRI-BOLD signal in isoflurane-anesthetized rats. fMRI-BOLD measurements at 9.4 T were made during normal and exsanguinated condition previously known to alter cerebral blood flow (CBF) fluctuations in anesthetized rats. fMRI signal time series were low pass filtered and studied by spectral analysis. During normal conditions, baseline mean arterial pressure (MAP) was 110+/-10 mm Hg and low-frequency fluctuations in BOLD signal were observed in the frequency range of 0.01 to 0.125 Hz. Following blood withdrawal (exsanguination), MAP decreased to 68+/-7 mm Hg, resulting in an increase in the amplitude of the low-frequency fluctuations in BOLD signal time series and an increase in power at several frequencies between 0.01 and 0.125 Hz. Spatially, the BOLD fluctuations were confined to the cortex and thalamus spanning both hemispheres with sparse presence in the caudate putamen and hippocampus during both normal and exsanguinated states. Spatial distribution of the low-frequency fluctuations in BOLD signal, from cross-correlation analysis, indicates substantial inter-hemispheric synchrony similar to that observed in the conscious human brain. The behavior of the resting state BOLD signal fluctuations similar to CBF fluctuations during exsanguination indicates a myogenic dependence. Also, a high inter-hemispheric synchrony combined with different phase characteristics of the low-frequency BOLD fluctuations particularly in the hippocampus relative to the cortex emphasizes distinct functional networks.

MRI of spontaneous fluctuations after acute cerebral ischemia in nonhuman primates

PURPOSE

To study the spontaneous low-frequency blood oxygenation level-dependent (BOLD) functional MRI (fMRI) signal fluctuations during hyperacute focal cerebral ischemia.

MATERIALS AND METHODS

A stroke model in nonhuman primates (macaques) was used in this study. Spontaneous fluctuations were recorded using a series of gradient-recalled echo (GRE) echo-planar imaging (EPI) images. Fast Fourier transformation (FFT) was performed on the serial EPI data to calculate the frequency and magnitude of the spontaneous fluctuations. Diffusion tensor imaging (DTI) and perfusion-weighted imaging (PWI) were preformed to detect the ischemic lesion.

RESULTS

The frequency of these fluctuations decreased in the periinfarct tissue in the ipsilateral hemisphere, while their magnitude increased. This area of abnormal signal fluctuations often extended beyond the hyperacute diffusion/perfusion abnormality.

CONCLUSION

This study suggests that measurement of the spontaneous fMRI signal fluctuations provides different information than is available from diffusion/perfusion or T2-weighted MRI.

Large-amplitude, spatially correlated fluctuations in BOLD fMRI signals during extended rest and early sleep stages

A number of recent studies of human brain activity using blood-oxygen-level-dependent (BOLD) fMRI and EEG have reported the presence of spatiotemporal patterns of correlated activity in the absence of external stimuli. Although these patterns have been hypothesized to contain important information about brain architecture, little is known about their origin or about their relationship to active cognitive processes such as conscious awareness and monitoring of the environment. In this study, we have investigated the amplitude and spatiotemporal characteristics of resting-state activity patterns and their dependence on the subjects' alertness. For this purpose, BOLD fMRI was performed at 3.0 T on 12 normal subjects using a visual stimulation protocol, followed by a 27 min rest period, during which subjects were allowed to fall asleep. In subjects who were asleep at the end of the scan, we found (a) a higher amplitude of BOLD signal fluctuation during rest compared with subjects who were awake at the end of the scan; (b) spatially independent patterns of correlated activity that involve all of gray matter, including deep brain nuclei; (c) many patterns that were consistent across subjects; (d) that average percentage levels of fluctuation in visual cortex (VC) and whole brain were higher in subjects who were asleep (up to 1.71% and 1.16%, respectively) than in those who were awake (up to 1.15% and 0.96%) at the end of the scan and were comparable with those levels evoked by intense visual stimulation (up to 1.85% and 0.76% for two subject groups); (e) no confirmation of correlation, positive or negative, between thalamus and VC found in earlier studies. These findings suggest that resting-state activity continues during sleep and does not require active cognitive processes or conscious awareness.

Tissue specificity of nonlinear dynamics in baseline fMRI

In this work, recent advances in the field of nonlinear dynamics (NLD) were applied to fMRI data to examine the spatio-temporal properties of BOLD resting state fluctuations. Five human subjects were imaged during resting state (visual fixation) at 3T using single-shot gradient echo planar imaging (EPI). Respiration and cardiac signals were concurrently recorded for retrospectively removing fluctuations due to these physiologic activities. Patterns of singularity in the complex plane (PSC) and Lempel-Ziv complexity (LZ) were used to study the deterministic nonlinearity in resting state fMRI data. The results show that there is greater nonlinearity (higher PSC) and determinism (lower LZ) in gray matter compared to white matter and CSF. In addition, the removal of respiratory and cardiac pulsations decreases the nonlinearity and determinism but does not alter the relative difference between gray matter and white matter. Therefore, our results demonstrate that determinism and nonlinearity in the fMRI data are tissue-specific, suggesting that they reflect native physiologic and metabolic fluctuations and are not a result of physiologic artifacts due to respiration and cardiac pulsation.

Spontaneous slow hemodynamic oscillations are impaired in cerebral microangiopathy

Small-vessel disease or cerebral microangiopathy (CMA) is a common finding in elderly people. It is related to a variety of vascular risk factors and may finally lead to subcortical ischemic vascular dementia. Because vessel stiffness is increased, we hypothesized that slow spontaneous oscillations are reduced in cerebral hemodynamics. Accordingly, we examined spontaneous oscillations in the visual cortex of 13 patients suffering from CMA, and compared them with 14 age-matched controls. As an imaging method we applied functional near-infrared spectroscopy, because it is particularly sensitive to the microvasculature. Spontaneous low-frequency oscillations (LFOs) (0.07 to 0.12 Hz) were specifically impaired in CMA in contrast to spontaneous very-low-frequency oscillations (0.01 to 0.05 Hz), which remained unaltered. Vascular reagibility was reduced during visual stimulation. Interestingly, changes were tightly related to neuropsychological deficits, namely executive dysfunction. Vascular alterations had to be attributed mainly to the vascular risk factor arterial hypertension. Further, results suggest that the impairments might be, at least partly, reversed by medical treatment such as angiotensin-converting enzyme inhibitors/angiotensin II receptor blockers. Results indicate that functional near-infrared spectroscopy may detect changes in the microvasculature due to CMA, namely an impairment of spontaneous LFOs, and of vascular reagibility. Hence, CMA accelerates microvascular changes due to aging, leading to impairments of autoregulation.

Neurophysiological Architecture of Functional Magnetic Resonance Images of Human Brain

We investigated large-scale systems organization of the whole human brain using functional magnetic resonance imaging (fMRI) data acquired from healthy volunteers in a no-task or 'resting' state. Images were parcellated using a prior anatomical template, yielding regional mean time series for each of 90 regions (major cortical gyri and subcortical nuclei) in each subject. Significant pairwise functional connections, defined by the group mean inter-regional partial correlation matrix, were mostly either local and intrahemispheric or symmetrically interhemispheric. Low-frequency components in the time series subtended stronger inter-regional correlations than high-frequency components. Intrahemispheric connectivity was generally related to anatomical distance by an inverse square law; many symmetrical interhemispheric connections were stronger than predicted by the anatomical distance between bilaterally homologous regions. Strong interhemispheric connectivity was notably absent in data acquired from a single patient, minimally conscious following a brainstem lesion. Multivariate analysis by hierarchical clustering and multidimensional scaling consistently defined six major systems in healthy volunteers-- corresponding approximately to four neocortical lobes, medial temporal lobe and subcortical nuclei- - that could be further decomposed into anatomically and functionally plausible subsystems, e.g. dorsal and ventral divisions of occipital cortex. An undirected graph derived by thresholding the healthy group mean partial correlation matrix demonstrated local clustering or cliquishness of connectivity and short mean path length compatible with prior data on small world characteristics of non-human cortical anatomy. Functional MRI demonstrates a neurophysiological architecture of the normal human brain that is anatomically sensible, strongly symmetrical, disrupted by acute brain injury, subtended predominantly by low frequencies and consistent with a small world network topology.

Eigenvector-based spatial filtering for reduction of physiological interference in diffuse optical imaging

Diffuse optical imaging is an effective technique for noninvasive functional brain imaging. However, the measurements respond to systemic hemodynamic fluctuations caused by the cardiac cycle, respiration, and blood pressure, which may obscure or overwhelm the desired stimulus-evoked response. Previous work on this problem employed temporal filtering, estimation of systemic effects from background pixels, or modeling of interference signals with predefined basis functions, with some success. However, weak signals are still lost in the interference, and other complementary methods are desirable. We use the spatial behavior of measured baseline signals to identify the interference subspaces. We then project signals components in this subspace out of the stimulation data. In doing so, we assume that systemic interference components will be more global spatially, with higher energy, than the stimulus-evoked signals of interest. Thus, the eigenvectors corresponding to the largest eigenvalues of an appropriate correlation matrix form the basis for an interference subspace. By projecting the data onto the orthogonal nullspace of these eigenvectors, we can obtain more localized response, as reflected in improved contrast-to-noise ratio and correlation coefficient maps.

Spontaneous low-frequency oscillations decline in the aging brain

It is well known that aging leads to a degeneration of the vascular system. Hence, one may hypothesize that spontaneous oscillations decrease in the cerebral microvasculature with aging. Accordingly, the authors investigated the age dependency of spontaneous oscillations in the visual cortex during rest and functional activation. Functional near-infrared spectroscopy was used because it is particularly sensitive to the microvasculature. Visual stimulation led to an increase of oxyhemoglobin, total hemoglobin, and a decrease of deoxyhemoglobin, without any influence of age. Peaks of normalized power spectral density were detected for spontaneous low-frequency (0.07 to 0.11 Hz) and very-low-frequency (0.01 to 0.05 Hz) oscillations, with a higher amplitude for oxyhemoglobin than for deoxyhemoglobin. Spontaneous low-frequency oscillations of oxyhemoglobin and deoxyhemoglobin declined strongly with aging during both rest and visual stimulation. Reduction of spontaneous low-frequency oscillations might indicate a declining spontaneous activity in microvascular smooth muscle cells, in conjunction with an increased vessel stiffness with aging.

Regulation of oxygen consumption by vasomotion

OBJECTIVE

To describe a stochastic model of the variability and heterogeneity of blood flow through the microcirculation, and to show the ability of vasomotion to vary oxygen consumption at a steady blood flow.

METHODS

The description of vasomotion is based on whether each microvessel is open for blood flow or closed. Over a unit time period, let alpha be the probability that a given vessel is open and will remain open, beta be the probability that an open vessel will close, nu be the probability that a closed vessel will remain closed, and mu be the probability that a closed vessel will become open. Two main parameters that characterize such a scheme are: the fraction of open microvessels [n o= mu/(beta+mu)], and the rate of vasomotion defined as the rate of switching between open and closed microvessels (R=beta+mu). A model of O2 transport to tissues is based on the following assumptions: (a) the flux of O2 is due to passive diffusion, (b) the amount of O2 dissolved in tissue is negligible as compared with that contained in arterial blood, and (c) aerobic metabolism is proportional to the delivery of O2.

RESULTS

The stochastic model substantiates the possibility that vasomotion can control O2 consumption. The rate of vasomotion activity can change O2 consumption by 2-8-fold, depending on the fraction of open microvessels.

CONCLUSION

A stochastic description of blood flow through the microcirculation system demonstrates that vasomotion rate could be a factor influencing O2 consumption.

Resting fluctuations in arterial carbon dioxide induce significant low frequency variations in BOLD signal

Carbon dioxide is a potent cerebral vasodilator. We have identified a significant source of low-frequency variation in blood oxygen level-dependent (BOLD) magnetic resonance imaging (MRI) signal at 3 T arising from spontaneous fluctuations in arterial carbon dioxide level in volunteers at rest. Fluctuations in the partial pressure of end-tidal carbon dioxide (Pet(CO(2))) of +/-1.1 mm Hg in the frequency range 0-0.05 Hz were observed in a cohort of nine volunteers. Correlating with these fluctuations were significant generalized grey and white matter BOLD signal fluctuations. We observed a mean (+/-standard error) regression coefficient across the group of 0.110 +/- 0.033% BOLD signal change per mm Hg CO(2) for grey matter and 0.049 +/- 0.022% per mm Hg in white matter. Pet(CO(2))-related BOLD signal fluctuations showed regional differences across the grey matter, suggesting variability of the responsiveness to carbon dioxide at rest. Functional magnetic resonance imaging (fMRI) results were corroborated by transcranial Doppler (TCD) ultrasound measurements of the middle cerebral artery (MCA) blood velocity in a cohort of four volunteers. Significant Pet(CO(2))-correlated fluctuations in MCA blood velocity were observed with a lag of 6.3 +/- 1.2 s (mean +/- standard error) with respect to Pet(CO(2)) changes. This haemodynamic lag was adopted in the analysis of the BOLD signal. Doppler ultrasound suggests that a component of low-frequency BOLD signal fluctuations is mediated by CO(2)-induced changes in cerebral blood flow (CBF). These fluctuations are a source of physiological noise and a potentially important confounding factor in fMRI paradigms that modify breathing. However, they can also be used for mapping regional vascular responsiveness to CO(2).

Detection of fast neuronal signals in the motor cortex from functional near infrared spectroscopy measurements using independent component analysis

Fast changes, in the range of milliseconds, in the optical properties of cerebral tissue are associated with brain activity and can be detected using non-invasive near-infrared spectroscopy (NIRS). These changes are assumed to be caused by changes in the light scattering properties of the neuronal tissue. The aim of this study was to develop highly sensitive data analysis algorithms to detect this fast signal, which is small compared with other physiological signals. A frequency-domain tissue oximeter, whose laser diodes were intensity modulated at 110 MHz, was used. The amplitude, mean intensity and phase of the modulated optical signal were measured at a sample rate of 96 Hz. The probe, consisting of four crossed source detector pairs was placed above the motor cortex, contralateral to the hand performing a tapping exercise consisting of alternating rest and tapping periods of 20 s each. An adaptive filter was used to remove the arterial pulsatility from the optical signals. Independent component analysis allowed further separation of a signal component containing the fast signal. In nine out of 14 subjects, a significant fast neuronal signal related to the finger tapping was found in the intensity signals. In the phase signals, indications of the fast signal were found in only two subjects.

High field human imaging

This review article examines the state of knowledge regarding human imaging using MRI at high main magnetic field strengths. The article starts with a summary of the technical issues associated with magnetic field strengths in the range of 3-8 T, including magnet characteristics and the properties of radiofrequency magnetic fields, with special reference to sensitivity, power deposition, and homogeneity. The published data on tissue-water relaxation times in the brain is tabulated and the implications for contrast and pulse sequence implementation is elucidated. The behavior of the major fast imaging sequences, fast low angle shot (FLASH), rapid acquisition with relaxation enhancement (RARE), and echo planar imaging (EPI), is examined in this context. A number of anatomical images from 3 T systems are presented as examples. Particular attention is given to various forms of vascular imaging, namely, time of flight angiography, venography, and arterial spin labeling. The most complex changes in contrast with main magnetic field strength are in activation studies utilizing the blood oxygen level dependent mechanism, which are examined in detail. Improvements in spatial specificity are emphasized, particularly in conjunction with spin-echo imaging. The article concludes with a discussion of the current status and the potential impact of technical developments such as parallel imaging.

Differential activation of ion channels by inositol 1,4,5-trisphosphate (IP3)- and ryanodine-sensitive calcium stores in rat basilar artery vasomotion

Spontaneous, rhythmical contractions, or vasomotion, can be recorded from cerebral vessels under both normal physiological and pathophysiological conditions. Using electrophysiology to study changes in membrane potential, the ratiometric calcium indicator Fura-2 AM to study changes in [Ca(2+)](i) in both the arterial wall and in individual smooth muscle cells (SMCs), and video microscopy to study changes in vessel diameter, we have investigated the cellular mechanisms underlying vasomotion in the juvenile rat basilar artery. During vasomotion, rhythmical oscillations in both membrane potential and [Ca(2+)](i) were found to precede rhythmical contractions. Nifedipine depolarized SMCs and abolished rhythmical contractions and depolarizations. [Ca(2+)](i) oscillations in the arterial wall became reduced and irregular, while [Ca(2+)](i) oscillations in adjacent SMCs were no longer synchronized. BAPTA-AM, thapsigargin and U73122 hyperpolarized SMCs, relaxed the vessel, decreased basal calcium levels and abolished vasomotion. Chloride substitution abolished rhythmical activity, depolarized SMCs, increased basal calcium levels and constricted the vessel, while niflumic acid and DIDS abolished vasomotion. Ryanodine, charybdotoxin and TRAM-34, but not iberiotoxin, 4-aminopyridine or apamin, each depolarized SMCs and increased the frequency of rhythmical depolarizations and [Ca(2+)](i) oscillations. We conclude that vasomotion in the basilar artery depends on the release of intracellular calcium from IP(3) (inositol 1,4,5,-trisphosphate)-sensitive stores which activates calcium-dependent chloride channels to depolarize SMCs. Depolarization in turn activates voltage-dependent calcium channels, synchronizing contractions of adjacent cells through influx of extracellular calcium. Subsequent calcium-induced calcium release from ryanodine-sensitive stores activates an intermediate conductance potassium channel, hyperpolarizing the SMCs and providing a negative feedback pathway for regeneration of the contractile cycle.

Sound repetition rate in the human auditory pathway: representations in the waveshape and amplitude of fMRI activation

Sound repetition rate plays an important role in stream segregation, temporal pattern recognition, and the perception of successive sounds as either distinct or fused. This study was aimed at elucidating the neural coding of repetition rate and its perceptual correlates. We investigated the representations of rate in the auditory pathway of human listeners using functional magnetic resonance imaging (fMRI), an indicator of population neural activity. Stimuli were trains of noise bursts presented at rates ranging from low (1-2/s; each burst is perceptually distinct) to high (35/s; individual bursts are not distinguishable). There was a systematic change in the form of fMRI response rate-dependencies from midbrain to thalamus to cortex. In the inferior colliculus, response amplitude increased with increasing rate while response waveshape remained unchanged and sustained. In the medial geniculate body, increasing rate produced an increase in amplitude and a moderate change in waveshape at higher rates (from sustained to one showing a moderate peak just after train onset). In auditory cortex (Heschl's gyrus and the superior temporal gyrus), amplitude changed somewhat with rate, but a far more striking change occurred in response waveshape-low rates elicited a sustained response, whereas high rates elicited an unusual phasic response that included prominent peaks just after train onset and offset. The shift in cortical response waveshape from sustained to phasic with increasing rate corresponds to a perceptual shift from individually resolved bursts to fused bursts forming a continuous (but modulated) percept. Thus at high rates, a train forms a single perceptual "event," the onset and offset of which are delimited by the on and off peaks of phasic cortical responses. While auditory cortex showed a clear, qualitative correlation between perception and response waveshape, the medial geniculate body showed less correlation (since there was less change in waveshape with rate), and the inferior colliculus showed no correlation at all. Overall, our results suggest a population neural representation of the beginning and the end of distinct perceptual events that is weak or absent in the inferior colliculus, begins to emerge in the medial geniculate body, and is robust in auditory cortex.

Voltage independence of vasomotion in isolated irideal arterioles of the rat

The cellular mechanisms underlying vasomotion of irideal arterioles from juvenile rats have been studied using electrophysiological methods, ratiometric calcium measurements and video microscopy. Vasomotion was not affected by removal of the endothelium. Spontaneous contractions were preceded by spontaneous depolarizations. Both were abolished by the intracellular calcium chelator, BAPTA AM (20 microM), but not by ryanodine (10 microM), suggesting a dependence on the cyclical release of calcium from intracellular stores, other than those operated by ryanodine receptors. Oscillations were little changed when the membrane potential of short segments of arteriole was either depolarized or hyperpolarized. When the segments were voltage clamped, oscillating inward currents were recorded, indicating that the changes in membrane potential were voltage independent. Vasomotion was preceded by intracellular calcium oscillations and both were abolished by inhibitors of phospholipase C (U73122, 10 microM), phospholipase A(2) (AACOCF(3), 30 microM) and protein kinase C (chelerythrine chloride, 5 microM, and myristoylated protein kinase C peptide, 10 microM). Inhibition of vasomotion by the dual lipoxygenase and cyclo-oxygenase inhibitor, NDGA (10 microM), the lipoxygenase inhibitor, ETI (1 microM) but not by the cyclo-oxygenase inhibitors, aspirin (10 microM) and indomethacin (10 microM), or the cytochrome P450 inhibitor 17-ODYA (10 microM), suggested an involvement of the lipoxygenase pathway. The observations suggest that vasomotion of iris arterioles is voltage independent and results from the cyclical release of calcium from IP(3)-sensitive stores which are activated by cross talk between the phospholipase C and phospholipase A(2) pathways in vascular smooth muscle.