Left ventricular blood flow kinetic energy after myocardial infarction - insights from 4D flow cardiovascular magnetic resonance

Background Myocardial infarction (MI) leads to complex changes in left ventricular (LV) haemodynamics that are linked to clinical outcomes. We hypothesize that LV blood flow kinetic energy (KE) is altered in MI and is associated with LV function and infarct characteristics. This study aimed to investigate the intra-cavity LV blood flow KE in controls and MI patients, using cardiovascular magnetic resonance (CMR) four-dimensional (4D) flow assessment. Methods Forty-eight patients with MI (acute-22; chronic-26) and 20 age/gender-matched healthy controls underwent CMR which included cines and whole-heart 4D flow. Patients also received late gadolinium enhancement imaging for infarct assessment. LV blood flow KE parameters were indexed to LV end-diastolic volume and include: averaged LV, minimal, systolic, diastolic, peak E-wave and peak A-wave KEiEDV. In addition, we investigated the in-plane proportion of LV KE (%) and the time difference (TD) to peak E-wave KE propagation from base to mid-ventricle was computed. Association of LV blood flow KE parameters to LV function and infarct size were investigated in all groups. Results LV KEiEDV was higher in controls than in MI patients (8.5 ± 3 μJ/ml versus 6.5 ± 3 μJ/ml, P = 0.02). Additionally, systolic, minimal and diastolic peak E-wave KEiEDV were lower in MI (P < 0.05). In logistic-regression analysis, systolic KEiEDV (Beta = − 0.24, P < 0.01) demonstrated the strongest association with the presence of MI. In multiple-regression analysis, infarct size was most strongly associated with in-plane KE (r = 0.5, Beta = 1.1, P < 0.01). In patients with preserved LV ejection fraction (EF), minimal and in-plane KEiEDV were reduced (P < 0.05) and time difference to peak E-wave KE propagation during diastole increased (P < 0.05) when compared to controls with normal EF. Conclusions Reduction in LV systolic function results in reduction in systolic flow KEiEDV. Infarct size is independently associated with the proportion of in-plane LV KE. Degree of LV impairment is associated with TD of peak E-wave KE. In patient with preserved EF post MI, LV blood flow KE mapping demonstrated significant changes in the in-plane KE, the minimal KEiEDV and the TD. These three blood flow KE parameters may offer novel methods to identify and describe this patient population. Electronic supplementary material The online version of this article (10.1186/s12968-018-0483-6) contains supplementary material, which is available to authorized users.


Background
Myocardial infarction (MI) is one of the leading causes of death and disability worldwide. The effects of acute MI on left ventricular (LV) haemodynamics [1] and the prognostic importance of infarct size and LV ejection fraction (EF) following MI are well established [2]. Large infarcts are associated with an increased LV wall stress secondary to increased pre-load on the LV [3,4]. This and subsequent adverse effects on LV afterload play a crucial role in infarct expansion and dilatation leading to adverse LV remodelling [5]. The assessment of LV haemodynamics and flow provides incremental prognostic evaluation in patients with MI [6]. Two-dimensional (2D) Doppler echocardiography is the mainstay for the non-invasive assessment of LV haemodynamics. However, Doppler echocardiography provides unidirectional velocity information and does not inform about three-dimensional intra-cavity flow [7].
The kinetic energy (KE) of the blood represents a fundamental component of work performed by the heart which results in the movement of the blood [8]. Four-dimensional flow (4D flow) cardiovascular magnetic resonance (CMR) imaging allows semi-automatic quantification of intra-cavity LV flow KE parameters in three dimensions (3D) [9][10][11][12][13][14][15]. This further allows to compute several parameters of KE at different time points in the cardiac cycle (Fig. 1). Previous studies have demonstrated LV systolic average KE is higher in patients with heart failure compared to healthy controls, but lower when indexing to heart size [16]. Another study has shown that the KE of the portion of flow that passes directly through the LV in a single cardiac cycle or 'direct flow' diminishes with progressive LV dilatation [17]. Recently, Wong et al. demonstrated that in healthy individuals, early peak diastolic KE progressively decreases with age, whereas systolic peaks remain constant [18]. The effects of MI on 3D LV blood flow KE have been described by Kanski et al. in a study which recruited MI patients as a sub-group [16]. However, more focused validation studies with larger MI patient population are needed to investigate if there is association of mechanical function, infarct characteristics and LV blood flow KE. It also remains unclear if significant differences in LV KE are seen in acute versus chronic MI.
Hence, the aims of this study were as follows: a) quantify dynamic parameters of LV blood flow KE in patients with acute and chronic MI and age and gender matched healthy controls; b) investigate the reproducibility of LV blood flow KE parameters; c) investigate if there are significant changes in LV blood flow KE in different sub-groups of left-ventricular ejection fraction; d) investigate the association of LV KE to infarct size and if LV KE is altered in acute versus chronic MI patients.

Study population
The study design is detailed in Fig. 2. Patients with acute MI patients were prospectively identified following an admission with ST-elevation myocardial infarction (STEMI) at Leeds. Patients with chronic MI were prospectively identified from the cardiology outpatient clinics in Leeds. Age and gender matched healthy controls were prospectively recruited from two centres: Leeds, UK and Leiden, The Netherlands.
The inclusion criteria for patients with acute MI were: patients with first time acute ST elevation MI (STEMI) treated by primary percutaneous coronary intervention (PPCI) within 12-h of onset of chest pain. Acute STEMI was defined as per the current international guidelines [19]. Acute STEMI patients were scheduled for CMR imaging within 72-h of indexed presentation. The inclusion criteria for patients with chronic MI were: previous history of MI and the presence of ischaemic scar on late gadolinium enhancement (LGE) imaging. Exclusion criteria for all patients included: non-ischaemic cardiomyopathy, severe Fig. 1 Illustration of all the left ventricular blood flow kinetic energy (KE) parameters investigated in this study in a healthy control valvular heart disease, estimated glomerular filtration rate < 30 ml/min/1.73m 2 , haemodynamic instability (Killip class III/IV requiring on-going intravenous diuretic therapy [20]) and contraindications to CMR imaging.
Healthy control subjects had no history or symptoms of cardiovascular disease, were not on cardiovascular or other relevant medication and had no contraindications to CMR.

CMR examination
All healthy control subjects and patients underwent CMR imaging on identical 1.5 T systems at both sites (Ingenia, Philips Healthcare, Best, The Netherlands); with a 28-channel flexible torso coil and digitization of the CMR signal in the receiver coil.

CMR protocol and image acquisition
The CMR protocol was as follows: 1. Survey images 2. The following cines were defined using survey: vertical long-axis, horizontal long-axis, 3-chamber (LVoutflow tract-views), and the LV volume contiguous short axis stack. All cines were acquired with a balanced steady-state free precession (bSSFP), single-slice breath-hold sequence.  [21].

4D flow error corrections and quality checks
The effects of concomitant gradient terms were compensated using Maxwell correction methods by the CMR scanner. Remaining background errors were corrected by the local phase correction (LPC) filter on the CMR scanner performed in two-dimensional way-slice by slice. The LPC is a magnitude-weighted spatial low pass filter; pixels that are expected to be part of the static background are used with a higher weight than noisy background pixels or pixels that are expected to contain flow to determine the local phase offset. LPC uses surrounding tissue to determine "static" areas [24,25]. All 3D phase contrast data sets were investigated for phase aliasing artefacts. If present then phase unwrapping was performed as per previously published guidelines on phase-contrast methods [26]. Additionally, any spatial misalignment of 4D flow data to cine imaging was corrected before any flow analysis was performed. This was done by visualizing streamlines in 4-chamber view at peak systole and repositioning them over descending aorta. Similar checks were done during diastole in 4-chamber and 2-chamber views for peak mitral inflow streamlines.

Image analysis
All images were analysed by PG (3-years-experience in advanced CMR techniques, PG did the blinded tissue characterisation for infarct location and size for all MI cases), SC (1-year-experience in advanced CMR techniques, SC did the blinded volumetric assessment) and RVDG (> 5 years experience in advanced CMR techniques, RVDG did the blinded LV flow KE mapping). Images were evaluated offline using research software (MASS; Version 2016EXP, Leiden University Medical Center, Leiden, The Netherlands). LV volumes, stroke volume, cardiac index, EF and infarct location on LGE imaging were determined according to standard methods. Acute and chronic infarct size were determined by the full-width at half-maximum semi-automated technique [27].

Kinetic energy mapping
For calculation of LV blood flow KE parameters, the LV volumetric mesh was resliced into short-axis sections of 2 mm thickness and pixel spacing equal to the original reconstructed pixel size of the short-axis cine acquisition (1.0-1.2 mm). This high-resolution LV mesh is construc ted by representing the mesh in cylinder coordinates. The LV radius for a given angle and LV level is derived by linear interpolation. Correction for translational and rotational misalignment between the short-axis cine and the 4D Flow CMR acquisition was performed using automated image registration as previously described [28].
For each volumetric element (voxel) the KE was computed as KE ¼ 1 2 ρ blood :V voxel :v 2 , with ρ blood being the density of blood (1.06 g/cm 3 ), V voxel the voxel volume and v the velocity magnitude. For each phase, the total KE within the LV was obtained by summation of the KE of every voxel. Time-resolved kinetic energy curves were generated to derive physiologically relevant parameters.
LV blood flow KE parameters which were derived for the complete cardiac cycle are described in Table 1/ Fig. 1.
The novel KE parameters studied in this study are detailed below:

In-plane KE
We define the in-plane KE as the sum of all kinetic energy in the x-y direction, in the short-axis LV from base to apex. In this study, the in-plane KE is represented as a percentage of the total LV KE, this negates the issue with volume normalisation. Most of the blood flow in the LV should be in the through-plane directionboth in systole (towards aortic valve) and during diastole (towards apex from the mitral valve) and an increase in the proportion of in-plane flow may be associated with pathological flow. Hence, in this study, the in-plane KE parameter was computed mainly to better understand the in-plane flow dynamics within the LV cavity.

Time difference (TD)
We also computed the time difference (TD) to peak early mitral inflow velocity (E-wave) from the base of the LV to mid-ventricle. This transit time or TD should be higher if the mitral valve propagation velocity (Vp), as measured by M-Mode echocardiography is lower [29]. Hence, the transit time of the peak KE from base to mid-ventricle, described as the TD in this study, may represent delayed filling.

Minimal KE
This parameter simply represents the minimal KE of the LV flow at any time point during the complete cardiac cycle. We computed this parameter mainly to understand how KE is preserved in health versus disease. The minimal KE of the LV is likely to be happen when the haemodynamic forces are minimal, and LV is not moving. A rise in the minimal KE may be influenced by heterogenous haemodynamic forces described in the LV fluid dynamics model secondary to dyssynchronous or late systolic contractility seen in infarcted segments [30].
To allow comparison between patients and healthy subjects, all KE parameters were normalized to the LV end-diastolic volume (EDV) and accordingly reported in μJ/ml. In sub-analysis, we also normalised the KE parameters to the stroke volume to develop insight into KE spent per unit of stroke volume.

Intra-/inter-observer variability
For inter-observer variability, SC and RVDG contoured the short-axis LV cine volumetric stack in 10 random healthy control subjects and 10 random MI patients, blinded to each other's analysis. Automated KE parameters were again generated using the endocardial contours from second observer. For intra-observer variability, SC re-analysed the LV short-axis cines for the same 10 healthy control subjects and 10 MI patients MI after 3 months. Akin to inter-observer variability, automated KE parameters were generated using the new endocardial contours by the same observer using the time-resolved methods previously described [31].

Statistical analysis
Statistical analysis was performed using SPSS® Statistics 21.0 (International Business Machines, Armonk, New York, USA). Continuous measurements are presented as mean ± standard deviation. Quantitative flow imaging parameters expected to be non-parametric were presented as median and inter-quartile ranges (IQR). Demographic comparisons were performed with an independent samples t-test. Comparisons between two groups with non-parametric data were made using Mann-Whitney U test. Intra−/inter-observer reliability tests were done by inter-class correlation coefficient. In different categories of LVEF, post hoc analysis was done by Kruskal-Wallis H test. Association of infarct size to KE parameters was done by Spearman's rank correlation coefficient test. In multivariate analysis, a forward-conditional method was used for regression and parameters with statistical significance from one-way analysis (p < 0.05) were chosen for multi-variate analysis. A p-value < 0.05 was considered statistically significant.

Demographic characteristics
Of the 55 patients with acute and chronic MI patients initially identified and recruited, 48 patients completed the full study protocol (Fig. 2). Controls and patients were matched for age, (52 ± 11 years versus 57 ± 11 years, P = 0.1), gender and body surface area ( Table 2).

Baseline CMR data
Heart rate was comparable between MI patients and healthy controls (65 ± 12 bpm versus 64 ± 13 bpm, P = 0.75). Most baseline volumetric CMR parameters were significantly higher in MI patients than healthy controls ( Table 2). Infarct size was similar among patients with acute and chronic MI (24.3 ± 14% vs. 19.6 ± 11%, P = 0.21). The majority of patients had anterior MI (n = 35) versus inferior/posterior MI (n = 13). The 4D flow CMR acquisition time was 8 ± 2 min.

Global LV flow KE in MI
Average LV flow KEi EDV was significantly higher in healthy controls than in MI patients (8.5 ± 3 μJ/ml versus 6.5 ± 3 μJ/ml, P = 0.02) ( Table 2) (Fig. 3). Similarly, LV flow minimal KEi EDV was significantly lower in patients. The proportion of in-plane KE was not significantly different in healthy controls and MI patients. Diastolic peak E-wave KEi EDV was significantly lower in MI patients compared to the healthy controls (P < 0.05). Time differences to peak E-wave KEi EDV were significantly higher in MI patients compared to healthy controls (11 ± 31msecs versus 34 ± 29msecs, P < 0.01) (Fig. 4).
In logistic regression analysis of all the significantly different KEi EDV parameters in MI patients (Table 2), LV blood flow systolic KEi EDV demonstrated the strongest association with the presence of MI (Beta = − 0.24, standard error = 0.08, P < 0.01).

Association with cardiac haemodynamic parameters
Kinetic energy parameters that were significantly associated with heart rate included the following: LV, minimal and diastolic KEi EDV (P < 0.05) ( Table 3). For these parameters, between EF group analysis was adjusted to the heart rate.
Kinetic energy parameters that were significantly associated with stroke volume included the following: systolic and E-wave KEi EDV , the proportion of in-plane KE and TD (P < 0.05) ( Table 3). In sub-group analysis, we did not identify significant association between stroke volume and systolic KEi EDV in the controls (R = 0.31, P = 0.17). This association was mainly driven by variation of stroke volume in the MI cohort (R = 0.39, P < 0.01). However, KE parameters associated with cardiac index included the following: LV, systolic, diastolic, E-wave KEi EDV and in-plane KE (P < 0.05).

LV systolic impairment and flow KE
Comparison of LV flow KE parameters (raw, normalised to end-diastolic volume and normalised to stroke volume) with the degree of LV systolic impairment are detailed in Table 4. Several LV flow KEi EDV parameters demonstrated significant differences between groups of patients, defined by the degree of LV systolic impairment (P-value< 0.05). These included: minimal, systolic, in-plane KEi EDV parameters and TD to peak E-wave KE propagation. The proportion of in-plane KE and minimal KEi EDV parameters were significantly reduced inpreservedpEF patients when compared to normal EF healthy controls but were increased compared with controls in patients with worsening EF (Figs. 5, 6c and d). In exploratory analysis, there was a modest association between the proportion of in-plane KE and minimal KEi EDV parameters (R = 0.35, P < 0.01). The systolic KEi EDV decreased significantly with worsening EF (Fig. 6a). TD to peak E-wave KE propagation demonstrated a linear increasing trend with worsening EF (P < 0.01).

