|
|
||||||||
1Hydraulics Laboratory, Institute Biomedical Technology, Ghent University, Gent, Belgium; Division of Cardiology, 2Department of Medicine, and 3Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, Maryland; 4Millar Instruments, Houston, Texas; 5Laboratory of Hemodynamics and Cardiovascular Technology, Swiss Federal Institute of Technology (Écóle Polytechnique Fédérale de Lausanne), Lausanne, Switzerland; and 6Laboratory for Physiology, ICaR-VU, Vrije Universiteit Medical Center, Amsterdam, Netherlands
Submitted 13 May 2004 ; accepted in final form 19 October 2004
| ABSTRACT |
|---|
|
|
|---|
hemodynamics; augmentation index; windkessel; arterial stiffness; compliance; ventriculovascular coupling
Arterial Zin, derived from simultaneously measured aortic pressure (Pao) and flow, however, remains the most complete description of global arterial function (14). In recent work, Reddy et al. (20) measured aortic Zin in the mouse using Doppler ultrasound for flow velocity measurement, whereas pressure was measured invasively using human intracoronary pressure sensor. They demonstrated that the mouse input impedance spectrum is similar to that in other mammals and exhibits the same changes with aging as reported in humans (14, 16). The technology, however, does not allow for simultaneous assessment of cardiac parameters.
In many disease states, alterations in both cardiac and vascular properties occur to maintain blood pressure and flow, thereby making it difficult to independently assess relative changes in each system (23). Therefore, simultaneous acquisition of all components of cardiovascular function is a prerequisite for a meaningful interpretation of individual changes. In this study, we use a high-fidelity 1.4-Fr custom-made, dual-pressure sensor mouse conductance system to simultaneously assess Pao and intraventricular pressure and volume. Given the fact that aortic flow can be estimated as the derivative of ventricular volume during ejection (with the assumption of intact mitral valve function), it is possible to simultaneously acquire cardiac function and arterial Zin using a single catheter. Demonstrating the feasibility of this technique of Zin measurement is the first goal of the study.
Zin is a mathematically complex property that encompasses mechanical characteristics of the cardiovascular system and the effects of wave propagation and reflection. Its interpretation may be facilitated by the use of lumped parameter "windkessel" models (39). These models consist of a limited number of parameters [such as total systemic vascular resistance, aortic characteristic impedance (Z0), and total arterial compliance], which represent biophysical characteristics of the arterial system. By matching measured impedance to the impedance of the model, the Zin is translated in terms of model parameter values that correspond to physiological properties, thereby providing a convenient way to describe the arterial system. The second goal of this study is therefore to fit the measured data to a four-element windkessel model (27) and to provide reference values for the lumped model parameters in the mouse. In the past, we have also done research on methods to assess total arterial compliance, i.e., the area method, the ratio of SV to pulse pressure (PP), and their relation to the PP method (24, 26). We make use of the available data sets to apply these methods and verify whether our findings in large animals also apply to the mouse.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The chest was entered through an anterior thoracotomy by visualization under a dissecting microscope, and a small apical stab was made at the left ventricular (LV) apex, leaving the pericardium as intact as possible. A custom-made, four-electrode conductance catheter with a dual-pressure sensor (Millar Instruments, Houston, TX) was then advanced retrogradely into the LV along the cardiac longitudinal axis with the distal tip (containing a pressure sensor) in the aortic root and the proximal electrode just within the endocardial wall of the LV apex. With the catheter fixed in place, the animal was turned over on its side. Care was taken to maintain catheter position by online visualization of the shape and position of the pressure-volume (PV) loops. A limited lateral thoracotomy was performed, and the descending aorta was dissected free from the spinal column just above the level of the diaphragm. A flow probe (1RB, Transonic; Ithaca, NY) was placed around the thoracic aorta and filled with conducting gel, and the measured flow served for calibration of the conductance catheter (gain factor
) using regressions of SV from thoracic flow and SV from conductance catheter on a beat-by-beat basis during inferior vena cava (IVC) occlusion for each mouse. This provides an estimate of the gain over a large loading range instead of a single estimate at steady state (7). Hypertonic saline infusion was performed for the assessment of parallel conductance (Vc) by bolus injection of 510 µl 35% saline.
Data were measured during steady-state conditions (baseline, BL) and during transient preload decline, achieved by manual compression of the IVC (VCO). Data were also recorded in a subgroup of eight animals during transient afterload increase, achieved by progressive occlusion of the aorta (AOO). To facilitate this manipulation, a 32-gauge stainless steel wire was deployed in the shape of an "L" and slipped underneath the esophagus as support when the thoracic aorta was manually compressed just above the diaphragm, which progressively occluded the aorta until complete occlusion was obtained. All signals were digitized at 2 kHz and stored to disk for subsequent analysis.
Hemodynamic analysis. The VCO procedure yielded 1020 successive cardiac cycles over the ensuing 2 s. We approximated the ESPVR as linear, fitting a linear regression line through the end-systolic points of typically 1525 cycles, selected from the moment that LV pressure (PLV) started to decline until achievement of a new steady state. An iterative method programmed in Matlab (The Mathworks; Natick, MA) was used to identify the end-systolic points, where elastance E(t) was first calculated as PLV/(VLV V0) (where VLV is LV volume) with an initial value for V0 = 0. The points in the PV plane corresponding to maximal E(t) for each cycle were identified as the end-systolic points. Linear regression analysis on these points yielded a first estimate of Emax and a new estimate of V0. This procedure was repeated with the resulting V0 until successive values for V0 (and Emax) did not differ by >0.1%, typically within three to four iterations. We also assessed the slope and intercept of the ESPVR during transient aortic occlusion in the subgroup of eight animals. The period with transient data was shorter than in VCO and yielded 5 to 10 cycles in transient conditions.
Steady-state measurements, obtained during a 1- to 2-s data sequence, contained 1020 cycles and were averaged to construct a representative steady-state beat with Pao, PLV, and VLV tracings. Aortic flow (Qao) was obtained as the time derivative (dVLV/dt) of VLV during the ventricular ejection period: Qao = dVLV/dt, as illustrated in Fig. 1. Because the numerical differentiation of VLV may introduce noise into the Qao signal, ventricular volume data were first filtered using a Savitsky-Golay filter of the third-order and 15-sample frame width (Matlab, The Mathworks), which implies that smoothing is obtained by fitting a third-order polynomial over a 15-sample point window (1 cardiac cycle typically contains 200 samples) sliding over the data. We chose this type of filter because it is better in preserving the pertinent high-frequency components of a signal than more traditional low-pass filters (17). The Savitsky-Golay filter does not have a clear cut-off frequency above which the signal is filtered, but the power spectra demonstrated that the additional smoothing effect of the filter (on top of smoothing effect of averaging the data) was only apparent for frequencies above 200 Hz. When applied by itself to the raw data with the parameters used, the Savitsky-Golay filter reduced the power of the frequencies above 100 Hz.
|
For each animal, Pao and Qao were decomposed into a series of sinusoidal harmonics using a discrete Fourier transform (Matlab, The Mathworks), transforming them from the time domain into the frequency domain. The pressure and flow wave are thus considered as simultaneously present sinusoidal waves (harmonics) with frequencies that are natural multiples of the heart frequency, and the arterial system can then be studied at a specific frequency. It is common to represent these harmonics in their complex form by their amplitude (or modulus) and phase angle. For a given frequency, Zin is then calculated as the ratio of the pressure and flow harmonics at that frequency, and we considered the steady (or DC component) and the first 10 harmonics. Zin is also a complex number, represented in the frequency domain, with the amplitude of Zin for a given harmonic frequency being the ratio of the amplitudes of the same pressure and flow harmonic frequency, and the phase angle the difference between pressure and flow phase angle of the same harmonics. The ratio of mean Pao and Qao (DC component of Zin) is systemic vascular resistance (R). Aortic Z0 was estimated 1) in the frequency domain (Z0-FD) as the average modulus of the fifth to tenth harmonic and 2) in the time domain (Z0-TD) by calculating the slope to the Qao-Pao relation in the early systolic upstroke, where this relation was linear (6). Pao was separated into a forward (Pf) and backward (Pb) running components using the linear wave separation method (38), and the reflection coefficient was calculated as the ratio of the amplitude of Pb and Pf. We also calculated the augmentation index (AIx) as a measure of arterial wave reflection intensity (16). To further allow comparison of our data with parameters commonly used in arterial function analysis, we estimated pulse-wave velocity (PWV) from Z0-FD as PWV = Z0-FDA/
, assuming density of blood
= 1,060 kg/m3 and (constant) aortic cross-sectional area A = 1.13 mm2 (diameter 1.2 mm) (8). Because Z0 represents the characteristics of the proximal aorta, the calculated PWV applies to the proximal part of the aorta.
The averaged steady-state data at baseline, VCO, and AOO were subsequently fitted to a four-element lumped parameter windkessel model consisting of total peripheral resistance, total arterial compliance (C), total inertance (L), and aortic Z0. Fitting was performed using a Matlab nonlinear least-squares fitting algorithm (Gauss-Newton method) with default settings for convergence and tolerance criteria. Fitted Pao was calculated from the response of the four-element windkessel model to the measured Qao, which was applied as an input into the model (25). Simultaneous fitting of all four parameters with this routine introduced large uncertainty in the estimation of C and L, with results depending on initial values. Therefore, these two parameters were fitted in an iterative way. First, C was given a fixed arbitrary starting value (0.5 µl/mmHg), and R, Z0, and L were assessed. Subsequently, R, Z0, and C were fitted with fixed L, using the value from the first iteration. This alternating scheme was repeated until values for the four parameters no longer changed, which occurred maximally after two repetitions. Confidence intervals of 95% of the four parameters were calculated.
The difference between measured (Pao) and fitted (Pfit) pressure was quantified by the root mean square error (RMSE) calculated as
![]() |
![]() |
Total arterial compliance was further assessed using earlier described methods, including the PP method (CPPM) (26), the area method (Carea) (12), and the ratio of SV and PP (CSV/PP) (3).
Statistical analysis. The statistical analysis related to the four-element windkessel model fitting was described earlier in MATERIALS AND METHODS. Paired t-tests were applied to assess the effect of VCO or AOO on arterial mechanical parameters. Linear regression and Pearson correlation analysis were used to quantify linear relations between parameters. To assess the effect of VCO or AOO on Zin modulus or phase, one-way repeated measures ANOVA was performed. When the ANOVA reached statistical significance (P < 0.05), post hoc analysis was done at each harmonic frequency using Bonferroni correction to assess the statistical significance level of the measured difference. All statistics were performed using SPSS (version 11.5, SPSS, Chicago, IL).
| RESULTS |
|---|
|
|
|---|
|
|
60 Hz).
|
|
Our results demonstrate that the mouse Zin spectrum is well described by the four-element windkessel model. For the baseline and VCO occlusion data, the average RMSE and %AD, indicating the difference between measured and fitted Pao, are 1.68 ± 0.60 mmHg (range 0.873.53 mmHg) and 1.77 ± 0.99% (range 0.726.16%), respectively. These numbers increase somewhat with AOO, the average RMSE difference being 7.81 ± 3.45 mmHg (range 4.9215.53 mmHg), and %AD 4.86 ± 1.70% (range 3.208.42%). The values of the four-element windkessel parameters and the 95% confidence intervals are given in Table 3.
|
At baseline, effective Ea is 4.46 ± 1.62 mmHg/µl, so that Ea/Emax = 1.44 ± 0.43. This ratio remains unchanged with VCO (Ea/Emax = 1.36 ± 0.58, P = 0.47) and AOO (Ea/Emax = 1.21 ± 0.37, P = 0.50).
| DISCUSSION |
|---|
|
|
|---|
In humans, the progressive stiffening of the arterial tree with age or in cardiovascular disease (hypertension, diabetes) has an impact on arterial hemodynamics. In terms of pressure and flow-wave phenomena, the increase in PWV (which may double from about 5 m/s to over 10 m/s from 20 to 80 yr) (1) leads to an early return of reflected wave, an increase in pressure augmentation [from negative values for individuals younger than 30 yr to values in the range of 40% or higher in older subjects (5)], and an increase in SBP (18, 37). In the "lumped parameter" view, the decreased compliance of the arterial system leads to an increase in PP (19, 28), presenting itself as an increase in SBP, often reaching values designated as cut-off values for hypertension (>140 mmHg), and a decrease in DBP (<80 mmHg). This phenomenon is called "isolated systolic hypertension."
Whereas Zin was measured earlier in larger rodents such as rats (13, 15), the assessment of arterial properties in mouse models is fairly new and made possible through the further miniaturization of (pressure) sensor technology and new developments in ultrasound technology. Recently, Reddy et al. (20) studied the effects of aging on arterial stiffness and reported data on aortic Zin, PWV, and AIx in adult and old mice. Velocities were measured with ultrasound and pressures with the RADI-wire system commonly used in interventional cardiology. Reddy et al. (20) reported significant anticipated age-related changes in SBP (from 88 to 116 mmHg), PP (from 29 to 42 mmHg), PWV (from 2.86 to 4.16 m/s), and AIx (from 14 to 37%) in adult and old mice (20). The minimum in Zin shifted from the second to the fourth harmonic in the older mice, further demonstrating the feasibility of measuring age-related changes in arterial function with the applied methodology. Note, however, that the absolute values of blood pressure, PP, and PWV (in both adult and old mice) are somewhat lower than one would expect in conscious animals, possibly due to specifics of the experimental preparation and the anesthetics used, which may also be responsible for the relatively low heart rate (400 beats/min) in their setting. Given the high (nonlinear) dependency of arterial mechanical parameters on blood pressure (11), the reported numbers on PWV and Z0 may therefore be an underestimation of their values in the old, conscious animal.
In our model, aortic banding induces a pressure-modulated increase in arterial stiffness (Z0; decrease in total arterial compliance), increasing SBP, PP, and AIx up to 177 mmHg, 80 mmHg, and 146%, respectively. Although our model does not account for changes in aortic structure and function related to aging or cardiovascular disease, these values do resemble what is observed in older subjects and in isolated systolic hypertension. When we make some assumptions on the dimensions of the mouse aorta [i.e., assuming a constant diameter of 1.2 mm (8)], it is possible to further quantitatively compare some of our data to what was found by Reddy et al. In particular, the Z0 values reported by Reddy et al. (20) would vary from 0.020 (adult) to 0.031 (old) mmHg·µl1·s, which is somewhat lower than what we found, but which may again be attributed to the lower pressures in their setting. Converting Z0 data into the more intuitively interpretable PWV, the increase in Z0 with AOO is equivalent with an increase in PWV from 3.50 to 5.95 m/s. Given the fact that the proximal aortic cross-sectional area is expected to increase with AOO, the estimated value of 5.95 m/s is even expected to underestimate PWV during AOO. Again, these values are close to what is found in human pathophysiological conditions associated with increased vessel stiffening. These data at least suggest that aortic banding may be an interesting mechanical model to study (chronic) effects of increased afterload in the mouse, and that the conductance catheter-based methodology for assessing arterial function parameter is sensitive enough to detect pressure-modulated increase of arterial stiffness.
There are also some important methodological differences between our study and that of Reddy et al. (20). We approximated aortic flow by differentiating an averaged and filtered ventricular volume tracing, operations that have an effect on the (high) frequency contents of the signal, and hence potentially on some of our results. Analysis of the power spectrum of the volume and flow data demonstrated that averaging of the time series (construction of ensemble average), an operation that is common in the analysis of hemodynamic data, smoothed the data by reducing the high-frequency contents above 100 Hz (greater than the tenth harmonic). As explained before, the Savitsky-Golay filter smoothing reduced the power only for frequencies higher than 100 Hz. We are thus confident that the followed procedure preserved the information contained in the first 10 harmonics (the harmonic range used to derive Zin). Also, the methodology did not seem to have an impact on the assessment of Z0, which depends on the high-frequency components of pressure and flow and is therefore susceptible to filter settings. We used both a time and frequency domain method to assess Z0, with acceptable agreement between both methods as expressed by the linear regression and Bland-Altman analysis (Fig. 4). The mean difference and limits of the agreement between the two methods are 0.00087 and 0.00911 per 0.0108 mmHg·µl1·s, respectively, and paired t-test indicates no difference (P = 0.26) between the two approaches. Although this is not a fool-proof validation of the methodology (which would require simultaneous measurement of velocity or flow), it is an indication that the frequency contents of the flow signal is sufficiently high for the calculation of impedance spectra. Reddy et al. (20) calculated impedance from flow velocities in the aorta measured with ultrasound equipment, which certainly has the advantage of being noninvasive. On the other hand, the sampling volume covers most of the aortic cross-section, which has the advantage that all velocities within the volume are being sampled but which also has the drawback that there is no information on the precise spatial distribution of these captured velocities. As such, one needs to make assumptions about the shape of the velocity profile (flat vs. parabolic), which may not be easily predicable in complex bending, three-dimensional geometries such as the aorta and which may be modified by physiological or pathological alterations in SV or heart rate.
|
It is common, especially in studies on heart-arterial interaction, to characterize the arterial system in terms of the effective Ea (4, 10), whereas the ventricle is characterized by the slope of the end-systolic PV relation Emax. The link between the Ea/Emax framework and ventricular mechanoenergetics further favored the widespread use of Ea/Emax as the heart-arterial coupling parameter (29, 30). Theoretical work (2), which was subsequently confirmed in large mammals (33, 34), has shown that Ea/Emax is 0.5 to 1 in the normal heart, with the left ventricle operating close to optimal efficiency or stroke work, respectively. In heart failure, coupling deteriorates, resulting in progressive increase of Ea/Emax. In our data set, average Ea/Emax was 1.44 ± 0.43 at baseline, which is higher than the value of 0.49 reported earlier by Georgakopoulos et al. (7). On the other hand, the values for Ea and Emax are similar to recent findings of Reyes et al. (21) who reported Ea and Emax values of 5.9 and 3.1 mmHg/µl, respectively (Ea/Emax about 1.9) (21). We attribute the difference between these and our previous data at least in part to the animal preparation, which we consider to be much more physiological in this setting. Since these initial experiments, the preparation has undergone substantial modifications, including the use of different anesthesia, fluid supplementation, optimized ventilation parameters, and minimal surgery. This has subsequently resulted in increased systolic pressure, end-diastolic volume, ejection fraction, and maximal first derivative of pressure with a simultaneous decrease in heart rate compared with our first study.
There are some limitations inherent to our study. First, the outcome of the study, in terms of absolute values for the different cardiovascular parameters, highly depends on appropriate calibration of the conductance catheter, which was based on an aortic flow probe around the thoracic aorta. It is, however, technically very difficult to place a flow probe on the ascending aorta of the mouse with a conductance catheter in the heart so as not to distort the signals. The values that we found for Emax and Ea are significantly lower than those previously reported by our group (7) (using earlier prototype conductance catheters in other animal preparation settings) but they are in good agreement with recent data of Reyes et al. (21) obtained with a dual-frequency conductance catheter system. Second, it has been demonstrated that, especially in conditions of enhanced contractility (35), the ESPVR may be better approximated by a nonlinear (e.g., quadratic) curve, particularly in small mammals such as the rat (22). The data presented in Fig. 2 demonstrate that the ESPVR is also curvilinear in the mouse. However, given the fact that the focus of the paper is in quantification of arterial function with eventual extension to ventriculovascular coupling using the Ea/Emax framework (where linear ESVPR relation is assumed), we chose to use the simple linear approximation of the ESPVR, using all transient beats during AOO and VCO to avoid bias by solely selecting beats at high or low pressures. Third, experiments were conducted in an open thorax setting, which affects intrapleural pressure and thus transmural pressure across the wall of the thoracic aorta. This may have an effect on pressure-dependent mechanical properties, and thus on Zin, PWV, and AIx. On the other hand, we did experiments with the conductance catheter placed into the heart through the right carotid and recorded PV data before and after the diaphragmatic incision and found no changes in the values. Fourth, although data were measured in mice of different strains, the groups are too small for subgroup analysis of arterial system properties in a statistically reliable way. Finally, in terms of practical applicability of our technique, its invasive character remains a drawback, making the technique less suitable for follow-up measurements in chronic experiments. In these settings, ultrasound technology may be more suitable, especially now that it has been shown that diameter distension of superficial arteries can be measured in the mouse (9). The morphology of these diameter distension traces closely resembles intra-arterial pressure morphology and may offer possibilities to obtain central pressure in a noninvasive way.
In conclusion, we have shown that a dual-pressure conductance catheter can be used to assess both cardiac and arterial parameters. Because no additional methodology is required when PV studies are performed in mice, this technique provides for complete characterization of the cardiovascular system. This should be useful in the assessment of genetic models of heart failure and potential therapies.
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Maloyan, H. Osinska, J. Lammerding, R. T. Lee, O. H. Cingolani, D. A. Kass, J. N. Lorenz, and J. Robbins Biochemical and Mechanical Dysfunction in a Mouse Model of Desmin-Related Myopathy Circ. Res., April 24, 2009; 104(8): 1021 - 1028. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Shioura, D. L. Geenen, and P. H. Goldspink Sex-related changes in cardiac function following myocardial infarction in mice Am J Physiol Regulatory Integrative Comp Physiol, August 1, 2008; 295(2): R528 - R534. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Pokreisz, G. Marsboom, and S. Janssens Pressure overload-induced right ventricular dysfunction and remodelling in experimental pulmonary hypertension: the right heart revisited Eur. Heart J. Suppl., December 1, 2007; 9(suppl_H): H75 - H84. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Shioura, D. L. Geenen, and P. H. Goldspink Assessment of cardiac function with the pressure-volume conductance system following myocardial infarction in mice Am J Physiol Heart Circ Physiol, November 1, 2007; 293(5): H2870 - H2877. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Marsh, P. C. Powell, A. Agarwal, L. J. Dell'Italia, and J. C. Chatham Cardiovascular dysfunction in Zucker obese and Zucker diabetic fatty rats: role of hydronephrosis Am J Physiol Heart Circ Physiol, July 1, 2007; 293(1): H292 - H298. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Wu and S. J. Kovacs Frequency-based analysis of the early rapid filling pressure-flow relation elucidates diastolic efficiency mechanisms Am J Physiol Heart Circ Physiol, December 1, 2006; 291(6): H2942 - H2949. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Jegger, R. da Silva, X. Jeanrenaud, M. Nasratullah, H. Tevaearai, L. K. von Segesser, P. Segers, V. Gaillard, J. Atkinson, I. Lartaud, et al. Ventricular-arterial coupling in a rat model of reduced arterial compliance provoked by hypervitaminosis D and nicotine Am J Physiol Heart Circ Physiol, October 1, 2006; 291(4): H1942 - H1951. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. E. Claessens, D. Georgakopoulos, M. Afanasyeva, S. J. Vermeersch, H. D. Millar, N. Stergiopulos, N. Westerhof, P. R. Verdonck, and P. Segers Nonlinear isochrones in murine left ventricular pressure-volume loops: how well does the time-varying elastance concept hold? Am J Physiol Heart Circ Physiol, April 1, 2006; 290(4): H1474 - H1483. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |