Magnetic resonance microscopic imaging based on high-order intermolecular multiple-quantum coherences

Most imaging studies using intermolecular multiple-quantum coherences (iMQCs) have focused on the two-spin dipolar interactions--zero and double quantum coherences. Here, we report the results of various experimental studies to assess the feasibility of magnetic resonance microscopy with high-order iMQCs in model systems at 7 and 14 T. Experimental results demonstrated that the iMQC microscopic images with high coherence orders are readily observable at high field and have unique contrast depending on the sample microstructure and coherence order.     

Simultaneous acquisition of multiple orders of intermolecular multiple-quantum coherence images

Recent studies have demonstrated the ability to detect images based on intermolecular multiple-quantum coherences (iMQCs) that correspond to flipping of two or more separated spins simultaneously, as opposed to conventional magnetic resonance where only one spin is flipped at a time. Until now, iMQC imaging has only acquired one coherence signal per pulse sequence. Here we report a new sequence that successfully detects five orders of coherence (2, 1, 0, −1, and −2-quantum coherence images) in one pulse sequence, with each signal having its full intensity. The simultaneous acquisition highlights substantial contrast differences between conventional and iMQC images, and between the different types of iMQC images.     

Effects of mild hypoxic hypoxia on poststimulus undershoot of blood-oxygenation-level-dependent fMRI signal in the human visual cortex

Characteristics of the blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal poststimulus undershoot in the visual cortex were studied at varying levels of arterial blood oxygen saturation (Ysat). Undershoot with an amplitude of -0.6+/-0.2% appeared after positive BOLD response (+1.7+/-0.5%) under control conditions. Cerebral blood volume (CBV), as determined with vascular-space-occupancy-dependent fMRI, increased by 26-43% during the positive BOLD peak, but the CBV proceeded at baseline level during the BOLD poststimulus undershoot. Mild hypoxic hypoxia (Ysat ranging from 0.82 to 0.89) had no effect on the amplitude or duration of poststimulus undershoot in activated BOLD pixels. Hypoxia did not influence CBV during the BOLD poststimulus undershoot. In contrast, the positive BOLD signal at the level of all activated pixels was smaller in hypoxia than in normoxia. The present results show that the BOLD poststimulus undershoot is not influenced by curtailed oxygen availability and that, during the undershoot, CBV is not different from resting state.     

Simultaneous acquisition of multiple orders of intermolecular multiple-quantum coherence images in vivo

Until recently, NMR imaging with intermolecular multiple-quantum coherences (iMQCs) has been based on the acquisition of a single echo. In vivo studies of iMQC image contrast would greatly benefit from a method that could acquire several orders of quantum coherence during the same acquisition. This would enable comparison of the image contrast for various orders and eliminate image coregistration problems between scans. It has previously been demonstrated that multiple orders of iMQC images can be simultaneously acquired of a simple phantom. Here, we examine the technique and its effect on biological tissue, both in vivo and in vitro. First, we establish the effectiveness of the iMQC sequence in vivo using earthworms as specimens. We then further show that the multi-CRAZED sequence enhances detection of next generation (nanoparticle) contrast agents on excised tumor tissue.     

Effects of hypoglycemia on human brain activation measured with fMRI

Functional magnetic resonance imaging (fMRI) was used to measure the effects of acute hypoglycemia caused by passive sensory stimulation on brain activation. Visual stimulation was used to generate blood-oxygen-level-dependent (BOLD) contrast, which was monitored during hyperinsulinemic hypoglycemic and euglycemic clamp studies. Hypoglycemia (50 +/- 1 mg glucose/dl) decreased the fMRI signal relative to euglycemia in 10 healthy human subjects: the fractional signal change was reduced by 28 +/- 12% (P < .05). These changes were reversed when euglycemia was restored. These data provide a basis of comparison for studies that quantify hypoglycemia-related changes in fMRI activity during cognitive tasks based on visual stimuli and demonstrate that variations in blood glucose levels may modulate BOLD signals in the healthy brain.     

Functional magnetic resonance imaging reference phantom

Functional magnetic resonance imaging (fMRI) is widely used to pinpoint active brain areas. Changes in neuronal activity modulate the local blood oxygenation level, and the associated modulation of the magnetic field homogeneity can be detected with magnetic resonance imaging. Thus, the blood oxygenation level-dependent (BOLD) fMRI indirectly measures neuronal activity. Similar modulation of magnetic field homogeneity was here elicited by other means to generate a BOLD-like change in a new phantom constructed to provide reference activations during fMRI. Magnetic inhomogeneities were produced by applying current to coils located near the phantom containing 1.5 ml of Gd-doped water. The signal-to-noise ratio of the images, produced by gradient-recalled echo-planar imaging, varied between 104 and 107 at a selected voxel when the field was and was not inhomogenized, respectively. The contrast of signals between homogeneous and inhomogeneous conditions was generally stable, except in 3% of time points. During the periods of greatest deviations an observable change would have been detected in a simultaneously measured BOLD signal. Such changes could result from the imaging method or occur through glitches in hardware or alterations in the measurement environment. With identical measurement setups, the phantom could allow comparing intersession or intersubject brain activations.     

