Nuclear magnetic resonance signal from flowing nuclei in rapid imaging using gradient echoes

Real-time quantification of in vivo cerebrospinal fluid velocity using the functional magnetic resonance imaging inflow effect

In vivo estimation of cerebrospinal fluid (CSF) velocity is crucial for understanding the glymphatic system and its potential role in neurodegenerative disorders such as Alzheimer's disease and Parkinson's disease. Current cardiac or respiratory-gated approaches, such as 4D flow magnetic resonance imaging (MRI), cannot capture CSF movement in real time because of limited temporal resolution and, in addition, deteriorate in accuracy at low fluid velocities. Other techniques like real-time phase-contrast-MRI or time-spatial labeling inversion pulse are not limited by temporal averaging but have limited availability, even in research settings. This study aims to quantify the inflow effect of dynamic CSF motion on functional MRI (fMRI) for in vivo, real-time measurement of CSF flow velocity. We considered linear and nonlinear models of velocity waveforms and empirically fit them to fMRI data from a controlled flow experiment. To assess the utility of this methodology in human data, CSF flow velocities were computed from fMRI data acquired in eight healthy volunteers. Breath-holding regimens were used to amplify CSF flow oscillations. Our experimental flow study revealed that CSF velocity is nonlinearly related to inflow effect-mediated signal increase and well estimated using an extension of a previous nonlinear framework. Using this relationship, we recovered velocity from in vivo fMRI signal, demonstrating the potential of our approach for estimating CSF flow velocity in the human brain. This novel method could serve as an alternative approach to quantifying slow flow velocities in real time, such as CSF flow in the ventricular system, thereby providing valuable insights into the glymphatic system's function and its implications for neurological disorders.

Real-time Quantification of in vivo cerebrospinal fluid velocity using fMRI inflow effect

In vivo estimation of cerebrospinal fluid (CSF) velocity is crucial for understanding the glymphatic system and its potential role in neurodegenerative disorders such as Alzheimer's disease and Parkinson's disease. Current cardiac or respiratory gated approaches, such as 4D flow MRI, cannot capture CSF movement in real time due to limited temporal resolution and in addition deteriorate in accuracy at low fluid velocities. Other techniques like real-time PC-MRI or time-spatial labeling inversion pulse are not limited by temporal averaging but have limited availability even in research settings. This study aims to quantify the inflow effect of dynamic CSF motion on functional magnetic resonance imaging (fMRI) for in vivo, real-time measurement of CSF flow velocity. We considered linear and nonlinear models of velocity waveforms and empirically fit them to fMRI data from a controlled flow experiment. To assess the utility of this methodology in human data, CSF flow velocities were computed from fMRI data acquired in eight healthy volunteers. Breath holding regimens were used to amplify CSF flow oscillations. Our experimental flow study revealed that CSF velocity is nonlinearly related to inflow effect-mediated signal increase and well estimated using an extension of a previous nonlinear framework. Using this relationship, we recovered velocity from in vivo fMRI signal, demonstrating the potential of our approach for estimating CSF flow velocity in the human brain. This novel method could serve as an alternative approach to quantifying slow flow velocities in real time, such as CSF flow in the ventricular system, thereby providing valuable insights into the glymphatic system's function and its implications for neurological disorders.

Measuring Arterial Pulsatility With Dynamic Inflow Magnitude Contrast

The pulsatility of blood flow through cerebral arteries is clinically important, as it is intrinsically associated with cerebrovascular health. In this study we outline a new MRI approach to measuring the real-time pulsatile flow in cerebral arteries, which is based on the inflow phenomenon associated with fast gradient-recalled-echo acquisitions. Unlike traditional phase-contrast techniques, this new method, which we dub dynamic inflow magnitude contrast (DIMAC), does not require velocity-encoding gradients as sensitivity to flow velocity is derived purely from the inflow effect. We achieved this using a highly accelerated single slice EPI acquisition with a very short TR (15 ms) and a 90° flip angle, thus maximizing inflow contrast. We simulate the spoiled GRE signal in the presence of large arteries and perform a sensitivity analysis. The sensitivity analysis demonstrates that in the regime of high inflow contrast, DIMAC shows much greater sensitivity to flow velocity over blood volume changes. We support this theoretical prediction with in-vivo data collected in two separate experiments designed to demonstrate the utility of the DIMAC signal contrast. We perform a hypercapnia challenge experiment in order to experimentally modulate arterial tone within subjects, and thus modulate the arterial pulsatile flow waveform. We also perform a thigh-cuff release challenge, designed to induce a transient drop in blood pressure, and demonstrate that the continuous DIMAC signal captures the complex transient change in the pulsatile and non-pulsatile components of flow. In summary, this study proposes a new role for a well-established source of MR image contrast and demonstrates its potential for measuring both steady-state and dynamic changes in arterial tone.

Assessment of intravoxel incoherent motion MRI with an artificial capillary network: analysis of biexponential and phase-distribution models

Purpose

To systematically analyze intravoxel incoherent motion (IVIM) MRI in a perfusable capillary phantom closely matching the geometry of capillary beds in vivo and to compare the validity of the biexponential pseudo‐diffusion and the recently introduced phase‐distribution IVIM model.

Methods

IVIM‐MRI was performed at 12 different flow rates (0.2⋯2.4mL/min) in a capillary phantom using 4 different DW‐MRI sequences (2 with monopolar and 2 with flow‐compensated diffusion‐gradient schemes, with up to 16b values between 0 and 800s/mm2). Resulting parameters from the assessed IVIM models were compared to results from optical microscopy.

Results

The acquired data were best described by a static and a flowing compartment modeled by the phase‐distribution approach. The estimated signal fraction f of the flowing compartment stayed approximately constant over the applied flow rates, with an average of f=0.451±0.023 in excellent agreement with optical microscopy (f=0.454±0.002). The estimated average particle flow speeds v=0.25⋯2.7mm/s showed a highly significant linear correlation to the applied flow. The estimated capillary segment length of approximately 189um agreed well with optical microscopy measurements. Using the biexponential model, the signal fraction f was substantially underestimated and displayed a strong dependence on the applied flow rate.

Conclusion

The constructed phantom facilitated the detailed investigation of IVIM‐MRI methods. The results demonstrate that the phase‐distribution method is capable of accurately characterizing fluid flow inside a capillary network. Parameters estimated using the biexponential model, specifically the perfusion fraction f, showed a substantial bias because the model assumptions were not met by the underlying flow pattern.

Cine-MR Imaging of Aqueductal CSF Flow in Normal Pressure Hydrocephalus Syndrome before and after CSF Shunt

Reproducibility of the aqueductal CSF signal intensity on a gradient echo cine-MR sequence exploiting through plane inflow enhancement was tested in 11 patients with normal or dilated ventricles. Seven patients with normal pressure hydrocephalus (NPH) syndrome were investigated with the sequence before and after CSF shunting. Two patients exhibiting central flow void within a hyperintense aqueductal CSF improved after surgery and the flow void disappeared after shunting. One patient with increased maximum and minimum aqueductal CSF signal as compared to 18 healthy controls also improved and the aqueductal CSF signal was considerably decreased after shunting. Three patients with aqueductal CSF values similar to those in the controls did not improve, notwithstanding their maximum aqueductal CSF signals decreasing slightly after shunting. No appreciable aqueductal CSF flow related enhancement consistent with non-communicating hydrocephalus was found in the last NPH patient who improved after surgery. Cine-MR with inflow technique yields a reproducible evaluation of flow-related aqueductal CSF signal changes which might help in identifying shunt responsive NPH patients. These are likely to be those with hyperdynamic aqueductal CSF or aqueductal obstruction.

Mr Imaging, Flow and Motion

The present work is intended as a nonmathematical review of the role of flow and motion in nuclear magnetic resonance (MR) imaging. A historical review of MR flow measurement techniques is given, followed by a short overview of flow models in vitro and in vivo. The theory behind the influence of motion on the modulus and phase MR signal information is discussed and effects such as washin/washout, flow-induced signal void, phase offset, and phase dispersion are defined. A simple approach to the concept of MR angiography is given, and methods for quantitative flow measurements such as the phase mapping technique, are surveyed. Aspects of the measurement of diffusion and mirocirculation are given, and finally, an overview of the role of MR flow imaging in present and future clinical application is given.

Chemical shift-induced phase errors in phase-contrast MRI

Phase-contrast magnetic resonance imaging is subject to numerous sources of error, which decrease clinical confidence in the reported measures. This work outlines how stationary perivascular fat can impart a significant chemical shift induced phase-contrast magnetic resonance imaging measurement error using computational simulations, in vitro, and in vivo experiments. This chemical shift error does not subtract in phase difference processing, but can be minimized with proper parameter selection. The chemical shift induced phase errors largely depend on both the receiver bandwidth and the TE. Both theory and an in vivo comparison of the maximum difference in net forward flow between vessels with and without perivascular fat indicated that the effects of chemically shifted perivascular fat are minimized by the use of high bandwidth (814 Hz/px) and an in-phase TE (high BW-TE(IN)). In healthy volunteers (N = 10) high BW-TE(IN) significantly improves intrapatient net forward flow agreement compared with low bandwidth (401 Hz/px) and a mid-phase TE as indicated by significantly decreased measurement biases and limits of agreement for the ascending aorta (1.8 ± 0.5 mL vs. 6.4 ± 2.8 mL, P = 0.01), main pulmonary artery (2.0 ± 0.9 mL vs. 11.9 ± 5.8 mL, P = 0.04), the left pulmonary artery (1.3 ± 0.9 mL vs. 5.4 ± 2.5 mL, P = 0.003), and all vessels (1.7 ± 0.8 mL vs. 7.2 ± 4.4 mL, P = 0.001).

Inflow effects on functional MRI

Blood inflow from the upstream has contribution or contamination to the blood oxygen level-dependent (BOLD) functional signal both in its magnitude and time courses. During neuronal activations, regional blood flow velocity increases which results in increased fMRI signals near the macrovasculatures. The inflow effects are dependent on RF pulse history, slice geometry, flow velocity, blood relaxation times and imaging parameters. In general, the effect is stronger with more T(1) weighting in the signal, e.g. by using a short repetition time and a large flip angle. This article reviews the basic principle of the inflow effects, its appearances in conventional GRE, fast spin-echo (FSE) and echo-planar imaging (EPI) acquisitions, methods for separating the inflow from the BOLD effect as well as the interplay between imaging parameters and other physiological factors with the inflow effects in fMRI. Based on theoretical derivation and human experiments, the inflow effects have been shown to contribute significantly in conventional GRE but negligible in FSE acquisitions. For gradient-echo EPI experiments, the blood inflow could modulate both amplitude and the temporal information of the fMRI signal, depending on the imaging parameters and settings.

Neurophysiologic correlates of fMRI in human motor cortex

The neurophysiological underpinnings of functional magnetic resonance imaging (fMRI) are not well understood. To understand the relationship between the fMRI blood oxygen level dependent (BOLD) signal and neurophysiology across large areas of cortex, we compared task related BOLD change during simple finger movement to brain surface electric potentials measured on a similar spatial scale using electrocorticography (ECoG). We found that spectral power increases in high frequencies (65-95 Hz), which have been related to local neuronal activity, colocalized with spatially focal BOLD peaks on primary sensorimotor areas. Independent of high frequencies, decreases in low frequency rhythms (<30 Hz), thought to reflect an aspect of cortical-subcortical interaction, colocalized with weaker BOLD signal increase. A spatial regression analysis showed that there was a direct correlation between the amplitude of the task induced BOLD change on different areas of primary sensorimotor cortex and the amplitude of the high frequency change. Low frequency change explained an additional, different part of the spatial BOLD variance. Together, these spectral power changes explained a significant 36% of the spatial variance in the BOLD signal change (R(2) = 0.36). These results suggest that BOLD signal change is largely induced by two separate neurophysiological mechanisms, one being spatially focal neuronal processing and the other spatially distributed low frequency rhythms.

Dynamic magnetic resonance imaging of cerebral blood flow using arterial spin labeling

Modern functional neuroimaging techniques, including positron emission tomography, optical imaging of intrinsic signals, and magnetic resonance imaging (MRI) rely on a tight coupling between neural activity and cerebral blood flow (CBF) to visualize brain activity using CBF as a surrogate marker. Because the spatial and temporal resolution of neuroimaging modalities is ultimately determined by the spatial and temporal specificity of the underlying hemodynamic signals, characterization of the spatial and temporal profiles of the hemodynamic response to focal brain stimulation is of paramount importance for the correct interpretation and quantification of functional data. The ability to properly measure and quantify CBF with MRI is a major determinant of progress in our understanding of brain function. We review the dynamic arterial spin labeling (DASL) method to measure CBF and the CBF functional response with high temporal resolution.

T1 measurement of flowing blood and arterial input function determination for quantitative 3D T1-weighted DCE-MRI

PURPOSE

To propose a simple, accurate method for measuring T(1) in flowing blood and the arterial input function (AIF), and to evaluate the impact on dynamic contrast-enhanced MRI (DCE-MRI) quantification of pharmacokinetic parameters.

MATERIALS AND METHODS

A total of 10 rabbits were scanned at 1.5 Tesla and administered a bolus of Gadomer. Preinjection T(1) and AIF measurements were acquired in the iliac arteries using a rapid three-dimensional (3D) spoiled gradient recalled echo (SPGR) approach. Correction was made for imperfect B(1) fields, in-flow, and partial volume effects. DCE-MRI parameters blood volume (v(b)) and endothelial transfer constant (K(trans)) in resting skeletal muscle were estimated from pharmacokinetic analysis using individually measured AIFs. Literature comparisons were made to assess accuracy.

RESULTS

Blood T(1) was more accurate and precise after correction for B(1) and partial volume errors (1267 +/- 72 msec). Measured AIFs followed reported biexponential decay characteristics for Gadomer clearance in rabbits. Parameters v(b) (2.47 +/- 0.65%) and K(trans) (3.6 +/- 1.0 x 10(-3) minute(-1)) derived from AIFs based on corrected blood T(1)s were more reproducible and in better agreement with literature values.

CONCLUSION

The proposed method enables accurate in vivo blood T(1) and AIF measurements and can be easily implemented in a range of DCE-MRI applications to improve both the accuracy and reproducibility of pharmacokinetic parameters.

Rapid quantitation of cardiovascular flow using slice-selective fourier velocity encoding with spiral readouts

Accurate flow visualization and quantitation is important for the assessment of many cardiovascular conditions such as valvular stenosis and regurgitation. Phase contrast based methods experience partial volume artifacts when flow is highly localized, complex and/or turbulent. Fourier velocity encoding (FVE) avoids such problems by resolving the full velocity distribution within each voxel. This work proposes the use of slice selective FVE with spiral readouts to acquire fully localized velocity distributions in a short breath-hold. Scan-plane prescription is performed using classic protocols, and an automatic algorithm is used for in-plane localization of the flow. Time and spatially-resolved aortic valve velocity distributions with 26-msec temporal resolution and 25 cm/sec velocity resolution over a 600 cm/sec field-of-view were acquired in a 12-heartbeat breath-hold. In carotid studies, scan time was extended to achieve higher spatial resolution. The method was demonstrated in healthy volunteers and patients, and the results compared qualitatively well with Doppler ultrasound. Acquisition time could be reduced to 7 heartbeats (a 42% reduction) using partial Fourier reconstruction along the velocity dimension.

Measurements of tissue T1 spin-lattice relaxation time and discrimination of large draining veins using transient EPI data sets in BOLD-weighted fMRI acquisitions

The signal intensity during the dynamic approach to the equilibrium state of longitudinal magnetization is a function of sequence parameters, such as repetition time and flip angle, and depends on tissue characteristics, including longitudinal relaxation time of stationary tissue and the rate of blood inflow. A method is presented to extract information from data acquired during the transient state prior to T1 equilibrium using echo-planar acquisitions in T2*-weighted functional magnetic resonance imaging (fMRI) experiments. A voxel in a single slice acquisition is assumed to contain either stationary tissue or large vessels with flowing blood. Models are presented to characterize longitudinal magnetization relaxation of heterogeneous stationary tissue and blood inflow. The data were fitted to theoretical models for longitudinal relaxation of stationary tissue and inflowing blood assuming no residual signal prior to each RF excitation. Parameters were estimated at 3 T for each model using least squares estimation. A goodness-of-fit criterion was applied to exclude voxels that have transient data that does not fit the selected (best fit) model. Voxels that best fit the inflow model, measured at various TR and flip angles, were assumed to contain large draining veins and were excluded from functional maps. Histogram analysis of T1 distributions for activated voxels in a visual paradigm demonstrated the distributions are centered at T1 values of gray matter with tails at both sides of the center due to partial voluming of gray matter with white matter and CSF respectively. The mean gray matter volume fraction in activated voxels was about 0.9. The results indicate that transient data sets can provide additional information that is useful for both localization and characterization of the functionally relevant BOLD response.

A comparison of T2*-weighted magnitude and phase imaging for measuring the arterial input function in the rat aorta following intravenous injection of gadolinium contrast agent

The arterial input function (AIF) is important for quantitative MR imaging perfusion experiments employing Gd contrast agents. This study compared the accuracy of T(2)*-weighted magnitude and phase imaging for noninvasive measurement of the AIF in the rat aorta. Twenty-eight in vivo experiments were performed involving simultaneous arterial blood sampling and MR imaging following Gd injection. In vitro experiments were also performed to confirm the in vivo results. At 1.89 T and TE=3 ms, the relationship between changes in 1/T(2)* in blood (estimated from MR signal magnitude) and Gd concentration ([Gd]) was measured to be approximately 19 s(-1) mM(-1), while that between phase and [Gd] was approximately 0.19 rad mM(-1). Both of these values are consistent with previously published results. The in vivo phase data had approximately half as much scatter with respect to [Gd] than the in vivo magnitude data (r(2)=.34 vs. r(2)=.17, respectively). This is likely due to the fact that the estimated change in 1/T(2)* is more sensitive than the phase to a variety of factors such as partial volume effects and T(1) weighting. Therefore, this study indicates that phase imaging may be a preferred method for measuring the AIF in the rat aorta compared to T(2)*-weighted magnitude imaging.

Dynamic imaging of perfusion and oxygenation by functional magnetic resonance imaging

Cerebral blood flow can be measured with magnetic resonance imaging (MRI) by arterial spin labeling techniques, where magnetic labeling of flowing spins in arterial blood water functions as the endogenous tracer upon mixing with the unlabeled stationary spins of tissue water. The consequence is that the apparent longitudinal relaxation time (T1) of tissue water is attenuated. A modified functional MRI scheme for dynamic CBF measurement is proposed that depends on extraction of T1 weighting from the blood oxygenation level-dependent (BOLD) image contrast, because the functional MRI signal also has an intrinsic T1 weighting that can be altered by variations of the excitation flip angle. In the alpha-chloralose-anesthetized rat model at 7T, the authors show that the stimulation-induced BOLD signal change measured with two different flip angles can be combined to obtain a T1-weighted MRI signal, reflecting the magnitude of the CBF change, which can be deconvolved to obtain dynamic changes in CBF. The deconvolution of the T1-weighted MRI signal, which is a necessary step for accurate reflection of the dynamic changes in CBF, was made possible by a transfer function obtained from parallel laser-Doppler flowmetry experiments. For all stimulus durations (ranging from 4 to 32 seconds), the peak CBF response measured by MRI after the deconvolution was reached at 4.5 +/- 1.0 seconds, which is in good agreement with (present and prior) laser-Doppler measurements. Because the low flip angle data can also provide dynamic changes of the conventional BOLD image contrast, this method can be used for simultaneous imaging of CBF and BOLD dynamics.

Time-of-flight quantitative measurements of blood flow in mouse hindlimbs

PURPOSE

To evaluate the feasibility of using time-of-flight (TOF) imaging to directly measure hindlimb blood flow in a mouse model of peripheral vascular disease.

MATERIALS AND METHODS

Four tubes were imaged simultaneously (diameters = 0.39 mm, 0.59 mm, and two at 1.46 mm) with a 1.0 mM copper sulfate solution for 19 flow velocities. In vivo measurements were performed in the hindlimbs of three mouse strains-C57BL/6 (N = 5), BALB/c (N = 5), and 129S2/Sv (N = 5)-three weeks after femoral artery ligation with a calibration standard.

RESULTS

The flow phantom showed that the intensity was linear (r2 = 0.92) over the pertinent blood flow velocities in the mouse hindlimbs. Measurements of the blood flow in the distal hindlimbs in different strains of mice (combination of both the venous and arterial flows) were obtained 21 days after right-sided femoral artery occlusion. The results showed that under similar conditions of anesthesia and temperature, SV129 mice on the nonligated side had the highest flows (0.50 +/- 0.07 mL/minute), followed by C57BL/6 (0.28 +/- 0.04 mL/minute) and BALB/c (0.23 +/- 0.05 mL/minute), P < 0.02. The ligated side measurements (SV129, 0.31 +/- 0.05 mL/minute (P = 0.02); C57BL/6, 0.21 +/- 0.02 mL/minute (P = 0.13); and BALB/c, 0.12 +/- 0.02 mL/minute (P= 0.06)) showed a trend in BALB/c and C57BL/6 and significant differences in SV129 for incomplete recovery three weeks after surgery, compared to the nonligated side.

CONCLUSION

Two-dimensional TOF imaging permits quantitative in vivo measurements of hindlimb blood flow in a mouse model of peripheral vascular disease without the need of contrast injection, offering advantages of serial imaging not limited by tissue penetration.

Contrast-enhanced magnetic resonance angiography: development and optimization of techniques for paramagnetic and hyperpolarized contrast media

Contrast-enhanced magnetic resonance angiography (CE-MRA) is a diagnostic method for imaging of vascular structures based on nuclear magnetic resonance. Vascular enhancement is achieved by injection of a contrast medium (CM). Studies were performed using two different types of CM: conventional paramagnetic CM, and a new type of CM based on hyperpolarized (HP) nuclei. The effects of varying CM concentration with time during image acquisition were studied by means of computer simulations using two different models. It was shown that a rapid concentration variation during encoding of the central parts of k-space could result in signal loss and severe image artifacts. The results were confirmed qualitatively with phantom experiments. A postprocessing method was developed to address problems with simultaneous enhancement of arteries and veins in CE-MRA of the lower extremities. The method was based on the difference in flow-induced phase in the two vessel types. Evaluation of the method was performed with flow phantom measurements and with CE-MRA in two volunteers using standard pulse sequences. The flow-induced phase in the vessels of interest was sufficient to distinguish arteries from veins in the superior-inferior direction. Using this method, the venous enhancement could be extinguished. The possibility of using HP nuclei as CM for CE-MRA was evaluated. Signal expressions for a flow of HP CM imaged with a gradient echo sequence were derived. These signal expressions were confirmed in phantom experiments using HP 129Xe dissolved in ethanol. Studies were also performed with a new CM based on HP 13C. The CM had very long relaxation times (T1, in vivo/T2, in vivo approximately 38/1.3 s). The long relaxation times were utilized in imaging with a fully balanced steady-state free precession pulse sequence (trueFISP), where the optimal flip angle was found to be 180 degrees. CE-MRA with the 13C-based CM in rats resulted in images with high vascular SNR (approximately 500). CE-MRA is a useful clinical tool for diagnosing vascular disease. With the development of new contrast media, based on hyperpolarized nuclei for example, there is a potential for further improvement in the signal levels that can be achieved, enabling a standard of imaging of vessels that is not possible today.

Gradient echo imaging of flowing hyperpolarized nuclei: theory and phantom studies on 129Xe dissolved in ethanol

The influence of flip angle and flow velocity on the signal intensity achieved when imaging a hyperpolarized substance with a spoiled gradient echo sequence was investigated. The study was performed both theoretically and experimentally using hyperpolarized xenon dissolved in ethanol. Analytical expressions regarding the optimal flip angle with respect to signal and the corresponding signal level are presented and comparisons with thermally polarized substances are made. Both experimentally and theoretically, the optimal flip angle was found to increase with increasing flow velocity. Numerical calculations showed that the velocity dependence of the signal differs between the cases of hyperpolarized and thermally polarized substances.

Inversion profiles of adiabatic inversion pulses for flowing spins: the effects on labeling efficiency and labeling accuracy in perfusion imaging with pulsed arterial spin-labeling

The inversion profile of adiabatic inversion pulses is essential to the accuracy of perfusion measurement with pulsed arterial spin-labeling (ASL). In this paper, the inversion profiles for flowing spins were investigated using a numerical solution of the modified Bloch equations including a term for moving spins. Inversion profiles for spins flowing at a constant or varying velocity were examined for hyperbolic secant (HS) and frequency-offset corrected inversion (FOCI) pulses. Distortions of the inversion profiles were found for both pulses with spins flowing within physiological velocity range. The effects of the distorted profiles on labeling efficiency and labeling accuracy in the application of pulsed ASL perfusion imaging were analyzed. These effects should be taken into account in ASL techniques, in order to obtain robust and accurate perfusion measurements.

MAGNETIC RESONANCE SIGNAL LOSS IN TURBULENT SHEAR FLOW

NMR signal loss due to turbulent shear flow is discussed, and a general expression for the phase fluctuation is derived. In the presence of flow shear, the velocity fluctuation perpendicular to the direction of magnetic gradient and the Reynolds stress can cause loss of MR signal Most of signal loss results from the boundary layer, where the flow shear is strong in turbulent pipe flaw, Half the signal loss within the mixing layer distal to a moderate stenosis is caused by the velocity fluctuation in the direction of magnetic gradient, while the remaining results from the velocity, fluctuation perpendicular to the magnetic gradient. The use of eddy diffusivity for the description of signal dephasing in a spin echo sequence is also addressed; A positive, constant eddy diffusivity can not describe the temporal change of phase fluctuation correctly.

Flow quantification using low-spatial-resolution and low-velocity-resolution velocity images

In this report, a flow-quantification method using Fourier velocity encoding (FVE) with limited spatial and velocity resolution is presented. The total flow rate in a vessel corresponds to the first moment of the velocity histogram of spins in the vessel, whereas the spin density of flowing spins is the normalization constant. Because the measured histogram using FVE is distorted by RF saturation effects, the RF saturation effects are first estimated and then accurately compensated by acquiring five velocity-encoded images. The spatial resolution in each image can be relatively low because all stationary spins vanish in the resultant flow map. In a phantom study, the errors in measured flow rates were within +/-10% even when the pixel size was greater than the vessel size. This method was also successfully applied to measure flow in the femoral artery. In general, this method constitutes a basis for analyzing multiple velocity-encoded images and is particularly useful for quantifying slow flow or flow in small vessels. Magn Reson Med 42:682-690, 1999.

Pseudo-continuous arterial spin labeling technique for measuring CBF dynamics with high temporal resolution

Cerebral blood flow (CBF) can be measured noninvasively with nuclear magnetic resonance (NMR) by using arterial water as an endogenous perfusion tracer. However, the arterial spin labeling (ASL) techniques suffer from poor temporal resolution due to the need to wait for the exchange of labeled arterial spins with tissue spins to produce contrast. In this work, a new ASL technique is introduced, which allows the measurement of CBF dynamics with high temporal and spatial resolution. This novel method was used in rats to determine the dynamics of CBF changes elicited by somatosensory stimulation with a temporal resolution of 108 ms. The onset time of the CBF response was 0.6 +/- 0.4 sec (mean +/- SD) after onset of stimulation (n = 10). The peak response was observed 4.4 +/- 3.7 sec (mean +/- SD) after stimulation began. These results are in excellent agreement with previous data obtained with invasive techniques, such as laser-Doppler flowmetry and hydrogen clearance, and suggest the appropriateness of this novel technique to probe CBF dynamics in functional and pathological studies with high temporal and spatial resolution. Magn Reson Med 42:425-429, 1999.

Magnetic resonance angiography with hyperpolarized 129Xe dissolved in a lipid emulsion

Hyperpolarized (HP) 129Xe can be dissolved in biologically compatible lipid emulsions while maintaining sufficient polarization for in vivo vascular imaging. For xenon in Intralipid 30%, in vitro spectroscopy at 2 T yielded a chemical shift of 197 +/- 1 ppm with reference to xenon gas, a spin-lattice relaxation time T1 = 25.3 +/- 2.1 sec, and a T2* time constant of 37 +/- 5 msec. Angiograms of the abdominal and pelvic veins in the rat obtained with 129Xe MRI after intravenous injection of HP 129Xe/Intralipid 30% into the tail demonstrated signal-to-noise ratios between 8 and 29. An analysis of the inflow effect on time-of-flight images of two segments of the inferior vena cava yielded additional information. The mean blood flow velocity was 34.7 +/- 1.0 mm/sec between the junction of the caudal veins and the kidneys and 13.3 +/- 0.8 mm/sec at the position of the diaphragm. The mean volume flow rates in these segments were 7.2 +/- 3.4 ml/min and 11.0 +/- 2.8 ml/min, respectively. Intravenous delivery of HP 129Xe dissolved in a carrier may lead to novel biomedical applications of laser-polarized gases.

Effects of slow flow on slice profile and NMR signal in fast imaging sequences

A computer program has been developed to evaluate the selective-slice profiles obtained in the steady state for fast gradient-echo imaging. Both spoiled and refocused gradient-echo pulse sequences have been considered. By numerically solving the Bloch equations modified for the effects of flow, for a three-dimensional volume of spins, for realistic RF excitations and linear gradient combinations, the program permits the combined effects of flow and imaging variables on the magnetization slice profile to be assessed quantitatively. We have found that the gradient pattern in gradient-echo pulse sequences is a significant factor for determining the steady-state slice profiles and the strength of the NMR signal from the flowing spins.

NMR signal from flowing nuclei in fast gradient-echo pulse sequences with refocusing

A theoretical description of the nuclear magnetic resonance (NMR) signal from flowing nuclei in refocused gradient-echo pulse sequences, both with continuous- and alternating-phase pulse trains, has been developed. Both laminar and plug flow models have been considered and formulae have been derived that relate mean signal intensity to flip angle, pulse sequence repetition interval (TR), and flow velocity. The degree of signal enhancement or reduction in various conditions of flow and pulse sequences depends on the precise phase relationships between the residual transverse magnetization and each radio-frequency (RF) pulse.

Pulse sequences for steady-state saturation of flowing spins

