Comparison between magnetic resonance phase contrast imaging and transcranial Doppler ultrasound with regard to blood flow velocity in intracranial arteries: work in progress

Effect of MRI acquisition parameters on accuracy and precision of phase-contrast measurements in a small-lumen vessel phantom

BACKGROUND

Phase-contrast magnetic resonance imaging (PC-MRI) quantifies blood flow and velocity noninvasively. Challenges arise in neurovascular disorders due to small vessels. We evaluated the impact of voxel size, number of signal averages (NSA), and velocity encoding (VENC) on PC-MRI measurement accuracy and precision in a small-lumen vessel phantom.

METHODS

We constructed an in vitro model with a constant flow rate using a 2.2-mm inner diameter plastic tube. A reservoir with a weight scale and timer was used as standard reference. Gradient-echo T1 weighted PC-MRI sequence was performed on a 3-T scanner with varying voxel size (2.5, 5.0, 7.5 mm3), NSA (1, 2, 3), and VENC (200, 300, 400 cm/s). We repeated measurements nine times per setting, calculating mean flow rate, maximum velocity, and least detectable difference (LDD).

RESULTS

PC-MRI flow measurements were higher than standard reference values (mean ranging from 7.3 to 9.5 mL/s compared with 6.6 mL/s). Decreased voxel size improved accuracy, reducing flow rate measurements from 9.5 to 7.3 mL/s. The LDD for flow rate and velocity varied between 1 and 5%. The LDD for flow rate decreased with increased voxel size and NSA (p = 0.033 and 0.042). The LDD for velocity decreased with increased voxel size (p < 10-16). No change was observed when VENC varied.

CONCLUSIONS

PC-MRI overestimated flow. However, it has high precision in a small-vessel phantom with constant flow rate. Improved accuracy was obtained with increasing spatial resolution (smaller voxels). Improved precision was obtained with increasing signal-to-noise ratio (larger voxels and/or higher NSA).

RELEVANCE STATEMENT

Phase-contrast MRI is clinically used in large vessels. To further investigate the possibility of using phase-contrast MRI for smaller intracranial vessels in neurovascular disorders, we need to understand how acquisition parameters affect phase-contrast MRI-measured flow rate and velocity in small vessels.

KEY POINTS

• PC-MRI measures flow and velocity in a small lumen phantom with high precision but overestimates flow rate. • The precision of PC-MRI measurements matches the precision of standard reference for flow rate measurements. • Optimizing PC-MRI settings can enhance accuracy and precision in flow rate and velocity measurements.

Comparison of computational fluid dynamics with transcranial Doppler ultrasound in response to physiological stimuli

Cerebrovascular haemodynamics are sensitive to multiple physiological stimuli that require synergistic response to maintain adequate perfusion. Understanding haemodynamic changes within cerebral arteries is important to inform how the brain regulates perfusion; however, methods for direct measurement of cerebral haemodynamics in these environments are challenging. The aim of this study was to assess velocity waveform metrics obtained using transcranial Doppler (TCD) with flow-conserving subject-specific three-dimensional (3D) simulations using computational fluid dynamics (CFD). Twelve healthy participants underwent head and neck imaging with 3 T magnetic resonance angiography. Velocity waveforms in the middle cerebral artery were measured with TCD ultrasound, while diameter and velocity were measured using duplex ultrasound in the internal carotid and vertebral arteries to calculate incoming cerebral flow at rest, during hypercapnia and exercise. CFD simulations were developed for each condition, with velocity waveform metrics extracted in the same insonation region as TCD. Exposure to stimuli induced significant changes in cardiorespiratory measures across all participants. Measured absolute TCD velocities were significantly higher than those calculated from CFD (P range < 0.001-0.004), and these data were not correlated across conditions (r range 0.030-0.377, P range 0.227-0.925). However, relative changes in systolic and time-averaged velocity from resting levels exhibited significant positive correlations when the distinct techniques were compared (r range 0.577-0.770, P range 0.003-0.049). Our data indicate that while absolute measures of cerebral velocity differ between TCD and 3D CFD simulation, physiological changes from resting levels in systolic and time-averaged velocity are significantly correlated between techniques.

Quantitative time-of-flight MR angiography for simultaneous luminal and hemodynamic evaluation of the intracranial arteries

PURPOSE

To report a quantitative time-of-flight (qTOF) MRA technique for simultaneous luminal and hemodynamic evaluation of the intracranial arteries.

METHODS

Implemented using a thin overlapping slab 3D stack-of-stars based 3-echo FLASH readout, qTOF was tested in a flow phantom and for imaging the intracranial arteries of 10 human subjects at 3 Tesla. Display of the intracranial arteries with qTOF was compared to resolution-matched and scan time-matched standard Cartesian 3D time-of-flight (TOF) MRA, whereas quantification of mean blood flow velocity with qTOF, done using a computer vision-based inter-echo image analysis procedure, was compared to 3D phase contrast MRA. Arterial-to-background contrast-to-noise ratio was measured, and intraclass correlation coefficient was used to evaluate agreement of flow velocities.

RESULTS

For resolution-matched protocols of similar scan time, qTOF portrayed the intracranial arteries with good morphological correlation with standard Cartesian TOF, and both techniques provided superior contrast-to-noise ratio and arterial delineation compared to phase contrast (20.6 ± 3.0 and 37.8 ± 8.7 vs. 11.5 ± 2.2, P < .001, both comparisons). With respect to phase contrast, qTOF showed excellent agreement for measuring mean flow velocity in the flow phantom (intraclass correlation coefficient = 0.981, P < .001) and good agreement in the intracranial arteries (intraclass correlation coefficient = 0.700, P < .001). Stack-of-stars data sampling used with qTOF eliminated oblique in-plane flow misregistration artifacts that were seen with standard Cartesian TOF.

CONCLUSION

qTOF is a new 3D MRA technique for simultaneous luminal and hemodynamic evaluation of the intracranial arteries that provides significantly greater contrast-to-noise ratio efficiency than phase contrast and eliminates misregistration artifacts from oblique in-plane blood flow that occur with standard 3D TOF.

A cross-sectional, case-control study of intracranial arterial wall thickness and complete blood count measures in sickle cell disease

In sickle cell disease (SCD), cerebral oxygen delivery is dependent on the cerebral vasculature's ability to increase blood flow and volume through relaxation of the smooth muscle that lines intracranial arteries. We hypothesised that anaemia extent and/or circulating markers of inflammation lead to concentric macrovascular arterial wall thickening, visible on intracranial vessel wall magnetic resonance imaging (VW-MRI). Adult and pediatric SCD (n = 69; age = 19.9 ± 8.6 years) participants and age- and sex-matched control participants (n = 38; age = 22.2 ± 8.9 years) underwent 3-Tesla VW-MRI; two raters measured basilar and bilateral supraclinoid internal carotid artery (ICA) wall thickness independently. Mean wall thickness was compared with demographic, cerebrovascular and haematological variables. Mean vessel wall thickness was elevated (P < 0·001) in SCD (1·07 ± 0·19 mm) compared to controls (0·97 ± 0·07 mm) after controlling for age and sex. Vessel wall thickness was higher in participants on chronic transfusions (P = 0·013). No significant relationship between vessel wall thickness and flow velocity, haematocrit, white blood cell count or platelet count was observed; however, trends (P < 0·10) for wall thickness increasing with decreasing haematocrit and increasing white blood cell count were noted. Findings are discussed in the context of how anaemia and circulating inflammatory markers may impact arterial wall morphology.

Cerebral Blood Flow Response to Simulated Hypovolemia in Essential Hypertension: A Magnetic Resonance Imaging Study

Supplemental Digital Content is available in the text.

Hypertension is associated with raised cerebral vascular resistance and cerebrovascular remodeling. It is currently unclear whether the cerebral circulation can maintain cerebral blood flow (CBF) during reductions in cardiac output (CO) in hypertensive patients thereby avoiding hypoperfusion of the brain. We hypothesized that hypertension would impair the ability to effectively regulate CBF during simulated hypovolemia. In the present study, 39 participants (13 normotensive, 13 controlled, and 13 uncontrolled hypertensives; mean age±SD, 55±10 years) underwent lower body negative pressure (LBNP) at −20, −40, and −50 mmHg to decrease central blood volume. Phase-contrast MR angiography was used to measure flow in the basilar and internal carotid arteries, as well as the ascending aorta. CBF and CO decreased during LBNP (P<0.0001). Heart rate increased during LBNP, reaching significance at −50 mmHg (P<0.0001). There was no change in mean arterial pressure during LBNP (P=0.3). All participants showed similar reductions in CBF (P=0.3, between groups) and CO (P=0.7, between groups) during LBNP. There was no difference in resting CBF between the groups (P=0.36). In summary, during reductions in CO induced by hypovolemic stress, mean arterial pressure is maintained but CBF declines indicating that CBF is dependent on CO in middle-aged normotensive and hypertensive volunteers. Hypertension is not associated with impairments in the CBF response to reduced CO.

Analysis of the time-velocity curve in phase-contrast magnetic resonance imaging: a phantom study

The aim of this study was to analyze the characteristics of time-velocity curve acquired by phase-contrast magnetic resonance imaging (PC-MRI) using an in-vitro flow model as a reference for hemodynamic studies. The time- velocity curves of the PC-MRI were compared with Doppler ultrasonography (US) and also compared with those obtained in the electromagnetic flowmeter. The correlation between techniques was analyzed using an electromagnetic flowmeter as a reference standard; the maximum, minimum, and average velocities, full-width at half-maximum (FWHM), and ascending gradient (AG) were measured from time-velocity curves. The correlations between an electromagnetic flowmeter and the respective measurement technique for the PC-MRI and Doppler US were found to be high (mean R² > 0.9, p < 0.05). These results indicate that these measurement techniques are useful for measuring blood flow information and reflect actual flow. The PC-MRI was the best fit for the minimum velocity and FWHM, and the maximum velocity and AG were the best fit for Doppler US. The PC-MRI showed lower maximum velocity value and higher minimum velocity value than Doppler US. Therefore, PC-MRI demonstrates more obtuse time-velocity curve than Doppler US. In addition, the time- velocity curve of PC-MRI could be calibrated by introducing formulae that can convert each measurement value to a reference standard value within a 10% error. The PC-MRI can be used to estimate the Doppler US using this formula.

Normal ranges and test-retest reproducibility of flow and velocity parameters in intracranial arteries measured with phase-contrast magnetic resonance imaging

INTRODUCTION

The purpose of the present study was to investigate normal ranges and test-retest reproducibility of phase-contrast MRI (PC-MRI)-measured flow and velocity parameters in intracranial arteries.

METHODS

Highest flow (HF), lowest flow (LF), peak systolic velocity (PSV), and end diastolic velocity (EDV) were measured at two dates in the anterior (ACA), middle (MCA), and posterior (PCA) cerebral arteries of 30 healthy volunteers using two-dimensional PC-MRI at 3 T. Least detectable difference (LDD) was calculated.

RESULTS

In the left ACA, HF was (mean (range, LDD)) 126 ml/min (36-312, 59 %), LF 61 ml/min (0-156, 101 %), PSV 64 cm/s (32-141, 67 %), and EDV 35 cm/s (18-55, 42 %); in the right ACA, HF was 154 ml/min (42-246, 49 %), LF 77 ml/min (0-156, 131 %), PSV 75 cm/s (26-161, 82 %), and EDV 39 cm/s (7-59, 67 %). In the left MCA, HF was 235 ml/min (126-372, 35 %), LF 116 ml/min (42-186, 48 %), PSV 90 cm/s (55-183, 39 %), and EDV 46 cm/s (20-66, 28 %); in the right MCA, HF was 238 ml/min (162-342, 44 %), LF 120 ml/min (72-216, 48 %), PSV 88 cm/s (55-141, 35 %), and EDV 45 cm/s (26-67, 23 %). In the left PCA, HF was 108 ml/min (42-168, 54 %), LF 53 ml/min (18-108, 64 %), PSV 50 cm/s (24-77, 63 %), and EDV 28 cm/s (14-40, 45 %); in the right PCA, HF was 98 ml/min (30-162, 49 %), LF 49 ml/min (12-84, 55 %), PSV 47 cm/s (27-88, 59 %), and EDV 27 cm/s (16-41, 45 %).

CONCLUSION

PC-MRI-measured flow and velocity parameters in the main intracranial arteries have large normal ranges. Reproducibility is highest in MCA.