Synchronized detection of minute electrical currents with MRI using Lorentz effect imaging

The blood oxygenation level-dependent (BOLD) effect is the most commonly used contrast mechanism in functional magnetic resonance imaging (fMRI), due to its relatively high spatial resolution and sensitivity. However, the ability of BOLD fMRI to accurately localize neuronal activation in space and time is limited by the inherent hemodynamic modulation. There is hence a need to develop alternative MRI methods that can directly image neuroelectric activity, thereby achieving both a high temporal resolution and spatial specificity as compared to conventional BOLD fMRI. In this paper, we extend the Lorentz effect imaging technique, which can detect spatially incoherent yet temporally synchronized minute electrical activity in a strong magnetic field, and demonstrate its feasibility for imaging randomly oriented electrical currents on the order of microamperes with a temporal resolution on the order of milliseconds in gel phantoms. This constitutes a promising step towards its application to direct imaging of neuroelectric activity in vivo, which has the same order of current density and temporal synchrony.     

High-resolution absorptive intermolecular multiple-quantum coherence NMR spectroscopy under inhomogeneous fields

Intermolecular multiple-quantum coherence (iMQC) is capable of improving NMR spectral resolution using a 2D shearing manipulation method. A pulse sequence termed CT-iDH, which combines intermolecular double-quantum filter (iDQF) with a modified constant-time (CT) scheme, is designed to achieve fast acquisition of high-resolution intermolecular zero-quantum coherences (iZQCs) and intermolecular double-quantum coherences (iDQCs) spectra without strong coupling artifacts. Furthermore, double-absorption lineshapes are first realized in 2D intermolecular multi-quantum coherences (iMQCs) spectra under inhomogeneous fields through a combination of iZQC and iDQC signals to double the resolution without loss of sensitivity. Theoretically the spectral linewidth can be further reduced by half compared to original iMQC high-resolution spectra. Several experiments were performed to test the feasibility of the new method and the improvements are evaluated quantitatively. The study suggests potential applications for in vivo spectroscopy.     

人脑初级视皮层内事件相关型视觉刺激fMRI BOLD响应特性初步研究

基于脑血氧水平依赖(BOLD)对比的功能磁共振成像（fMRI）方法是研究脑功能活动的一种重要的无损伤探测手段，BOLD的机制及它与脑神经活动关系一直是国际上十分活跃的研究领域，尤其是在实验研究方面. 视觉刺激所引起的脑的初级视皮层(V1) 的BOLD响应的时间特性已有较多的研究，但是在这些研究中，有的结论认为BOLD对视觉刺激的响应是线性的，有的结论却是相反的. 我们采用事件关联型核磁共振功能成像（ER-fMRI）方法，研究不同的短暂视觉刺激持续时间下的BOLD响应，得到了视觉刺激下的脑激活图. 同时找出了视觉刺激时间分别是1、2、3、4、5和6 s的V1区的BOLD响应曲线，并初步显示出BOLD响应与刺激持续时间的线性关系.     

To smooth or not to smooth? ROC analysis of perfusion fMRI data

Blood oxygenation level dependent (BOLD) contrast has been widely used for visualizing regional neural activation. Temporal filtering and parameter estimation algorithms are generally used to account for the intrinsic temporal autocorrelation present in BOLD data. Arterial spin labeling perfusion imaging is an emerging methodology for visualizing regional brain function both at rest and during activation. Perfusion contrast manifests different noise properties compared with BOLD contrast, represented by the even distribution of noise power and spatial coherence across the frequency spectrum. Consequently, different strategies are expected to be employed in the statistical analysis of functional magnetic resonance imaging (fMRI) data based on perfusion contrast. In this study, the effect of different analysis methods upon signal detection efficacy, as assessed by receiver operator characteristic (ROC) measures, was examined for perfusion fMRI data. Simulated foci of neural activity of varying amplitude and spatial extent were added to resting perfusion data, and the accuracy of each analysis was evaluated by comparing the results with the known distribution of pseudo-activation. In contrast to the BOLD fMRI, temporal smoothing or filtering reduces the power of perfusion fMRI data analyses whereas spatial smoothing is beneficial to the efficacy of analyses.     

A kernel machine-based fMRI physiological noise removal method

Functional magnetic resonance imaging (fMRI) technique with blood oxygenation level dependent (BOLD) contrast is a powerful tool for noninvasive mapping of brain function under task and resting states. The removal of cardiac- and respiration-induced physiological noise in fMRI data has been a significant challenge as fMRI studies seek to achieve higher spatial resolutions and characterize more subtle neuronal changes. The low temporal sampling rate of most multi-slice fMRI experiments often causes aliasing of physiological noise into the frequency range of BOLD activation signal. In addition, changes of heartbeat and respiration patterns also generate physiological fluctuations that have similar frequencies with BOLD activation. Most existing physiological noise-removal methods either place restrictive limitations on image acquisition or utilize filtering or regression based post-processing algorithms, which cannot distinguish the frequency-overlapping BOLD activation and the physiological noise. In this work, we address the challenge of physiological noise removal via the kernel machine technique, where a nonlinear kernel machine technique, kernel principal component analysis, is used with a specifically identified kernel function to differentiate BOLD signal from the physiological noise of the frequency. The proposed method was evaluated in human fMRI data acquired from multiple task-related and resting state fMRI experiments. A comparison study was also performed with an existing adaptive filtering method. The results indicate that the proposed method can effectively identify and reduce the physiological noise in fMRI data. The comparison study shows that the proposed method can provide comparable or better noise removal performance than the adaptive filtering approach.     

Comparison between end-tidal CO₂ and respiration volume per time for detecting BOLD signal fluctuations during paced hyperventilation

Respiratory motion and capnometry monitoring were performed during blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) of the brain while a series of paced hyperventilation tasks were performed that caused significant hypocapnia. Respiration volume per time (RVT) and end-tidal carbon dioxide (ETCO(2)) were determined and compared for their ability to explain BOLD contrast changes in the data. A 35% decrease in ETCO(2) was observed along with corresponding changes in RVT. A best-fit ETCO(2) response function, with an average initial peak delay time of 12 s, was empirically determined. ETCO(2) data convolved with this response function was more strongly and prevalently correlated to BOLD signal changes than RVT data convolved with the corresponding respiration response function. The results suggest that ETCO(2) better models BOLD signal fluctuations in fMRI experiments with significant transient hypocapnia. This is due to hysteresis in the ETCO(2) response when moving from hypocapnia to normocapnia, compared to moving from normocapnia to hypocapnia.     

Assessment of spatial BOLD sensitivity variations in fMRI using gradient-echo field maps

Clinical blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is becoming increasingly valuable in, e.g., presurgical planning, but the commonly used gradient-echo echo-planar imaging (GE-EPI) technique is sometimes hampered by macroscopic field inhomogeneities. This can affect the degree of signal change that will occur in the GE-EPI images as a response to neural activation and the subsequent blood oxygenation changes, i.e., the BOLD sensitivity (BS). In this study, quantitative BS maps were calculated directly from gradient-echo field maps obtainable on most clinical scanners. In order to validate the accuracy of the calculated BS-maps, known shim gradients were applied and field maps and GE-EPI images of a phantom were acquired. Measured GE-EPI image intensity was then compared with the calculated (predicted) image intensity (pII) which was obtained from the field maps using theoretical expressions for image-intensity loss. The validated expressions for pII were used to calculate the corresponding predicted BOLD sensitivity (pBS) maps in healthy volunteers. Since the field map is assumed to be valid throughout an entire fMRI experiment, the influence of subject motion on the pBS maps was also assessed. To demonstrate the usefulness of such maps, pBS was investigated for clinically important functional areas including hippocampus, Broca's area and primary motor cortex. A systematic left/right pBS difference was observed in Broca's area and in the hippocampus, most likely due to magnetic field inhomogeneity of the particular MRI-system used in this study. For all subjects, the hippocampus showed pBS values above unity with a clear anterior–posterior gradient and with an abrupt drop to zero pBS in the anterior parts of hippocampus. It is concluded that GE field maps can be used to accurately predict BOLD sensitivity and that this parameter is useful to assess spatial variations which will influence fMRI experiments.     

Integration of EEG source imaging and fMRI during continuous viewing of natural movies

Electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are noninvasive neuroimaging tools which can be used to measure brain activity with excellent temporal and spatial resolution, respectively. By combining the neural and hemodynamic recordings from these modalities, we can gain better insight into how and where the brain processes complex stimuli, which may be especially useful in patients with different neural diseases. However, due to their vastly different spatial and temporal resolutions, the integration of EEG and fMRI recordings is not always straightforward. One fundamental obstacle has been that paradigms used for EEG experiments usually rely on event-related paradigms, while fMRI is not limited in this regard. Therefore, here we ask whether one can reliably localize stimulus-driven EEG activity using the continuously varying feature intensities occurring in natural movie stimuli presented over relatively long periods of time. Specifically, we asked whether stimulus-driven aspects in the EEG signal would be co-localized with the corresponding stimulus-driven BOLD signal during free viewing of a movie. Secondly, we wanted to integrate the EEG signal directly with the BOLD signal, by estimating the underlying impulse response function (IRF) that relates the BOLD signal to the underlying current density in the primary visual area (V1). We made sequential fMRI and 64-channel EEG recordings in seven subjects who passively watched 2-min-long segments of a James Bond movie. To analyze EEG data in this natural setting, we developed a method based on independent component analysis (ICA) to reject EEG artifacts due to blinks, subject movement, etc., in a way unbiased by human judgment. We then calculated the EEG source strength of this artifact-free data at each time point of the movie within the entire brain volume using low-resolution electromagnetic tomography (LORETA). This provided for every voxel in the brain (i.e., in 3D space) an estimate of the current density at every time point. We then carried out a correlation between the time series of visual contrast changes in the movie with that of EEG voxels. We found the most significant correlations in visual area V1, just as seen in previous fMRI studies (Bartels A, Zeki, S, Logothetis NK. Natural vision reveals regional specialization to local motion and to contrast-invariant, global flow in the human brain. Cereb Cortex 2008;18(3):705–717), but on the time scale of milliseconds rather than of seconds. To obtain an estimate of how the EEG signal relates to the BOLD signal, we calculated the IRF between the BOLD signal and the estimated current density in area V1. We found that this IRF was very similar to that observed using combined intracortical recordings and fMRI experiments in nonhuman primates. Taken together, these findings open a new approach to noninvasive mapping of the brain. It allows, firstly, the localization of feature-selective brain areas during natural viewing conditions with the temporal resolution of EEG. Secondly, it provides a tool to assess EEG/BOLD transfer functions during processing of more natural stimuli. This is especially useful in combined EEG/fMRI experiments, where one can now potentially study neural-hemodynamic relationships across the whole brain volume in a noninvasive manner.     

Direct measurement of oxygen extraction with fMRI using 6% CO2 inhalation

The blood-oxygenation-level-dependent (BOLD) signal is an indirect hemodynamic signal that is sensitive to cerebral blood flow (CBF), cerebral blood volume (CBV) and cerebral metabolic rate of oxygen. Therefore, the BOLD signal amplitude and dynamics cannot be interpreted unambiguously without additional physiological measurements, and thus, there remains a need for a functional magnetic resonance imaging (fMRI) signal, which is more closely related to the underlying neuronal activity. In this study, we measured CBF with continuous arterial spin labeling, CBV with an exogenous contrast agent and BOLD combined with intracortical electrophysiological recording in the primary visual cortex of the anesthetized monkey. During inhalation of 6% CO2, it was observed that CBF and CBV are not further increased by a visual stimulus, although baseline CBF for 6% CO2 is below the maximal value of CBF. In contrast, the electrophysiological response to the stimulation was found to be preserved during hypercapnia. As a consequence, the simultaneously measured BOLD signal responds negatively to a visual stimulation for 6% CO2 inhalation in the same voxels responding positively during normocapnia. These observations suggest that the fMRI response to a sensory stimulus for 6% CO2 inhalation occurs in the absence of a hemodynamic response, and it therefore directly reflects oxygen extraction into the tissue.     

BOLD signal and vessel dynamics: a hierarchical cluster analysis

The aim of the present study was to analyze blood oxygenation level-dependent (BOLD) signal variation during an apnea-based task in order to assess the capability of a functional magnetic resonance imaging (fMRI) procedure to estimate cerebral vascular dynamic effects. We measured BOLD contrast by hierarchical cluster analysis in healthy subjects undergoing an fMRI experiment, in which the task paradigm was one phase of inspirational apnea (IA). By processing the time courses of the fMRI experiment, analysis was performed only on a subclass of all the possible signal patterns; basically, root mean square and absolute variation differences have been calculated. Considering the baseline value obtained by computing the mean value of the initial rest period as reference, particular voxels showed relative important variations during the IA task and during the recovery phase following the IA. We focused our interest on the signal response of voxels that would correspond mainly to white and gray matter regions and that also may be affected by the proximity of large venous vessels. The results are presented as maps of space-temporal distribution of time series variations with two levels of hierarchical clustering among voxels with low to high initial response. Furthermore, we have presented a clustering of the signal response delay, conducting to a partition and identification of specified brain sites.     