It is useful to be able to suppress the NMR signal from spins in a flowing fluid, for example for "black-blood" visualization of blood vessels in vivo, for the suppression of flow artifacts, and for the estimation of tissue perfusion by continuous labeling of inflowing arterial spins. This work considers the flow of fluid through a region in which it is subjected to a train of saturation pulses. Computer simulations and in vitro measurements show that a train of equal-duration spoiler pulses produces less effective suppression than does a train of pulses of geometrically increasing duration. It is shown analytically that a long train of ideal equal-duration spoiler pulses converts initial magnetization (0, 0, M0) into a combination of longitudinal and transverse magnetization equal to 0. 29 (-M0, 0, M0) and is therefore unsatisfactory for continuous saturation.

On the characteristics of functional magnetic resonance imaging of the brain

In this review we discuss various recent topics that characterize functional magnetic resonance imaging (fMRI). These topics include a brief description of MRI image acquisition, how to cope with noise or signal fluctuation, the basis of fMRI signal changes, and the relation of MRI signal to neuronal events. Several observations of fMRI that show good correlation to the neurofunction are referred to. Temporal characteristics of fMRI signals and examples of how the feature of real time measurement is utilized are then described. The question of spatial resolution of fMRI, which must be dictated by the vascular structure serving the functional system, is discussed based on various fMRI observations. Finally, the advantage of fMRI mapping is shown in a few examples. Reviewing the vast number of recent fMRI application that have now been reported is beyond the scope of this article.

Fuzzy clustering of gradient-echo functional MRI in the human visual cortex. Part II: quantification

Fuzzy cluster analysis (FCA) is a new exploratory method for analyzing fMRI data. Using simulated functional MRI (fMRI) data, the performance of FCA, as implemented in the software package Evident, was tested and a quantitative comparison with correlation analysis is presented. Furthermore, the fMRI model fit allows separation and quantification of flow and blood oxygen level dependent (BOLD) contributions in the human visual cortex. In gradient-recalled echo fMRI at 1.5 T (TR = 60 ms, TE = 42 ms, radiofrequency excitation flip angle [theta] = 10 degrees-60 degrees) total signal enhancement in the human visual cortex, ie, flow-enhanced BOLD plus inflow contributions, on average varies from 5% to 10% in or close to the visual cortex (average cerebral blood volume [CBV] = 4%) and from 100% to 20% in areas containing medium-sized vessels (ie, average CBV = 12% per voxel), respectively. Inflow enhancement, however, is restricted to intravascular space (= CBV) and increases with increasing radiofrequency (RF) flip angle, whereas BOLD contributions may be obtained from a region up to three times larger and, applying an unspoiled gradient-echo (GRE) sequence, also show a flip angle dependency with a minimum at approximately 30 degrees. This result suggests that a localized hemodynamic response from the microvasculature at 1.5 T may be extracted via fuzzy clustering. In summary, fuzzy clustering of fMRI data, as realized in the Evident software, is a robust and efficient method to (a) separate functional brain activation from noise or other sources resulting in time-dependent signal changes as proven by simulated fMRI data analysis and in vivo data from the visual cortex, and (b) allows separation of different levels of activation even if the temporal pattern is indistinguishable. Combining fuzzy cluster separation of brain activation with appropriate model calculations allows quantification of flow and (flow-enhanced) BOLD contributions in areas with different vascularization.

Quantification of signal changes in gradient recalled echo FMRI

Understanding and quantifying the various contributions to functional magnetic resonance imaging (FMRI) signal changes in activated cortical areas is paramount for a clinical application of brain mapping by FMRI. Therefore, all significant contributions to FMRI signal changes, both extra- and intravascular, from macrovessels down to the capillary network, should be taken into account. We present a gradient-recalled-echo FMRI model based on in-flow effects described by the Bloch equations, adding susceptibility effects empirically via T2* differences measured in vitro in human blood samples. Results of these calculations (by systematically varying alpha, echo time (TE), repetition time (TR), as well as blood velocity and T2* upon stimulation) may be used to (a) simulate functional MRI experiments with different measurement protocols and (b) estimate realistic values for important anatomical and physiological details that influence local signal changes in FMRI (i.e., size and distribution of vessels, effective relaxation times of blood, etc.). The excellent agreement between our model calculations and experimental results from conventional gradient recalled echo fMRI in vivo suggests a significant contribution from very slow flow and oxygenation changes, predominantly in small vessels (vblood = 1-4 mm/s). The actual contribution of T1- and T2-related effects is strongly dependent on sequence design and actual sequence parameters used. Thus, the model simulations presented may also be used to optimize measurement protocols for investigating various neurophysiological phenomena.

Steady state effects in fast gradient echo magnetic resonance imaging

Fourier methodology is applied to analyze steady state effects in fast gradient echo imaging. Simulations and MR imaging experiments on phantoms demonstrate the effectiveness of existing schemes under different experimental conditions, and modifications are introduced which result in reduced sensitivity to slow object motion as compared to conventional phase modulation schemes. In addition, phase modulation schemes are introduced which more than double temporal signal stability in human brain scans, at flip angles well above the Ernst angle and TR << T1, T2. Application of phase modulation to the measurement of higher order signals shows increased sensitivity, particularly for the highest order signals.

Magnetization and diffusion effects in NMR imaging of hyperpolarized substances

The special magnetization characteristics of hyperpolarized noble gases have led to an interest in using these agents for new MRI applications. In this note, the magnetization effects and NMR signal dependence of two noble gases, 3He and l29Xe, are modeled across a range of gradient-echo imaging parameters. Pulse-sequence analysis shows a wide variation in optimum flip angles between imaging of gas (e.g., 3He or 129Xe) in air spaces (e.g., trachea and lung) and in blood vessels. To optimize imaging of the air spaces, it is also necessary to reduce the otherwise substantial signal losses from diffusion effects by increasing voxel size. The possibility of using hyperpolarized 129Xe for functional MRI (fMRI) is discussed in view of the results from the blood flow analysis. The short-lived nature of the hyperpolarization opens up new possibilities, as well as new technical challenges, in its potential application as a blood-flow tracer.

Fast 3D functional magnetic resonance imaging at 1.5 T with spiral acquisition

A new method to perform rapid 3D fMRI in human brain is introduced and evaluated in normal subjects, on a standard clinical scanner at 1.5 Tesla. The method combines a highly stable gradient echo technique with a spiral scan method, to detect brain activation related changes in blood oxygenation with high sensitivity. A motor activation paradigm with a duration of less than 5 min, performed on 10 subjects, consistently showed significant changes in signal intensity in the area of the motor cortex. In all subjects, these changes survived high statistical thresholds.

Quantitative assessment of blood inflow effects in functional MRI signals

Functional MRI (fMRI) signal dependence on changes in blood flow velocities were analyzed for both conventional and echo-planar (EPI) gradient-echo pulse sequences. As the flow velocity increases, the fMRI signal increases monotonically in spoiled gradient-echo sequences, while the fMRI signal may increase or decrease in conventional refocused gradient-echo sequences. A larger flip angle generates a larger inflow contribution to the fMRI signal. For conventional gradient-echo sequences, the inflow contribution to the fMRI images is dominated by the cortical draining veins, while its effect on capillaries is generally small and may be negligible in the spoiled sequences. For EPI gradient-echo sequences, the contribution from inflow effects is relatively small, as compared with the blood oxygen level-dependent (BOLD) contribution, to the fMRI signal, not only for capillaries but also for the cortical draining veins.

Measurement of the arterial concentration of Gd-DTPA using MRI: a step toward quantitative perfusion imaging

A noninvasive method using an inversion recovery turbo-FLASH for dynamic measurement of the arterial input function represented by the bolus passage of Gd-DTPA in the descending aorta is presented, and the results are compared with the input function obtained by arterial blood samples. A good accordance between the two input functions was found, indicating that it is possible to measure the input function to the myocardium using MRI. A variation between the two concentration curves of 5% at upslope, 2.7% at peak point, and < 7% at downslope was found. The study also indicates that a short inversion time < 250 ms has to be used to ensure correct measurement of peak concentration.

Determination of the optimal imaging parameters of the RODEO pulse sequence by computer simulation

A computer program has been developed for evaluating the NMR signal response of various imaging parameters and its efficiency in fat suppression of the RODEO (ROtating Delivery of Excitation Off-resonance) pulse sequence. Both spoiled and refocused RODEO pulse sequences have been considered. By numerically solving the Bloch equation modified for a three-dimensional volume of spins, for realistic RF excitations and Lorentzian distribution of the frequency spectrum for both fat and water, the program permits the imaging contrast and fat suppression of the RODEO pulse sequences to be assessed quantitatively. We have found that excellent fat suppression can be achieved by choosing appropriate imaging parameters. Imaging contrast for different tissues can be enhanced by using longer repetition time (TR) in the spoiled scheme. The complex pattern of NMR signal response and imaging contrast has been observed in the refocused scheme.

Left ventricular radial tagging acquisition using gradient-recalled-echo techniques: sequence optimization

Myocardial tagging with magnetic resonance (MR) imaging offers unique possibilities for noninvasive left ventricular (LV) strain analysis. True three-dimensional strain analysis can be achieved with tags implemented in cardiac short axis and long axis images. Spin-echo (SE) techniques have been used for these studies. However, this approach is time-consuming: images at different phases of the cardiac cycle have to be obtained in successive measurements and hence the total number of measurements equals the number of time frames. Moreover, the images are often degraded by flow and motion artifacts. The purpose of this study was to optimize a faster and more robust MR tagging sequence for use on a clinical whole-body 1 T MR system with optimal persistence of the tags during the entire cardiac cycle. The tagging pulses were implemented in gradient-recalled-echo (GRE) sequences and compared to SE-based acquisitions. The effects of the use of flow-compensating gradients, the excitation angles, and the angles of the saturation pulses have been studied with MR signal simulations and in comparative measurements in volunteers. GRE acquisitions with flow-compensating gradients are robust techniques for myocardial tagging acquisitions. Use of optimized flip angles and saturation pulses can significantly improve delineation of the tag and can be used up to at least 700 ms after the R-wave. Therefore, LV tagging with GRE acquisitions using optimized MR parameters is a robust and promising technique.

Applications of magnetic resonance imaging in food science

The physical and chemical changes that occur in foods during growth, harvest, processing, storage, preparation, and consumption are often very difficult to measure and quantify. Magnetic resonance imaging (MRI) is a pioneering technology, originally developed in the medical field, that is now being used in a large number of disciplines to study a wide variety of materials and processes. In food science, MRI techniques allow the interior of foods to be imaged noninvasively and nondestructively. These images can then be quantified to yield information about several processes and material properties, such as mass and heat transfer, fat and ice crystallization, gelation, water mobidity, composition and volume changes, food stability and maturation, flow behavior, and temperature. This article introduces the fundamental principles of MRI, presents some of the recent advances in MRI technology, and reviews some of the current applications of MRI in food science research.

Decomposition of inflow and blood oxygen level-dependent (BOLD) effects with dual-echo spiral gradient-recalled echo (GRE) fMRI

Image contrast with gradient-recalled echo sequences (GRE) used for fMRI can have both blood oxygen level-dependent (BOLD) and inflow components, and the latter is often undesirable. A dual-echo technique can be used to differentiate these mechanisms, because modulation of signal from inflow is common to both echoes, whereas susceptibility and diffusion-related signal losses are larger in the second echo. An efficient dual-echo interleaved spiral sequence was developed for use with a conventional scanner. It uses a k-space trajectory that spirals out from the origin while the first echo is collected, then spirals back in while collecting the second echo. Decomposition of the data provides separate images of the inflow and T2-weighted components. Results demonstrate the decomposition with phantom experiments and with photic stimulation in normal volunteers.

Fast spin-echo characteristics of visual stimulation-induced signal changes in the human brain

A fast spin-echo (FSE) technique used in a conventional MR imaging scanner has been successfully developed for obtaining functional MR images with high spatial resolution and multiple slices. Our preliminary visual stimulation studies using the FSE technique show that the nuclear MR signal increases by 2.6% during activation in the primary visual cortex. These results provide evidence that the diffusion of tissue water molecules plays a key role in determining functional MR signal amplitude. Because the FSE functional MR imaging signal is extremely sensitive to microvascular (brain capillaries) hemodynamics, the FSE technique can be a powerful tool for studying the neuronal activity of the human brain.

Investigation of coronary vessels in microscopic dimensions by two- and three-dimensional NMR microscopic imaging in the isolated rat heart. Visualization of vasoactive effects of endothelin 1

