Background

Echocardiography is used widely to investigate cardiac function in patients on mechanical ventilation due to acute respiratory failure, cardiogenic pulmonary edema, or other critical illnesses, or for general anesthesia [1]. The echocardiographic parameters best suited for the evaluation of left or right ventricular dysfunction remain unclear. This issue has been given high priority in research agendas for critical care echocardiography [1, 2]. Diastolic dysfunction (diagnosed by echocardiography with tissue Doppler and mitral flow pattern analyses) is a frequent cause of weaning failure [3,4,5,6], has been termed “understudied” [7], and is used to guide therapeutic decision making [5, 8]. Questions remain as to whether Doppler-derived parameters truly reflect intrinsic diastolic cardiac properties or represent preload- and afterload-dependent filling phenomena [9, 10]. The positive intrathoracic pressure applied by mechanical ventilation leads to an overall preload reduction and an increase in right ventricular and decrease in left ventricular afterload [11,12,13,14]. Positive end-expiratory pressure (PEEP) influences the echocardiographic assessment of diastolic function in anesthetized and critically ill patients [15,16,17], questioning the validity of echo-Doppler assessments under mechanical ventilation. For the investigation of diastolic dysfunction independently of cardiac loading conditions, invasive pressure–volume-loop analysis with conductance catheters is the gold standard. Invasive animal studies have provided conflicting results [18, 19], and we are not aware of any reporting of invasively measured data from humans under mechanical ventilation. As the pleural space is the effective working environment of the heart, its transmural pressure may be approximated by subtracting the esophageal pressure from the intracavitary pressure [11, 20, 21].

The widespread use of positive pressure ventilation in anesthesia and critical care and the load dependency of Doppler-derived parameters for diastolic function underscore the need for the validation of echocardiographic diastology parameters measured under mechanical ventilation. With this study, we aimed to elucidate intrinsic ventricular properties under PEEP and to obtain an assembled picture of diastolic left ventricular function together with Doppler-derived filling parameters. We hypothesized that mitral inflow and tissue Doppler parameters would reflect altered cardiac loading conditions attributable to PEEP, rather than changes in myocardial properties.

Materials and methods

Participants

This prospective, controlled single-center observational study conducted at the University Hospital Bern was approved by the Ethics Committee of the Canton of Bern (KEK 104/14). Patients provided written informed consent. The report follows the STROBE guidelines. Adult patients with no known heart, pulmonary, renal, or esophageal disease scheduled for elective coronary angiography between June 2015 and January 2018 were included prospectively, unless diastolic dysfunction > grade 1 was present on transthoracic echocardiography.

Protocol

A pulmonary-artery balloon catheter was placed in the right pulmonary artery. Cardiac output was measured by thermodilution and the Fick method [22]. Left-ventricular pressure volume loops were obtained with a 7 French combined pressure-conductance catheter system (Sentron Europe BV, Roden, Netherlands and INCA; CDLeycom, Hengelo, Netherlands), calibrated as described previously [23]. An Amplatzer sizing balloon (Abbott Medical, The corporate village, Zaventeem, Belgium) was placed in the inferior vena cava for intermittent reduction of the preload (Fig. 1). Comprehensive transthoracic echocardiography (Vivid E9 and E95; GE Medical Systems, Glattbrugg, Switzerland) was performed during the pressure–volume-loop assessment. Data were stored for blinded offline analysis (EchoPac software, version 08; GE Medical Systems).

Fig. 1
figure 1

Exemplary loop family from a caval occlusion maneuver with the INCA device, recorded at a PEEP level of 10 cmH2O. The loop family corresponds to the echocardiography shown in Fig. 5

The assessment was performed with the participants breathing spontaneously, and with continuous positive airway pressure (CPAP) applied via a tight face mask; this mode of ventilation maintains continuous PEEP. Pressures of 5 and 10 cmH2O were applied in randomized order (Servo-i; Maquet, Solna, Sweden). An air-filled esophageal balloon (Nutrivent; Sidam s.r.l., Mirandola, Italy) [24, 25] was placed behind the left atrium under fluoroscopic guidance, with occlusion testing to confirm proper positioning [26, 27]. This balloon served as a surrogate for pleural pressure, which is needed for the calculation of transmural cardiac pressures. Pulmonary-artery and left-ventricular blood gases were taken at each stage of the assessment (spontaneous breathing, 5 cmH2O PEEP, and 10 cmH2O PEEP), together with noninvasive blood pressure, pulmonary artery pressure, and left-ventricular pressure–volume loops during normal respiration and expiratory holds in parallel with the echocardiographic assessment. At each stage, the cardiac preload was reduced temporarily by occlusion of the inferior vena cava to obtain pressure–volume loop families for end-systolic and end-diastolic pressure–volume relationships (ESPVR and EDPVR, respectively, Fig. 1). All invasive measurements were taken during expiratory hold to avoid additional pressure swings from respiration. Additional details are provided in Additional file 1.

Data collection, processing, and analysis

Pressure data were entered into the dedicated acquisition system with standard pressure transducers, and subsequently extracted to Matlab (Release 2018b; MathWorks, Natick, Massachusetts, USA). The isovolumetric relaxation time constant (τ) was recalculated from heartbeats during expiratory hold using exponential (τGlantz) and logistic (τlogistic) models for intraventricular and transmural ventricular pressures by subtraction of the esophageal pressure [28]. Diastole lengths were resampled to a standardized diastole length. For the intraventricular and transmural pressures, end-diastolic pressure–volume relationships were calculated using an exponential curve (P = CeβV) fit to the diastolic pressure–volume points to determine the chamber stiffness constant (β) [29] and the pressure at an end-diastolic volume of 0 ml (= C). V equals chamber volume. Additional methodological details are provided in Additional file 1.

Statistical analysis

The normality of the data distributions was assessed using the Shapiro–Wilk test and QQ plots. Comparisons among groups were made by one-way repeated-measures or Friedman’s analysis of variance, as appropriate, with post-hoc Bonferroni correction. ESPVRs were assessed using linear mixed-effects regression, with systolic pressure serving as the dependent variable. EDPVRs were assessed using generalized linear mixed models with a log link and gamma distribution for diastolic pressure. In all models, the volume, PEEP group, and interaction between them served as independent variables. Additionally, a random intercept and slope for volume were introduced at the participant level to allow for heterogeneity in the pressure–volume relationship among participants. For better interpretability, systolic and diastolic volumes were centered at their mean values. The analyses were performed with SPSS (version 21; IBM Corporation, Armonk, NY, USA) and Stata (version 16; StataCorp, College Station, TX, USA) software, and two-tailed p values < 0.05 were considered to be significant.

Results

Sample characteristics

30 cardiac patients (aged 18–60 years) consented to study participation. Five could not undergo assessment due to capacity constraints in the catheter laboratory. The analyses were performed with data from 25 patients (7 women, 18 men, Table 1). The sample for steady-state assessment comprised 1298 pressure–volume loops. After the exclusion of occlusion maneuver artifacts, 1598 and 1052 loops were entered into linear and generalized linear mixed-effects models for the assessment of ESPVRs and EDPVRs, respectively. The study procedures were performed without complication.

Table 1 Characteristics of the study population

The thermodilution cardiac output decreased slightly with increasing PEEP at stable stroke volumes due to a slower heart rate with higher PEEPs (Table 2). The pulmonary artery occlusion pressure and the esophageal pressure (as a surrogate of pleural pressure) increased with PEEP. Echocardiography showed a stable mean ejection fraction of 65% ± 5% with unchanging chamber volumes (Table 2).

Table 2 Baseline hemodynamic data

Invasive assessment and modeling properties

Systolic function and ventriculo-aortic coupling

The ejection fraction, end-systolic pressure, the pressure change over time (dP/dtmax, reflecting passive filling) and ESPVR slope (reflecting end-systolic elastance) did not alter with the PEEP. The aortic elastance did not change with increasing PEEP. The ratio of end-systolic to aortic elastance, an indicator of ventriculo-aortic coupling, remained unchanged with changes in PEEP, as did the stroke work and preload recruitable stroke work (Table 3). The linear mixed model revealed small, but significant, changes in the ESVPR and transmural ESPVR slopes (except for the transmural ESPVR at 5 cmH2O PEEP; Table 4). The changes in the intercepts (except for the transmural ESPVR at 10 cmH2O PEEP) were also small but significant.

Table 3 Systolic impedance catheter data and ventriculo-aortic coupling
Table 4 Linear mixed-effect model with interactions for ESPVR

Relaxation

The application of PEEP significantly increased τ values for intracavitary ventricular pressures, indicating slower early relaxation of the left ventricle. This effect could be observed independently of the decay model used, and was sustained for the transmural cardiac pressures (Fig. 2, Table 5). When τ was corrected for the diastole length, the effect disappeared for the transmural pressure decay.

Fig. 2
figure 2

Mean left (a) and left transmural (b) ventricular pressures following an exponential decay model during isovolumetric relaxation, with 95% confidence intervals as dashed lines, resampled for 35 samples (140 ms). All data from all patients are included. PEEP, positive end-expiratory pressure

Table 5 Diastolic impedance catheter and modeled pressure decay

Passive filling properties

Early and mid-diastolic ventricular pressures and simultaneous pulmonary-artery occlusion pressures were significantly higher with PEEP (5 and 10 cmH2O) than without PEEP, with no difference between PEEP levels (Table 6, Fig. 3a). This effect disappeared when transmural pressures were calculated (Table 6, Fig. 2b). The instantaneous transmural trans-mitral pressure gradient remained unchanged with PEEP (Fig. 2).

Table 6 Diastolic left heart pressure properties
Fig. 3
figure 3

a Left ventricular filling curves represented by diastolic ventricular pressure and PAOP, with 95% confidence intervals as dashed lines. The data were resampled for a filling period of 125 samples (500 ms). The dotted lines for early, mid- and late diastolic data points indicate the data presented in Table 4. b Transmural left ventricular filling curves Means with 95% confidence intervals as dashed lines are presented. The data were resampled to a filling period of 125 samples (500 ms). The dotted lines for early, mid- and late diastolic data points. PAOP pulmonary-artery occlusal pressure, LVP left ventricular pressure, PEEP positive end-expiratory pressure

Positions of the EDPVR and chamber stiffness

The generalized linear mixed-effects model revealed small, but significant, changes in the natural logarithms of the intercept changes in the absolute and transmural EDPVRs. The resulting intercept pressures from the inverse function (for intracavitary pressures: no PEEP, 9.57 mmHg; 5 cmH2O PEEP, 10.31 mmHg; 10 cmH2O PEEP, 12.57 mmHg; for transmural pressures, no PEEP 4.73 mmHg, 5 cmH2O PEEP, 3.82 mmHg; 10 cmH2O PEEP, 3.68 mmHg) indicate a downward shift of the EDPVR for the transmural pressure, with a resolution of between-PEEP differences. The change in slope, representing the chamber stiffness, lost significance when corrected for transmural pressure (Table 7). The shift in position and change in stiffness with increasing PEEP disappeared when transmural pressures were considered (Fig. 4).

Table 7 Generalized linear mixed-effects model for EDPVR
Fig. 4
figure 4

a End-systolic pressure–volume relationships (ESPVRs) for 0, 5, and 10 cmH2O PEEP, according to the mixed linear model equations in Table 4. The equation values are reported at the mean centered volume (50.4 ml), indicated by the dotted line. b End-diastolic pressure–volume relationships (EDPVRs), according to the mixed-effect generalized linear model equations in Table 7. The equation values are reported at the mean centered volume (136.2 ml), indicated by the dotted line. PEEP positive end-expiratory pressure

Echocardiographic diastole evaluation

Trans-mitral flow velocities

The peak early and late trans-mitral flow velocities (E and A) decreased significantly with increasing PEEP, with no change in the E/A ratio. The duration of the E wave increased in absolute terms, but remained constant relative to the duration of diastole and the full cardiac cycle. The duration of the A wave and the deceleration time remained constant (Table 5, Fig. 5).

Fig. 5
figure 5

Exemplary echocardiography tracings at a PEEP level of 10 cm H2O. The frames that correspond to the PV-Loop family in Fig. 1 are shown. The left panel shows a transmitral pulsed wave Doppler trace, the right panel a tissue Doppler tracing

Tissue Doppler properties of the mitral annulus

The peak early septal tissue velocity (Eʹ) decreased with PEEP. The lateral Eʹ and late tissue velocities (Aʹ) were less affected. The isovolumetric relaxation time increased with PEEP (Table 8).

Table 8 Echocardiographic diastology parameters

Discussion

The main finding of this study is that the increase in intrathoracic pressure caused by mechanical ventilation had no intrinsic influence on cardiac function, assessed invasively by the gold standard (impedance catheter). Inotropy, represented by the ESPVR, increased slightly with PEEP. The invasive assessment of diastolic function revealed a small influence of PEEP on active relaxation, with an apparent increase in τ, upward shift of the EDPVR, and decreasing chamber stiffness, indicating altered passive ventricular filling and elastic properties, but only when the pressures were referenced to atmosphere [11]. These changes in invasive measures disappeared with the measurement of transmural chamber pressures, which reflect the working condition of the left ventricle. These results stand in contrast to echocardiographic parameters, as modest PEEP influences the transmitral flow pattern and annular tissue Doppler properties.

These findings challenge current paradigms and must be placed in clinical and physiological contexts. The inaccuracy of transmitral indices for the assessment of diastolic dysfunction has long been discussed [9] 30. Comparative Doppler-conductance catheter studies with non-ventilated patients have shown that tissue Doppler indices may better reflect diastolic dysfunction than do transmitral flow velocities, and that increasing E/E’ ratios are the best indicators of diastolic dysfunction [30, 31]. Our data, support this notion. The load dependencies of transmitral E and A waves are recognized [17], albeit not in comparison with invasive assessment under positive-pressure ventilation. Juhl-Ohlsen and colleagues . [15, 16] demonstrated the influence of PEEP on E and A waves in anesthetized cardiac patients. PEEP increases the pleural pressure [20], explaining its well-established preload-reducing effect [11, 21, 25, 32]. Our echocardiographic findings are in line with these previous observations. The E/Eʹ ratio, strongly predictive of ventilator weaning failure [3], was not affected by PEEP in our population.

We observed apparent slowing of active relaxation, with increases in τ with the application of PEEP. Despite its significance in the intracavitary pressure assessment, the magnitude of this change appears to be clinically marginal. For its proper interpretation, two physiological phenomena need to be taken into account. The first is the heart rate, and thus the duration of diastole. Several studies have shown that τ declines with increasing heart rate and vice versa [33,34,35]. As the heart rate changes with PEEP, we normalized the diastole duration in a second model of τ to exclude influences of changing diastole length in this study. The second issue is the exposure of the heart to pleural, rather than atmospheric, pressure [11]. In routine cardiac catheterization and hemodynamic monitoring, intracavitary and intravascular pressures are measured with a zero reference to the atmosphere. The pericardial pressure, however, approximates the pleural pressure [36]. Cardiac working conditions are better represented by transmural rather than intracavitary pressures. The subtraction of the esophageal pressure from intracavitary pressures may serve as a sufficient approximation [11, 20, 36,37,38]. When this intrathoracic (i.e., pleural) pressure reference is taken into account, the apparent effect of PEEP on τ disappears. Thus, findings may reflect the use of an inappropriate pressure reference, rather than the occurrence of a true physiological phenomenon. The consistency of the findings obtained with various modeling approaches for τ reflects reliability, particularly as the logistic modeling of pressure decay is relatively resistant to respiration influences [28]. In conductance volumetry, unchanged ventricular volumes may exclude ventricular interdependence as a cause of impaired relaxation [39].

PEEP applies a mechanical constraint to the heart by increasing the pleural pressure and elevating the functional residual capacity of the lung, with compression of the cardiac fossa [40]. Increased pericardial pressure shifts the pressure–volume curve upward [41], and a similar shift in the EDPVR slope was observed in this study. This upward shift, and the apparent decrease in chamber stiffness, disappeared when transmural pressures were considered. These findings may thus be considered artifacts of the use of the atmospheric pressure reference. The clinically modest PEEP levels used in this study had no apparent effect on cardiac chamber volumes; shifts may differ with the application of greater pressure, particularly with invasive mechanical ventilation [42]. We may only speculate on the effects of higher PEEP, as the preload reducing effect together with mechanical constraint may dominate afterload reduction and resemble restrictive filling with signs of obstructive shock [11]. Such data would need to be gained from invasive ventilation because higher pressure levels on NIV are usually badly tolerated.

