Skip to main content

Stent interventions for pulmonary artery stenosis improve bi-ventricular flow efficiency in a swine model



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.


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.


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).


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.


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.


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).

Fig. 1

Representative 3D pulmonary artery reconstructions for sham, pulmonary artery (PA) stenosis (PAS) and stent interventions. The left pulmonary artery (LPA) diameter measurement locations are shown on the sham PAs

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 (, [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

$$KE = \frac{1}{2}\rho {\varvec{u}}^{2}$$
$${\varvec{\omega}} = \nabla \times {\varvec{u}}$$
$$\varepsilon = 2\mu \left( {\varvec{S}} \right)^{2}$$

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.

$$KE^{*} = \frac{KE}{{\rho CO^{2} EDV^{{ - \frac{4}{3}}} }} = \frac{KE}{{\rho \left( {HR*SV} \right)^{2} EDV^{{ - \frac{4}{3}}} }}$$
$${\varvec{\omega}}^{*} = \frac{{\varvec{\omega}}}{{CO*EDV^{ - 1} }} = \frac{{\varvec{\omega}}}{{HR*SV*EDV^{ - 1} }}$$
$$\varepsilon^{*} = \frac{\varepsilon }{{\mu CO^{2} EDV^{ - 2} }} = \frac{\varepsilon }{{\mu \left( {HR*SV} \right)^{2} EDV^{ - 2} }}$$

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].


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.


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.

Table 1 Hemodynamic measurements

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).

Table 2 Standard CMR and CT angiography

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).

Table 3 Ventricular strain analysis

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.

Fig. 2

Dimensional RV and LV flow parameters. a peak systolic and peak diastolic kinetic energy (KE), b vorticity and c energy dissipation. Both KE and energy dissipation were indexed by stroke volume. *p < 0.05

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.

Fig. 3

Time curves and bar graphs of non-dimensional flow parameters are shown for both the RV and the LV. Standard error is not shown on the time curve plots. Row 1: kinetic energy, Row 2: vorticity, Row 3: energy dissipation rate. Sys – systolic, Dia – diastolic. *p < 0.05

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.

Fig. 4

As expected, quadratic relationships are found between vorticity and energy dissipation rate for both (a). the idealized-theoretical vortex ring and (b). 4D flow CMR measurements of ventricular flow. c Percent energy dissipation from vorticity is compared between groups. *p < 0.05


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 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].


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.



Body surface area


Balanced steady state free precession


Congenital heart disease


Cardiac index


Cardiovascular magnetic resonance


Cardiac output


Computed tomography


End-diastolic volume


Ejection fraction


End-systolic volume


Interventricular septum


Kinetic energy


Left pulmonary artery


Left ventricle/left ventricular


Left ventricular ejection fraction


Pulmonary artery


Pulmonary artery stenosis


Pulmonary capillary wedge pressure


Right heart catheterization


Right pulmonary artery


Right ventricle/right ventricular


Right ventricular ejection fraction


Stroke volume


Tetralogy of Fallot


  1. 1.

    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.

    CAS  Article  Google Scholar 

  2. 2.

    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.

    CAS  Article  Google Scholar 

  3. 3.

    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.

    CAS  PubMed  Google Scholar 

  4. 4.

    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.

    Article  Google Scholar 

  5. 5.

    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.

    Article  Google Scholar 

  6. 6.

    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.

    CAS  Article  Google Scholar 

  7. 7.

    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.

    CAS  Article  Google Scholar 

  8. 8.

    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.

    CAS  Article  Google Scholar 

  9. 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.

    Article  Google Scholar 

  10. 10.

    Sierra-Galan LM, François CJ. Clinical applications of MRA 4D-flow. Curr Treat Options Cardiovasc Med. 2019;21(10):1–18.

    Article  Google Scholar 

  11. 11.

    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.

    Article  Google Scholar 

  12. 12.

    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.

    Article  Google Scholar 

  13. 13.

    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.

    Google Scholar 

  14. 14.

    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.

    CAS  Article  Google Scholar 

  15. 15.

    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.

    Article  Google Scholar 

  16. 16.

    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.

    Article  Google Scholar 

  17. 17.

    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.

    Article  Google Scholar 

  18. 18.

    Pasipoularides A. Diastolic filling vortex forces and cardiac adaptations: Probing the epigenetic nexus. Hellenic J Cardiol. 2012;53(6):458–69.

    PubMed  PubMed Central  Google Scholar 

  19. 19.

    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.

    CAS  Article  Google Scholar 

  20. 20.

    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.

    Article  Google Scholar 

  21. 21.

    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.

    CAS  PubMed  Google Scholar 

  22. 22.

    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.

    Article  Google Scholar 

  23. 23.

    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.

    Article  Google Scholar 

  24. 24.

    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.

    Article  Google Scholar 

  25. 25.

    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.

    Article  Google Scholar 

  26. 26.

    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.

    Article  Google Scholar 

  27. 27.

    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.

    Article  Google Scholar 

  28. 28.

    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.

    Article  Google Scholar 

  29. 29.

    Kaplanskii FB, Rudi YA. Evolution of a viscous vortex ring. Fluid Dyn. 2001;36(1):20–30.

    Article  Google Scholar 

  30. 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.

    Article  Google Scholar 

  31. 31.

    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.

    Article  Google Scholar 

  32. 32.

    Stoll VM, et al. Left ventricular flow analysis. Circ Cardiovasc Imaging. 2019;12(5):e008130.

    Article  Google Scholar 

  33. 33.

    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.

    Article  Google Scholar 

  34. 34.

    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.

    Article  Google Scholar 

  35. 35.

    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.

    Article  Google Scholar 

  36. 36.

    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.

    Article  Google Scholar 

  37. 37.

    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.

    Article  Google Scholar 

  38. 38.

    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.

    Article  Google Scholar 

  39. 39.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    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.

    Article  Google Scholar 

  41. 41.

    Wu J-Z, Zhou Y, Fan M. A note on kinetic energy, dissipation and enstrophy. Phys Fluids. 1999;11:503.

    CAS  Article  Google Scholar 

Download references


Not applicable.


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




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

Correspondence to Ryan J. Pewowaruk.

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.



Appendix 1

See Table 4

Table 4 P-values for all statistical comparisons

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

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).

Download citation


  • 4D flow MRI
  • Congenital heart disease
  • Pediatrics
  • Right ventricle