Skip to main content

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

Abstract

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.

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.

Fig. 1
figure 1

Illustration of all the left ventricular blood flow kinetic energy (KE) parameters investigated in this study in a healthy control

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.

Methods

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.

Fig. 2
figure 2

Study design. 4D = four-dimensional, 4DF = four-dimensional flow, CMR = cardiovascular magnetic resonance imaging, HV = healthy volunteers, KE = kinetic energy, LGE = late gadolinium enhancement imaging, LV = left ventricle, m = month, MI = myocardial infarction, CMR = cardiovascular magnetic resonance and SAX = short-axis

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 valvular heart disease, estimated glomerular filtration rate < 30 ml/min/1.73m2, 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. 1.

    Survey images

  2. 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. Typical parameters for bSSFP cine were as follows: SENSE factor 2, flip angle 60°, echo time (TE) 1.5 milliseconds, repetition time (TR) 3 milliseconds, field of view (FOV) 320–420 mm according to patient size, slice thickness 8 mm, and 30 phases per cardiac cycle.

  3. 3.

    Contrast injection (Magnevist: 0.2 mmol/kg) [only done in the patient population]

  4. 4.

    For whole heart 4D flow, field of view (FOV) was planned in trans-axial plane making sure whole heart was in FOV. 4D flow was done using fast field echo (FFE) pulse sequence (EPI based, 3D) with retrospective electrocardiogram (ECG)-triggering. The acquisition voxel size was kept as close as possible to 3 × 3 × 3 mm3. FOV and number of slices (i.e., the 3D volume) was adapted to the subject’s size. Scan parameters: TE 3.5 ms, TR 10 ms, flip angle 10°, FOV 400 mm, number of signal averages 1. VENC 150 cm/sec. Acceleration was achieved by EPI factor 5. Free breathing was allowed, and no respiratory motion compensation was used. Number of slices was 40, temporal resolution used was 40 ms and the number of reconstructed phases was set to 30. EPI sequence is detailed in the Additional file 1: Table S1 as previously validated [21].

  5. 5.

    Modified Look-Locker inversion recovery (MOLLI) to determine the T1-inversion time [only done in the patient population]

  6. 6.

    LGE imaging was done at 15-min from gadolinium-based contrast injection in MI patients only. LGE-imaging was done by phase sensitive inversion recovery (PSIR) spoiled gradient echo (GE) sequence as per previous published methods [22, 23]. PSIR sequence details are as follows: SENSE factor 1.7, typical TE/TR of 3.0/6.1 msec, flip angle of 25°, slice thickness of 10 mm and with Look-Locker scout determined T1-inversion time. [only done in the patient population]

4D flow data reconstruction method

4D flow encoding was performed by standard 4-point encoding. Online/offline 4D flow data quality assurance checks were done as per previous published literature [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].

The degree of LV systolic impairment was categorised as follows: normal EF in controls, normal LVEF. In patients with MI, preserved LVEF was ≥55%, mildly impaired LVEF 45–54%, moderately impaired LVEF 35–44% and severely impaired LVEF < 35%" In patients with MI.

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 constructed 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 \( \mathrm{KE}=\frac{1}{2}\ {\rho}_{blood}.{V}_{voxel}.{v}^2 \), with ρblood being the density of blood (1.06 g/cm3), Vvoxel 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.

Table 1 Description of all left ventricular (LV) kinetic energy (KE) flow parameters investigated in this study. In this study all primary KE parameters are indexed to LV end diastolic volume (EDV) and expressed as KEiEDV

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 direction – both 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 multi-variate 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.

Results

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

Table 2 Demographics and CMR parameters in healthy controls versus MI patients

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 KEiEDV 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 KEiEDV 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 KEiEDV was significantly lower in MI patients compared to the healthy controls (P < 0.05). Time differences to peak E-wave KEiEDV were significantly higher in MI patients compared to healthy controls (11 ± 31msecs versus 34 ± 29msecs, P < 0.01) (Fig. 4).

Fig. 3
figure 3

Left ventricular (LV) blood flow KE mapping. The short-axis whole LV blood flow KE maps represented here are during peak systole. The KE of blood flow in the out-flow tract is higher in the control versus the myocardial infarction (MI) patient

Fig. 4
figure 4

Case examples of KE curves in control and two MI patients with preserved LVEF and reduced LVEF. First panel: The total KE of the LV is reduced in MI patients. The minimal KE demonstrates an increase in MI patient with reduced LVEF. Second panel: Diastolic KE curves at different levels of the LV (base, mid-ventricular and apex). The diastolic KE curves demonstrate blunting of peak E-wave KE in the two patient groups. Time differences to peak E-wave KE propagation increase in each group (orange arrows between the two dotted black lines). EF = ejection fraction, nEF = normal ejection fraction, pEF = preserved ejection fraction, KE = kinetic energy, LV = left ventricle and MI = myocardial infarction

In logistic regression analysis of all the significantly different KEiEDV parameters in MI patients (Table 2), LV blood flow systolic KEiEDV 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 KEiEDV (P < 0.05) (Table 3). For these parameters, between EF group analysis was adjusted to the heart rate.

Table 3 Association of LV kinetic energy parameters to LV haemodynamic parameters

Kinetic energy parameters that were significantly associated with stroke volume included the following: systolic and E-wave KEiEDV, 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 KEiEDV 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 KEiEDV 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 KEiEDV 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 KEiEDV parameters and TD to peak E-wave KE propagation. The proportion of in-plane KE and minimal KEiEDV 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 KEiEDV parameters (R = 0.35, P < 0.01). The systolic KEiEDV 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).

Table 4 Post-hoc analysis of several LV KE flow parameters in different LV function sub-groups
Fig. 5
figure 5

In-plane velocity streamlines in two case examples from the study demonstrating higher in-plane (xy-axis) flow during diastasis in a patient with myocardial infarction and worse EF (Panel b) versus an MI patient with a smaller infarct and preserved LVEF (Panel a). The in-plane KE was computed as the sum of all KE in the x-y direction, in the short-axis left ventricle from base to apex. It is represented as a percentage of the total LV KE. KE = kinetic energy and EF = ejection fraction

Fig. 6
figure 6

Bar charts of left-ventricular (LV) flow kinetic energy parameters and their variations within ejection fraction (EF) sub-groups. The LV systolic KEiEDV (Panel a) demonstrated reduction with worsening LVEF versus TD (Panel b) which demonstrated significant increase with worsening LVEF. Note the U-shaped pattern of minimal (Panel c) and in-plane (Panel d) KE of the LV. (bars = median; error bars = inter-quartile range). KE = kinetic energy, LV = left ventricle, nEF = normal ejection fraction, pEF = preserved ejection fraction, Mild rEF = mildly reduced ejection fraction, Moderate rEF = moderately reduced EF and Severe rEF = severely reduced EF

Infarct size and flow KE

Infarct size was significantly associated with minimal, systolic and peak A-wave KEiEDV, in-plane KE and the TD for peak E-wave KE (Table 5). Infarct size was negatively associated with systolic KEiEDV (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).

Table 5 LV flow KEiEDV parameters and their association to the infarct size in patients with MI
Fig. 7
figure 7

Scatter plots of infarct size and its positive association to in-plane KE and negative association to LV systolic flow KE. KE = kinetic energy, LV = left ventricle and MI = myocardial infarction

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 KEiEDV between anterior and inferior/posterior MI.

Table 6 Kinetic energy distribution in patients with anterior and inferior infarcts. Units for all KE parameters are μJ/ml

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). KEiEDV parameters and LV EF were comparable in both groups.

Table 7 Comparison of study results in acute versus chronic myocardial infarction

Kinetic energy of blood flow per stroke volume

All KEiSV parameters (LV, minimal, systolic, diastolic, E-wave and A-wave) demonstrated significant rise between different groups of EF (P < 0.05) (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).

Fig. 8
figure 8

Comparison of systolic and diastolic kinetic energy per LV stroke volume in different categories of LV impairment. KE = kinetic energy, nEF = normal ejection fraction, pEF = preserved ejection fraction, Mild rEF = mildly reduced ejection fraction, Moderate rEF = moderately reduced EF, Severe rEF = severely reduced EF and SV = stroke volume

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 KEiEDV 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 KEiEDV parameters are reduced in patients with MI. In addition, decrease in systolic KEiEDV 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 KEiEDV, 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 KEiEDV, 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 KEiEDV 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 KEiEDV 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 KEiEDV 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 KEiEDV to decrease.

The TDs for peak E-wave KE propagation demonstrated a more linear trend to increase with worsening EF. The minimal KEiEDV 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 EF – a 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 KEiEDV, we speculate that the rise in the proportion of in-plane KE with worsening EF is associated with an increase in minimal KEiEDV. 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 KEiEDV 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 heterogenous 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 KEiEDV 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 quantification [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 KEiEDV and the TD. Further longitudinal studies are warranted to investigate the long-term clinical significance of mapping LV blood flow KE parameters.

Abbreviations

2D:

Two-dimensional

3D:

Three-dimensional

4D:

Four-dimensional

bSSFP:

Balanced steady state free precession

CMR:

Cardiovascular magnetic resonance

ECG:

Electrocardiogram

EF:

Ejection fraction

EPI:

Echo-planar Imaging acceleration

FOV:

Field of view

GE:

Gradient echo

HFpEF:

Heart failure with preserved ejection fraction

ICC:

Intra-class correlation coefficient

IQR:

Inter-quartile range

KE:

Kinetic energy

LGE:

Late gadolinium enhancement

LPC:

Local phase correction

LV:

Left ventricle/left ventricular

MI :

Myocardial infarction

MOLLI:

Modified Look Locker inversion recovery

PPCI:

Primary percutaneous coronary intervention

PSIR:

Phase sensitive inversion recovery

STEMI:

ST-elevation myocardial infarction

TD:

Time difference

TE:

Echo time

TR:

Repetition time

VP:

Propagation velocity

References

  1. Yamamoto K, Nishimura RA, Chaliki HP, Appleton CP, Holmes DR, Redfield MM. Determination of left ventricular filling pressure by Doppler echocardiography in patients with coronary artery disease: critical role of left ventricular systolic function. J Am Coll Cardiol. 1997;30:1819–26.

    Article  PubMed  CAS  Google Scholar 

  2. Hammermeister KE, DeRouen TA, Dodge HT. Variables predictive of survival in patients with coronary disease. Selection by univariate and multivariate analyses from the clinical, electrocardiographic, exercise, arteriographic, and quantitative angiographic evaluations. Circulation. 1979;59:421–30.

    Article  PubMed  CAS  Google Scholar 

  3. Aikawa Y, Rohde L, Plehn J, Greaves SC, Menapace F, Arnold JMO, et al. Regional wall stress predicts ventricular remodeling after anteroseptal myocardial infarction in the healing and early afterload reducing trial (HEART): an echocardiography-based structural analysis. Am Heart J. 2001;141:234–42.

    Article  PubMed  CAS  Google Scholar 

  4. Alter P, Koczulla AR, Nell C, Figiel JH, Vogelmeier CF, Rominger MB. Wall stress determines systolic and diastolic function - characteristics of heart failure. Int J Cardiol. 2016;202:685–93.

    Article  PubMed  Google Scholar 

  5. French BA, Kramer CM. Mechanisms of post-infarct left ventricular remodeling. Drug Discov Today Dis Mech. 2007;4:185–96.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Meta-Analysis Research Group in Echocardiography (MeRGE) AMI Collaborators, Møller JE, Whalley GA, Dini FL, Doughty RN, Gamble GD, et al. Independent Prognostic Importance of a Restrictive Left Ventricular Filling Pattern After Myocardial Infarction: An Individual Patient Meta-Analysis: Meta-Analysis Research Group in Echocardiography Acute Myocardial Infarction. Circulation. 2008;117:2591–8.

    Article  Google Scholar 

  7. Pearlman AS, Stevenson JG, Baker DW. Doppler echocardiography: applications, limitations and future directions. Am J Cardiol. 1980;46:1256–62.

    Article  PubMed  CAS  Google Scholar 

  8. PREC O, KATZ LN, Sennett L, Rosenman RH, Fishman AP, Hwang W. Determination of kinetic energy of the heart in man. Am J Phys. 1949;159:483–91.

    CAS  Google Scholar 

  9. Carlsson M, Töger J, Kanski M, Bloch KM, Ståhlberg F, Heiberg E, et al. Quantification and visualization of cardiovascular 4D velocity mapping accelerated with parallel imaging or k-t BLAST: head to head comparison and validation at 1.5 T and 3 T. J Cardiovasc Magn Reson. 2011;13:55.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Eriksson J, Bolger AF, Ebbers T, Carlhäll C-J. Four-dimensional blood flow-specific markers of LV dysfunction in dilated cardiomyopathy. Eur Heart J Cardiovasc Imaging. 2013;14:417–24.

    Article  PubMed  Google Scholar 

  11. Eriksson J, Dyverfeldt P, Engvall J, Bolger AF, Ebbers T, Carlhall CJ. Quantification of presystolic blood flow organization and energetics in the human left ventricle. AJP Hear Circ Physiol. 2011;300:H2135–41.

    Article  CAS  Google Scholar 

  12. Bolger AF, Heiberg E, Karlsson M, Wigström L, Engvall J, Sigfridsson A, et al. Transit of blood flow through the human left ventricle mapped by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2007;9:741–7.

    Article  PubMed  Google Scholar 

  13. Crandon S, Elbaz MSM, Westenberg JJM, van der Geest RJ, Plein S, Garg P. Clinical applications of intra-cardiac four-dimensional flow cardiovascular magnetic resonance: a systematic review. Int J Cardiol. 2017;249:486–93.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Eriksson J, Carlhäll CJ, Dyverfeldt P, Engvall J, Bolger AF, Ebbers T. Semi-automatic quantification of 4D left ventricular blood flow. J Cardiovasc Magn Reson. 2010;12:9.

    Article  PubMed  PubMed Central  Google Scholar 

  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. AJP Hear Circ Physiol. 2012;302:H893–900.

    Article  CAS  Google Scholar 

  16. Kanski M, Arvidsson PM, Töger J, Borgquist R, Heiberg E, Carlsson M, et al. Left ventricular fluid kinetic energy time curves in heart failure from cardiovascular magnetic resonance 4D flow data. J Cardiovasc Magn Reson. 2015;17:111.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Svalbring E, Fredriksson A, Eriksson J, Dyverfeldt P, Ebbers T, Bolger AF, et al. Altered Diastolic Flow Patterns and Kinetic Energy in Subtle Left Ventricular Remodeling and Dysfunction Detected by 4D Flow MRI. Aliseda A, editor. PLoS One. 2016;11:e0161391.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Wong J, Chabiniok R, deVecchi A, Dedieu N, Sammut E, Schaeffter T, et al. Age-related changes in intraventricular kinetic energy: a physiological or pathological adaptation? Am J Physiol Hear Circ Physiol. 2016;310:H747 LP–H755.

    Article  Google Scholar 

  19. Thygesen K, Alpert JS, Jaffe AS, Simoons ML, Chaitman BR, White HD, et al. Third universal definition of myocardial infarction. J Am Coll Cardiol. 2012;60:1581–98.

    Article  PubMed  Google Scholar 

  20. Killip T, Kimball J. Treatment of myocardial infarction in a coronary care unit. A two year experience with 250 patients. Am J Cardiol. 1967;20:457–64.

    Article  PubMed  Google Scholar 

  21. Garg P, Westenberg JJM, van den Boogaard PJ, Swoboda PP, Aziz R, Foley JRJ, et al. Comparison of fast acquisition strategies in whole-heart four-dimensional flow cardiac MR: two-center, 1.5 tesla, phantom and in vivo validation study. J Magn Reson Imaging. 2017;47:272–81.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Garg P, Kidambi A, Foley JRJ, Al MT, Ripley DP, Swoboda PP, et al. Ventricular longitudinal function is associated with microvascular obstruction and intramyocardial haemorrhage. Open Hear. 2016;3:e000337.

    Article  Google Scholar 

  23. Garg P, Kidambi A, Swoboda PP, Foley JRJ, Musa TA, Ripley DP, et al. The role of left ventricular deformation in the assessment of microvascular obstruction and intramyocardial haemorrhage. Int J Cardiovasc Imaging. 2017;33:361–70.

    Article  PubMed  Google Scholar 

  24. Walker PG, Cranney GB, Scheidegger MB, Waseleski G, Pohost GM, Yoganathan AP. Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. J Magn Reson Imaging. 1993;3:521–30.

  25. MacDonald ME, Forkert ND, Pike GB, Frayne R, Drangova M, Fenster A. Phase Error Correction in Time-Averaged 3D Phase Contrast Magnetic Resonance Imaging of the Cerebral Vasculature. Zhang H, editor. PLoS One. 2016;11:e0149930.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Lotz J, Meier C, Leppert A, Galanski M. Cardiovascular flow measurement with phase-contrast MR imaging: basic facts and implementation. Radiographics. 2002;22:651–71.

    Article  PubMed  Google Scholar 

  27. Garg P, Broadbent DA, Swoboda PP, Foley JRJ, Fent GJ, Musa TA, et al. Acute infarct extracellular volume mapping to quantify myocardial area at risk and chronic infarct size on cardiovascular magnetic resonance imaging. Circ Cardiovasc Imaging. 2017;10:e006182.

    Article  PubMed  Google Scholar 

  28. Elbaz MSM, van der Geest RJ, Calkoen EE, de Roos A, Lelieveldt BPF, Roest AAW, et al. Assessment of viscous energy loss and the association with three-dimensional vortex ring formation in left ventricular inflow: in vivo evaluation using four-dimensional flow MRI. Magn Reson Med. 2017;77(2):794-805.

  29. Brun P, Tribouilloy C, Duval AM, Iserin L, Meguira A, Pelle G, et al. Left ventricular flow propagation during early filling is related to wall relaxation: a color M-mode Doppler analysis. J Am Coll Cardiol. 1992;20:420–32.

    Article  PubMed  CAS  Google Scholar 

  30. Domenichini F, Pedrizzetti G. Hemodynamic forces in a model left ventricle. Phys Rev Fluids. 2016;1:083201.

    Article  Google Scholar 

  31. 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:821–8.

    Article  PubMed  Google Scholar 

  32. De Boeck BWL, Oh JK, Vandervoort PM, Vierendeels JA, van der Aa RPLM, Cramer M-JM. Colour M-mode velocity propagation: a glance at intra-ventricular pressure gradients and early diastolic ventricular performance. Eur J Heart Fail. 2005;7:19–28.

    Article  PubMed  Google Scholar 

  33. Greenberg NL, Vandervoort PM, Firstenberg MS, Garcia MJ, Thomas JD. Estimation of diastolic intraventricular pressure gradients by Doppler M-mode echocardiography. Am J Physiol Heart Circ Physiol. 2001;280:H2507–15.

    Article  PubMed  CAS  Google Scholar 

  34. Yotti R, Bermejo J, Antoranz JC, Desco MM, Cortina C, Rojo-Alvarez JL, et al. A noninvasive method for assessing impaired diastolic suction in patients with dilated cardiomyopathy. Circulation. 2005;112:2921–9.

    Article  PubMed  Google Scholar 

  35. Abe H, Caracciolo G, Kheradvar A, Pedrizzetti G, Khandheria BK, Narula J, et al. Contrast echocardiography for assessing left ventricular vortex strength in heart failure: a prospective cohort study. Eur Heart J Cardiovasc Imaging. 2013;14:1049–60.

    Article  PubMed  Google Scholar 

  36. Pedrizzetti G, La Canna G, Alfieri O, Tonti G. The vortex--an early predictor of cardiovascular outcome? Nat Rev Cardiol. Nature Publishing Group, a division of Macmillan Publishers Limited. All Rights Reserved. 2014;11:545–53.

    Article  PubMed  Google Scholar 

  37. Eriksson J, Zajac J, Alehagen U, Bolger AF, Ebbers T, Carlhäll C-J. Left ventricular hemodynamic forces as a marker of mechanical dyssynchrony in heart failure patients with left bundle branch block. Sci Rep. 2017;7:2971.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Eriksson J, Bolger AF, Ebbers T, Carlhall C-J. Assessment of left ventricular hemodynamic forces in healthy subjects and patients with dilated cardiomyopathy using 4D flow MRI. Physiol Rep. 2016;4(3).

  39. Arvidsson PM, Töger J, Carlsson M, Steding-Ehrenborg K, Pedrizzetti G, Heiberg E, et al. Left and right ventricular hemodynamic forces in healthy volunteers and elite athletes assessed with 4D flow magnetic resonance imaging. Am J Physiol Circ Physiol. 2017;312:H314–28.

    Article  Google Scholar 

  40. Pedrizzetti G, Arvidsson PM, Töger J, Borgquist R, Domenichini F, Arheden H, et al. On estimating intraventricular hemodynamic forces from endocardial dynamics: a comparative study with 4D flow MRI. J Biomech. 2017;60:203–10.

    Article  PubMed  Google Scholar 

  41. Kanski M, Töger J, Steding-Ehrenborg K, Xanthis C, Bloch KM, Heiberg E, et al. Whole-heart four-dimensional flow can be acquired with preserved quality without respiratory gating, facilitating clinical use: a head-to-head comparison. BMC Med Imaging. 2015;15:20.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Gavin Bainbridge, Caroline Richmond, Margaret Saysell and Petra Bijsterveld for their assistance in recruiting and collecting data for this study.

Dr. Hakan Arheden served as a JCMR Guest Editor for this manuscript.

Funding

This work was supported by the British Heart Foundation [FS/10/62/28409 to S.P.] and Dutch ZonMw [Project Number: 104003001 to J.W].

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

PG was involved in recruitment, CMR analysis, statistics and drafted the initial study manuscript. SV assisted in the recruitment and design of study. JJMW, MSME, MB, RN and SP developed the CMR protocol for the study. LAEB, PPS, JRJF, GJF, PGC and MECJH were major contributors in writing the manuscript and provided cross-cover for the study. SC did significant segmentation for the CMR studies and provided administrative support. RJvG contributed substantially in smooth running of the core lab at Leiden (LUMC) and analysis. AAM provided critical insight into heart failure associated mechanisms and how KE may provide mechanistic insight. RN, JPG, provided critical senior review of the final drafted manuscript. EDA and SP provided overall supervision for the study. All authors read and approved the final manuscript and agreed to be accountable for all aspects of the work.

Corresponding author

Correspondence to Pankaj Garg.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the National Research Ethics Service (12/YH/0169) in the UK and the institutional Medical Ethical Committee (P11.136) in Leiden. The study complied with the Declaration of Helsinki and all controls and patients gave written informed consent.

Consent for publication

All patients gave written informed consent before randomization.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Supplementary Material. (DOCX 30 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Garg, P., Crandon, S., Swoboda, P.P. et al. Left ventricular blood flow kinetic energy after myocardial infarction - insights from 4D flow cardiovascular magnetic resonance. J Cardiovasc Magn Reson 20, 61 (2018). https://doi.org/10.1186/s12968-018-0483-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12968-018-0483-6

Keywords