Infarct size and flow KE
Infarct size was significantly associated with minimal, systolic and peak A-wave KEi EDV , in-plane KE and the   (Table 5). Infarct size was negatively associated with systolic KEi EDV (rho = − 48, P = 0.01) and positively associated with the proportion of in-plane KE (rho = 0.5, P = 0.001) (Fig. 7). TD of peak E-wave KE propagation had a modest positive association with MI size (rho = 0.33, P = 0.02). In multiple regression analysis, in-plane KE was independently associated with infarct size (Beta = 1.1, P = 0.001) ( Table 5).

Anterior versus inferior/posterior infarction
Thirty-five patients (73%) had anterior MI and 13 (27%) had inferior/posterior MI. CMR baseline parameters were similar in both infarct locations (Table 6). Infarct size was significantly larger in patients with anterior MI (14.4 ± 12% vs 23 ± 18%; P < 0.01). Controlling for infarct size, no significant differences were noted for the baseline CMR parameters and the different parameters of KEi EDV between anterior and inferior/posterior MI.

Acute versus chronic MI
Baseline demographics were comparable between patients with either acute or chronic MI (Table 7). From the baseline CMR parameters, only stroke volume was significantly higher in chronic MI versus acute (37 ± 10 ml/m2 versus 44 ± 8 ml/m2, P-value = 0.02). KEi EDV parameters and LV EF were comparable in both groups.    Table 4). We noted that rise in KE in diastole per stroke volume over took systolic KE per stroke volume from the moderately impaired LV group. Hence, in moderately and severely impaired LV, there is more KE per stroke volume in diastole than systole (Fig. 8).

Intra-observer and inter-observer variability
Intra-observer and inter-observer variability results for global KE parameters are detailed in the online Additional file 1: Table S2. In summary, all global KE parameters demonstrated excellent intra-class correlation coefficient (ICC) (average 0.99, P-value> 0.9). In addition, the mean bias (%) for all KEi EDV parameters was 2 ± 9% for intra-observer tests and 3 ± 9% for inter-observer tests. TD to peak E-wave flow KE also demonstrated very high intra-class correlation coefficient (ICC = 0.94, 95% CI 0.88 to 0.97).

Discussion
In this study, we mapped different parameters of LV blood flow KE in MI patients and matched healthy control subjects. We found that the majority of LV blood flow KEi EDV parameters are reduced in patients with MI. In addition, decrease in systolic KEi EDV of the LV blood flow is independently associated with the presence of MI and the proportion of in-plane KE of the LV blood flow is independently associated with the infarct size. Another important finding of this study was that worsening degree of LV systolic function is associated TD of peak E-wave KE propagation. MI patients with preserved EF have significantly reduced minimal KEi EDV , the proportion of in-plane KE and increased TD of peak E-wave KE propagation when compared with healthy controls with nEF. Kinetic energy per stroke volume revealed higher systolic KE versus diastolic KE in normal EF, preserved EF and mildly reduced EF groups. However, patients with moderate and severely reduced EF demonstrated reversal and had higher diastolic KE per stroke volume than systolic KE. Finally, global LV blood flow KE parameters demonstrated good intra−/inter-observer variability.

Global LV flow KE in MI
In the present study, the global LV flow KE iEDV , which in essence is the global density of the velocity profile of intra-cavity blood flow, was significantly lower in patients with MI than controls. This is in keeping with results by Kanski et al., who found reduced indexed LV KE in patients with moderate-to-severe heart failure versus controls. (6.3 ± 2.2 μJ/ml vs 8.0 ± 2.1 μJ/ml, p = 0.025) [16]. Even though the diastolic parameter of LV KE was not significantly different in MI patients, the peak E-wave KE was significantly lower in MI patients. This would be consistent with early pre-clinical studies that have demonstrated that any condition that interferes with normal regional systolic function is expected to modify the pattern of the normal early diastolic intraventricular pressure gradients. The present study also showed that the presence of MI is most strongly associated with systolic KE iEDV of LV blood flow. As a consequence of MI and regional hypo-contractility, there is reduced regional systolic mechanical force imparted on the intra-cavity blood consequently due to lowered pressure gradient between the LV and the aorta, the overall thrust of systolic blood flow is reduced as demonstrated by drop in the systolic blood flow energetics. This is also supported by the finding that LV stroke volume was only associated to systolic KEi EDV in MI cohort and not in the controls.

Delayed filling of the LV
During diastolic filling of the LV in health, we already know through M-mode Doppler echocardiography that blood flow into the LV cavity from base to apex happens very fast with a mitral valve flow propagation velocity (Vp) of above 50 cm/sec [32]. Even though in this study we did not measure the Vp, we measured the transit time (or TD in this study) for the peak E-wave to travel from base to mid-ventricle and hence, it is plausible to conclude that if the Vp will decrease, the time delay for the peak inflow will increase. Accordingly, we observed an increase in the time differences for peak E-wave KE propagation in MI patients compared with controls. This delay reflects restrictive LV filling and possible normalisation of intra-ventricular pressure gradients. Intra-ventricular pressure gradients represent a suction force and have been attributed to energy consuming, active LV relaxation [33,34]. The significant increase in TDs in patients with MI and preserved EF versus controls with normal EF demonstrates that restrictive filling precedes any detectable systolic impairment. The TDs measured in this study offer physiological insight similar to mitral valve propagation velocity on echocardiography. Compared with mitral valve propagation velocity, the TDs measured here are semi-automatically computed and have low intra−/inter-observer variability. From the clinical perspective, this parameter may offer semi-automated and possible reliable assessment of the degree of LV diastolic impairment and the LV filling pressures. The reliability and the prognostic implications of using TD as a bio-imaging marker needs to be further tested in future studies.

Association of LV flow KE to the degree of systolic impairment
The patterns of changes between these LV flow KE parameters to the degree of LV systolic impairment were different. The systolic KEi EDV had a more linear trend to decrease with worsening EF. Decreased LV contractile function leads to impaired emptying and increased preload in subsequent cycles. Even if this causes stroke volumes to return to normal due to the Frank-Starling mechanism, the increased ventricular dimensions will cause the KEi EDV to decrease. The TDs for peak E-wave KE propagation demonstrated a more linear trend to increase with worsening EF. The minimal KEi EDV and the proportion of in-plane KE significantly reduced from normal EF in healthy controls to MI patients with preserved EF and then had a strong trend to increase with decreasing EFa finding that at first seems paradoxical. We speculate that the in-plane KE of LV flow is made up of several in-plane blood flow movements during the complete cardiac cycle. These include: the physiological systolic in-plane movement of blood towards the outflow track during LV contraction, in-plane flow within the vortex in early and late diastolic filling. The systolic and vortex associated in-plane flow are plausibly decreased in MI patients with preserved EF as they have subtle mechanical dysfunction with preserved LV volumes. This finding would be consistent with other studies demonstrating alterations in diastolic inflow vortex strength in heart failure patients with preserved EF (HFpEF) [35]. However, progressive LV impairment and dilatation also causes increased sphericity, which will in turn change flow conditions inside the cavity to a 'meta-stable' state with a large, swirling vortex which encompasses the majority of the LV [36]. This vortex flow includes transversal thrusts which will show up as in-plane KE and become increasingly prominent as the ventricle remodels further. In addition, given the fact our reference spatial plane is the atrioventricular valve, the angular differences between the inflow and outflow direction which follows from the non-parallel orientation of the left atrium and aorta may result in significant outflow in-plane component.
As the pattern of changes to the degree of LV impairment was similar for in-plane KE and minimal KEi EDV , we speculate that the rise in the proportion of in-plane KE with worsening EF is associated with an increase in minimal KEi EDV . This is possibly explained by the fact that progressive LV impairment leads to LV dilatation and reduced filling of the LV causing an overall increase in the in-plane flow, which results in higher minimal KEi EDV for the LV. This observation is further supported by the analysis of the KE parameters per stroke volume which demonstrate the all blood flow KE parameters rise significantly for each stroke volume with the degree of LV impairment. In addition, once the degree of LV impairment reaches the moderate zone, the systolic KE per stroke volume becomes lesser than the diastolic KE per stroke volume. This perhaps reflect the loss of mechanical push by LV in systole combined with restrictive early filling in diastole associated with high velocities.

Infarct characteristics and LV flow KE
Results from our study show that an increase in infarct size increases the in-plane KE of LV blood flow. Such pathological rise in the in-plane KE may exert hete rogenous haemodynamic forces on the LV wall, possibly contributing to more dilatation and increase in endothelial dysfunction in the endocardium [37][38][39][40]. In our study, anterior infarcts were significantly larger than inferior infarcts. When adjusted for MI size, MI location (anterior versus inferior) was not significantly associated with LV function or blood flow KE parameters. In a previous study done by Kanski et al., in sub-group of patients with MI, multiple stepwise linear regression analysis, averaged LV KE was associated with LV end-diastolic volume and not IS. Similarly, in this study averaged LV KEi EDV was not associated with IS (P = 0.09), however, in regression, the in-plane demonstrated independent association to the IS. This would indicate that some components of LV blood flow energetics are more influenced by infarct size than others and its associated mechanical dysfunction than its location.

Future clinical applications
We have demonstrated that using 4D flow that it is possible to detect subtle changes in cardiac function after MI even in patients with preserved EF. Longitudinal prospective studies are needed to establish whether LV flow KE has incremental prognostic impact over established markers such as EF, scar burden and microvascular obstruction. Studies are also needed to establish whether changes in intra-cardiac flow can be altered by medical therapy. After MI, LV remodelling leads to progressive LV cavity dilatation and subsequent heart failure. Using 4D flow it will be possible to investigate these mechanisms and assess the impact of existing and potentially novel pharmacological interventions on these processes.
To compute LV blood flow KE parameters, time-resolved LV volumetric endocardial contours are used [31]. Hence, the main influence on LV energetics variation is mainly secondary to endocardial contour delineation. If this is done consistently and methodically for the complete cardiac cycle as in this study and previous published literature, it results in not only highly reproducible LV volumetric assessment but also LV blood flow energetics.

Study limitations
Respiratory navigation was omitted for the 4D flow acquisition which could have influenced KE parameters. However, whole-heart 4D flow head-to-head comparison studies have also demonstrated that non-respiratory navigated acquisition of 4D flow is comparable to respiratory navigated acquisition for intra-cardiac KE quan tification [41]. In addition, a recent study validated a non-respiratory navigated 4D Flow EPI acceleration sequence for clinical use [21]. The temporal resolution of the 4D flow was 40 ms, which may affect the quality of KE and TD assessment. MI patients were prescribed 4D flow post contrast versus healthy controls who didn't receive any contrast. Post contrast 4D flow improves the signal-noise-ratio and may have introduced some bias when comparing controls versus patients. This study used the start of the cardiac cycle as defined by the scanner's electrocardiogram to calculate TDs. Scanner's electrocardiogram is not always correct, and this may result in less precise TD estimates. We used The LV geometry was defined by LV cine stack which was done using breath-hold technique while the 4D flow was done using free breathing. Hence, although spatial miss-registration was corrected for, other issues still remain including difference in heart rate and physiological conditions. This may have impact on the time-varying flow characteristics which could not be corrected for. Results from this study cannot be applied to patients with significant valvulopathy, cardiomyopathies and congenital heart disease.

Conclusions
Post MI, reduction in LV function results in reduction in LV blood flow KE and is most strongly associated with the systolic LV blood flow KE. Increase in infarct size results in proportionate increase in the in-plane KE of LV blood flow. Degree of LV impairment is associated with TD of peak E-wave KE. Even in patient with preserved EF post MI, LV blood flow KE mapping demonstrated significant changes in the in-plane KE, the minimal KEi EDV and the TD. Further longitudinal studies are warranted to investigate the long-term clinical significance of mapping LV blood flow KE parameters.