BACKGROUND

Nuclear magnetic resonance (NMR) imaging of macroscopic coronary vessels is rapidly advancing, whereas little attention has focused on development of NMR techniques for investigation of coronary microvessels. Such techniques would be of particular importance, since conventional methods to visualize coronary microvessels have specific limitations. The aim of our study was to develop two- and three-dimensional (2D and 3D) high-resolution imaging of coronary microvessels. Quantitative analysis of vessel size was performed in tomograms and applied to evaluate the vasoconstrictor effect of endothelin 1.

METHODS AND RESULTS

Angiographic imaging was performed on an 11.75-T magnet by 2D and 3D gradient-echo pulse sequences. In tomograms, the validity of this method in providing correct vessel size was tested by phantom experiments. Experiments were carried out in the isolated constant-pressure-perfused rat heart with continuous registration of coronary flow and left ventricular pressure. NMR pulse sequences were pressure-triggered in mid diastole. Four groups of hearts were studied. In group 1 (n = 20), 2D imaging perpendicular and parallel to the long axis of the heart was performed. Cross sections of vessels with diameter > 140 microns were clearly detectable. In group 2 (control, n = 5) and group 3 (n = 13), tomograms perpendicular to the long axis were obtained before and after administration of vehicle (group 2) and 200 pmol endothelin 1 bolus (group 3). Vehicle had no effect on vessel cross section. Endothelin 1, which decreased global coronary flow by 47%, reduced vessel cross section by 38 +/- 19%. A weak but, on average, significant inverse correlation between area of cross section and vessel size was found. In group 4 (n = 10), 3D imaging was performed in 7 normal hearts and 3 hearts with anterior myocardial infarction. A 3D image of the entire coronary artery tree was obtained, revealing excellent agreement with anatomic studies. In infarcted rat hearts, occlusion of the left coronary artery was demonstrated.

CONCLUSIONS

Visualization and quantification of coronary microvessels are feasible by NMR microscopy. NMR microscopy bears the potential of becoming a powerful tool for the investigation of the coronary microcirculation.

A complex-difference phase-contrast technique for measurement of volume flow rates

Magnetic resonance (MR) phase-difference methods work well for measuring volumetric flow rates when the vessel diameter is large compared with the in-plane voxel dimensions. For small vessels (eg, coronary arteries), partial-volume effects introduce substantial errors in the measured volume flow rate. To correctly measure flow rates through a voxel, both the fraction of the voxel containing moving spins and the phase shift imparted to those spins must be known. The authors propose a flow measurement method that combines information obtained with both the complex-difference and phase-difference processing techniques and thereby provides the fractional volume occupied by the moving spins and the phase of those spins. The complex-difference flow map method proposed results in improved accuracy of MR phase-contrast flow measurements in the presence of partial-volume effects.

Pulsation artifact in short TR MR imaging and angiography: exacerbation with signal averaging

Averaging the signals from more than one excitation per phase-encoding view increases the signal-to-noise ratio and, in conventional spin-echo magnetic resonance imaging, reduces most motion artifacts. To determine the effects of signal averaging on two-dimensional gradient-echo images, acquisitions with different TRs and with no averaging versus multiple-signal averaging were compared in a pulsatile flow phantom and the human abdominal aorta. Intraview (each view repeated before changing the phase-encoding value) and interview (obtaining all views sequentially and then repeating the entire set) averaging methods were used. Pulsation artifacts were present on all images of the flow phantom and the aorta. Intraview signal averaging, the method most commonly used, exacerbated rather than ameliorated pulsation artifacts with short TR sequences. Pulsation artifacts on two-dimensional images obtained with a short TR can be minimized by completing the acquisition as rapidly as possible, avoiding signal averaging. If signal averaging is used for short TR images, it should be interview averaging.

Turbine flow sensor for volume-flow rate verification in MR

A turbine flow sensor for MR flow experiments has been evaluated using reference volume-flow rate measurements obtained using an electromagnetic (EM) flow meter measurements and simultaneous phase contrast (PC) MR acquisitions. After calibration, the device was found to have accuracy (compared with the EM flow meter), linearity, and precision of better than +/- 1%, +/- 3.5%, 3.5%, respectively, in constant flow mode (0 to 30 ml s-1). The frequency response of the flow sensor was flat (within +/- 10%) up to 13.9 Hz. Volume-flow rate measurements on constant and simulated physiologic flow waveforms were in close agreement with both the electromagnetic (EM) flow meter and the gated MR PC estimates.

Pulsatile flow artifacts in fast magnetization-prepared sequences

Fast magnetization-prepared magnetic resonance imaging sequences allow clinical acquisitions in about 1 second, with the preparation phase providing the desired contrast. Pulsatile flow artifacts, although reduced by rapid acquisition, can degrade image quality. The authors explore the causes of aortic pulsatile flow artifacts in inversion-recovery-prepared acquisitions of the abdomen, taking into consideration various parameters. The flow signal within an 8-mm-thick section was simulated and subsequently Fourier transformed to determine the location and extent of flow artifacts. Results of simulations were validated with abdominal images of human subjects. Recording all encodings within one cardiac cycle reduced pulsatile flow artifacts in nonsegmented acquisitions with sequential phase-encoding order, regardless of the location of magnetization preparation within the cardiac cycle. In segmented acquisitions, however, the sequential order always increased flow artifacts. To reduce the artifacts in short TI acquisitions, the magnetization should be prepared during diastole. In clinical acquisitions, flow artifacts were further reduced by modifying the phase-encoding scheme.

Potential pitfalls of functional MRI using conventional gradient-recalled echo techniques

The conventional gradient-recalled echo technique, FLASH, has widely been used for functional MRI. FLASH results at 4 T with short TEs of 10-20 ms mimic those at 1.5 T with TEs of 25-50 ms or longer. Under these conditions, large venous vessels dominate the activated area; however, the use of longer TEs at 4 T reveals activation in gray matter areas as well as large vessels. Inflow effects of large vessels can be greatly reduced with centric-reordering of phase-encoding steps and inter-image delay. Finger and toe movement paradigms show that functional activation maps are consistent with classical somatotopic maps, and are specific to the tasks. Navigator-based motion correction generates functional maps with larger activation areas by reducing physiological noise.