The slope of the ESPVR at 10 cm H2O PEEP increased, with a positive (intracavitary) or no (transmural) relevant increase in the pressure intercept, in this study. These findings are consistent with an increased inotropic state of the left ventricle [29]. At 5 cmH2O PEEP, the increase in slope was less prominent and of borderline significance. We have no mechanistically convincing explanation for this unexpected result. Data from animal studies suggest that PEEP negatively affects the coronary blood flow [43], which would have the opposite effect as we observed. The Anrep effect, or slow force response to increased stretch, could explain the increased inotropy [44, 45], but the afterload seems to be constant in our context with stable ventriculo-arterial coupling across PEEP levels. Verification of the direct positive inotropic effects of PEEP would provide an additional mechanistic explanation of the beneficial effects of positive pressure ventilation for acute cardiac failure, beyond the classic concept of preload and afterload reduction.

Limitations

Several limitations of this study must be considered. First, the study participants were cardiac patients with no apparent diastolic dysfunction at the time of enrollment, as verified by echocardiographic grading and low τ values on spontaneous breathing. Whether similar results would be obtained in patients with such dysfunction or those under invasive mechanical ventilation in an intensive care unit setting warrants further investigation. The advantage of conducting the study with patients without diastolic dysfunction is that we could demonstrate the appearance of such dysfunction with echocardiography. We can only speculate about effects of PEEP in patients with preexisting diastolic function. These effects may be less unidirectional, since the mechanical constraint caused by PEEP may be counterbalanced by its afterload reducing effect in the clinical context of heart failure with preserved ejection fraction or obstructive sleep apnea. Clinical studies indicate long term amelioration of diastolic dysfunction with positive pressure ventilation [46, 47].

Second, our echocardiographic finding of increasing diastolic dysfunction severity with increasing airway pressure may be clinically elusive. Its importance lies not in its magnitude, but in its difference from the conductance-catheter gold standard. Third, voluntary control of respiration in an awake person under an invasive hemodynamic study is demanding. We took meticulous care to instruct the participants in performing expiratory holds and monitored their respiratory drives with esophageal pressure swings. An attending intensivist supervised all maneuvers, and the esophageal pressure swings were screened for inadvertent Valsalva maneuvers. Still, unintended respiratory variation cannot be excluded. Thus, the use of τ logistic, which is least influenced by respiration [28], was of particular value. Fourth, an esophageal balloon is the only feasible pleural pressure surrogate in an awake patient and its use is recommended [38], although this technique may not directly reflect pericardial pressure. We have experience in esophageal balloon use [20, 25], and calibrated the balloon using an accepted technique in this study [37, 38]. As each patient served as his/her own control in the randomized crossover setting, the esophageal pressure was a valid relative reference, although absolute values might have been slightly inaccurate. Fifth, calibration of the PV loops required a two-step approach. The stroke volume was calibrated via thermodilution at every PEEP level. We followed a standardized approach for thermodilution by averaging the three closest cardiac output measurements from five saline injections [48]. The PV loops were positioned on the volume axis by calculating the end-diastolic volume from a baseline ejection fraction measurement (obtained by echocardiography or with hypertonic saline injection). The variability of the PV loop position through the ejection fraction may have contributed to shifts in intercept and slope, but was partially accounted for with the use of mixed-effects models for the ESPVR and EDPVR. Last, our sample size could not be based on a power calculation, because of the exploratory nature of this study and the lack of preexisting data.

Conclusion

We conclude that dynamic echocardiographic filling parameters reflect changing loading conditions, rather than diastolic function, under the application of positive airway pressure. Furthermore, invasive assessments should be referenced to intrathoracic pressure conditions, rather than to atmospheric pressure. The role of echocardiography in the detection of diastolic dysfunction in mechanically ventilated patients needs to be reassessed.