Inhomogeneity-free heteronuclear iMQC

Intermolecular dipolar interactions between proton and carbon spins can be used to indirectly detect carbon spectra with high sensitivity. In this communication, we present a modified sequence that, in addition to the high sensitivity of heteronuclear intermolecular multiple quantum coherence (iMQC) experiments, retains the line narrowing capability characteristic of homonuclear zero-quantum coherences. We demonstrate that this sequence can be used to obtain high resolution (13)C spectra in the presence of magnetic field inhomogeneities, both for thermal and hyperpolarized samples, and discuss applications to water-hyperpolarized carbon imaging.     

Component structure of event-related fMRI responses in the different neurovascular compartments

In most functional magnetic resonance imaging (fMRI) studies, brain activity is localized by observing changes in the blood oxygenation level-dependent (BOLD) signal that are believed to arise from capillaries, venules and veins in and around the active neuronal population. However, the contribution from veins can be relatively far downstream from active neurons, thereby limiting the ability of BOLD imaging methods to precisely pinpoint neural generators. Hemodynamic measures based on apparent diffusion coefficients (ADCs) have recently been used to identify more upstream functional blood flow changes in the capillaries, arterioles and arteries. In particular, we recently showed that, due to the complementary vascular sensitivities of ADC and BOLD signals, the voxels conjointly activated by both measures may identify the capillary networks of the active neuronal areas. In this study, we first used simultaneously acquired ADC and BOLD functional imaging signals to identify brain voxels activated by ADC only, by both ADC and BOLD and by BOLD only, thereby delineating voxels relatively dominated by the arterial, capillary, and draining venous neurovascular compartments, respectively. We then examined the event-related fMRI BOLD responses in each of these delineated neurovascular compartments, hypothesizing that their event-related responses would show different temporal componentries. In the regions activated by both the BOLD and ADC contrasts, but not in the BOLD-only areas, we observed an initial transient signal reduction (an initial dip), consistent with the local production of deoxyhemoglobin by the active neuronal population. In addition, the BOLD-ADC overlap areas and the BOLD-only areas showed a clear poststimulus undershoot, whereas the compartment activated by only ADC did not show this component. These results indicate that using ADC contrast in conjunction with BOLD imaging can help delineate the various neurovascular compartments, improve the localization of active neural populations, and provide insight into the physiological mechanisms underlying the hemodynamic signals.     

Functional magnetic resonance imaging based on SEEP contrast: response function and anatomical specificity

Functional magnetic resonance imaging based on a non-BOLD-contrast mechanism, which we have termed "SEEP" (Signal Enhancement by Extravascular water Protons), has previously been demonstrated. Here the reproducibility of areas of activity identified with both SEEP and BOLD contrast is assessed in duplicate experiments in healthy volunteers, with relatively high resolution (1.6 mm) image data at 1.5 T. These areas of activity are equally well localized to voxels containing primarily gray matter with the two contrast mechanisms. As in previous studies, areas of SEEP activity are observed to be immediately adjacent to areas of BOLD activity, with very little overlap. The response functions were estimated for both SEEP and BOLD contrast, and are observed to be distinct. The peak SEEP response is observed to lag the BOLD response by approximately 1 s and to decay more slowly with no poststimulus undershoot. Average BOLD signal changes (GE-EPI, TE=50 ms) were observed to be 3.4+/-1.2% (mean+/-S.D.), whereas SEEP signal changes (SE-EPI, TE=23 ms) were 1.9+/-0.5%, consistent with previous studies carried out at 0.2 and 3 T. These observations provide further support for the existence of a non-BOLD-contrast mechanism for fMRI, based on changes in extravascular spin density.     

Neonatal auditory activation detected by functional magnetic resonance imaging

The objective of this study was to detect auditory cortical activation in non-sedated neonates employing functional magnetic resonance imaging (fMRI). Using echo-planar functional brain imaging, subjects were presented with a frequency-modulated pure tone; the BOLD signal response was mapped in 5 mm-thick slices running parallel to the superior temporal gyrus. Twenty healthy neonates (13 term, 7 preterm) at term and 4 adult control subjects. Blood oxygen level-dependent (BOLD) signal in response to auditory stimulus was detected in all 4 adults and in 14 of the 20 neonates. FMRI studies of adult subjects demonstrated increased signal in the superior temporal regions during auditory stimulation. In contrast, signal decreases were detected during auditory stimulation in 9 of 14 newborns with BOLD response. fMRI can be used to detect brain activation with auditory stimulation in human infants.