(Circulation. 2005;112:1771-1779.)
© 2005 American Heart Association, Inc.
Imaging |
From the Department of Cardiology (R.Y., J.B., C.C., C.A., M.M., M.A.G.-F.), Unit of Experimental Medicine and Surgery (M.M.D.), and Department of Cardiovascular Surgery (H.R.-A.), Hospital General Universitario Gregorio Marañón, Madrid; Department of Mathematical Physics and Fluids, Facultad de Ciencias, Universidad Nacional de Educación a Distancia, Distancia (J.C.A.); and Department of Signal Theory and Communications, Universidad Carlos III de Madrid, Madrid (J.L.R.-A.), Spain.
Correspondence to Javier Bermejo, MD, PhD, Laboratory of Echocardiography, Department of Cardiology, Hospital General Universitario Gregorio Marañón, Dr Esquerdo 46, 28007 Madrid, Spain. E-mail javbermejo{at}jet.es
Received June 16, 2004; revision received February 10, 2005; accepted June 1, 2005.
| Abstract |
|---|
|
|
|---|
Methods and Results CDMM images and pressure-volume (conductance) signals were simultaneously acquired in 9 minipigs undergoing pharmacological interventions and acute ischemia. Inertial, convective, and total EIVPD curves were calculated from CDMM recordings. Peak EIVPD closely correlated with indices of systolic function based on the pressure-volume relationship: peak elastance (within-animal R=0.98; between-animals R=0.99), preload recruitable stroke work (within-animal R=0.81; between-animals R=0.86), and peak of the first derivative of pressure corrected for end-diastolic volume (within-animal R=0.88; between-animals R=0.91). The correlation of peak inertial EIVPD with these indices was also high (all R>0.75). Load dependence of EIVPDs was studied in another 5 animals in which consecutive beats obtained during load manipulation were analyzed. During caval occlusion (40% EDV reduction), dP/dtmax, ejection fraction, and stroke volume significantly changed, whereas peak EIVPD remained constant. Aortic occlusion (40% peak LV pressure increase) significantly modified dP/dtmax, ejection fraction, and stroke volume; a nearly significant trend toward decreasing peak EIVPD was observed (P=0.06), whereas inertial EIVPD was unchanged (P=0.6). EIVPD beat-to-beat and interobserver variabilities were 2±12% and 5±11%, respectively.
Conclusions Doppler-derived EIVPDs provide quantitative, reproducible, and relatively load-independent indices of global systolic chamber function that correlate closely with currently available reference methods.
Key Words: hemodynamics systole echocardiography imaging pressure
| Introduction |
|---|
|
|
|---|
See p 1684
Intraventricular pressure gradients were first measured with micromanometers >50 years ago.4 Generated by the active force of the contracting ventricle,5,6 the instantaneous ejection intraventricular pressure difference (EIVPD) between the LV apex and the outflow tract (LVOT) has shown to be related to the inotropic state.79 In the absence of outflow obstruction, the EIVPD reaches its peak early during systole and is fundamentally caused by impulsive forces. It is worth noting that the diagnostic value of the ventricular impulse was the basis for developing well-established indices of systolic function such as the first time derivative of pressure (dP/dtmax) and flow acceleration.10 Theoretical evidence further supports that peak EIVPD may be a sensitive index of systolic chamber function.11
We have recently developed and validated a new method to obtain noninvasive measurements of EIVPDs by postprocessing color Doppler M-mode (CDMM) echocardiograms.12,13 The present study was designed to study the utility of EIVPD as a noninvasive index of systolic chamber function in vivo with a high-fidelity pressure and conductance instrumentation animal setup. The specific objectives were (1) to correlate Doppler-derived EIVPDs with invasive reference indices of systolic function in a range of inotropic states induced by pharmacological interventions and acute ischemic cardiomyopathy (study 1), (2) to assess the modification of EIVPDs caused by acute load changes (study 2), and (3) to assess the potential error caused by scan line misalignment on the estimation of EIVPDs (study 3).
| Methods |
|---|
|
|
|---|
Animals were preanesthetized with ketamine and xylazine and mechanically ventilated. Complete anesthesia and relaxation were maintained by propofol infusion (0.2 mg · kg1 · min1) and by repetitive boluses of pentobarbital (15 mg/kg IV plus 5 mg · kg1 · 15 min1) and pancuronium (0.2 mg · kg1 · 15 min1). Through the right carotid artery, a 5F pigtail 12-pole multielectrode combination conductance-pressure catheter (Millar Instruments or CD-Leycom, Zoetermeer) was placed in the LV and connected to a dual-field conductance processor (Sigma 5DF, CD-Leycom). Adequate position of the electrodes was confirmed by examination of the segmental volume signals. The volume signal was calibrated with high-frequency harmonic 2D echocardiographic images (biplane Simpsons method).16
Study 1: Comparison With Invasive Indices of Systolic Chamber Function
Nine pigs (54±10 kg; range, 35 to 65 kg) underwent median sternotomy without the pericardium being opened, and the heart was cradled. A snare was placed around the inferior vena cava for preload manipulation. Animals were studied at baseline, during dobutamine (n=9; 1 to 10 µg · kg1 · min1) and esmolol (n=4; 25 to 200 µg · kg1 · min1) infusions, and after left main coronary microembolization of polystyrene microspheres (n=4; mean diameter, 45 µm; Polysciences).17 Instead of using fixed pre-established doses, we titrated perfusions to obtain a broad range of maximal elastance (Emax) measured online during each experiment.
In each state, simultaneous CDMM images and signals were acquired during end-expiratory apnea, followed immediately by transient caval occlusion while invasive data storage was maintained. This acquisition process was repeated 3 times for each state, with periods of waiting for stabilization of >5 minutes. The average of Doppler-derived EIVPD measurements of the 3 consecutive beats recorded just before caval occlusion was compared with the pressure-volume relationship indices obtained during the preload reduction maneuver performed immediately thereafter.
Study 2: Load Dependence
In 5 pigs, a small subxyphoid aperture was created in the abdominal wall to introduce the ultrasound probe, opening neither the peritoneum nor the thorax. This approach provides a suitable 5-chamber view in these animals. Occlusion balloons were placed in the inferior vena cava (right atrium junction) and in the descending thoracic aorta via the femoral vessels. Because the inotropic state may modify load dependence,18 animals were studied at baseline and during infusion of dobutamine and esmolol (same dose range as study 1). For each inotropic steady state, CDMM and invasive data were simultaneously acquired immediately before and during 10 seconds of caval and aortic occlusion. Within each hemodynamic run, measurements of each individual beat recorded during balloon occlusion were compared with the average of 3 consecutive beats recorded at baseline. Data were obtained at end-expiration apnea without the ultrasound probe or the Doppler cursor being moved.
Echocardiographic Image Acquisition and Analysis
Broadband 2.0- to 4.0-MHz and 2.5- to 5-MHz transducers were used on a Sequoia C-256 system (Siemens AG). CDMM images were obtained from 5-chamber views as described.13 The method used for digital image processing to obtain EIVPDs has been reported12 and validated elsewhere.13 Other authors have also demonstrated the accuracy of a similar approach to measure intraventricular pressure differences in diastole and its applicability in the clinical setting.19,20 Briefly, if the M-mode cursor closely approximates a flow streamline, the spatiotemporal velocity distribution of a discrete blood sample is provided by the value of its corresponding pixel color: v(s,t), where v represents velocity, s represents the linear dimension of the streamline, and t is time. Thus, the color-Doppler M-mode recording provides the data necessary to solve Eulers momentum equation: (
p/
s)=
· [(
v/
t)+v(
v/
s)], where p is pressure and
is blood density. The first and second terms on the right side of the equation account for inertial and convective acceleration, respectively. Once pressure gradient maps are obtained, total, inertial, and convective EIVPDs are calculated by spatial integration between the apex and the LVOT (Figure 1). Instead of a fixed distance between locations, these 2 positions are traced in each image based on the grayscale layer and the pressure-gradient overlay. Peak and time to peak values of each of the 3 EIVPD curves were measured constrained to the ejection period, and all temporal values are reported with the QRS onset as reference. Intraobserver, interobserver, and beat-to-beat variabilities of EIVPD measurements (20 unselected patients referred for a conventional echocardiographic examination; independently and blindly acquired CDMM recordings) were 0.1±0.5 mm Hg (2±9%), 0.2±0.6 mm Hg (5±11%), and 0.1±0.7 mm Hg (2±12%), respectively.
|
Invasive Hemodynamic Measurements
All signals were digitized at 1000 Hz. By cross correlation of a synchronicity signal connected to the ultrasound scanner and the signal acquisition system, CDMM and invasive hemodynamic measurements were obtained for the same beat (Figure 1).13
Three indices of systolic function were derived from the pressure-volume loops obtained during caval occlusion: (1) Emax, defined as the slope of the end-systolic pressure-volume relationship;3 (2) preload recruitable stroke work (PRSW), defined as the slope of the relationship between stroke work and end-diastolic volume (EDV);21 and (3) the slope of the relationship between dP/dtmax and EDV (dP/dtmax/EDV).22 The iterative linear regression method was used to calculate Emax.23 Beats were selected after a 3- to 4-mm Hg decrease in end-systolic LV pressure, and a mean of 10 subsequent cycles was processed (6.5±3.0 seconds) to minimize the effect of autonomic reflexes. The end-systolic pressure-volume relationship was linear at all hemodynamic states (all R>0.91). Intraobserver, interobserver (10 randomly selected runs), and consecutive run-to-run (32 hemodynamic states) variabilities of Emax were 0.11±0.49, 0.23±0.51, and 0.14±0.37 mm Hg/mL, respectively.
Statistical Analysis
Linear mixed-effects (LME; S-Plus version 2000) models were used for analysis because repeated-measures designs were heavily unbalanced. Best linear unbiased estimators are reported for these models. Best linear unbiased estimators account for the mean expected values for fixed effects once the source of variation resulting from random effects is omitted. Significant LME models were followed by simulation contrasts.24 For study 1, EIVPD was compared with Emax, PRSW, and dP/dtmax/EDV (fixed effects). For this purpose, 2 sets of models were calculated using a random animal effect, allowing either for different intercepts and slopes or only for different intercepts of regression lines within each animal. Within-animal correlation coefficients25 and between-animals weighted regression models26 were also computed. Pooled data obtained at baseline and during pharmacological interventions were categorized in 3 levels of systolic function based on Emax tertiles. Hemodynamic differences between categories (fixed effects) were then assessed by LME models (random animal effect) and all-pair contrasts. For study 2, the relative change from preocclusion values of EDV and peak systolic pressure was calculated for all occlusion beats and categorized in 5% intervals. Then, LME models were used to study the effects of the load percentage-change category (fixed effect) on the hemodynamic variables, nesting the hemodynamic run within the animal random effect. Simulation contrasts were calculated against baseline values. A value of P<0.05 was considered significant.
Scan Line Misalignment Error Analysis (Study 3)
The method used to estimate EIVPDs assumes that flow follows a linear streamline and that this streamline is coaxially interrogated to obtain the CDMM image. These assumptions may not always be met in clinical practice and can lead to significant bias, particularly because outflow is known to be skewed and axis asymmetric. Therefore, the potential error related to streamline linearization and nonoptimal scan line positioning was addressed in vivo. Using a biplane scanner (Philips Sonos 7500) and a matrix transducer, we acquired 36 color Doppler cine loops of LV outflow at end expiration every 5° in a healthy volunteer (see the expanded Methods section in the Data Supplement, found online at http://circ.ahajournals.org/cgi/content/full/CIRCULATIONAHA.104.485128/DC1). Loops were postprocessed to render a full 5D (3D+time+velocity) outflow data set.27 The core position of the jet was calculated for all flow cross sections obtained in the long axis from the apex to the LVOT at peak flow rate.28 The ideal scan line orientation was estimated as the biplane linear projection of this core position. Against the reference EIVPD calculated from this line, errors related to angle and parallel displacement were measured every 2.5° and 0.625 mm, respectively. Errors are reported as the mean±SD of the 360° spatial rotation at each misaligned position.
| Results |
|---|
|
|
|---|
|
|
|
|
|
Load Dependence of EIVPDs
Caval occlusion remarkably decreased EDV and end-diastolic LV pressure without increasing heart rate (Table 3). Stroke volume, ejection fraction, and dP/dtmax significantly changed during preload manipulation. However, peak EIVPD remained constant. When EDV decreased >10%, peak inertial EIVPD increased by 15% to 23%. During aortic occlusion, peak LV pressure increased rapidly up to a relative change of 40%. Ejection fraction and stroke volume decreased, whereas dP/dtmax increased significantly (Table 4). Although nonsignificant, a trend toward a reduction of peak EIVPD at high afterload (P=0.06) was observed. Peak inertial EIVPD was unchanged by afterload interventions.
|
|
Scan Line Misalignment Error
The ejection jet core followed a predominantly linear trajectory in the frontal and sagittal planes of the LVOT (linear regression coefficients of their respective projections >0.9). Consequently, the error related to the linear trajectory assumption was only 0.15 mm Hg. The error in EIVPD resulting from scan line misalignment is represented in Figure 5. As shown, an angle misalignment of ±35° (in any spatial direction) was associated with an error of 0.2 to 0.6 mm Hg (95% confidence limits of agreement), whereas a ±0.5-cm coaxial displacement from centerline (representing a 50% uncertainty of the LVOT diameter size) induced an error of 0.3 to 0.7 mm Hg.
|
| Discussion |
|---|
|
|
|---|
Intraventricular Pressure Gradients and Systolic Function
Regional pressure gradients inside the LV take place during ejection in normal hearts. EIVPDs, the consequence of the resistance of intraventricular blood to inertial and convective acceleration, represent the force applied by the contracting ventricle per unit volume of accelerated blood.5,6 The magnitude of EIVPDs has been shown to increase with exercise8 and adrenergic stimulation9 and to decrease with ß-blocking agents in invasive studies.7 Furthermore, peak EIVPD has shown a greater sensitivity and a more linear pattern than dP/dtmax for changes in the inotropic state in a computer fluid-dynamic simulation study.11 Applying the same noninvasive method used in the present study, we have previously demonstrated that peak EIVPD is depressed in patients with dilated cardiomyopathy and is highly sensitive to inotropic drug interventions.13
Load Dependence
The nature of load dependence in vivo is often complex and nonlinear; thus, several indices can be significantly affected by preload and afterload yet can be relatively independent of load over a specified range.15,18,31 Even methods based on the pressure-volume relationship can be slightly influenced by load in the intact circulation.32 In our study, peak EIVPD showed to be relatively load independent. Importantly, the magnitude of the load interventions performed in our study was higher than those explored in previous studies assessing the stability of Doppler-derived indices of systolic function.30,33
A nearly significant trend toward partial afterload dependence of total EIVPD was found. Because of the competitiveness between intrinsic and extrinsic components of systolic load,5 the intraventricular pressure drop increases when aortic pressure is reduced and vice versa. This interaction between afterload components is the manifestation of the inverse force-velocity relationship of contraction and has been demonstrated in situations of predominantly convective gradients such as aortic valve stenosis.34 However, without outflow obstruction, EIVPDs are mediated mostly by inertial forces5,13 and behave differently; as shown in the present study, inertial EIVPD is unaffected by afterload (Table 4). The fact that inertial EIVPD reaches its peak very early during ejection (Table 1) probably explains why this index is less afterload dependent than other ejection phase indices. Furthermore, the close correlation observed between inertial EIVPD and Emax conforms to the well-known relationship between flow acceleration and inotropic state.11,33
Clinical Implications
Ejection fraction is the pivotal index of systolic function in clinical practice. Although the clinical utility of ejection fraction is well established, this index may be unreliable because of inaccuracy in the estimation of LV volumes and, most importantly, its load dependence.31 Ejection fraction is particularly modified by load when contractility is impaired,18 an issue that limits its applicability under extreme or changing load conditions. In these situations, Doppler-derived EIVPDs could provide a more reliable measurement of systolic chamber function.
Because the Doppler signal is less sensitive to noise than cross-sectional echocardiography, EIVPDs are obtainable even in patients in whom poor endocardial definition precludes an accurate estimation of LV volumes. The present study demonstrates that the measurements of EIVPDs are relatively accurate as long as scan line misalignment is <30° and displacement from the LVOT center is <0.5 mm. Because CDMM acquisition is always guided by 2D echocardiography, these constraints generally can be met during clinical examination, as suggested by our reproducibility results.
Study Limitations
The present study sought to characterize LV systolic function through the analysis of chamber pump parameters. Assessing myocardial function is a different issue that involves the analysis of aspects such as fiber strain and sarcomere length. To correlate indices of chamber function with true myocardial function, it needs to be assumed that instantaneous modifications of chamber volume directly follow changes in fiber length. However, in the presence of abnormal ventricular geometry, this may not always be the case. In fact, chamber dynamics are known to overestimate myocardial function when LV wall thickness is increased.1 Therefore, the performance of EIVPDs in the presence of abnormal chamber geometry resulting from hypertrophy or regional wall motion abnormalities deserves further investigation. In addition, because the sizes of animals we used for our correlation study were relatively homogeneous, the value of EIVPDs in a wider range of LV volumes should be explored. In dilated ventricles, the greater disproportion between the LV chamber and the LVOT hypothetically could increase the convective component of EIVPD. However, we have previously shown that, in patients with dilated cardiomyopathy, reduced ejection velocities cause convective forces to be lower than in normal hearts; in fact, total EIVPDs are smaller.13 Additionally, the small sample size used for studies 1 and 2 limits the precision of statistical estimates and reduces power to detect modest and moderate fixed effects.
For study 3, we used 3D color Doppler to obtain the volumetric data set, therefore constraining data acquisition to an axial interrogation of flow velocity. Phase-contrast cardiac magnetic resonance can provide the full 3D vectorial velocity information noninvasively but at a lower spatial, temporal, and velocity resolution than color Doppler. Although ignoring the other 2 velocity components could theoretically lead to underestimation of the true EIVPD, several pieces of evidence support that the error related to an axial interrogation of ejection flow is negligible. First, magnetic resonance studies have confirmed that ejection flow follows a mostly 1D trajectory without vorticity.35 Second, although LV diastolic filling flow is known to also follow a slightly curvilinear trajectory,35 a similar 1D simplification has been validated for LV filling.19,20,36 Third, the excellent accuracy of Eulers equation for calculating EIVPDs would not stand if significant nonaxial components of flow velocity were present during ejection.13 Another limitation of study 3 is that this type of reliability analyses is obviously model dependent. Consequently, misalignment errors could be higher in ventricles with geometrical variations such as localized septal hypertrophy. In situations like this one in which 2D color Doppler imaging shows that flow follows a tightly curved trajectory, magnetic resonance may be an alternative for measuring EIVPDs.37
Conclusions
Doppler-derived EIVPDs provide quantitative, reproducible, noninvasive, and relatively load-independent indices of LV systolic chamber function that closely correlate with reference methods. This new method should improve the echocardiographic assessment of LV systolic chamber function in the clinical setting.
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
| Footnotes |
|---|
The online-only Data Supplement can be found at http://circ.ahajournals.org/cgi/content/full/112/12/1771/DC1.
Related Article:
Circulation 2005 112: 1684-1686.
This article has been cited by other articles:
![]() |
C. Cortina, J. Bermejo, R. Yotti, M. M. Desco, D. Rodriguez-Perez, J. C. Antoranz, J. L. Rojo-Alvarez, D. Garcia, M. A. Garcia-Fernandez, and F. Fernandez-Aviles Noninvasive Assessment of the Right Ventricular Filling Pressure Gradient Circulation, August 28, 2007; 116(9): 1015 - 1023. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Thomas and Z. B. Popovic Assessment of Left Ventricular Function by Cardiac Ultrasound J. Am. Coll. Cardiol., November 21, 2006; 48(10): 2012 - 2025. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. B. Popovic, K. E. Richards, N. L. Greenberg, A. Rovner, J. Drinko, Y. Cheng, M. S. Penn, K. Fukamachi, N. Mal, B. D. Levine, et al. Scaling of diastolic intraventricular pressure gradients is related to filling time duration Am J Physiol Heart Circ Physiol, August 1, 2006; 291(2): H762 - H769. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Yotti, J. Bermejo, J. C. Antoranz, M. M. Desco, C. Cortina, J. L. Rojo-Alvarez, C. Allue, L. Martin, M. Moreno, J. A. Serrano, et al. A Noninvasive Method for Assessing Impaired Diastolic Suction in Patients With Dilated Cardiomyopathy Circulation, November 8, 2005; 112(19): 2921 - 2929. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Thomas and Z. B. Popovic Intraventricular Pressure Differences: A New Window Into Cardiac Function Circulation, September 20, 2005; 112(12): 1684 - 1686. [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2005 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |