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.     

Functional MRI using steady-state arterial water labeling

Flow-sensitive functional MRI (fMRI) was performed using steady-state arterial water labeling (SS-AWL). Arterial water labeling was accomplished by flow induced adiabatic fast passage. The signal intensity of the visual cortex in arterial water labeled images decreased by ?1.4% during visual stimulation of the brain. Acquisition of arterial water unlabeled and labeled images allows measurement of relative cerebral blood flow increase during brain activation. During visual stimulation, cerebral blood flow in the visual cortex increased by 17 to 35% as measured by SS-AWL. Quantitation of brain activation in terms of a physiological parameter using SS-AWL will facilitate comparative fMRI studies under different conditions.     

Quantitative measurement of high flow velocities by a spin echo MR technique.

A new method of flow measurement using a spin echo (SE) technique has been developed on the basis of the flow effect that at high velocities signal intensity decreases linearly with increasing flow velocity. Flow velocity is calculated from the signal intensity ratio of the flowing material in two images with the same imaging parameters but different echo times. The linear relationship between the signal intensity and flow velocity was examined with a steady flow phantom. When assessed with steady flows in the phantom, flow velocities calculated by this method were in good agreement with velocities measured by a flow meter. This method was used with ECG gating to measure the blood flow of the right common carotid artery of a healthy volunteer. The measured peak flow velocity and the pattern of flow velocities during systole correlated well with the results obtained by Doppler ultrasound.     

BOLD and CBV-weighted functional magnetic resonance imaging of the rat somatosensory system.

A multislice spin echo EPI sequence was used to obtain functional MR images of the entire rat brain with blood oxygenation level dependent (BOLD) and cerebral blood volume (CBV) contrast at 11.7 T. Maps of activation incidence were created by warping each image to the Paxinos rat brain atlas and marking the extent of the activated area. Incidence maps for BOLD and CBV were similar, but activation in draining veins was more prominent in the BOLD images than in the CBV images. Cerebellar activation was observed along the surface in BOLD images, but in deeper regions in the CBV images. Both effects may be explained by increased signal dropout and distortion in the EPI images after administration of the ferumoxtran-10 contrast agent for CBV fMRI. CBV-weighted incidence maps were also created for 10, 20, and 30 mg Fe/kg doses of ferumoxtran-10. The magnitude of the average percentage change during stimulation increased from 4.9% with the 10 mg Fe/kg dose to 8.7% with the 30-mg Fe/kg dose. Incidence of activation followed a similar trend.     

Balanced phase-contrast steady-state free precession (PC-SSFP): a novel technique for velocity encoding by gradient inversion.

A technique for measuring velocity is presented that combines cine phase contrast (PC) MRI and balanced steady-state free precession (SSFP) imaging, and is thus termed PC-SSFP. Flow encoding was performed without the introduction of additional velocity encoding gradients in order to keep the repetition time (TR) as short as in typical SSFP imaging sequences. Sensitivity to through-plane velocities was instead established by inverting (i.e., negating) all gradients along the slice-select direction. Velocity sensitivity (VENC) could be adjusted by altering the first moments of the slice-select gradients. Disturbances of the SSFP steady state were avoided by acquiring different flow echoes in consecutively (i.e., sequentially) executed scans, each over several cardiac cycles, using separate steady-state preparation periods. A comparison of phantom measurements with those from established 2D-cine-PC MRI demonstrated excellent correlation between both modalities. In examinations of volunteers, PC-SSFP exhibited a higher intrinsic signal-to-noise ratio (SNR) and consequently low phase noise in measured velocities compared to conventional PC scans. An additional benefit of PC-SSFP is that it relies less on in-flow-dependent signal enhancement, and thus yields more uniform SNRs and better depictions of vessel geometry throughout the whole cardiac cycle in structures with slow and/or pulsatile flow.     

Extravascular proton-density changes as a non-BOLD component of contrast in fMRI of the human spinal cord.

The fractional signal intensity change (Delta S/S) observed during activation in T(2)-weighted fMRI of the spinal cord has previously been shown to depend linearly on the echo time (TE) but to have a positive value of roughly 2.5% extrapolated to zero TE. In this study we investigated the origin of this finding by measuring the Delta S/S in spinal fMRI with very short TEs. Our results demonstrate that the Delta S/S does not approach zero, but has a value as high as 3.3% at TE = 11 ms. At TEs > 33 ms we observed the linear relationship between Delta S/S and TE as in previous studies. These data demonstrate that there is a non-BOLD contribution to signal changes observed in spinal fMRI. We hypothesize that this contribution is a local proton density increase due to increased water exudation from capillaries with increased blood flow during neuronal activation, and term this effect "signal enhancement by extravascular protons" (SEEP).     

Blood flow: magnetic resonance imaging

The appearance of flowing fluid has been evaluated in several clinical situations using a flow phantom, computer simulation, and clinical magnetic resonance (MR) images. Unsaturated protons just entering the imaging volume can emit a strong signal relative to the partially saturated adjacent tissue ("flow-related enhancement"). Slow laminar flow in veins can be distinguished on the basis of a stronger second echo due to rephasing effects ("even echo rephasing"). Synchronization of the cardiac cycle and the MR pulsing sequence produces increased signal in sections acquired during diastole ("diastolic pseudogating"). Intraluminal signal is shown to decrease as velocity is increased ("high velocity signal loss"). Onset of turbulence causes further loss of signal. Direction of flow oblique to the imaging plane can be predicted on the basis of decreased upstream and increased downstream signal.     

High-resolution MR cisternography of the cerebellopontine angle: 2D versus 3D fast spin-echo sequences.

BACKGROUND AND PURPOSE: The clinical usefulness of MR cisternography of the cerebellopontine angle, applying 2D or 3D fast spin-echo sequences, has been reported recently. Our purpose was to investigate the cause of signal loss in CSF in the prepontine or cerebellopontine angle cistern on 2D FSE MR images and to compare the cisternographic effects of 2D and 3D FSE sequences. METHODS: Preliminary experiments were performed in four volunteers to assess the causes of signal loss. Initially, using a 2D cardiac-gated cine phase-contrast method with a velocity encoding value of 6 cm/s, we measured the velocity and flow pattern of CSF. Comparisons were made to assess the effects of intravoxel dephasing, amplitude of the section-selecting gradient, echo time (TE), and section thickness. Four healthy subjects and 13 patients with ear symptoms were examined, and multisection 3-mm-thick 2D images and 30-mm-slab, 1-mm-section 3D images were compared qualitatively and quantitatively. Then, 3D MR cisternography was performed in 400 patients with ear symptoms, and qualitative evaluation was performed. RESULTS: In volunteers, the average peak velocity of CSF was 1.2 cm/s. With TE = 250, CSF may move an average of 3 mm, and can be washed out of a 3-mm-thick 2D section volume. The CSF signal relative to that of a water phantom decreased gradually as TE increased on single-section 3-mm-thick 2D images. The CSF signal relative to that of the water phantom increased gradually as section thickness increased. No significant differences were noted in intravoxel dephasing and amplitude of the section-selecting gradient. The contrast-to-noise ratio (CNR) between CSF and the cerebellar peduncle, and the visibility of the cranial nerves and vertebrobasilar artery were significantly improved on 3D images in 17 subjects. In images from 400 patients, no significant signal loss in the cistern was observed using 3D FSE. CONCLUSION: CSF signal loss in thin-section 2D MR cisternography is mainly attributable to the wash-out phenomenon. 3D acquisition can reduce this phenomenon and provide thinner sections. The scan time for 3D acquisition is not excessive when a long echo train length and half-Fourier imaging are used. MR cisternography should be performed using a 3D acquisition.     

Improving the temporal resolution of functional MR imaging using keyhole techniques

Using a keyhole technique, it is shown that the data acquisition rate of gradient-echo imaging for functional MRI (fMRI) studies can be increased substantially. The resulting enhancement of the temporal resolution of fMRIs was accomplished without modifying the hardware of a conventional MRI system. High spatial resolution fMRI images were first collected with conventional full k-space acquisition and image reconstruction. Using the same data set, simulation reconstruction using the keyhole principle and zero-padding were performed for comparison with the full k-space reconstruction. No significant changes were found for fMRI images generated from the keyhole technique with a data sharing profile of 50% of the k-space. As k-space data sharing profiles increased to 75 and 87.5%, the keyhole fMRI images began to show only modest changes in activation intensity and area compared with the standard images. In contrast, zero-padding fMRI images produced a significant disparity both in activation intensity and area relative to the truly high-resolution fMRI images. The keyhole technique's ability to retain the intensity and area of fMRI information, while substantially reducing acquisition time, makes it a promising method for fMRI studies.     

