The effect of physiological noise in phase functional magnetic resonance imaging: from blood oxygen level-dependent effects to direct detection of neuronal currents

Recently, the possibility to use both magnitude and phase image sets for the statistical evaluation of fMRI has been proposed, with the prospective of increasing both statistical power and the spatial specificity. In the present work, several issues that affect the spatial and temporal stability in fMRI phase time series in the presence of physiologic noise processes are reviewed, discussed and illustrated by experiments performed at 3 T. The observed phase value is a fingerprint of the underlying voxel averaged magnetic field variations. Those related to physiological processes can be considered static or dynamic in relation to the temporal scale of a 2D acquisition and will play out on different spatial scales as well: globally across the entire images slice, and locally depending on the constituents and their relative fractions inside the MRI voxel. The 'static' respiration-induced effects lead to magneto-mechanic scan-to-scan variations in the global magnetic field but may also contribute to local BOLD fluctuations due to respiration-related variations in arterial carbon dioxide. Likewise, the 'dynamic' cardiac-related effects will lead to global susceptibility effects caused by pulsatile motion of the brain as well as local blood pressure-related changes in BOLD and changes in blood flow velocity. Finally, subject motion may lead to variations in both local and global tissue susceptibility that will be especially pronounced close to air cavities. Since dissimilar manifestations of physiological processes can be expected in phase and in magnitude images, a direct relationship between phase and magnitude scan-to-scan fluctuations cannot be assumed a priori. Therefore three different models were defined for the phase stability, each dependent on the relation between phase and magnitude variations and the best will depend on the underlying noise processes. By experiments on healthy volunteers at rest, we showed that phase stability depends on the type of post-processing and can be improved by reducing the low-frequency respiration-induced mechano-magnetic effects. Although the manifestations of physiological noise were in general more pronounced in phase than in magnitude images, due to phase wraps and global Bo effects, we suggest that a phase stability similar to that found in magnitude could theoretically be achieved by adequate correction methods. Moreover, as suggested by our experimental data regarding BOLD-related phase effects, phase stability could even supersede magnitude stability in voxels covering dense microvascular networks with BOLD-related fluctuations as the dominant noise contributor. In the interest of the quality of both BOLD-based and nc-MRI methods, future studies are required to find alternative methods that can improve phase stability, designed to match the temporal and spatial scale of the underlying neuronal activity.     

Neural and vascular variability and the fMRI-BOLD response in normal aging

Neural, vascular and structural variables contributing to the blood oxygen level-dependent (BOLD) signal response variability were investigated in younger and older humans. Twelve younger healthy human subjects (six male and six female; mean age: 24 years; range: 19–27 years) and 12 older healthy subjects (five male and seven female; mean age: 58 years; range: 55–71 years) with no history of head trauma and neurological disease were scanned. Functional magnetic resonance imaging measurements using the BOLD contrast were made when participants performed a motor, cognitive or a breath hold (BH) task. Activation volume and the BOLD response amplitude were estimated for the younger and older at both group and subject levels. Mean activation volume was reduced by 45%, 40% and 38% in the elderly group during the motor, cognitive and BH tasks, respectively, compared to the younger. Reduction in activation volume was substantially higher compared to the reduction in the gray matter volume of 14% in the older compared to the younger. A significantly larger variability in the intersubject BOLD signal change occurred during the motor task, compared to the cognitive task. BH-induced BOLD signal change between subjects was significantly less-variable in the motor task-activated areas in the younger compared to older whereas such a difference between age groups was not observed during the cognitive task. Hemodynamic scaling using the BH signal substantially reduced the BOLD signal variability during the motor task compared to the cognitive task. The results indicate that the origin of the BOLD signal variability between subjects was predominantly vascular during the motor task while being principally a consequence of neural variability during the cognitive task. Thus, in addition to gray matter differences, the type of task performed can have different vascular variability weighting that can influence age-related differences in brain functional response.     

Comparison of analytical strategies for EEG-correlated fMRI data in patients with epilepsy

The simultaneous recording of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) can be used to localize interictal epileptiform discharges (IEDs). Previous studies have reported varying degrees of concordance of EEG-fMRI with electroclinical findings. The aim of the present study is to evaluate to what extent this variability is determined by the analytical strategy or by the properties of the EEG data. For that purpose, 42 IED sets obtained in 29 patients with epilepsy were reanalyzed using a finite impulse response approach, which estimates the hemodynamic response function (HRF) from the data and allows non-causal effects. Cardiac effects were treated as additional confounders in the model. This approach was compared to the classical approach assuming a fixed HRF for each voxel in the brain. The performance of each method was assessed by comparing the fMRI results to the EEG focus. The flexible model revealed more significantly activated voxels, which resulted in more activated brain regions concordant with the EEG focus (26 vs. 16). Correction for cardiac effects improved the results in 7 out of the 42 data sets. Furthermore, design theory for event-related experiments was applied in order to determine the influence of the number of IEDs and their temporal distribution on the success of an experiment. It appeared that this success is highly dependent upon the number of IEDs present during the recording and less on their temporal spacing. We conclude that the outcome of EEG-fMRI can be improved by using an optimized analytical strategy, but also depends on the number of IEDs occurring during the recording.     

A new automatic phase mask filter for high-resolution brain venography at 3 T: theoretical background and experimental validation

To improve vessel contrast in high-resolution susceptibility-based brain venography, an automatic phase contrast enhancing procedure is proposed, based on a new phase mask filter suitable for maximizing contrast of venous MR signals. The effectiveness of the new approach was assessed both on digital phantoms and on acquired MR human brain images, and then compared with venographic results of phase masking methods in recent literature. The digital phantom consisted of a simulated MR dataset with given signal-to-noise ratios (SNRs), while real human data were collected by scanning healthy volunteers with a 3.0-T MR system and a 3D gradient echo pulse sequence. The new phase mask (NM) was more effective than the conventional mask (CM) both on the digital phantoms and on the acquired MR images. A quantitative comparison based on phantom venograms indicates how this phase enhancement can lead to a significant increase in the contrast-to-noise ratio (CNR) for all considered phase values as well as for all vessel sizes of clinical interest. Likewise, the in vivo brain venograms reveal a better depiction of the smallest venous vessels and the enhancement of many details undetectable in conventional venograms.     

In contrast to BOLD: signal enhancement by extravascular water protons as an alternative mechanism of endogenous fMRI signal change

Despite the popularity and widespread application of functional magnetic resonance imaging (fMRI) in recent years, the physiological bases of signal change are not yet fully understood. Blood oxygen level-dependant (BOLD) contrast — attributed to local changes in blood flow and oxygenation, and therefore magnetic susceptibility — has become the most prevalent means of functional neuroimaging. However, at short echo times, spin-echo sequences show considerable deviations from the BOLD model, implying a second, non-BOLD component of signal change. This has been dubbed “signal enhancement by extravascular water protons” (SEEP) and is proposed to result from proton-density changes associated with cellular swelling. Given that such changes are independent of magnetic susceptibility, SEEP may offer new and improved opportunities for carrying out fMRI in regions with close proximity to air–tissue and/or bone–tissue interfaces (e.g., the prefrontal cortex and spinal cord), as well as regions close to large blood vessels, which may not be ideally suited for BOLD imaging. However, because of the interdisciplinary nature of the literature, there has yet to be a thorough synthesis, tying together the various and sometimes disparate aspects of SEEP theory. As such, we aim to provide a concise yet comprehensive overview of SEEP, including recent and compelling evidence for its validity, its current applications and its future relevance to the rapidly expanding field of functional neuroimaging. Before presenting the evidence for a non-BOLD component of endogenous functional contrast, and to enable a more critical review for the nonexpert reader, we begin by reviewing the fundamental principles underlying BOLD theory.     

A preliminary study of blood-oxygen-level-dependent MRI in patients with chronic kidney disease

Background

Blood-oxygen-level-dependent (BOLD) magnetic resonance imaging (MRI) can provide regional measurements of oxygen content using deoxyhemoglobin paramagnetic characteristics. The apparent relaxation rate or R2*(=1/T2*) can be determined from the slope of log (intensity) versus echo time and is directly proportional to the tissue content of deoxyhemoglobin. Thus, as the level of deoxyhemoglobin increases, T2* will decrease, leading to an increase in R2*. Chronic kidney disease (CKD) can affect oxygenation levels in renal parenchyma, which influences the clinical course of the disease. The goal of this study was to detect and assess renal oxygenation levels in CKD using BOLD MRI.

Methods

Fifteen healthy subjects and 11 patients with CKD underwent a renal scan using multigradient-recalled-echo sequence with eight echoes. R2* (1/s) of the renal cortex and medulla was measured on BOLD images. Of the 11 patients, nine had biopsy-proven chronic glomerulonephritis, and two had a similar diagnosis based on clinical symptoms and investigations.

Results

Mean medullary R2* (MR2*) and cortex R2* (CR2*) levels were significantly higher in patients (22 kidneys, MR2*=24.79±4.84 s⁻¹, CR2*=18.97±2.72 s⁻¹) than in controls (30 kidneys, MR2*=19.98±1.19 s⁻¹, CR2*=16.03±1.23 s⁻¹) (P<.01), and MR2* was increased more than CR2*. Medullary to cortical R2* ratios (MCR2*) of patients were significantly increased when compared with those of controls (P<.01). In the patient group, estimated glomerular filtration rate levels were greater than or equal to 60 ml/min/1.73 m² in six patients (12 kidneys), whose MR2* and CR2* were also significantly higher than those of controls (P<.01). Serum creatinine levels were normal in seven patients (14 kidneys), whose MR2*, CR2* and MCR2* were also higher than those of controls (P<.01).

Conclusions

BOLD MRI can be used to evaluate changes in renal oxygenation in CKD, suggesting that it has the potential to be an excellent noninvasive tool for the evaluation of renal function.     

Functional MRI using super-resolved spatiotemporal encoding

Recently, new ultrafast imaging sequences such as rapid acquisition by sequential excitation and refocusing (RASER) and hybrid spatiotemporal encoding (SPEN) magnetic resonance imaging (MRI) have been proposed, in which the phase encoding of conventional echo planar imaging (EPI) is replaced with a SPEN. In contrast to EPI, SPEN provides significantly higher immunity to frequency heterogeneities including those caused by B₀ inhomogeneities and chemical shift offsets. Utilizing the inherent robustness of SPEN, it was previously shown that RASER can be used to successfully perform functional MRI (fMRI) experiments in the orbitofrontal cortex — a task which is challenging using EPI due to strong magnetic susceptibility variation near the air-filled sinuses. Despite this superior performance, systematic analyses have shown that, in its initial implementation, the use of SPEN was penalized by lower signal-to-noise ratio (SNR) and higher radiofrequency power deposition as compared to EPI-based methods. A recently developed reconstruction algorithm based on super-resolution principles is able to alleviate both of these shortcomings; the use of this algorithm is hereby explored within an fMRI context. Specifically, a series of fMRI measurements on the human visual cortex confirmed that the super-resolution algorithm retains the statistical significance of the blood oxygenation level dependent (BOLD) response, while significantly reducing the power deposition associated with SPEN and restoring the SNR to levels that are comparable with those of EPI.     

A nonlinear BOLD model accounting for refractory effect by applying the longitudinal relaxation in NMR to the linear BOLD model

A mathematical model to regress the nonlinear blood oxygen level-dependent (BOLD) fMRI signal has been developed by incorporating the refractory effect into the linear BOLD model of the biphasic gamma variate function. The refractory effect was modeled as a relaxation of two separate BOLD capacities corresponding to the biphasic components of the BOLD signal in analogy with longitudinal relaxation of magnetization in NMR. When tested with the published fMRI data of finger tapping, the nonlinear BOLD model with the refractory effect reproduced the nonlinear BOLD effects such as reduced poststimulus undershoot and saddle pattern in a prolonged stimulation as well as the reduced BOLD signal for repetitive stimulation.     

BOLD functional MRI in disease and pharmacological studies: room for improvement?

In the past decade the use of blood oxygen level-dependent (BOLD) fMRI to investigate the effect of diseases and pharmacological agents on brain activity has increased greatly. BOLD fMRI does not measure neural activity directly, but relies on a cascade of physiological events linking neural activity to the generation of MRI signal. However, most of the disease and pharmacological studies performed so far have interpreted changes in BOLD fMRI as "brain activation," ignoring the potential confounds that can arise through drug- or disease-induced modulation of events downstream of the neural activity. This issue is especially serious in diseases (like multiple sclerosis, brain tumours and stroke) and drugs (like anaesthetics or those with a vascular action) that are known to influence these physiological events. Here we provide evidence that, to extract meaningful information on brain activity in patient and pharmacological BOLD fMRI studies, it is important to identify, characterise and possibly correct these influences that potentially confound the results. We suggest a series of experimental measures to improve the interpretability of BOLD fMRI studies. We have ranked these according to their potential information and current practical feasibility. First-line, necessary improvements consist of (1) the inclusion of one or more control tasks, and (2) the recording of physiological parameters during scanning and subsequent correction of possible between-group differences. Second-line, highly recommended important aim to make the results of a patient or drug BOLD study more interpretable and include the assessment of (1) baseline brain perfusion, (2) vascular reactivity, (3) the inclusion of stimulus-related perfusion fMRI and (4) the recording of electrophysiological responses to the stimulus of interest. Finally, third-line, desirable improvements consist of the inclusion of (1) simultaneous EEG-fMRI, (2) cerebral blood volume and (3) rate of metabolic oxygen consumption measurements and, when relevant, (4) animal studies investigating signalling between neural cells and blood vessels.     

Hemodynamic scaling of fMRI-BOLD signal: validation of low-frequency spectral amplitude as a scalability factor

Functional magnetic resonance imaging blood-oxygenation-level-dependent (fMRI-BOLD) signal representing neural activity may be optimized by discriminating MR signal components related to neural activity and those related to intrinsic properties of the cortical vasculature. The objective of this study was to reduce the hemodynamic change independent of neural activity to obtain a scaled fMRI-BOLD response using two factors, namely, low-frequency spectral amplitude (LFSA) and breath-hold amplitude (BHA). Ten subjects (age range, 22–38 years) were scanned during four task conditions: (a) rest while breathing room air, (b) bilateral finger tapping while breathing room air, (c) rest during a partial inspirational breath-hold, and (d) rest during moderate hypercapnia (breathing 5% CO₂, 20% O₂ and 75% N₂). In all subjects who breathed 5% CO₂, regions with significant BOLD response during breath-hold correlated significantly with the percent signal increase during 5% CO₂ inhalation. Finger-tapping-induced responses in the motor cortex were diminished to a similar extent after scaling using either LFSA or BHA. Inter- and intrasubject variation in the amplitude of the BOLD signal response reduced after hemodynamic scaling using LFSA or BHA. The results validated the hemodynamic amplitude scaling using LFSA with the earlier established BHA. LFSA free from motor-task contamination can be used to calibrate the fMRI-BOLD response in lieu of BHA or hypercapnia to minimize intra- and intersubject variation arising from vascular anatomy and vasodilative capacity.