- Open Access
- Open Peer Review
Area of the pressure-strain loop during ejection as non-invasive index of left ventricular performance: a population study
Cardiovascular Ultrasound volume 17, Article number: 15 (2019)
Previous studies highlighted the usefulness of integrating left ventricular (LV) deformation (strain) and hemodynamic parameters to quantify LV performance. In a population sample, we investigated the anthropometric and clinical determinants of a novel non-invasive index of LV systolic performance derived from simultaneous registration of LV strain and brachial pressure waveforms.
Three hundred fifty-six randomly recruited subjects (44.7% women; mean age, 53.9 years; 47.5% hypertensive) underwent echocardiographic and arterial data acquisition. We constructed pressure-strain loops from simultaneously recorded two-dimensional LV strain curves and brachial pressure waveforms obtained by finger applanation tonometry. We defined the area of this pressure-strain loop during ejection as LV ejection work density (EWD). We reported effect sizes as EWD changes associated with a 1-SD increase in covariables.
In multivariable-adjusted analyses, higher EWD was associated with age, female sex and presence of hypertension (P ≤ 0.0084). In both men and women, EWD increased independently with augmentation pressure (effect size: + 59.1 Pa), central pulse pressure (+ 65.7 Pa) and pulse wave velocity (+ 44.8 Pa; P ≤ 0.0006). In men, EWD decreased with relative wall thickness (− 29.9 Pa) and increased with LV ejection fraction (+ 23.9 Pa; P ≤ 0.040). In women, EWD increased with left atrial (+ 76.2 Pa) and LV end-diastolic (+ 43.8 Pa) volume indexes and with E/e’ ratio (+ 51.1 Pa; P ≤ 0.026).
Older age, female sex and hypertension were associated with higher EWD. Integration of the LV pressure-strain loop during ejection might be a useful tool to non-invasively evaluate sex-specific and interdependent effects of preload and afterload on LV myocardial performance.
The performance of the left ventricle (LV) is determined by its intrinsic contractility (inotropy), by the tension on the LV wall at end-diastole (preload) and by the load against which the LV needs to eject blood (afterload) . If myocardial and vascular properties are matched, the heart can adequately respond to changes in pre- and afterload and regulate cardiac output and blood pressure (BP). However, the adaptive response to increased afterload appears impaired in patients with chronic hypertension  and with symptomatic heart failure [2, 3].
In line, community-based studies demonstrated the detrimental impact of increased afterload on LV structure and function [4,5,6,7]. Indeed, a long-term increased afterload and, consequently, a chronically increased cardiac performance lead to adverse LV maladaptation and increased LV oxygen requirements . In particular, older women appear susceptible to the detrimental effects of increased pulsatile load on LV diastolic function [4, 6, 7], which might be explained by the higher aortic stiffness, enhanced LV systolic performance, higher preload sensitivity and lower LV compliance in women as in men .
Echocardiographic techniques such as two-dimensional (2D) speckle tracking allow quantification of the relative myocardial deformation (i.e. strain) . Furthermore, LV strain and arterial function can be assessed simultaneously in a non-invasive way . In fact, recent developments in construction and quantification of pressure-strain loops illustrates the pressing need for a non-invasive, clinical tool that integrates both LV deformation and its loading conditions [2, 12]. Several indexes derived from these simultaneous recordings were suggested to reflect the interaction between the heart and the arteries (i.e. ventricular-arterial coupling, VAC) [2, 12]. For instance, we previously constructed LV pressure-strain loops from simultaneously recorded LV strain curves and carotid pressure waveforms, and defined the area of these loops during ejection as the LV ejection work density (EWD) .
In fact, EWD might better reflect LV systolic performance than peak LV deformation, as it integrates the instantaneous deformation as well as the instantaneous afterload against which the myocardium has to shorten [2, 11]. However, we should better understand this LV performance index before evaluating its utility in clinical decision-making. Therefore, we investigated in a population sample the anthropometric and clinical determinants of EWD. We also explored the relationship of EWD with indexes reflecting LV structure, LV diastolic function and arterial stiffness.
Materials and methods
The Flemish Study on Environment, Genes and Health Outcomes (FLEMENGHO) was approved by the Ethics Committee of the University of Leuven. We randomly recruited a family-based population sample in northern Belgium as described before . All subjects provided written informed consent. Between 2011 and 2016, we performed a VAC protocol including simultaneous echocardiography and finger applanation tonometry in 405 participants. We excluded 49 subjects from statistical analysis because of a history of myocardial infarction or ischemic heart disease (n = 14), atrial fibrillation (n = 4) or symptomatic heart failure (n = 1), or because of insufficient quality of the echocardiograms (n = 6) or finger pressure waves (n = 24). In total, we thus analysed data from 356 participants.
Echocardiography and arterial phenotyping were performed after the subject had rested for at least 15 min in supine position. Details on the echocardiography, arterial phenotyping and other measurements are provided as Data Supplement.
Data acquisition - Briefly, one experienced physician (T.K.) did the ultrasound examination using a Vivid E9 (GE Vingmed, Norway) interfaced with a 2.5- to 3.5-MHz phased-array probe, in accordance to recommendations  and previous reports .
During echocardiography, we recorded continuous finger pressure waves at the subject’s right middle-finger using a Finometer Pro (Finapres Medical Systems, The Netherlands), which were converted to brachial pressure curves by a validated transfer function implemented in the Finometer software .
Off-line analysis - One observer (TK) analyzed the echocardiograms using EchoPac software (GE Vingmed). End-diastolic LV dimensions were used to calculate relative wall thickness (RWT) and LV mass. We measured transmitral peak early (E) and late (A) diastolic velocities as well as early (e’) and late (a’) diastolic peak mitral annular velocities at 4 acquisition sites (septal, lateral, inferior and posterior). The E/e’ ratio, a surrogate of LV filling pressure, was transmitral E peak divided by e’ averaged from the 4 acquisition sites. Two observers (N.C and T.K.) measured LV global longitudinal strain (LS) using myocardial speckle-tracking software (Q-analysis, GE Vingmed) as described before . We used peak systolic, mid-wall global LS for statistical analysis. Additional file 1: Figure S1 shows the relative inter-observer variability of individual segmental LS values.
VAC analysis – We post-processed the simultaneously recorded LV deformation curves and pressure waves using a custom Matlab algorithm (The MathWorks, Inc., USA). The performance of the software was validated in 50 subjects using intermediate and final quality checkpoints and is available upon request from the corresponding author (Additional file 1: Figure S2). First, we constructed pressure-strain loops by plotting the calibrated brachial pressure wave against the global LS curve (Fig. 1). Next, we calculated EWD as the area of the pressure-strain loop during LV ejection, representing the cumulative work density on the muscle that instantaneously shortens a given amount (i.e. change in strain) against an instantaneous pressure (Fig. 1) . We averaged EWD from apical 4 and 2 chamber views (Additional file 1: Figure S3). We additionally included the average EWD standardized by RWT in sensitivity analyses to better account for LV wall stress.
We recorded carotid, femoral and radial arterial waveforms using a SPC-301 micromanometer (Millar Instruments Inc., USA) linked to a computer running SphygmoCor software (AtCor Medical Pty. Ltd., Australia). Pulse waves were calibrated by the supine brachial BP measured immediately before tonometry. From radial signals, SphygmoCor software constructed the aortic (central) pulse wave using a validated generalized transfer function. Central pulse pressure (PP) was central systolic minus diastolic pressure. Augmentation pressure (AP) was the pressure difference between the first and second shoulder of the central waveform. In 157 men and 117 women, we measured aortic pulse wave velocity (PWV), the non-invasive gold standard of arterial stiffness, as the carotid-femoral distance divided by the carotid-femoral pressure transit time .
Conventional BP was the average of 5 auscultatory readings obtained with the subject in seated position. Hypertension was defined as a BP of at least 140 mmHg systolic or 90 mmHg diastolic or the use of antihypertensive drugs. Participants on antihypertensive therapy were defined as having either uncontrolled (BP > 140/90) or controlled (BP ≤ 140/90) hypertension. Diabetes mellitus was determined by self-report, a fasting glucose level of at least 126 mg/dL, or the use of antidiabetic agents.
We used SAS software v9.4 (SAS Institute, Cary, NC) for database management and statistical analysis. We compared means and proportions using a large sample z-test and χ2-test, respectively. Statistical significance was a two-sided P value < 0.05. Using mixed models, we assessed multivariable-adjusted associations of EWD with anthropometric and clinical characteristics, hemodynamics, arterial stiffness and echocardiographic indexes of LA volume index (LAVi) and LV structure and function. All models were adjusted for age, sex, heart rate and body height and weight, and accounted for family clusters modelled as a random effect. We reported multivariable-adjusted regression coefficients for EWD per 1-SD increase in the covariable. We repeated regression analyses for men and women separately.
Characteristics of participants
The 356 participants (44.7% women) included 169 (47.5%) hypertensive subjects, of whom 92 (54.4%) were on antihypertensive drug treatment. Mean age was 53.9 ± 13.7 years. Tables 1 and 2 list the clinical, arterial and echocardiographic characteristics of the study cohort by sex. EWD averaged 635.8 ± 179.6 Pa in men and 782.3 ± 224.7 Pa in women (P < 0.0001; Table 2).
Anthropometric and clinical determinants of EWD
EWD increased significantly with age in both unadjusted and multivariable-adjusted analyses (P < 0.0001; Fig. 2a). As shown in Fig. 2b, the pressure-strain loop extended greater along the pressure axis with only little fluctuation in peak global LS in older subjects as compared to younger ones. We observed similar age-related changes in pressure-strain loop and EWD in both men and women (Additional file 1: Figure S4). Of note, EWD was higher in postmenopausal than in premenopausal women (P < 0.0001; Additional file 1: Figure S5).
On average, after adjustment for age, heart rate and body size, the pressure-strain loop area during ejection was in its entirety broader in women than in men (Fig. 2c). Hence, averaged EWD was significantly higher in women as compared to men (P < 0.0001; Fig. 2c). Furthermore, adjusted EWD was significantly greater in participants whose hypertension was either untreated (P = 0.018) or uncontrolled (P = 0.0036) as compared to subjects with controlled hypertension (Fig. 3a). Of note, EWD did not differ between normotensives and subjects with controlled hypertension (P = 0.49; Fig. 3a).
After full adjustment, EWD decreased with the use of β-blockers (− 65.4 Pa; P = 0.046), but was not related to the use of ACE-inhibitors/ARBs, calcium channel blockers or diuretics (P ≥ 0.51). In addition, after accounting for sex and age, EWD was significantly higher in obese than in non-obese women (P ≤ 0.043; Fig. 3b). In both unadjusted and fully-adjusted analysis, EWD was not associated with heart rate (P ≥ 0.23) or with smoking (P ≥ 0.22).
Our findings were consistent for the 4 and 2 chamber EWDs analyzed separately (data not shown).
Associations of EWD with arterial hemodynamics and stiffness
Table 3 presents the overall and sex-specific multivariable-adjusted estimates (95% CI) for EWD associated with a 1-SD increase in hemodynamic and arterial indexes.
As expected, in multivariable-adjusted analyses, higher EWD was independently related to higher systolic (effect size: + 95.2 Pa) and diastolic BP (+ 37.2 Pa; P ≤ 0.0003) as well as to higher brachial PP (+ 82.1 Pa) and mean arterial pressure (MAP; + 69.6 Pa; P < 0.0001). Moreover, after full adjustment, EWD increased with higher AP (+ 59.1 Pa), central PP (+ 65.7 Pa) and PWV (+ 44.8 Pa; P ≤ 0.0006; Table 3; Additional file 1: Figure S6). We also observed that EWD increased stronger with higher brachial PP and MAP in women than in men (Table 3; Pint ≤ 0.043).
We confirmed the associations between EWD and arterial characteristics in a sensitivity analysis excluding subjects on antihypertensive treatment (Additional file 1: Table S1). In addition, EWD standardized by RWT increased with all brachial and central BP components (P ≤ 0.011), but not with PWV (P = 0.13; Additional file 1: Table S2).
Associations between EWD and LV structure and function
We determined the overall and sex-specific multivariable-adjusted estimates (95% CI) for EWD associated with a 1-SD increase in echocardiographic indexes reflecting LA and LV geometry (Table 4) and LV diastolic function (Table 5).
EWD in relation to LA and LV geometry
In all subjects, EWD significantly increased with LAVi (+ 41.2 Pa), end-diastolic volume index (EDVi; + 23.0 Pa) and LV stroke volume (+ 27.9 Pa) after full adjustment (P ≤ 0.039 for all; Table 4). EWD decreased significantly with increased LV posterior wall thickness (− 34.3 Pa) and RWT (− 29.9 Pa) only in men (P ≤ 0.013; Table 4; Fig. 4). In contrast, EWD increased significantly with higher LV mass (+ 46.1 Pa), LAVi (+ 76.2 Pa) and EDVi (+ 43.8 Pa) only in women (P ≤ 0.035; Table 4; Fig. 4).
We confirmed these associations when excluding subjects receiving antihypertensive drug therapy (Additional file 1: Table S3).
Similarly, when standardizing by RWT, EWD increased with higher LAVi, EDVi and end-systolic volume index (ESVi) and stroke volume (P ≤ 0.012; Additional file 1: Table S4).
EWD in relation to LV function
Elevated EWD correlated independently to greater early (+ 55.3 Pa) and late (+ 47.2 Pa) diastolic transmitral peak velocities (P ≤ 0.0005; Table 5). After full adjustment, EWD related directly with the E/e’ ratio (+ 32.1 Pa; P = 0.0062; Table 5). However, in sex-specific analyses, this relationship was only observed in women (Fig. 4). EWD correlated directly with peak global LS in both men and women (+ 65.4 Pa; P < 0.0001).
We confirmed these associations when excluding subjects on antihypertensive drug treatment (Additional file 1: Table S5). Standardized by RWT, EWD increased independently with transmitral velocities, E/A ratio and e’ peak (P ≤ 0.027 for all; Additional file 1: Table S4).
Here, we investigated the anthropometric and clinical determinants of a novel index of LV myocardial performance as derived from non-invasive, simultaneous pressure-strain recordings. In addition, we explored its relation to arterial stiffness and LV function and structure. The key findings of our study were as follows: (i) older age, female sex and hypertensive status are associated with higher EWD; (ii) in both men and women, EWD increased with steady and pulsatile BP components and arterial stiffness; and (iii) EWD decreased independently with LV wall thickness in men, yet increased with higher LAVi and LV filling pressure (E/e’ ratio) in women only.
LV performance is determined by its intrinsic contractility (inotropy), by the tension on the LV myocardium at end-diastole (preload) and by the load against which the LV needs to contract to eject blood (afterload) . Within this context, the work density at which the LV generates a stroke volume against a given afterload might be represented by the area of the LV pressure-strain loop during ejection (i.e. EWD). Indeed, EWD integrates the instantaneous LV deformation against an instantaneous pressure and might thus reflect LV systolic performance [2, 11].
Numerous population studies demonstrated the detrimental impact of a chronically increased afterload on LV structure and function [4,5,6,7]. Indeed, at increased afterload, the heart needs to generate greater force to preserve stroke volume, thereby increasing its energy expenditure and, eventually, promoting cardiac dysfunction and remodeling. Along these lines, an experimental study showed that the peak rate of changes in LV pressure (dP/dt), an invasive index of contractility, was 51% greater in hypertensive than in normotensive rabbits . We previously observed that the higher arterial load in subjects with hypertension matched with higher LV systolic stiffness and enhanced LV myocardial performance . In result, EWD was 24% higher in hypertensive subjects as compared to normotensives . In line with a small case-control study in 74 patients , we confirmed that hypertensive subjects had a higher pressure-strain area than normotensive participants if untreated or if their hypertension was uncontrolled, even after accounting for important confounders such as age, sex and body size. In contrast, EWD was not different between normotensives and effectively controlled hypertensives. Evidently, the observational and cross-sectional nature of our study did not allow to infer true causality between controlled hypertension and normalization of EWD. Moreover, EWD increased with PWV, reflecting aortic stiffness.
Previous population studies demonstrated that LV diastolic function, particularly in women, is sensitive to the detrimental effects of increased central pulse pressure and arterial stiffness [4, 6, 7]. This sex-dependent vulnerability of LV diastolic function to cardiac afterload especially in older women might be explained by the higher aortic pulsatile load and stiffness [5, 7] and the enhanced LV systolic performance [9, 19] in women as compared to men. Indeed, women have on average smaller heart chamber volumes, even after standardization for body size , with higher LV ejection fraction and global LS , and a steeper slope of the stroke work-EDV relationship . Of note, the passive LV diastolic elastance, a major determinant of LV filling pressure, also appears higher in women than in men . In line, we observed that women had higher pulsatile load and LV LS and, therefore, overall higher EWD than men. Moreover, EWD increased significantly with echocardiographic surrogates of LV filling pressure (E/e’ ratio)  and preload indexes like LAVi and EDVi in women only . As such, EWD might reflect the sex-dependent interplay between preload, afterload and LV systolic performance.
Besides undergoing changes in function, the LV progressively remodels in response to chronically elevated LV afterload [4, 22]. Indeed, the LV walls thicken during chronic pressure overload in order to normalize LV wall stress . In fact, the LV tends towards concentric remodeling to cope with the progressive stiffening of the large arteries during life . In our study, yet only observed in men, EWD decreased with LV concentric remodeling (higher RWT).
Recent interest in constructing and evaluating pressure-strain loops supports the need for a non-invasive clinical tool that integrates LV deformation and loading. Indeed, inspired by Russel and colleagues , GE Healthcare recently added a feature to their echocardiographic post-processing software (EchoPAC) to construct pressure-strain loops from 2D LV strain and an estimated LV pressure curve. The LV pressure-strain loop area derived from such estimated loops correlated strongly with invasive measurements . Yet, Hubert et al. observed substantial shortcomings in the estimation of the pressure curve, particularly at greater pressures and loop areas . In contrast, in our study, we utilized the true area of subject-specific pressure waveforms at the specific time of strain recording. As such, our approach does not only account for differences in the shape of pressure waves between individuals, but also for temporal variability in pressure magnitude within individuals. However, validation of our approach to invasive methods is required.
The present study must be interpreted in context of its limitations and strengths. First, echocardiography is prone to measurement errors. However, one experienced observer recorded all echocardiographic images using a standardized protocol. Moreover, images were post-processed by two experienced observers with good reproducibility. Second, we derived EWD by simultaneous finger tonometry and 2D speckle tracking to overcome the technical challenges linked to simultaneous carotid tonometry and Tissue Doppler Imaging . Of notice, EWD values derived from our 356 participants corresponded to those obtained by carotid tonometry and TDI in a random cohort of 148 participants . Finally, EWD requires further validation in outcome and experimental studies.
Older age, female sex and hypertensive status are associated with higher EWD. Integration of the LV pressure-strain loop during ejection might be a useful tool to non-invasively evaluate sex-specific and interdependent effects of preload and afterload on LV myocardial performance. Future studies should investigate the utility of pressure-strain loops in clinical decision-making, particularly in subjects at risk for heart failure.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author (T.K.) upon reasonable request.
End-diastolic volume index
End-systolic volume index
Ejection work density
Flemish Study on Environment, Genes and Health Outcomes
Left atrial volume index
Left ventricular mass index
Pulse wave velocity
Spinale FG. Assessment of cardiac function--basic principles and approaches. Compr Physiol. 2015;5(4):1911–46.
Kuznetsova T, D’hooge J, Kloch-Badelek M, Sakiewicz W, Thijs L, Staessen JA. Impact of hypertension on ventricular-arterial coupling and regional myocardial work at rest and during isometric exercise. J Am Soc Echocardiogr. 2012;25(8):882–90.
Borlaug BA. The pathophysiology of heart failure with preserved ejection fraction. Nat Rev Cardiol. 2014;11(9):507–15.
Cauwenberghs N, Knez J, D’hooge J, Thijs L, Yang W-Y, Wei F-F, et al. Longitudinal changes in LV structure and diastolic function in relation to arterial properties in general population. JACC Cardiovasc Imaging. 2017;10(11):1307–16.
Cauwenberghs N, Knez J, Tikhonoff V, D’hooge J, Kloch-Badelek M, Thijs L, et al. Doppler indexes of left ventricular systolic and diastolic function in relation to the arterial stiffness in a general population. J Hypertens. 2016;34(4):762–71.
Shim CY, Park S, Choi D, Yang WI, Cho IJ, Choi EY, et al. Sex differences in central hemodynamics and their relationship to left ventricular diastolic function. J Am Coll Cardiol. 2011;57(10):1226–33.
Coutinho T, Borlaug BA, Pellikka PA, Turner ST, Kullo IJ. Sex differences in arterial stiffness and ventricular-arterial interactions. J Am Coll Cardiol. 2013;61(1):96–103.
Kawaguchi M, Hay I, Fetics B, Kass DA. Combined ventricular systolic and arterial stiffening in patients with heart failure and preserved ejection fraction: implications for systolic and diastolic reserve limitations. Circulation. 2003;107(5):714–20.
Borlaug BA. Sex, load, and relaxation: are women more susceptible to load-dependent diastolic dysfunction? J Am Coll Cardiol. 2011;57(10):1234–6.
Geyer H, Caracciolo G, Abe H, Wilansky S, Carerj S, Gentile F, et al. Assessment of myocardial mechanics using speckle tracking echocardiography: fundamentals and clinical applications. J Am Soc Echocardiogr. 2010;23(4):351–5.
Urheim S, Rabben SI, Skulstad H, Lyseggen E, Ihlen H, Smiseth OA. Regional myocardial work by strain Doppler echocardiography and LV pressure: a new method for quantifying myocardial function. Am J Physiol Heart Circ Physiol. 2005;288(5):H2375–80.
Hubert A, Le Rolle V, Leclercq C, Galli E, Samset E, Casset C, et al. Estimation of myocardial work from pressure-strain loops analysis: an experimental evaluation. Eur Heart J Cardiovasc Imaging. 2018;19(12):1372–9.
Gottdiener JS, Bednarz J, Devereux RB, Gardin J, Klain A, Manning WJ, et al. American Society of Echocardiography recommendations for use of echocardiography in clinical trials. J Am Soc Echocardiogr. 2004;17(10):1086–119.
Schutte AE, Huisman HW, van Rooyen JM, Malan NT, Schutte R. Validation of the Finometer device for measurement of blood pressure in black women. J Hum Hypertens. 2004;18(2):79–84.
Kuznetsova T, Cauwenberghs N, Knez J, Yang W-Y, Herbots L, D’hooge J, et al. Additive prognostic value of left ventricular systolic dysfunction in a population-based cohort. Circ Cardiovasc Imaging. 2016;9(7).
Vlachopoulos C, Xaplanteris P, Aboyans V, Brodmann M, Cifkova R, Cosentino F, et al. The role of vascular biomarkers for primary and secondary prevention. A position paper from the European Society of Cardiology Working Group on peripheral circulation: endorsed by the Association for Research into arterial structure and physiology. Atherosclerosis. 2015;241(2):507–32.
Aylward PE, McRitchie RJ, Chalmers JP, West MJ. Baroreflex control of myocardial contractility in conscious normotensive and renal hypertensive rabbits. Hypertension. 1983;5(6):916–26.
Chan J, Edwards NFA, Khandheria BK, Shiino K, Sabapathy S, Anderson B, et al. A new approach to assess myocardial work by non-invasive left ventricular pressure-strain relations in hypertension and dilated cardiomyopathy. Eur Heart J Cardiovasc Imaging. 2019;20(1):31–9.
Hayward CS, Kalnins WV, Kelly RP. Gender-related differences in left ventricular chamber function. Cardiovasc Res. 2001;49(2):340–50.
Ommen SR, Nishimura RA, Appleton CP, Miller FA, Oh JK, Redfield MM, et al. Clinical utility of Doppler echocardiography and tissue Doppler imaging in the estimation of left ventricular filling pressures: a comparative simultaneous Doppler-catheterization study. Circulation. 2000;102(15):1788–94.
Slama M. J.L. T. Assessment of cardiac preload and volume responsiveness using echocardiography. In: Vincent JL, editor. Intensive Care Medicine. 1st ed: Springer; 2003. p. 491–8.
Zhang Z, Tendulkar A, Sun K, Saloner DA, Wallace AW, Ge L, et al. Comparison of the young-Laplace law and finite element based calculation of ventricular wall stress: implications for postinfarct and surgical ventricular remodeling. Ann Thorac Surg. 2011;91(1):150–6.
Russell K, Eriksen M, Aaberge L, Wilhelmsen N, Skulstad H, Remme EW, et al. A novel clinical method for quantification of regional left ventricular pressure-strain loop area: a non-invasive index of myocardial work. Eur Heart J. 2012;33(6):724–33.
The Research Unit Hypertension and Cardiovascular Epidemiology (Leuven, Belgium) received grants from the Flanders Research Foundation (FWO grants G.0880.13; 11Z0916N and G0C5319N).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the University of Leuven. All participants provided written informed consent prior to the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.