|
|
||||||||
Eur J Cardiothorac Surg 2005;27:462-467
© 2005 Elsevier Science NL
Division of Image Processing, Departments of Radiology, Cardiology and Cardio-thoracic Surgery, Leiden University Medical Center, 2333ZA Leiden, The Netherlands
Received 31 August 2004; received in revised form 12 November 2004; accepted 22 November 2004.
* Corresponding author. Tel.: +31 71 526 4846; fax: +31 71 526 6801. (E-mail: j.j.m.westenberg{at}lumc.nl).
| Abstract |
|---|
|
|
|---|
Key Words: Magnetic resonance imaging Mitral valve Regurgitation Flow quantification Image processing
| 1. Introduction |
|---|
|
|
|---|
Besides measurement restrictions, geometrical models are applied to deduct physiological parameters describing the cardiac condition from the echo Doppler measurements. These geometrical models are not generally applicable to all individual subjects.
Furthermore, imaging the region of interest in TTE can be complicated by the acoustic attenuation from structures inside the thorax (lungs, subcutaneous tissues and ribs) and the heart (prosthetic valve construction materials and calcification). TEE is less limiting, but this technique is semi-invasive and patients often experience this as a discomforting examination.
Magnetic resonance imaging (MRI) is a non-invasive and patient-friendly technique, which is applied for the determination of global and regional left ventricular anatomy and function [13]. As a three-dimensional (3D) imaging technique, volumetric measurements are performed directly and do not rely on model assumptions. Velocity-encoded cine MRI provides quantitative information on moving spins and can be applied to determine intra-ventricular blood flow [14,15]. These qualifications make MRI a very suitable modality for quantifying trans-valvular blood flow.
Previous attempts have all encountered the motion of the heart during contraction and relaxation as the main obstacle towards direct acquisition of the trans-valvular flow. The mitral annulus moves 1216mm towards the apex during contraction [16]. Kayser et al. showed that correction for this through-plane motion is indispensable [17]. This correction is performed retrospectively. But since, the flow quantification is performed in a single, fixed plane during the complete cardiac cycle, the location of the data acquisition is not identical to the plane of the MV during the whole cardiac cycle. This implies that flow quantification does not represent the true trans-valvular flow.
An acquisition method with a motion-adapted acquisition plane is described by Kozerke et al. [18]. The accurate positioning of the acquisition plane enables measurement of the true valvular flow during the whole cycle. This technique requires dedicated scanning software not widely available on all commercial MR scanners, and the prospective triggering that is necessarily used implies inaccurate flow quantification at a regurgitant mitral valve. With prospective triggering, the final part of the RR-interval is not acquired and therefore, the diastolic flow quantification will not be complete.
In this study, a new acquisition method is introduced using velocity-encoded MRI to obtain the velocity vector field of the intra-ventricular blood flow during the complete cardiac cycle. In this 3D velocity vector field, the MV-plane is indicated retrospectively. The vector components of the flow velocity going through the MV-plane will be reconstructed which results in the trans-valvular flow velocity. Regurgitation will be evident as the flow fraction passing the valve during systole, flowing from the left ventricle into the left atrium. In 10 patients with severe mitral valve regurgitation and selected for mitral valve repair, the MV-flow will be quantified and the regurgitant volume will be determined.
| 2. Material and methods |
|---|
|
|
|---|
2.2. MRI acquisition
MRI was performed on a 1.5 T scanner (ACS-NT15 Gyroscan with the Powertrack 6000 gradient system; Philips Medical Systems, Best, The Netherlands). The intra-ventricular blood flow was sampled with a radial stack of six acquisition planes parallel to the long-axis of the left ventricle (Fig. 1). The specific 3-directional velocity-encoded MRI [19] scanning parameters and the image processing methods to obtain the 3D intra-ventricular flow velocity vector field from the MRI data were described previously [20]. In each consecutive imaging plane, velocity values were acquired in three directions, with a maximum velocity sensitivity value of 100cm/s in each direction. The reconstructed velocity values in three directions represent the velocity vector components. Slice thickness of the imaging planes was 8mm, Field-of-view=370mm (60% rectangular), scan matrix=128x102. The flip angle
=10°, TR/TE=5.8/3.5; two signal averages were used to increase the signal-to-noise ratio. For each subject, the intra-ventricular flow velocity vector field was obtained in 30 phases during the cardiac cycle (obtained with retrospective cardiac synchronisation). The spatial resolution was 2.89x2.89x8.0mm3. Typical examination time was 912min, depending on the subject's heart rate.
|
|
2.4. Validation
For each subject, the MV-flow volume determined with the new method described above (from hereon noted as the 3-dir MV-flow) was validated by comparing it with the flow volume measured with MRI at the ascending aorta (described by van der Geest et al. [21]), referred to as the AO-flow volume. The same MRI acquisition protocol was used as for the MV-flow, but now with a single imaging plane positioned perpendicular to the ascending aorta. Only velocity values in through-plane direction, perpendicular to the acquisition plane, were considered. The aortic flow was determined using the FLOW software (Medis, Leiden, The Netherlands) [21].
Furthermore, for each subject, the 3-dir MV-flow volume was compared with the conventional one-directional MV-flow volume measurement with through-plane motion correction [17]. The same imaging protocol was used, with the single imaging plane positioned at the MV during end-systole. Again, only the velocity values in through-plane direction were considered. The through-plane velocity of the myocardium was assessed in a region-of-interest, drawn inside the myocardium for each of the phases, by measuring the mean velocity inside this region-of-interest. This velocity value was subtracted from the velocity measured at the MV annulus. The MV-flow volume measured with this method is referred to as the 1-dir MV-flow volume.
For ten volunteers and ten patients, the MV-flow volume was determined with the 3-dir MV-flow method, as well as with the 1-dir MV-flow method, and compared with the AO-flow. The regurgitant flow volume in the patients was determined by calculating the flow volume during systole going from the left ventricle into the left atrium. The regurgitant flow fraction, a measure for the severity of regurgitation, was determined by the ratio between the regurgitant flow volume during systole and the inflow volume through the MV during diastole.
2.5. Statistical analysis
Correlations between the AO-flow volume measurement and both the conventional 1-dir MV-flow volume and the new 3-dir MV-flow volume were tested with the Pearson correlation coefficient rP, under the assumption that the parameters were distributed normally. This assumption was tested with KolmogorovSmirnov tests. The approach described by Bland and Altman [22] was followed to study the differences between both MV-flow volume measurements and the AO-flow volume. The mean differences were determined and the statistical significance of these differences was tested by a paired-samples t-test.
| 3. Results |
|---|
|
|
|---|
|
Next, the new method for measuring the MV-flow volume was performed on the ten patients selected for mitral valve repair. The regurgitant flow fractions were between 3 and 30% (mean±standard deviation=16±8%) in these patients. The results of the MV-flow measurements are presented in Fig. 4A. Again, KolmogorovSmirnov tests were performed and prove that the flow values can be considered as distributed normally (P=0.98 for 1-dir MV-flow, P=0.91 for 3-dir MV-flow, P=0.68 for AO-flow). The correlation between the 1-dir MV-flow volume and the AO-flow volume is not statistically significant: (rP=0.08, P=0.82), whereas the correlation between the 3-dir MV-flow volume and the AO-flow volume is statistically significant: (rP=0.90, P<0.01).
|
| 4. Discussion |
|---|
|
|
|---|
In this study, a new MRI method is used to provide information on the exact regurgitant flow fraction. MRI is a reliable non-invasive imaging modality, readily applied for studying global and regional left ventricular anatomy and function [13]. Velocity-encoded MRI enables intra-ventricular flow quantification. Kayser et al. used single-slice 1-directional velocity-encoded MRI for measuring the flow at heart valves [16]. From their data, through-plane motion correction was proven to be necessary. But this is not the only condition for acquiring accurate transvalvular blood flow. Adaptation of the acquisition plane to the motion of the valve of interest needs to be performed, as is illustrated by the data of our study. In the volunteers, a statistically significant over-estimation of 25±22ml of the MV-flow volume was found when performing the conventional 1-dir MV-flow method. Furthermore, there was no correlation with the AO-flow volume. Also for the patient data, the same results were found. The difference between 1-dir MV-flow volume and AO-flow volume (19±28ml) was not statistically significant (P=0.06) due to the large standard deviation. The large spread in data also indicates the inaccuracy and unreliability of this method.
The new 3-dir MV-flow method provides small, non-significant differences compared with the AO-flow volume and good agreement between the flow measured at both locations. In a previous study, this 3-dir MV-flow method was subjected to an intra- and inter-observer variation study [20]. A high reproducibility was found, excellent correlation between repeated analyses and the differences between analyses were very small and not statistically significant.
When measuring the MV-flow volume and using the AO-flow volume as a standard of reference, the flow to the coronaries has to be taken into account. The AO-flow was acquired at a level distal to the branches of the coronaries. This implies that the AO-flow volume measured at this location will be smaller than the MV-flow volume. In our study, this small systematic difference was not detectable. The contribution of the flow to the coronaries, only 0.5% of the cardiac output [24], was too low to be detected by this method.
The new 3-dir MV-flow method enables accurate and reproducible quantification of the true trans-valvular MV-flow in a patient-friendly and easy-to-use manner. The three-dimensional velocity vector field of the intra-ventricular blood flow enables quantification of asymmetric regurgitant flow jets, a limitation of TEE which can miss asymmetric jets due to prolapse of the commissurae.
The new method can be applied on most commercially available MRI scanners. The acquisition of the radial stack of 6 imaging planes takes 912min and depends on the patient's heart rate. Scout images, standard 2- and 4-chamber views and a set of short-axis images are needed to plan the radial stack of slices. This process takes approximately 1520min of examination time. The manual post processing is currently very laborious and needs to be automated before this method can be introduced in daily clinical practice. The reconstruction of the velocity vector field and the manual contour drawing, all needed for the flow analysis, can take up to 1h per patient. Automating the manual processing (i.e. the indication of the mitral valve in the long-axis images as well as the contour drawing in the reconstructed trans-valvular velocity images) will result in a clinically applicable method.
| Appendix A. Conference discussion |
|---|
|
|
|---|
Dr Westenberg: Yes, I think that is very much possible. At this point this method is being used for scientific purposes, not clinically generally applied yet, but I think the possibilities are certainly there to do this in a clinical routine.
Dr P. Kolh (Liege, Belgium): Have you compared and validated your results with an invasive method such as, for example, the conductance catheter?
The Leiden group, Baan and Steendijk for example, has extensively published on the use of the conductance catheter to generate pressure-volume loops. So I wonder if you have any plan to validate your measurements with this technique?
Dr Westenberg: At this point that is not being done but we are indeed going into that direction by measuring flows with the same technique in the aorta and compare that with intravascular measured flows, that's right.
Dr R. De Simone (Heidelberg, Germany): Actually the quantification of mitral regurgitation with MRI is nothing new; there are many groups working on this. I have two questions for you. First, you stated that your measurement, the positioning of the mitral valve, was taken manually. Do you think that you can define manually very well the complex geometry of the mitral orifice? You showed also in your data that the orifice changed during systole. This is the first question.
The second question, you just used six planes, and we saw that the velocities between the planes were interpolated. How accurate are your measurements? You know that in our experience we use 90 slices in order to define the complex geometry of the regurgitant jets and we were able to show the data. And just a comment. By the way, echocardiography is not based on a model assumption. It is based on real measurements
Dr Westenberg: First I want to come back to that last remark. Indeed, echocardiography is a real measurement, too, but to quantify the amount of regurgitation, for instance, there has to be applied, for the PISA method, for instance, some model assumptions.
Now coming back to your questions, indeed, we are manually indicating the mitral valvular plane, but as you can see in the movies that I showed here, it is very evident to find the mitral valvular annulus in these images.
And for your second question, we are using six acquisition planes. That is because there is some time constriction involved. It is indeed true if you use more acquisition planes that you have to do less interpolation. We use triangular interpolation between these six planes and more measurements will be more accurate. But as you can see from the results that I have presented, there are no large statistical over or under estimations when we are doing it this way with six planes.
Dr R. Frater (Bronxville, New York, USA): I am sure you recognize that what people are worrying about is whether this isn't truly accurate. It is an extraordinarily beautiful demonstration of how to get the data and so on. I would suggest really strongly you go back to 30, 40 years ago, put a flow probe above the mitral valve, and you can quantify it absolutely precisely and then see if that correlates with these beautiful images and calculations that you have been doing. It is the gold standard, after all, an electromagnetic flow probe over the mitral valve.
Dr Westenberg: Yes, thank you, but we didn't go that invasive with this study.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. A. Marsan, J. J.M. Westenberg, C. Ypenburg, V. Delgado, R. J. van Bommel, S. D. Roes, G. Nucifora, R. J. van der Geest, A. de Roos, J. C. Reiber, et al. Quantification of Functional Mitral Regurgitation by Real-Time 3D Echocardiography: Comparison With 3D Velocity-Encoded Cardiac Magnetic Resonance J. Am. Coll. Cardiol. Img., November 1, 2009; 2(11): 1245 - 1252. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. B. Srichai, R. P. Lim, S. Wong, and V. S. Lee Cardiovascular Applications of Phase-Contrast MRI Am. J. Roentgenol., March 1, 2009; 192(3): 662 - 675. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. M. Westenberg, S. D. Roes, N. Ajmone Marsan, N. M. J. Binnendijk, J. Doornbos, J. J. Bax, J. H. C. Reiber, A. de Roos, and R. J. van der Geest Mitral Valve and Tricuspid Valve Blood Flow: Accurate Quantification with 3D Velocity-encoded MR Imaging with Retrospective Valve Tracking Radiology, December 1, 2008; 249(3): 792 - 800. [Abstract] [Full Text] [PDF] |
||||
![]() |
Authors/Task Force Members, A. Vahanian, H. Baumgartner, J. Bax, E. Butchart, R. Dion, G. Filippatos, F. Flachskampf, R. Hall, B. Iung, et al. Guidelines on the management of valvular heart disease: The Task Force on the Management of Valvular Heart Disease of the European Society of Cardiology Eur. Heart J., January 26, 2007; (2007) ehl428v1. [Full Text] [PDF] |
||||
![]() |
R. De Simone, I. Wolf, S. Mottl-Link, R. Hoda, B. Mikhail, F.-U. Sack, H.-P. Meinzer, and S. Hagl A clinical study of annular geometry and dynamics in patients with ischemic mitral regurgitation: new insights into asymmetrical ring annuloplasty Eur. J. Cardiothorac. Surg., March 1, 2006; 29(3): 355 - 361. [Abstract] [Full Text] [PDF] |
||||
![]() |
T A M Kaandorp, H J Lamb, E E van der Wall, A de Roos, and J J Bax Cardiovascular MR to access myocardial viability in chronic ischaemic LV dysfunction Heart, October 1, 2005; 91(10): 1359 - 1365. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ANN THORAC SURG | ASIAN CARDIOVASC THORAC ANN | EUR J CARDIOTHORAC SURG |
| J THORAC CARDIOVASC SURG | ICVTS | ALL CTSNet JOURNALS |