Theory and validation of magnetic resonance fluid motion estimation using intensity flow data

Background

Motion tracking based on spatial-temporal radio-frequency signals from the pixel representation of magnetic resonance (MR) imaging of a non-stationary fluid is able to provide two dimensional vector field maps. This supports the underlying fundamentals of magnetic resonance fluid motion estimation and generates a new methodology for flow measurement that is based on registration of nuclear signals from moving hydrogen nuclei in fluid. However, there is a need to validate the computational aspect of the approach by using velocity flow field data that we will assume as the true reference information or ground truth.

Methodology/Principal Findings

In this study, we create flow vectors based on an ideal analytical vortex, and generate artificial signal-motion image data to verify our computational approach. The analytical and computed flow fields are compared to provide an error estimate of our methodology. The comparison shows that the fluid motion estimation approach using simulated MR data is accurate and robust enough for flow field mapping. To verify our methodology, we have tested the computational configuration on magnetic resonance images of cardiac blood and proved that the theory of magnetic resonance fluid motion estimation can be applicable practically.

Conclusions/Significance

The results of this work will allow us to progress further in the investigation of fluid motion prediction based on imaging modalities that do not require velocity encoding. This article describes a novel theory of motion estimation based on magnetic resonating blood, which may be directly applied to cardiac flow imaging.

Cerebral hemodynamics and white matter hyperintensities in CADASIL

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a hereditary small-vessel disease caused by mutations in the NOTCH3 gene on chromosome 19. On magnetic resonance imaging (MRI), subcortical white matter hyperintensities and lacunar infarcts are visualized. It is unknown whether a decrease in cerebral blood flow or cerebrovascular reactivity is primarily responsible for the development of white matter hyperintensities and lacunar infarcts. The authors used phase-contrast MRI in 40 NOTCH3 mutation carriers (mean age 45 +/- 10 years) and 22 nonmutated family members (mean age 39 +/- 12 years), to assess baseline total cerebral blood flow (TCBF) and cerebrovascular reactivity after acetazolamide. Mean baseline TCBF was significantly decreased in NOTCH3 mutation carriers. In young subjects, baseline TCBF was significantly lower than in nonmutation carriers (mean difference 124 mL/min). Furthermore, baseline TCBF did not differ significantly between mutation carriers with minimal and mutation carriers with moderate or severe white matter hyperintensities. No significant difference in mean cerebrovascular reactivity was found between mutation carriers and nonmutation carriers. This study suggests that a decrease in baseline TCBF in NOTCH3 mutation carriers precedes the development of white matter hyperintensities.

Quantification of blood flow in the carotid arteries: comparison of Doppler ultrasound and three different phase-contrast magnetic resonance imaging sequences

RATIONALE AND OBJECTIVES

To compare blood flow velocities in the carotid arteries measured with three different magnetic resonance (MR) phase-contrast imaging techniques and with percutaneous Doppler ultrasound.

METHODS

Fourteen healthy male volunteers with a mean age of 33 +/- 3.8 years were studied. Ultrasound and MR phase velocity mapping of both common carotid arteries (n = 28) was performed within 5 hours. A two-dimensional fast low-angle shot sequence with retrospective cardiac gating, a sequence with prospective cardiac triggering, and a breath-hold sequence with prospective cardiac triggering were used. Resistance indexes and pulsatility indexes were calculated for all modalities.

RESULTS

The comparison of flow velocities obtained with ultrasound and the different MR techniques led to a moderate correlation of the retrospective gated and prospective triggered MR techniques (eg, r = 0.73 for maximum systolic velocity). The worst correlation was found between the breath-hold technique and retrospective cardiac gating (eg, r = 0.004 for pulsatility index). There was a weak correlation of all three MR sequences compared with ultrasound (r = 0.19-0.60)

CONCLUSIONS

A moderate correlation was found between velocities and indexes measured with the prospective cardiac-triggered phase-contrast MR technique and the retrospective cardiac-gated phase-contrast MR technique. A weak correlation was found between the three different MR techniques and ultrasound, as well as between the breath-hold prospective cardiac-triggered MR sequence and both of the other MR sequences. The influence of temporal and spatial resolution on MR phase-contrast velocity mapping was confirmed.