Volumetric assessment of kinetic energy and vorticity in the pulmonary artery: alteration of flow hemodynamics in patients with repaired tetralogy of Fallot using 4D flow MRI
Julio Garcia1, Silvia Hidalgo Tobon2,3, Guadalupe Sagaon Rojas4, Benito de Celis Alonso5, Manuel Obregon2, Porfirio Ibanez2, Julio Erdmenger6, and Pilar Dies-Suarez2

1Radiology, Northwestern University, Chicago, IL, United States, 2Investigacion en Imagen y Resonancia Magnetica Nuclear, Hospital Infantil de Mexico Federico Gomez, Mexico City, Mexico, 3Physics, Universidad Autonoma Metropolitana, Mexico, Mexico, 4Physics, Universidad Autonoma Metropolitana, Mexico City, Mexico, 5Faculty of Physics and Mathematics, Benemérita Universidad Autónoma de Puebla, Puebla, Mexico, 6Pediatric Cardiology, Hospital Infantil de Mexico Federico Gomez, Mexico City, Mexico

Synopsis

Flow alterations in the pulmonary artery (PA) of patients with repaired tetralogy of Fallot (rTOF) may be link with elevated kinetic energy (KE). 4D flow MRI allows for the non-invasive volumetric assessment of flow hemodynamics, vorticity, and KE in patients with rTOF in the pulmonary (PA). Thus, the aim was to investigate the impact of flow alterations in the PA and its association with KE and vorticity.

Purpose:

Flow alterations in the pulmonary artery (PA) of patients with repaired tetralogy of Fallot (rTOF) may be link with the expending of kinetic energy (KE)1. Furthermore, complex flow patterns (vortical and helical) may be associated with the energy dissipation within the PA2 and can be evaluated by the mean of flow vorticity3. Three-dimensional time-resolved phase contrast MRI (4D flow) allows for the non-invasive volumetric assessment of flow hemodynamics, vorticity, and KE in patients with rTOF in the PA. Thus, the aim was to investigate the impact of flow alterations in the PA and its association with KE and vorticity.

Methods:

15 pediatric patients with rTF (age=9±6 yrs, 6 females) underwent aortic 4D flow MRI as part of an IRB-approved protocol. 4D flow MRI4 was performed at 1.5T (Philips, Achieva, Best, The Netherlands) with full 3D coverage of the thoracic aorta and PA (spatial resolution=2.5×2.1×3.2 mm3; temporal resolution=40-50 ms) using prospective PPU and respiratory navigator gating. Pulse sequence parameters were as follows: 1.5 T scan parameters ranged from TE/TR=2.3–3.4/4.8–6.6 ms, flip angle α=15°, Venc= 1.5-5 m/s, and a field of view of 340–400×200–300 mm. 4D flow dataset pre-processing5 include: eddy-current correction, flow aliasing, and calculation of 3D phase contrast angiography (3D PC-MRA). A 3D segmentation of the PA (Figure 1B) was obtained from the 3D PC-MRA using Matlab (The Mathworks, Natick, MA, USA) and was used to mask the PA velocity field and for flow pattern visualization in the PA. Masked velocity field was used to calculate KE (KE=1/2×rho×v2, were rho is the blood density = 1.06 g/mL and v the velocity field), and vorticity (ω= curl (v)). Maximum intensity projections (MIPs) were calculated for flow velocity, KE, and vorticity. Flow measurements were perform the main pulmonary artery (MPA), right pulmonary artery (RPA), and left pulmonary artery (LPA) using manually located analysis planes (Figure 1C). From each plane (MPA, RPA, LPA) the pulmonic lumen was segmented and used to extract following flow hemodynamic parameters: peak velocity (PV), maximal flow (Qmax), mean flow (Qmean). Volumetric median of KE was used to divide the evaluated cohort in two groups: 1) low KE and 2) elevated KE. The association of kinetic energy at MPA, RPA, and LPA with other flow parameters were assessed by Pearson’s correlation. Comparison between low and elevated KE groups was performed by Mann-Whiteney test.

Results:

Maximal, mean KE in the pulmonary artery showed a global correlation with PV (r=0.47, p=0.008; r=0.38, p<0.037), Qmax (r=0.49, p<0.005; r=0.45, p<0.014), and Qmean (r=0.49, p<0.006; r=0.44, p<0.015). Both maximal and mean KE were mainly originated from the RPA where associations with PV (r=0.87, p=0.001; r=0.84, p<0.002), Qmax (r=0.77, p<0.01; r=0.75, p<0.013), and Qmean (r=0.69, p=0.028; r=0.69, p=0.027) were more important. Maximal KE was 59% higher in the MPA than in the RPA, as well as mean KE with 33% increment. Flow distribution was the major contributor to these correlations and increment in the RPA for KE. Velocity, KE, and vorticity MIPs (Figure 1E) allowed to identify regions proximal to pulmonary bifurcation with elevated KE, and vortical flow in the dominant direction of the flow. PA volumetric mean and median KE were associated with volumetric mean vorticity, r=0.78 (p<0.001) and r=0.44 (p<0.001) respectively. When comparing low and elevated KE significant differences were found for volumetric mean KE (0.029±0.019 mJ vs. 0.047±0.022 mJ, p=0.02), median KE (0.041±0.012 mJ vs. 0.07±0.02 mJ, p<0.001), and mean vorticity (0.032±0.008 1/s vs. 0.037±0.006 1/s, p<0.04).

Discussion and Conclusions:

In this pilot study, maximal and mean KE in the RPA was associated with flow hemodynamic parameters, whereas KE in the MPA and LPA were not. This observation was explain by flow distribution within the PA and the regions (proximal PA bifurcation) where elevated energy dissipation occurs. A large cohort study is needed to evaluate the clinical usefulness of KE to survey patients with rTOF.

Acknowledgements

No acknowledgement found.

References

1. Jeong D, Anagnostopoulos P, Roldan-Alzate A, et al. Ventricular kinetic energy may provide a novel noninvasive way to assess ventricular performance in patients with repaired tetralogy of Fallot. J Thorac Cardiovasc Surg 2015; 149: 1339-47.

2. Francois CJ, Srinivasan S, Schiebler ML, et al. 4D cardiovascular magnetic resonance velocityu mapping of alterations of right heart flow patterns and in main pulmonary artery hemodynamics in tetralogy of Fallot. J Cardiovasc Magn Reson 2012; 14:16.

3. J. Garcia, E. Larose, P. Pibarot, L. Kadem. On the Evaluation of Vorticity using Cardiovascular Magnetic Resonance Velocity Measurements. J Biomech Eng 2013, 135(12), 12450. IF: 2.0. PMID: 24026138.

4. Markl M, Harloff A, Bley TA, et al. Time-resolved 3D MR velocity mapping at 3T: improved navigator-gated assessment of vascular anatomy and blood flow. J Magn Reson Imaging 2007; 25: 824-31.

5. Block J, Kreher BW, Henning J, Markl M. Optimized pre-processing of time-resolved 2D and 3D phase contrast MRI data. In. 15th Sci Meet Int Soc Magn Reson Med 2007:3138.

Figures

Figure 1. Data pre-processing. Panel A shows an example of acquired datasets. Panel B shows and example of pulmanary artery segmentation. Panel C shows the analysis planes locations and pathline flow visualization. Panel D shows an example of maximum intensity projections for flow velocity, kinetic energy and vorticity.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2719