Non‐contrast‐enhanced vascular magnetic resonance imaging using flow‐dependent preparation with subtraction

Recent concerns over contrast agent safety have encouraged new developments in non‐contrast‐enhanced vascular imaging techniques. This work investigates the potential for imaging both arteries and veins with vascular anatomy by nonenhanced static subtraction angiography (VANESSA), a method using controllable flow suppression together with subtraction of bright‐ and dark‐blood images. The lower legs of eight healthy volunteers and three patients were imaged using a modified motion‐sensitized driven equilibrium preparation, with three‐dimensional balanced steady‐state free precession readout. The vascular signal decreased with increasing motion‐suppression gradient amplitude, and was suppressed when the velocity‐encoding parameter was (approximately) less than the measured flow velocity. Selected pairs of images were subtracted to depict vessels with either fast flow (e.g. arteries), slow flow (e.g. veins), or both. Several methodological modifications improved image quality and reduced the background signal from static tissues. Subjectively assessed image quality in volunteers was rated as excellent for 56/64 arterial segments, and good or excellent for 35/64 veins. In conclusion, VANESSA enables rapid non‐contrast‐enhanced imaging of arteries and veins, combining information on both morphology and flow. This study demonstrates good technical performance in volunteers and evaluation in patients with vascular disease is warranted. Magn Reson Med, 2012. © 2011 Wiley Periodicals, Inc.     

Rapid mapping of flow velocity using a new PARSE method.

A new method for flow velocity mapping is presented here. Instead of the conventional approach of employing two images (velocity sensitive and control) to generate velocity information, in the new method one determines the velocity directly from a single-shot acquisition by solving an inverse problem. This technique is a variant of single-shot parameter assessment by retrieval from signal encoding (SS-PARSE). The results of simulation and phantom studies show strong agreement with the actual velocities. The prototype method can measure velocities in the range of -50 to 50 cm/s, which is roughly appropriate for future applications in dynamic blood flow measurement in carotid arteries.     

中文语言活动区功能磁共振研究

目的　利用BOLD功能磁共振 ,初步探索国人说汉语时大脑语言相关功能区。方法　对 2 2名正常中国人 ,说中文句子状态下进行BOLD磁共振脑功能成像 ,探索其大脑语言活动皮层相关功能区。结果　说中文时激活的脑区包括 :两侧运动区 ,左右侧额下回 ,左右侧颞上回 ,左侧岛叶及左右侧小脑半球。结论　BOLD功能磁共振可显示中文大脑活动相关皮层 ,是无创性研究人类语言的有力武器     

Quantification of blood velocity and flow rates in the uterine vessels using echo planar imaging at 0.5 Tesla

Purpose:

To measure uterine artery and vein blood velocity and flow rate profiles using MRI during normal pregnancy.

Materials and Methods:

A two‐shot velocity magnitude‐encoded echo planar imaging (EPI) sequence is used at a magnetic field 0.5T. Data analysis procedures, necessary to overcome problems associated with low signal to noise ratio (SNR), and a spatial resolution comparable to the vessel size were used.

Results:

The measured blood flow values averaged over nine volunteers for the mean velocity are 5.33 and 3.97 cm/s and for the unilateral flow rate are 203 and 274 mL/min (for the arteries and veins respectively). Values for the flow rate are consistent with ultrasound Doppler studies. Arterial velocity measurements are more pulsatile than venous ones and validation calculations performed on average velocity values would suggest that the nature of blood flow in the uterine vessels is laminar.

Conclusion:

This study presents the first report of noninvasive quantitative measurements of uterine artery and vein blood velocity and flow rate profiles using MRI during normal pregnancy. Consistent and reproducible measurements have been obtained by subject specific sequence optimization and data analysis procedures. J. Magn. Reson. Imaging 2010;31:921–927. ©2010 Wiley‐Liss, Inc.     

Lateralization of functional magnetic resonance imaging (fMRI) activation in the auditory pathway of patients with lateralized tinnitus

