(Circulation. 2004;109:1987-1993.)
© 2004 American Heart Association, Inc.
Clinical Investigation and Reports |
From the Clinic for Congenital Heart Disease and Institute for Magnetic Resonance Imaging (H.K., P. Barth, H.E., A.P., H.M., P. Beerbaum), Heart- and Diabetescenter Northrhine-Westfalia, Ruhr-University of Bochum, Germany, and Philips Medical Systems (J.G., R.H.), Best, The Netherlands.
Correspondence to Philipp Beerbaum, MD, Clinic for Congenital Heart Disease, Heart- and Diabetescenter Northrhine-Westfalia, Ruhr-University of Bochum, Georgstraße 11, D-32545 Bad Oeynhausen, Germany. E-mail pbeerbaum{at}hdz-nrw.de
Received November 4, 2003; revision received January 27, 2004; accepted February 4, 2004.
| Abstract |
|---|
|
|
|---|
Methods and Results In 14 pediatric patients (mean age 5.2±2.0 years) with cardiac left-to-right shunt, pulmonary (Qp) and aortic (Qs) flow rates were determined by nontriggered free-breathing real-time PC-MRI with single-shot echo-planar imaging combined with sensitivity encoding, which yielded 25 phase images per second at 2.7x2.7-mm in-plane resolution (field of view 30x34 cm2). Over a 9.5-second period that included 2 to 5 respiratory cycles, 16.6±2.6 subsequent stroke volumes (range 13 to 22) were acquired in each vessel. Results were compared with conventional retrospectively ECG-gated PC-MRI. Mean Qp/Qs by conventional PC-MRI was 1.91±0.64, and it was 1.94±0.68 (mean±SD) by real-time PC-MRI. For blood flow rate through pulmonary artery and aorta, we found differences of 2% to 3% (Bland-Altman analysis), with lower limits of agreement of 11% to 13% (mean2 SD) and upper limits of 18% to 19% (mean+2 SD), which demonstrated good agreement between both methods. Mean difference for Qp/Qs was 1%, with limits of agreement ranging between 18% and 22% (mean±2 SD). High repeatability but some flow overestimation was observed in vitro (pulsatile flow phantom) with real-time PC-MRI, whereas conventional PC-MRI was accurate. Beat-to-beat stroke-volume variation was 6.1±2.3% in vivo and 3.7±0.3% in vitro.
Conclusions Beat-to-beat quantification of pulmonary and aortic flows and hence left-to-right shunt within a few seconds is reliable by nontriggered real-time PC-MRI with echo-planar imaging and sensitivity encoding. Good spatial/temporal resolution and a large field of view may render the sequence valuable for multiple applications in congenital heart disease.
Key Words: magnetic resonance imaging blood flow heart defects, congenital pediatrics shunts
| Introduction |
|---|
|
|
|---|
A nontriggered free-breathing PC-MRI technique with image acquisition in real time at frame rates of
25/s (ie, comparable to an echocardiographic frame rate) may overcome these limitations and provide more physiological information, such as the influence of respiration on flow volumes. Current MR techniques to obtain flow data in real time, however, are limited either in spatial5 or temporal resolution6 or use a small field of view (FOV),7 thus largely limiting slice angulation and/or examination of larger subjects because of backfolding.
We evaluated a novel real-time magnetic resonance (MR) phase-contrast pulse sequence using echo-planar imaging (EPI) in combination with sensitivity encoding (SENSE).8 This flow sequence provides both high spatial and high temporal resolution, as is appropriate for examination of great thoracic vessels in pediatric patients, but also allows for a large FOV, as is necessary for larger subjects. The purpose of the study was to validate the method in pediatric patients with cardiac left-to-right shunt. We chose this clinical setting because conventional PC-MRI, which was used as a reference for the real-time sequence in this investigation, has been evaluated recently in these patients for in vitro accuracy and in vivo reliability.9,10
| Methods |
|---|
|
|
|---|
Each patient underwent MRI to measure through-plane flow in the ascending aorta and pulmonary artery. A standard PC-MRI measurement was repeated twice in each location. The mean value of both measurements served as reference and was compared with mean values of stroke volumes sampled by nontriggered free-breathing real-time single-shot echo-planar PC-MRI combined with SENSE (Table 1), over a data acquisition period of 9.5 seconds to include more than 2 respiratory cycles (Figure 1). This was done to provide several subsequent and independently sampled stroke volumes, to determine overall stroke-volume variability with respiration in a free-breathing acquisition. Patients were sedated with midazolam and thiopental intravenously, and blood pressure, oxygen saturation, heart rate, and respiratory rate were monitored.
|
|
In Vitro Study Design
A polyvinyl chloride tube (luminal diameter 12.7 mm; wall thickness 2.4 mm) was embedded in a glass tube (diameter 100 mm, length 340 mm) filled with copper sulfate solution (3 mmol). To approximate blood viscosity,11 the flow phantom was filled with a 60:40 mixture of distilled water and glycerol. The flow phantom was placed parallel to the bore in the isocenter of the magnet. The tubing was connected to a 60-mL artificial left ventricular system (Medos HIA VAD-System II, Medos Medizintechnik AG) to deliver pulsatile flow. Flow pump and MR scanner were triggered simultaneously by an external ECG simulator. With a constant heart rate of 84 bpm, 6 different flow rates (ie, 497, 854, 1127, 1401, 1827, and 2219 mL/min) were achieved by appropriate adjustment of the systolic and diastolic pressures of the ventricular assist device. Flow rate was determined 3 times before and after the MR experiments by stopwatch/graduated cylinder to obtain reference ("true") values. All in vitro MR flow measurements were performed as through-plane measurements with a velocity-encoding value of 150 cm/s and were repeated 3 times.
MR Imaging Technique
All examinations were performed on a 1.5-T MR scanner (Philips, ACS-NT, R8, maximum gradient performance 30 mT/m, slew rate 150 T · m1 · s1). The body coil was used for signal transmission, and a 5-element cardiac phased-array surface coil was used for signal detection. Plan scans in 3 orthogonal planes used ECG-triggered segmented-k-space turbo-field echo to allow double-oblique slice positioning for through-plane flow measurements in the pulmonary artery and ascending aorta. A free-breathing SENSE reference scan (<1 minute of scan time) was performed to obtain sensitivity information necessary for application of the SENSE technique.10,12 Both standard and real-time PC-MRI acquisitions (Table 1) used a phase-correction algorithm provided by the manufacturer.
Real-Time and Standard PC-MRI Quantitative Flow Sequences
We designed a nontriggered, free-breathing, real-time PC-MRI sequence with blipped gradient-switching approach (EPI) combined with SENSE (Table 1). A matrix size of 112x128 and an FOV of 300x340 mm2 resulted in a symmetrical in-plane resolution of 2.7x2.7 mm2. To increase temporal resolution and shorten echo-train length, a SENSE-reduction factor of 4 combined with half-Fourier technique was applied, which reduced the number of actually acquired lines of k-space per repetition time (TR) to 19 (EPI factor). Interleaved velocity encoding was performed by application of a bipolar gradient after the slice-selection gradient (Figure 2). Thus, a complete phase-difference flow image was obtained every 39 ms, ie,
25 flow images per second. An excitation angle of 40° was selected to provide sufficient contrast between moving and stationary spins.
|
As a reference, we used a standard T1-weighted cine phase-contrast gradient-echo sequence as described previously,10 using vector-ECG (VCG) for retrospective gating (Table 1). Both real-time and reference PC-MRI used a higher velocity-encoded value of 200 cm/s versus the in vitro value of 150 cm/s, to avoid phase wrap that may occur in the presence of a large left-to-right shunt in the pulmonary artery.
MR Data Analysis
All flow data were evaluated on an external workstation with a computer algorithm for semiautomatic vessel border detection (P. Barth) as described previously.9 To accelerate postprocessing for real-time PC-MRI, regions of interest were defined over 1 complete heart cycle (magnitude images) and copied to the remaining heart cycles, which left only a few frames for manual adaptation. Simultaneously recorded physiological data (VCG, respiration curves) were used for determination of stroke-volume variability during real-time PC-MRI measurements. Along with the real-time flow measurements (performed over a time period of
9.5 seconds), a VCG curve was recorded. During postprocessing, both curves were aligned. A flow volume that was sampled between 2 consecutive R-wave peaks was defined as a single stroke volume.
In the standard and real-time PC-MRI flow measurements of 1 subject (aged 4.5 years), we determined the signal-to-noise ratio (SNR) in the pulmonary artery and aorta (maximum flow, systole) in magnitude images,13
|
|
and the velocity-to-noise ratio (VNR) in phase images14
|
|
where SImean is mean signal intensity, SDnoise is the SD of noise, v is velocity, Venc is the velocity-encoding value, and SDv is the SD of velocity noise.
Statistical Analysis
All flow values are given as mean values and SD. In vitro flow measurements were compared by linear regression analysis, percent deviation, and paired Students t test. A probability value <0.05 was considered significant. Bland-Altman analysis15 was used to compare conventional with real-time PC-MRI and to determine interobserver variation. Relative stroke-volume variations assessed by real-time PC-MRI flow measurements both in vivo and in vitro were calculated from their mean±SD values.
| Results |
|---|
|
|
|---|
|
Conventional Versus Real-Time PC-MRI
The mean Qp/Qs (pulmonary [Qp] and aortic [Qs] flow rate) ratios determined by conventional and real-time PC-MRI were 1.91±0.64 (range 1.17 to 3.33) and 1.94±0.68 (range 1.05 to 3.17), respectively. Bland-Altman analysis was applied to the log-transformed data because differences increased linearly with mean stroke volumes.15 Estimation of precision of the limits of agreement (defined as mean±2 SD) was based on calculation of 95% CIs. A mean value of 1.0 after antilog transformation is expected in the case of no difference between 2 tested methods.15
When stroke volumes obtained with conventional and real-time PC-MRI were compared, a negligible difference of 3% was present in the pulmonary artery and 2% in the aorta, as reflected by a mean difference value of 1.03 and 1.02, respectively (Table 2). Upper and lower limits of agreement were 1.18 and 0.89 (pulmonary artery) and 1.19 and 0.87 (aorta). Thus, real-time flow measurements may differ from conventional PC-MRI by 18% above and 11% below for the pulmonary artery and by 19% above and 13% below for the aorta. For the flow ratio Qp/Qs, we found a mean difference of 1% (mean value 1.01), with upper and lower limits of agreement of 1.22 and 0.82, respectively.
|
Repeatability and Stroke Volume Variation
Excellent precision with a mean variation of only 2.3±1.9% was found in vivo between the 2 conventional flow measurements used as reference for the real-time MR technique. Real-time flow measurements were performed over 9.5 seconds to include more than 2 respiratory cycles. Depending on the heart rate, 16.6±2.6 (mean±SD; range 13 to 22) truly repeated stroke volumes were available from each measurement (Figure 1). On average, stroke volumes varied from each other by 6.1±2.3% (range 2.6% to 11.8%) in the 27 real-time data sets under sedated free-breathing conditions, with no significant difference between aorta and pulmonary artery.
Interobserver Variability
In 7 randomly selected patients, all stroke volumes in the ascending aorta and pulmonary artery (14 real-time and 28 conventional scans) were evaluated independently and in a blinded manner by 2 experienced observers (H.K. and Kerstin Berkemeier, RN). The mean difference for the averaged pulmonary stroke volume acquired with real-time PC-MRI was 5% (mean 0.95), with upper and lower limits of agreement of 1.09 (95% CI 1.00 to 1.19) and 0.83 (95% CI 0.76 to 0.91). For the averaged aortic stroke volume, the difference was 4% (mean 0.96), with upper and lower limits of agreement of 1.06 (95% CI 0.99 to 1.13) and 0.87 (95% CI 0.82 to 0.93). Negligible interobserver differences were found for conventional PC-MRI in the pulmonary artery, with a mean difference of 1% and upper and lower limits of agreement of 1.07 (95% CI 1.03 to 1.11) and 0.96 (95% CI 0.93 to 0.99), as well as in the aorta, with a mean difference of 1% and upper and lower limits of agreement of 1.05 (95% CI 1.01 to 1.08) and 0.94 (95% CI 0.91 to 0.98), respectively.
SNR and VNR
SNR and VNR values13,14 are expressed as percentages of corresponding values of standard PC-MRI (which equals 100%). With real-time PC-MRI, SNR in the magnitude images was reduced to 13% (6%) in the aorta (pulmonary artery), whereas VNR in the phase images was only reduced to 93% (55%) of standard PC-MRI (Figure 4).
|
In Vitro Flow Measurements
Repeatability (ie, precision) of the manually obtained ("true") flow values determined by stopwatch and graduated cylinder was excellent (error <1%), as was repeatability of the conventional PC-MRI measurements (error <2.2%). Pulsatile flow-phantom measurements with conventional and real-time PC-MRI correlated well with true flow, as reflected by the correlation coefficients of 0.999 for both techniques. The y-intercepts were negligible. Although the slope of 1.01 for conventional PC-MRI reflects close identity with "true" flow values (y=0.01+1.01x; relative differences 0.4% to 4.7%; P=NS), a moderate overestimation (systematic error) was observed for real-time PC-MRI (y=0.03+1.26x; relative differences 12.8% to 28.1%; P<0.006). During the 9.5-second data sampling period, 12 to 13 flow curves were registered. On average, stroke volumes differed by 3.7±0.3%.
| Discussion |
|---|
|
|
|---|
More recently, novel nontriggered real-time MR flow techniques have been introduced to save time and to overcome limitations of the conventional method.58 With real-time MR flow as opposed to the conventional techniques, no respiratory gating or ECG triggering is needed, and as many stroke volumes or flow volumes as needed may be acquired in a row within a few seconds. This makes real-time MR flow robust against motion artifacts, arrhythmia, or improper ECG signal detection, and repeated measurements at multiple sites may be possible within a short time.58
Furthermore, flow velocity mapping in real time may provide unique quantitative hemodynamic information that is not available by any other imaging technique. This was recently demonstrated by Hjortdal et al7 using real-time MR flow to determine the effects of respiration and exercise on blood flow in total cavopulmonary connection (ie, Fontan-type circulation) in patients with a single-ventricle physiology.
Current Obstacles in Use of Real-Time MR Flow
A more widespread use of real-time MR flow in clinical routine, however, is hampered by important technical limitations of currently available techniques. Hjortdal et al7 used a real-time phase-contrast segmented k-space fast-field-echo MR technique with echo-planar readout. Although they were able to achieve adequate temporal and spatial resolution (phase-difference images every
50 ms, pixel size 3.4x2.1 mm), this was only possible with a very small FOV of 9.0x13.6 cm.7 Therefore, owing to limitations from backfolding, the sequence is not suited for double angulation of slice position as is usually mandatory in the MR evaluation of congenital heart disease or for examination of larger subjects. Moreover, validation data for their MR flow method are lacking.7 Klein et al6 used a comparable MR sequence with a larger FOV in their adult population, but at the expense of considerably reduced temporal resolution (phase images every 124 ms, ie, 8 images/s). Eichenberger et al5 used real-time phase-contrast single-shot EPI. They achieved good temporal resolution (phase intervals of 37 ms, ie, 27 images/s) at the expense of limited spatial resolution (5.0x5.0 mm in-plane), which is likely to introduce systematic error from partial volume effects in pediatric patients.16
Real-Time MR Flow With EPI and SENSE
Clinical Considerations
In this context, we present the first real-time MR flow quantification technique that combines both improved spatial (in-plane resolution 2.7x2.7 mm) and temporal (phase interval 39 ms) resolution, which is considered adequate for flow measurements in the great thoracic vessels in pediatric patients.9 The sequence also allows for a relatively large FOV and thus for more flexible applications in both pediatric and adult congenital heart disease. This is achieved by the use of single-shot EPI and a high SENSE-reduction factor of 4.
We validated this real-time MR flow sequence in comparison with standard retrospectively ECG-gated PC-MRI and demonstrated the reliability in the quantification of pulmonary and aortic flow rates and hence left-to-right shunt in pediatric congenital heart disease. Although a stroke volume can be assessed in a single heart beat8 by application of this real-time technique, scan time was extended to
9.5 seconds to include more than 2 respiratory cycles and thus allow comparison of the mean stroke volumes with standard (reference) PC-MRI, by which an averaged stroke volume is acquired over multiple heart beats (Figure 3). There was almost no difference, and there was limited scatter of measurement results for aortic and pulmonary flow and for the ratio QP/QS compared with a standard PC-MRI protocol with a comparable spatial resolution.
Another result of the present study is the small variation (6.1±2.3%) of stroke volumes with respiration in both great arteries in sedated pediatric patients. In light of this, scan time may be shortened to only 1 to 2 seconds (instead of 9.5 seconds) to sample 1 to 3 stroke volumes. In contrast, a 10% decrease of aortic flow during inspiration was described in healthy adults.17 This demonstrates that the impact of respiration on flow in thoracic vessels must be considered when novel techniques for flow quantification are evaluated. For example, it is unknown at present whether and to what extent flow rate in the great thoracic vessels is attenuated by breath-holding. This question can now be addressed with real-time MR flow. Another advantage of real-time over conventional gated MR flow is that cardiac arrhythmia should no longer preclude safe and reliable flow measurements. It is now possible to analyze such flow data during postprocessing, because VCG curves displaying the arrhythmia and flow curves can be aligned. This enables the selection of usable flow curves for diagnostics in that, for example, a ventricular ectopic beat can be identified and such flow values rejected. In contrast, such flow data are usually nondiagnostic when conventional ECG-gated PC-MRI is used because of errors in phase reordering. This is important given that cardiac arrhythmia is frequently encountered in adult patients with congenital heart disease, who benefit most from noninvasive follow-up by MRI because Doppler echocardiography is often hampered by limited acoustic windows. Real-time MR flow provides a novel method to assess hemodynamics in this subgroup of patients.
The available degree of freedom in slice angulation and the good spatial/temporal resolution may extend the clinical and research applications of real-time MR quantitative flow in both pediatric and adult cardiovascular disease. It may be possible, for instance, to quantitatively determine immediate changes in the pulmonary flow rate in response to pharmacologically induced vasodilation in patients with pulmonary hypertension resulting from increased pulmonary arteriolar resistance. This may be especially useful in an interventional MR setting, allowing simultaneous pressure recording and flow assessment to determine rapid changes in pulmonary vascular resistance.18
Technical Considerations
In the SENSE approach, an array of multiple, simultaneously operated receiver coils for signal acquisition (eg, cardiac 5-element phased-array surface coil) is used to spatially encode the MR signal as a function of the sensitivity of each coil element.12 This allows the number of phase-encoding steps to be reduced by increasing the distance of readout lines in k-space (ie, parallel imaging) and thus accelerates acquisition, but at the expense of a lower SNR.12 Compared with standard PC-MRI, a significant loss of SNR in real-time PC-MRI magnitude and phase images is present (Figure 4). This not only relates to the SNR penalty bounded by square root to the number of samples acquired but also contains noise from sensitivity maps as propagated in the SENSE reconstruction process, especially at higher SENSE reduction factors, such as those used in the present real-time sequence.12 However, Conturo and Smith19 demonstrated that phase images have an inherently better dynamic range than magnitude images that result in a signal-to-noise advantage in clinical flow imaging. They proposed that useful phase-angle information can be extracted even when magnitude images are very noisy, such as with rapid data collection techniques.19 This may explain why VNR apparently was not a major limitation with real-time PC-MRI. Moreover, surprisingly, no obvious degrading effects of residual eddy currents (high-gradient switching rate with EPI) and concomitant fields (high-gradient amplitude) were obvious in the in vivo measurement results of the present study.
The main advantage of combining single-shot EPI with SENSE is to drastically reduce echo-train length while preserving high imaging speed.8 To achieve this goal, we used a high SENSE reduction factor of 4. This may be beneficial for 3 reasons. First, a higher SNR is achievable when SENSE is used because EPI data are collected during the early part of T2* decay (which contains more signal). In contrast, when single-shot EPI is used without SENSE (ie, long echo-trains) data collection occurs in part beyond T2* decay (which contains little signal). Therefore, some of the signal-to-noise penalty from the use of higher SENSE reduction factors is compensated. Second, phase errors/spatial distortions that result from off-resonance effects such as magnetic susceptibility and chemical shift20 or that are caused by local B0-field inhomogeneities21,22 are reduced with a shorter echo-train. Third, with a shorter echo-train, considerably shorter effective echo times are possible, thus diminishing the influence of turbulent flow. This may improve flow accuracy and image quality.
To avoid severe deterioration of image quality from backfolding artifacts (due to the SENSE reconstruction process), no further foldover should occur during planning of the desired FOV. Therefore, we chose a relatively large FOV of 30x34 cm2 to enable the use of a standardized real-time PC-MRI protocol at different body sizes and slice angulations with respect to geometric properties and gradient performance. As a result, however, the impact of SENSE reconstruction may have been variable, depending on angulation and body size as related to the FOV.
Although smaller flip angles are generally thought to reduce systematic error in quantitative flow assessment from partial volume effects by decreasing the signal magnitude difference between flowing blood and surrounding tissue,23 a more inflow-dependent acquisition with a flip angle of
=40° was considered important to enable clear vessel border detection in the corresponding magnitude images, thus avoiding major bias from displacement of regions of interest.
In Vitro Performance
Although we found high precision (repeatability) in vitro for real-time PC-MRI, some flow volume overestimation was present in the flow-phantom experiments, most likely because of partial volume effects, chemical shift, and susceptibility artifacts (eg, tube wall),16 for which single-shot EPI is known to be sensitive.20 The use of large (3-cm diameter)8 or very thin-walled (0.1 mm)23 tubes and nonpulsatile low-flow conditions8 has been shown to improve in vitro accuracy but does not resemble the pediatric measurement environment. Therefore, assessment of the accuracy of this real-time PC-MRI technique remains difficult, because both in vitro and in vivo, there is a lack of a true "gold standard."
Study Limitations
No conclusions from the present data are applicable to children with valvular stenosis and turbulent blood flow. The borderline in-plane resolution of both real-time and reference PC-MRI may render flow volume measurements in smaller vessels at risk for flow underestimation due to partial volume effects. Only relative changes should be emphasized in SNR/VNR measurements, because absolute values may be influenced by interpatient variations and background filtering in the reconstructed images, which reduce the noise component. Although real-time magnitude and phase-image quality were acceptable, skilled observers and computer-assisted vessel border detection software were needed. However, some region-of-interest displacement may have occurred in heart phases with slow flow, thus contributing to the small but statistically significant interobserver bias. In the in vitro experiments, the FOV was too large to ensure that 4-fold undersampling would lead to actual foldover in the single coil data, and therefore conclusions regarding the impact of SENSE acquisition and reconstruction on phase encoding of flow volume should be drawn with caution. An adopted FOV with a reduced matrix for a similar voxel size, however, would have been different from the real-time PC-MRI protocol used in vivo.
| Conclusions |
|---|
|
|
|---|
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
V. E. Hjortdal, T. D. Christensen, S. H. Larsen, K. Emmertsen, and E. M. Pedersen Caval Blood Flow During Supine Exercise in Normal and Fontan Patients Ann. Thorac. Surg., February 1, 2008; 85(2): 599 - 603. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2004 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |