- Research
- Open access
- Published:
Stent interventions for pulmonary artery stenosis improve bi-ventricular flow efficiency in a swine model
Journal of Cardiovascular Magnetic Resonance volume 23, Article number: 13 (2021)
Abstract
Background
Branch pulmonary artery (PA) stenosis (PAS) commonly occurs in patients with congenital heart disease (CHD). Prior studies have documented technical success and clinical outcomes of PA stent interventions for PAS but the impact of PA stent interventions on ventricular function is unknown. The objective of this study was to utilize 4D flow cardiovascular magnetic resonance (CMR) to better understand the impact of PAS and PA stenting on ventricular contraction and ventricular flow in a swine model of unilateral branch PA stenosis.
Methods
18 swine (4 sham, 4 untreated left PAS, 10 PAS stent intervention) underwent right heart catheterization and CMR at 20 weeks age (55 kg). CMR included ventricular strain analysis and 4D flow CMR.
Results
4D flow CMR measured inefficient right ventricular (RV) and left ventricular (LV) flow patterns in the PAS group (RV non-dimensional (n.d.) vorticity: sham 82 ± 47, PAS 120 ± 47; LV n.d. vorticity: sham 57 ± 5, PAS 78 ± 15 p < 0.01) despite the PAS group having normal heart rate, ejection fraction and end-diastolic volume. The intervention group demonstrated increased ejection fraction that resulted in more efficient ventricular flow compared to untreated PAS (RV n.d. vorticity: 59 ± 12 p < 0.01; LV n.d. vorticity: 41 ± 7 p < 0.001).
Conclusion
These results describe previously unknown consequences of PAS on ventricular function in an animal model of unilateral PA stenosis and show that PA stent interventions improve ventricular flow efficiency. This study also highlights the sensitivity of 4D flow CMR biomarkers to detect earlier ventricular dysfunction assisting in identification of patients who may benefit from PAS interventions.
Introduction
Branch pulmonary artery (PA) stenosis (PAS) commonly occurs in patients with complex congenital heart disease (CHD). Acute post-surgical PAS is associated with hemodynamic instability, prolonged in-hospital recovery and increased mortality [1,2,3,4,5]. Chronic PAS is associated with abnormal PA growth and pulmonary blood flow mal-distribution and can lead to PA hypertension [6, 7]. These vascular pathologies contribute to reduced exercise capacity and progressive right ventricular (RV) dysfunction thought to be related to chronic ventilation-perfusion mismatch and increased RV afterload [8]. For many patients the first line therapy for treatment of PAS includes catheter interventions with balloon angioplasty and intravascular stenting. Interventional outcome measures have long focused on the technical and anatomic success of the procedure, defined principally as an immediate post-procedure increase in PA size [9]. Thus far, there has been limited information detailing the impact of intravascular stenting on the ventricular response to these interventions.
In CHD patients with PAS standard measures of ventricular function are often normal until late in the disease process even with significant pulmonary blood flow mal-distribution [9] and indications for PAS interventions are not clearly defined. Ventricular flow analysis using four dimensional (4D) flow cardiovascular magnetic resonance (CMR; 4D flow CMR) is an increasingly utilized method to study ventricular function in patients with CHD before and after complex interventions [10,11,12,13,14]. Measurement of time-varying velocity of blood flow with 4D flow CMR can precisely define ventricular flow patterns and permit calculation of biomarkers of mechanical efficiency including kinetic energy and vorticity [15, 16]. Changes in these ventricular flow biomarkers are believed to represent early signs of cardiac dysfunction [17] as momentum lost through inefficient ventricular flow increases ventricular energy demand and contributes to adverse ventricular remodeling [18, 19]. We believe the advanced imaging capabilities of 4D flow CMR can detect early markers of ventricular dysfunction caused by PAS and may assist in identifying patients who could benefit from PAS interventions.
With evidence that 4D flow CMR detects early biomarkers of ventricular dysfunction, we performed this study applying CMR ventricular flow analysis in a swine model of treated and untreated unilateral branch pulmonary artery stenosis and compared the 4D flow data to conventional anatomic and hemodynamics measures. We hypothesized that untreated PAS is associated with abnormal and inefficient RV flow patterns that are quantifiable by 4D flow CMR and that PAS interventions with intravascular stenting lead to improvement in RV mechanics.
Methods
Experimental protocol
Isolated left PAS was created in neonatal piglets (n = 14, 5. 4 kg, 2 weeks of age) by suturing a short segment of 4.0 mm Gore-Tex tube graft around the proximal left PA (LPA) [20] with 4 additional piglets having sham surgery to serve as normal controls. Four swine served as untreated left PAS controls and 10 swine underwent stent interventions to treat the left PAS. Stent interventions occurred between 5 and 10 weeks of age with a possible redilation procedure to match the distal LPA diameter by 10 weeks of age. Further technical details for this animal model of left PAS and stent interventions are detailed in our previous publication [20]. All groups (sham: n = 4, left PAS n = 4, and intervention n = 10) had right heart catheterizations (RHC), computed tomography (CT) angiography and CMR at 20 weeks of age. Given the faster rate of cardiopulmonary development in swine compared to humans [21], this timeline is developmentally equivalent to performing stent interventions in a 3–5 year old human patient and then comparing outcomes in late adolescence or early adulthood. This study complied with all institutional and national requirements for the care and use of laboratory animals and the Institutional Animal Care and Use Committee of the University of Wisconsin reviewed and approved this protocol.
Hemodynamic measurements at 20-week catheterization
RHC was performed following induction (telazol and xylazine intramuscular) and anesthesia (isoflurane 2–5%) during mechanical ventilation. Venous access was obtained percutaneously and intravenous heparin was given. Beat-by-beat pressure waveforms were obtained in the right atrium, RV, main PA, right PA (RPA), LPA, and bilateral capillary wedge positions using end-hole fluid filled catheters attached to a pressure transducer.
Computed tomography angiography
Under the same anesthesia, CT angiography was performed using a 64-slice CT scanner (750 CT, General Electric Healthcare, Waukesha, Wisconsin, USA) with retrospective cardiac gating and intravenous Iopamidol-370 (Bracco Diagnostics Inc., Monroe Township, New Jersey, USA) for contrast enhancement. The LPA diameter was measured immediately after the PA bifurcation and at two more distal locations adjacent to the LPA first order branches as previously defined [20]. (See Fig. 1 for measurement locations and representative 3D PA reconstructions).
CMR protocol
Under the same anesthesia on a 3T system (Discovery MR 750, General Electric Healthcare) using a short axis balanced steady state free precession (bSSFP) CMR sequence and the under sampled-radial 4D flow CMR sequence PC-VIPR (Phase Contrast Vastly Under sampled Projection Imaging) [22, 23]. bSSFP CMR used a 50° flip angle, TR = 3.3 ms, TE = 1.3 ms, a voxel size of 1.4 × 1.4x8.0 mm and data was reconstructed to 20 time frames per cardiac cycle. 4D flow CMR acquisition parameters include TR = 6.3 ms, TE = 2.1 ms, an isotropic voxel size of 1.25 × 1.25x1.25 mm, VENC = 150 cm/sec and data was reconstructed to 20 time frames per cardiac cycle. 4D flow CMR was used to evaluate ventricular flow.
Cine CMR analysis
Ventricular volumes were segmented and ejection fraction (EF) and cardiac output (CO) were calculated and normalized to body surface area (BSA). Interventricular septum (IVS) curvature was calculated from the IVS radius at a mid-ventricular slice normalized by the cube root of LV end-diastolic volume (EDV).
Strain was calculated from short axis CMR using feature tracking in Segment 2.2 R6423 (http://segment.heiberg.se/, [24]). Global peak circumferential systolic strain, peak circumferential systolic strain rate and peak circumferential diastolic strain rate were calculated for both the left ventricle (LV) and the RV. One stenosis swine was excluded from strain analysis because of poor image quality.
4D flow CMR analysis
Lung perfusion was quantified from 4D flow CMR using flow visualization software (EnSight, Ansys, Canonsburg, Pennsylvania, USA). LPA flow rate was calculated by subtracting the RPA flow rate from the main PA flow rate as CMR signal is lost in a metal stent and distal to the stenosis there may be aliasing dephasing. A 4D flow CMR dataset was not obtained for one intervention swine.
The LV and RV were segmented from end diastolic cine frames. The end diastolic frame was chosen to capture the ventricles through the entire cardiac cycle. The use of a single ventricle mask has previously been validated against time resolved masks with a single ventricle mask mildly overestimating peak systolic kinetic energy [25]. Velocity field gradients were calculated using a centered finite difference in MATLAB (Mathworks, Natick, Massachusetts, USA). The density and viscosity of blood were assumed to be 1.06 g/cm3 and 0.04 poise respectively.
Using the velocity data from the 4D flow CMR acquisitions, kinetic energy (KE), vorticity (\({\varvec{\omega}}\)) and energy dissipation rate (\(\varepsilon\)) were calculated from the following equations
where \(\rho\) is blood density, u is the velocity vector, \(\nabla \times\) is = the curl operator, \(\mu\) is blood viscosity and S is the strain rate tensor (\({\varvec{S}}=\frac{1}{2}\left[\nabla {\varvec{u}}+\nabla {{\varvec{u}}}^{T}\right]\)). Flow parameters were reported similar to prior studies [26, 27]: KE was reported at peak systole and peak diastole, vorticity was averaged over the cardiac cycle and energy dissipation rate was summed over the cardiac cycle to calculate the total energy dissipation. Both KE and energy dissipation were normalized by stroke volume.
Ventricular flow parameters are typically either not normalized or normalized by stroke volume although heart rate and ventricular volume could also be different be subjects. To account for these differences, dimensionless KE, vorticity magnitude and energy dissipation rate were calculated from dimensional analysis with characteristic length and velocity scales calculated from CO and EDV.
where * refers to the dimensionless variable. Dimensional analysis is the proper fluid dynamics method to account for this variation in geometry and flow conditions.
Lastly, to understand the mechanism of ventricular energy dissipation we calculated the percentage of energy dissipation due to vorticity (\(100\%* {{\varvec{\omega}}}^{2}/2{{\varvec{S}}}^{2}\)) over the entire cardiac cycle to determine if energy dissipation was primarily due to vorticity or due to surface strain rate. As this is a novel biomarker for ventricular flow we additionally performed a “sanity check” by comparing 4D flow CMR measurements to the vorticity-energy dissipation rate relationship for an idealized-theoretical vortex ring. Theoretical vortex ring results were presented for a Reynold’s number of 900 which is both reasonable for ventricular flow [28] and the theoretical formulation has been experimentally validated for this Reynold’s number [29].
Statistics
All results were reported as mean ± standard deviation (SD). Statistical analysis was performed SPSS (Statistical Package for the Social Sciences, International Business Machines, Inc., Armonk, New York, USA). Data normality was assessed with a Shapiro-Wilkes test. In normally distributed data, statistical comparisons were made by ANOVA with Tukey’s honestly significant difference procedure for post hoc comparisons. When data were not normally distributed, comparisons between groups were made using the non-parametric Kruskal–Wallis test with a Bonferroni correction for post-hoc comparisons. The p-values for every comparison performed can be found in Appendix 1.
Results
Cardiac catheterization
Conventional hemodynamic measurements are summarized in Table 1. At 20 weeks, there were no differences in body weight, heart rate or RV/(LV + Septum) ratio. LPA banding created severe angiographic left PAS associated with a peak systolic main PA to distal LPA pressure gradient of 23 + 4 mmHg. The left PAS group had elevated mean right atrial, RV systolic, mean main PA pressure, main PA and right PA systolic pressures. The left PAS group also had decreased LPA systolic pressure.
Standard CMR and CT angiography
Results from standard CMR and CT angiography are shown in Table 2. Stent interventions significantly increased the size of the LPA at all measurement locations compared to the left PAS group (p < 0.01). Stent interventions did not increase the size of the proximal LPA (stented segment) or the first LPA measurement location (adjacent to the stent) to that of sham controls (p = 0.03 and p = 0.052). LPA diameters at the second, more distal, measurement location were similar between sham and stent interventions. Similarly, stenting increased left lung perfusion for the intervention group (42 ± 6%) compared to the left PAS group sixfold (7 ± 5%, p < 0.001) but did not reach perfusion levels of the sham control group (52 ± 9%, p = 0.07).
RV and LV stroke volume (SV) were not different between groups. The intervention group had smaller end systolic volumes (ESV) and EDV in both the RV and the LV compared to left PAS. The intervention group also had elevated RVEF versus left PAS (p = 0.03), trended towards increased RVEF versus sham (p = 0.06) and trended towards increased LVEF (p = 0.06) and increased cardiac index (CI) compared to the sham control group (p = 0.08).
Ventricular strain analysis
Ventricular strain measurements did not identify any statistically significant differences between groups (Table 3).
Dimensional ventricular flow analysis
KE, vorticity and energy dissipation indices are shown in Fig. 2. Stent interventions normalized all dimensional measures of ventricular mechanics to levels equal to sham controls. In the RV, the left PAS group trended towards increased systolic and diastolic KE, vorticity and energy dissipation although only systolic KE approached significance (p = 0.07 vs sham, p = 0.06 vs intervention). In the LV the left PAS group had increased systolic and diastolic KE, vorticity and energy dissipation.
Dimensionless ventricular flow analysis
Time curves of dimensionless KE, vorticity and energy dissipation are shown along with bar graphs of systolic and diastolic peak KE, average vorticity and average energy dissipation rate in Fig. 3. The time curves demonstrate that throughout the cardiac cycle, the left PAS group has significant abnormalities in all measures of ventricular mechanics in both ventricles and stent interventions normalized all levels to that of sham controls. The analysis of dimensionless flow parameters again confirms the benefits of a PA stent intervention for improving the efficiency of ventricular flow. Elevated RV systolic KE, diastolic KE, and elevated RV vorticity in the left PAS group were not statistically significant when analyzing dimensional parameters, but were statistically significant when comparing groups using dimensionless parameters.
Mechanism of energy dissipation
For an idealized vortex ring [29], vorticity and energy dissipation rate are quadratically related with 25% of energy dissipation being directly due to vorticity (Fig. 4a). For all groups and both ventricles, 4D flow CMR measurements also showed a strong quadratic relationship (adjusted R2 = 0.86) between vorticity and energy dissipation rate (Fig. 4b). This was theoretically expected and the percentage of energy dissipation due to vorticity (21 ± 2%) agreed well with the idealized vortex ring. In the RV no differences were found for the energy dissipation mechanism between groups (Fig. 4c). In the LV the percentage of energy dissipation due to vorticity was less in the stenosis group (p = 0.02 vs sham, p = 0.04 vs intervention) indicating increased energy dissipation from surface deformation (Fig. 4c). The mechanism of LV energy dissipation in the intervention group was not different than the sham control group.
Discussion
The results of this study reveal previously unknown effects of unilateral PAS on ventricular mechanics. In this animal model, 4D flow CMR detected abnormal and inefficient flow in both the RV and the LV associated with untreated PAS. Stent interventions to relieve PAS normalized both RV and LV KE, vorticity and energy dissipation. These results are clinically relevant as they demonstrate that in addition to decreasing RV afterload and improving PA blood flow distributions [10] successful PA stent interventions improve bi-ventricular flow dynamics that should contribute to the long-term myocardial health.
While alterations in cardiac flow patterns and parameters seem nearly ubiquitous in cardiovascular disease, few studies have shown the capability of CHD interventions to correct changes in intra cardiac flows. In a swine CHD model of pulmonary valve insufficiency and unilateral PAS Petit et al., demonstrated that relief of the PAS with stent interventions reduced pulmonary regurgitation and improved RV systolic function [30]. This frequently cited animal study is used to justify aggressive interventional management of PAS in patients with pulmonary valve insufficiency in hopes of slowing RV dilation and preserving ventricular function. Through application of more advanced 4D flow CMR techniques Sjoberg et al. [31] demonstrated a decrease in RV systolic and diastolic KE in adult Tetralogy of Fallot (TOF) patients following surgical pulmonary valve replacement that corresponds with a decrease in RV volumes (RV EDV, 164 ± 18–116 ± 18 ml/m2). This was the first CHD study to quantify KE during the entire cardiac cycle and following surgical pulmonary valve replacement KE was comparable to controls in both ventricles. While our left PAS model should have milder changes in RV mechanics than tetralogy of Fallot (TOF) patients with severe pulmonary regurgitation, the results of the current study similarly demonstrate that PA stent interventions normalized the ventricular flow measures of KE, vorticity and energy dissipation in the ventricles to levels similar to sham controls. This normalization in ventricular flow mechanics in combination with normal PA pressures documented at catheterization should reduce the energetic demands of the myocardium and we hypothesize this contributes to clinical improvements seen following successful PAS interventions [8,9,10]. Our previous swine PAS stenting study demonstrated that in response to dobutamine stress testing, the stent intervention animals increased cardiac output similar to sham controls while under the same stress conditions CO actually decreased in untreated PAS [20]. The proposed mechanisms responsible for clinical improvement in exercise capacity after relief of branch PAS include improved ventilation-perfusion matching, reduced RV afterload and improved stroke volume response to meet increased metabolic demands [8,9,10]. The current study would suggest that preserved ventricular flow mechanics play an additional role in clinical improvement following successful PAS interventions. Further study is warranted, including the effects of pulmonary regurgitation.
Studies of heart failure patients have found that exercise capacity is related to LV direct flow KE [32] and LV KE fluctuations [33]. In addition to exercise capacity, LV KE in heart failure patients is moderately associated with heart failure symptoms, serum markers of ventricular remodeling, and myocardial energetics [32]. While PA stent interventions improved ventricular flow dynamics, it is interesting that the intervention group remodeled to have higher EF compared to the sham and PAS controls. It appears that ventricular remodeling in the intervention group occurred to maintain normal ventricular flow and mechanical homeostasis for the myocardium. The effects of this remodeling seem beneficial (higher EF and improved flow efficiency) and consistent with improved exercise capacity in humans following relief of PAS.
4D flow CMR showed that the PAS group had abnormal and inefficient flow in the RV (Figs. 2 and 3) despite having similar heart rate, CI and EF values compared to sham controls. This finding was similar to 4D flow CMR results for TOF patients, who had RV KE inefficiency [34], elevated RV diastolic KE [31] and increased RV vortexes and vorticity [27, 35] compared to healthy subjects. Similarly, patients who have surgically repaired transposition of the great arteries have increased pulmonary artery vortices [36]. The changes in RV flow parameters are disease dependent. In swine pulmonary regurgitation increases RV KE and PAS decreases RV KE [37]. It is important to note that while Fernandes et al. found decreased RV KE in the pressure overloaded swine RV, the opposite of our finding, their PAS was surgically created at the level of the main PA and the resultant RV pressure overload was much greater [37].
We were surprised to find that LV flow was also altered in the PAS group when only the RV was directly affected by the stenosis and the degree of pressure elevation was mild. We would speculate These results seem similar to rTOF patients who despite having normal LV function exhibit decreased LV KE [31]. It has previously been shown that in pediatric PA hypertension patients, increased RV pressure decreases septal curvature which decreases LV torsion and ultimately LV function [38, 39]. We speculate that our results showing abnormal LV flow in the setting of mildly elevated RV afterload are early signs of the impact of septal mechanics on LV function. Again similar to our results following PAS stenting, a right side intervention, improved LV flow, pulmonary valve replacement in TOF patients improves LV KE compared to healthy subjects [31]. The detection of altered RV and LV flow parameters with only mild RV afterload and no changes in traditional functional strain parameters underscores the sensitivity of 4D flow CMR derived biomarkers to detect cardiac dysfunction.
Additional technical discussion of two points regarding ventricular flow analysis that are pertinent to this study (1. the mechanism of energy dissipation, 2. dimensional analysis) are included in Appendix 2.
Limitations
Limitations of this work include the use of an animal model of CHD in place of human subjects as myocardial and PA response to catheter interventions may be species specific. The short time course may limit the degree of acquired disease. In addition, only a small number of subjects were studied per group and there was likely not sufficient power to detect differences between groups for metrics that had greater variability (e.g. RV energy dissipation and RV strain). Myocardial strain analysis only quantified circumferential strain and future studies should consider longitudinal strain which is a more reproducible measure. All measurements occurred under general anesthesia which is known to decrease myocardial contractility in swine. We cannot rule out the possibility that greater RVEF in the intervention group could be a result of differential response to anesthesia rather than a true difference in myocardial function. This study also used only male swine so we were unable to assess for any potentially important sex differences [40].
Conclusion
In a swine model of PAS, stent interventions normalized RV and LV flow mechanics as measured by 4D flow CMR. Additionally, we identified significantly inefficient RV and LV flow associated with unilateral branch PAS, even in the setting of normal heart rate, EF and EDV index. The stent intervention group may have achieved normal ventricular flow by increased EF. These results provide evidence that PAS impacts ventricular function in CHD in subtle, yet likely important ways. Successful PA stent interventions improve ventricular flow efficiency and may promote long-term health of the ventricle. This study also highlights the sensitivity of 4D flow CMR biomarkers to detect ventricular dysfunction and treatment response in a pediatric setting.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- BSA:
-
Body surface area
- bSSFP:
-
Balanced steady state free precession
- CHD:
-
Congenital heart disease
- CI:
-
Cardiac index
- CMR:
-
Cardiovascular magnetic resonance
- CO:
-
Cardiac output
- CT:
-
Computed tomography
- EDV:
-
End-diastolic volume
- EF:
-
Ejection fraction
- ESV:
-
End-systolic volume
- IVS:
-
Interventricular septum
- KE:
-
Kinetic energy
- LPA:
-
Left pulmonary artery
- LV:
-
Left ventricle/left ventricular
- LVEF:
-
Left ventricular ejection fraction
- PA:
-
Pulmonary artery
- PAS:
-
Pulmonary artery stenosis
- PCWP:
-
Pulmonary capillary wedge pressure
- RHC:
-
Right heart catheterization
- RPA:
-
Right pulmonary artery
- RV:
-
Right ventricle/right ventricular
- RVEF:
-
Right ventricular ejection fraction
- SV:
-
Stroke volume
- TOF:
-
Tetralogy of Fallot
References
Hirsch JC, Mosca RS, Bove EL. Complete repair of tetralogy of fallot in the neonate results in the modern era. Ann Surg. 2000;232(4):508–14.
Parry AJ, Mcelhinney DB, Kung GC, Reddy VM, Brook MM, Hanley FL. Elective primary repair of acyanotic tetralogy of fallot in early infancy: overall outcome and impact on the pulmonary valve. J Am Coll Cardiol. 2000;36(7):2279–83.
Groh MA, et al. Repair of tetralogy of fallot in infancy. Effect of pulmonary artery size on outcome. Circulation. 1991;84(5 Suppl):III206–12.
Kim H, Sung SC, Chang YH, Lee HD, Park JA. Early and midterm outcomes of left pulmonary artery angioplasty using an anterior wall flap of the main pulmonary artery in tetralogy of fallot repair. J Thorac Cardiovasc Surg C. 2014;148(6):2597–601.
Wilder TJ, et al. Aggressive patch augmentation may reduce growth potential of hypoplastic branch pulmonary arteries after tetralogy of fallot repair. Ann Thorac Surg. 2016;101(3):996–1004.
Gatzoulis MA, et al. Risk factors for arrhythmia and sudden cardiac death late after repair of tetralogy of fallot: a multicentre study. Lancet. 2000;356(9234):975–81.
Harrison DA, et al. Sustained ventricular tachycardia in adult patients late after repair of tetralogy of fallot. J Am Coll Cardiol. 1997;30(5):1368–73.
Rhodes J, et al. Effect of pulmonary artery stenoses on the cardiopulmonary response to exercise following repair of tetralogy of fallot. Am J Cardiol. 1998;81(10):1217–9.
Hiremath G, et al. Balloon angioplasty and stenting for unilateral branch pulmonary artery stenosis improve exertional performance. JACC Cardiovasc Interv. 2019;12(3):289–97.
Sierra-Galan LM, François CJ. Clinical applications of MRA 4D-flow. Curr Treat Options Cardiovasc Med. 2019;21(10):1–18.
Zhong L, et al. Intracardiac 4D flow MRI in congenital heart disease: recommendations on behalf of the ISMRM flow & motion study group. J Magn Reson Imaging. 2019;50(3):jmri.26858.
Vasanawala SS, Hanneman K, Alley MT, Hsiao A. Congenital heart disease assessment with 4D flow MRI. J Magn Reson Imaging. 2015;42(4):870–86.
Isorni M-A, Tortigue M, Ben Moussa N, Hascoët S, Monnot S. 4D flow CMR analysis in repaired tetralogy of fallot: where we are and where we are going. Arch Cardiovasc Dis Suppl. 2018;10(3–4):285.
Isorni MA, et al. 4D flow MRI versus conventional 2D for measuring pulmonary flow after Tetralogy of Fallot repair. Int J Cardiol. 2020;300:132–6.
Carlsson M, Heiberg E, Toger J, Arheden H. Quantification of left and right ventricular kinetic energy using four-dimensional intracardiac magnetic resonance imaging flow measurements. Am J Physiol Hear Circ Physiol. 2012;302:893–900.
Töger J, et al. Vortex ring formation in the left ventricle of the heart: analysis by 4D flow MRI and lagrangian coherent structures. Ann Biomed Eng. 2012;40(12):2652–62.
Pedrizzetti G, La Canna G, Alfieri O, Tonti G. The vortex–an early predictor of cardiovascular outcome? Nat Rev Cardiol. 2014;11(9):545–53.
Pasipoularides A. Diastolic filling vortex forces and cardiac adaptations: Probing the epigenetic nexus. Hellenic J Cardiol. 2012;53(6):458–69.
Hove JR, Köster RW, Forouhar AS, Acevedo-Bolton G, Fraser SE, Gharib M. Intracardiac fluid forces are an essential epigenetic factor for embryonic cardiogenesis. Nature. 2003;421(6919):172–7.
Bates ML, et al. Consequences of an early catheter-based intervention on pulmonary artery growth and right ventricular myocardial function in a pig model of pulmonary artery stenosis. Catheter Cardiovasc Interv. 2018;92(1):78–87.
Rendas A, Branthwaite M, Reid L. Growth of pulmonary circulation in normal pig-structural analysis and cardiopulmonary function. J Appl Physiol Respir Environ Exerc Physiol. 1978;45(5):806–17.
Wieben O, Francois C, Reeder SB. Cardiac MRI of ischemic heart disease at 3 T: Potential and challenges. Eur J Radiol. 2008;65(1):15–28.
François CJ, et al. 4D cardiovascular magnetic resonance velocity mapping of alterations of right heart flow patterns and main pulmonary artery hemodynamics in tetralogy of Fallot. J Cardiovasc Magn Reson. 2012;14:16.
Heiberg E, Sjögren J, Ugander M, Carlsson M, Engblom H, Arheden H. Design and validation of segment—freely available software for cardiovascular image analysis. BMC Med Imaging. 2010;10(1):1–13.
Hussaini SF, Rutkowski DR, Roldán-Alzate A, François CJ. Left and right ventricular kinetic energy using time-resolved versus time-average ventricular volumes. J Magn Reson Imaging. 2017;45(3):821–8.
Rutkowski DR, Barton G, François CJ, Bartlett HL, Anagnostopoulos PV, Roldán-Alzate A. Analysis of cavopulmonary and cardiac flow characteristics in fontan patients: comparison with healthy volunteers. J Magn Reson Imaging. 2019;49(6):1786–99.
Hirtler D, Garcia J, Barker AJ, Geiger J. Assessment of intracardiac flow and vorticity in the right heart of patients after repair of tetralogy of Fallot by flow-sensitive 4D MRI. Eur Radiol. 2016;26(10):3598–607.
Vasudevan V, et al. Flow dynamics and energy efficiency of flow in the left ventricle during myocardial infarction. Biomech Model Mechanobiol. 2017;16(5):1503–17.
Kaplanskii FB, Rudi YA. Evolution of a viscous vortex ring. Fluid Dyn. 2001;36(1):20–30.
Petit CJ, et al. Relief of branch pulmonary artery stenosis reduces pulmonary valve insufficiency in a swine model. J Thorac Cardiovasc Surg. 2009;138(2):382–9.
Sjöberg P, et al. Disturbed left and right ventricular kinetic energy in patients with repaired tetralogy of fallot: pathophysiological insights using 4D-flow MRI. Eur Radiol. 2018;28(10):4066–76.
Stoll VM, et al. Left ventricular flow analysis. Circ Cardiovasc Imaging. 2019;12(5):e008130.
Kim IC, Hong GR, Pedrizzetti G, Shim CY, Kang SM, Chung N. Usefulness of left ventricular vortex flow analysis for predicting clinical outcomes in patients with chronic heart failure: a quantitative vorticity imaging study using contrast echocardiography. Ultrasound Med Biol. 2018;44(9):1951–9.
Jeong D, 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(5):1339–47.
Francois CJ, Srinivasan S, Landgraf BR, Niespodzany E, Wieben O, Frydrychowicz A. Analysis of right heart flow patterns in repaired Tetralogy of Fallot with 4D flow-sensitive MRI. J Cardiovasc Magn Reson. 2011;13(S1):1–2.
Geiger J, et al. Postoperative pulmonary and aortic 3D haemodynamics in patients after repair of transposition of the great arteries. Eur Radiol. 2014;24(1):200–8.
Fernandes JF, et al. Right ventricular energetics and power in pulmonary regurgitation vs. stenosis using four dimensional phase contrast magnetic resonance. Int J Cardiol. 2018;263:165–70.
Dufva MJ, Truong U, Tiwari P, Ivy DD, Shandas R, Kheyfets VO. Left ventricular torsion rate and the relation to right ventricular function in pediatric pulmonary arterial hypertension. Pulm Circ. 2018;8(3):204589401879135.
Kheyfets VO, Truong U, Ivy D, Shandas R. Structural and biomechanical adaptations of right ventricular remodeling-in pulmonary arterial hypertension-reduces left ventricular rotation during contraction: a computational study. J Biomech Eng. 2019. https://doi.org/10.1115/1.4042682.
Rutkowski DR, Barton GP, François CJ, Aggarwal N, Roldán-Alzate A. Sex differences in cardiac flow dynamics of healthy volunteers. Radiol Cardiothorac Imaging. 2020;2(1):e190058.
Wu J-Z, Zhou Y, Fan M. A note on kinetic energy, dissipation and enstrophy. Phys Fluids. 1999;11:503.
Acknowledgements
Not applicable.
Funding
This investigation was supported by the Clinical and Translational Science Award (CTSA) program, through the NIH National Center for Advancing Translational Sciences (NCATS), grant UL1TR002373 (AR, LL and CF), the NIH Ruth L. Kirschstein National Research Service Award T32 HL 007936 from the National Heart Lung and Blood Institute to the University of Wisconsin-Madison Cardiovascular Research Center (RP). This investigation was partially supported by T32AI007635 (GB) and the Summer Program in Undergraduate Urologic Research (CJ). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
Author information
Authors and Affiliations
Contributions
Processed medical imaging data (RP, GB, CJ, LL and AR), performed catheterization and stenting studies (LL), planned study and obtained funding (CF, LL and AR), wrote initial manuscript draft (RP), and made revisions to manuscript (all authors). All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study complied with all institutional and national requirements for the care and use of laboratory animals and the Institutional Animal Care and Use Committee of the University of Wisconsin reviewed and approved this protocol.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendices
Appendix 1
See Table 4
Appendix 2
In Appendix 2 we further expand our discussion to include two important technical points about ventricular flow analysis: 1. the mechanism of energy dissipation and 2. dimensional analysis. First, we calculated the percentage of energy dissipation directly due to vorticity as a novel metric of ventricular flow-structure interactions. Statistical correlations between flow parameters and ventricular kinematics [40] are useful but may only be associations and lack a direct physical basis. We propose a new metric by recognizing that viscous energy dissipation is due to two components, vorticity and surface strain rate [41]. The surface-strain rate component of energy dissipation is zero for flows with stationary walls and for these flows energy dissipation is solely due to vorticity. The percentage of energy dissipation directly due to vorticity is then a metric that will decrease if structural motion is excessively causing energy losses. We found percentage of energy dissipation from vorticity did not change in the right ventricle (RV) but in the left ventricle (LV) the pulmonary artery stenosis (PAS) group had a lower percentage of energy dissipation from vorticity. These results suggest that in PAS, the LV wall and mitral valve motion are excessively contributing to energy dissipation compared to the sham and intervention groups. The percentage of energy dissipation from vorticity could be a biomarker of the impact of ventricular wall and atrioventricular valve kinematics on ventricular flow efficiency in other diseases as well.
Second, we used dimensional analysis to calculate dimensionless ventricular flow parameters. It is important to normalize cardiovascular measurements given anatomical and functional differences between subjects, however 4D flow CMR derived metrics such as kinetic energy and vorticity are not systematically normalized. As ventricular flow analysis sits at a unique intersection between physiology and fluid dynamics, careful consideration must be given to the proper normalization of parameters. In this study calculating dimensionless parameters to account for subject variability in stroke volume, HR and ventricular volume improved the sensitivity of ventricular flow analysis to detect statistical differences in flow patterns between groups. Heart rate, stroke volume and ventricular volume can change in a variety of cardiovascular diseases and calculating dimensionless parameters may aid ventricular flow analysis in other diseases as well.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Pewowaruk, R.J., Barton, G.P., Johnson, C. et al. Stent interventions for pulmonary artery stenosis improve bi-ventricular flow efficiency in a swine model. J Cardiovasc Magn Reson 23, 13 (2021). https://doi.org/10.1186/s12968-021-00709-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12968-021-00709-4