Introduction Tinnitus is hypothesized to be an auditory phantom phenomenon resulting from spontaneous neuronal activity somewhere along the auditory pathway. We performed fMRI of the entire auditory pathway, including the inferior colliculus (IC), the medial geniculate body (MGB) and the auditory cortex (AC), in 42 patients with tinnitus and 10 healthy volunteers to assess lateralization of fMRI activation. Methods Subjects were scanned on a 3T MRI scanner. A T2*-weighted EPI silent gap sequence was used during the stimulation paradigm, which consisted of a blocked design of 12 epochs in which music presented binaurally through headphones, which was switched on and off for periods of 50 s. Using SPM2 software, single subject and group statistical parametric maps were calculated. Lateralization of activation was assessed qualitatively and quantitatively. Results Tinnitus was lateralized in 35 patients (83%, 13 right-sided and 22 left-sided). Significant signal change (P _corrected < 0.05) was found bilaterally in the primary and secondary AC, the IC and the MGB. Signal change was symmetrical in patients with bilateral tinnitus. In patients with lateralized tinnitus, fMRI activation was lateralized towards the side of perceived tinnitus in the primary AC and IC in patients with right-sided tinnitus, and in the MGB in patients with left-sided tinnitus. In healthy volunteers, activation in the primary AC was left-lateralized. Conclusion Our paradigm adequately visualized the auditory pathways in tinnitus patients. In lateralized tinnitus fMRI activation was also lateralized, supporting the hypothesis that tinnitus is an auditory phantom phenomenon.     

Turbulent Fluctuation Velocity: The Most Significant Determinant of Signal Loss in Stenotic Vessels

Studies of flow in a 90%-stenosis phantom were conducted to elucidate the parameters and mechanisms responsible for signal loss in MR angiographic images. The studies independently evaluated the effect of velocity, Reynolds number, turbulent fluctuation velocity, and turbulence intensity on the amount of post-stenotic signal loss. Results suggested that the magnitude of the turbulent fluctuation velocity, not merely the presence of turbulence or the intensity of turbulence, was the parameter that determined the extent of the signal loss. The study suggests that future flow phantom studies should be conducted with fluids having physiologic velocities and viscosities to obtain accurate levels of turbulent fluctuation velocities and hence reproduce results of in-vivo signal-loss patterns. The mechanism for signal loss is that the temporal and spatial variations of the turbulent fluctuation velocity cause a range of phases to be present within a voxel. Examination of the theoretical aspects of fluid turbulence suggest that shortening gradient durations and imaging during diastole may help reduce signal loss.     

Pseudolesion of the bile duct caused by flow effect: a diagnostic pitfall of MR cholangiopancreatography

OBJECTIVE: Our objective was to examine the influence of the shape of imaged structures and the velocity of flow on the appearance of flow artifacts seen on MR cholangiopancreatography (MRCP) in a phantom model. MATERIALS AND METHODS: Three types of phantoms representing the biliary system were constructed. The first phantom type was a straight tube; the second, a single tube in which the inlet and outlet diameters varied by a ratio of as much as 1:6; and the third, a tube that simulated a stricture in the biliary system and a gallstone. All experiments were repeated three times. RESULTS: We did not observe any flow artifacts in the experiments we performed with the straight tubes. A higher rate of flow resulted in decreased signal intensity in tubes simulating bile ducts; the decreased signal was most likely to be observed on images in which the speed of flow exceeded 5 mm/sec. Flow artifacts were seen only if the ratio between the inlet and outlet diameters was 1:4 or greater. Simulations of bile duct abnormalities--such as a 50% stricture or the presence of a gallstone--did not produce any flow artifacts. CONCLUSION: In our experiments, a flow artifact could be seen on images in which the ratio between the inlet and the outlet diameters in the phantom was equal to or greater than 1:4. This finding indicates that a flow artifact could be observed in dilated bile ducts on MRCP under clinical conditions. Knowing that a pseudo-filling defect can be caused by a flow artifact should help to prevent misinterpretation of MRCP images.     

Automatic model-based contour detection and blood flow quantification in small vessels with velocity encoded magnetic resonance imaging

RATIONALE AND OBJECTIVES: The quantitative assessment of blood flow in peripheral vessels from phase-contrast magnetic resonance imaging studies requires the accurate delineation of vessel contours in cross-sectional magnetic resonance images. The conventional manual segmentation approach is tedious, time-consuming, and leads to significant inter- and intraobserver variabilities. The aim of this study was to verify whether automatic model-based segmentation decreases these problems by fitting a model to the actual blood velocity profile. METHODS: In this study 2 new fully automatic methods (a static and a dynamic approach) were developed and compared with manual analyzes using phantom and in vivo studies of internal carotid and vertebral arteries in healthy volunteers. The automatic segmentation approaches were based on fitting a 3D parabolic velocity model to the actual velocity profiles. In the static method, the velocity profiles were averaged over the complete cardiac cycle, whereas the dynamic method takes into account the velocity data of each cardiac time bin individually. Materials consisted of the magnetic resonance imaging data from 3 straight phantom tubes and the blood velocity profiles of 8 volunteers. RESULTS: For the phantom studies, the automatic dynamic approach performed significantly better than the manual analysis (intraclass correlations [ICC] of 0.62-0.98 and 0.30-0.86, respectively). For the assessment of the total cerebral blood flow in the in vivo studies, the automatic static method performed significantly better than the manual 1 (ICC of 0.98-0.98 and 0.93-0.95, respectively). However, the automatic dynamic method was not significantly better than the manual 1 (ICC = 0.92-0.96) but had the advantage of providing additional parameters. CONCLUSION: Blood flow in magnetic resonance images of small vessels can be assessed accurately, rapidly, and fully automatically using model-based postprocessing by fitting a first approximation of the velocity profile to the actual flow data.     

Shared velocity encoding: a method to improve the temporal resolution of phase-contrast velocity measurements

Phase-contrast magnetic resonance imaging (PC-MRI) is used routinely to measure fluid and tissue velocity with a variety of clinical applications. Phase-contrast magnetic resonance imaging methods require acquisition of additional data to enable phase difference reconstruction, making real-time imaging problematic. Shared Velocity Encoding (SVE), a method devised to improve the effective temporal resolution of phase-contrast magnetic resonance imaging, was implemented in a real-time pulse sequence with segmented echo planar readout. The effect of SVE on peak velocity measurement was investigated in computer simulation, and peak velocities and total flow were measured in a flow phantom and in volunteers and compared with a conventional ECG-triggered, segmented k-space phase-contrast sequence as a reference standard. Computer simulation showed a 36% reduction in peak velocity error from 8.8 to 5.6% with SVE. A similar reduction of 40% in peak velocity error was shown in a pulsatile flow phantom. In the phantom and volunteers, volume flow did not differ significantly when measured with or without SVE. Peak velocity measurements made in the volunteers using SVE showed a higher concordance correlation (0.96) with the reference standard than non-SVE (0.87). The improvement in effective temporal resolution with SVE reconstruction has a positive impact on the precision and accuracy of real-time phase-contrast magnetic resonance imaging peak velocity measurements.     

The intravascular contribution to fmri signal change: monte carlo modeling and diffusion-weighted studies in vivo

Understanding the relationship between fMRI signal changes and activated cortex is paramount to successful mapping of neuronal activity. To this end, the relative extravascular and intravascular contribution to fMRI signal change from capillaries (localized), venules (less localized) and macrovessels (remote, draining veins) must be determined. In this work, the authors assessed both the extravascular and intravascular contribution to blood oxygenation level-dependent gradient echo signal change at 1.5 T by using a Monte Carlo model for susceptibility-based contrast in conjunction with a physiological model for neuronal activation-induced changes in oxygenation and vascular volume fraction. The authors compared our Model results with experimental fMRI signal changes with and without velocity sensitization via bipolar gradients to null the intravascular signal. The model and experimental results are in agreement and suggest that the intravascular spins account for the majority of fMRI signal change on T₂^*-weighted images at 1.5 T.     

The value of detection of flow voids between the uterus and the leiomyoma with MRI

Flow voids are occasionally seen between the uterus and the leiomyoma. This study was performed to determine the prevalence of this sign and its relation to tumor size and location, internal signal intensity of the leiomyomas, and the ease of detectability using different MR sequences. Also, to understand the pathologic causes for the phenomenon of flow voids. The MR images of the pelvis in 92 females with 359 uterine leiomyomas were analyzed. Flow voids located between the uterus and the leiomyomas were seen in 32 lesions on T1-weighted images and in 12 on T2-weighted images. Flow void was seen exclusively in leiomyomas of 3-cm diameter of larger. In five pedunculated subserous leiomyomas, flow voids were seen within the pedicle of the tumor. Large lieomyomas and leiomyomas with irregular signal intensity on T2-weighted images showed significantly greater frequency of signal void than other leiomyomas (P < .01). The sign was not seen in any other pelvic tumors. On pathologic evaluation, the flow voids were found to correspond to dilated feeding arteries located outside the capsule of the leiomyoma.