In vivo measurement of proximal pulmonary artery elastic modulus in the neonatal calf model of pulmonary hypertension: development and ex vivo validation

Kendall S. Hunter, Joseph A. Albietz, Po-Feng Lee, Craig J. Lanning, Steven R. Lammers, Stephen H. Hofmeister, Philip H. Kao, H. Jerry Qi, Kurt R. Stenmark, Robin Shandas


Developing clinical work suggests that vascular stiffening plays a role in the progression of pulmonary hypertension (PH), while recent studies in animal models of hypoxic PH have found significant proximal vascular stiffening in the diseased population. Here, we develop and validate a minimally invasive, clinically realizable method to estimate the local elastic modulus of the proximal pulmonary arteries from pressure-diameter (PD) data. PD measurements were made in the main pulmonary arteries of 16 calves; lumen diameter was assessed using color M-mode tissue Doppler imaging ultrasound, while pressure was measured via catheter. Two methods corresponding to thin-walled pressure vessel theory (“thin”) and Lame's equation for a thick-walled cylinder (“thick”) were used to approximate the artery elastic modulus from PD measurements. The harvested arteries were tested independently to determine their “true” ex vivo elastic modulus and stiffness. Both approximations displayed excellent correlation with ex vivo elastic modulus of the calf main pulmonary artery (thin r2 = 0.811; thick r2 = 0.844; both P < 0.01). Bland-Altman analysis indicated that the thick-walled approximation has better overall agreement with ex vivo modulus. The approximations displayed quantitatively distinct regression slopes that were statistically different (P = 0.02). The elastic modulus of the main pulmonary artery can be reasonably estimated from combined color M-mode tissue Doppler imaging ultrasound and catheter pressure measurements in calves. Such measurements may be a valuable tool in the diagnosis and treatment of human PH.

  • pulmonary hypertension
  • hypoxia
  • animal model
  • pulmonary vascular stiffness

pulmonary hypertension (PH) has been historically associated with only distal events, such as sustained vasoconstriction, occlusive vascular remodeling, and in situ thrombosis. These pathologies can result in increased pulmonary vascular resistance (PVR), a physiological parameter that is the principal diagnostic for PH in the clinical setting (3). PVR is traditionally thought to be the principal determinant of right ventricular (RV) afterload that, when chronically increased, can lead to right heart failure. However, recent evidence suggests that, despite its relative importance in PH diagnosis, PVR is only a modestly useful prognostic tool (2, 28). PVR quantifies only the static component of pulmonary hemodynamics to the neglect of the dynamic component, described by pulmonary vascular stiffness (PVS). PVS refers to an elastic pulmonary artery's innate resistance to deformation under a pressure load: it can change with the dimensions of vessel or due to changes in the vessel's inherent material properties.

Several recent clinical studies have suggested that global PVS is a better prognostic than PVR alone, in that elevations in PVS are associated with increased adult mortality (14, 29) and poorer pediatric outcomes (18). This association between vascular stiffening and disease severity is believed to be due to several factors: 1) stiffening reduces the windkessel effect of the elastic arteries, leading to greater flow inefficiency (33, 34); 2) stiffening increases RV afterload independent of PVR (30, 42); and 3) changes in the flow and pressure waveforms due to stiffening in both the systemic (39, 45) and pulmonary circulations (26, 27) have been shown to stimulate cellular signaling pathways in distal vessels and the lungs that exacerbate and perpetuate existent vascular pathologies.

Despite this increased appreciation that vascular stiffness plays an important role in vascular pathology, the contribution of stiffening of the great arteries [specifically the main (MPA), right, and left pulmonary arteries] to the pathophysiology of PH is not fully understood. The great arteries were recently found to chronically increase in stiffness in PH due predominantly to elastin-based extracellular matrix remodeling in an ex vivo study of severe PH in the neonatal calf (22). It has also been found that the stiffness of the proximal pulmonary vasculature of the mouse model increases in response to chronic hypoxia, with similar mechanical consequences (20, 21). It remains unclear whether significant proximal remodeling in the human produces similar mechanical consequences, and whether such stiffening is important to disease outcomes and/or pathogenesis.

To answer clinical questions about proximal stiffening, we sought to assess it in a minimally invasive way in a large animal model of severe PH (36). While the most comprehensive means of measuring mechanical properties, such as vascular stiffness, or elastic modulus, which quantifies only the inherent material properties, involves in vitro uniaxial or multiaxial stress-strain testing (1012, 22, 44), such a method is obviously not possible for serial monitoring or clinical studies. Methods do exist, however, to approximate mechanical properties from the in vivo response of an artery distending during systole; this response is measured by the artery's pressure-diameter (PD) curve. From PD curves and an estimate of thickness, previous work has obtained estimates of stiffness and elastic modulus (16, 23, 31), incremental elastic modulus (7), distensibility (1, 4, 24, 35), and compliance (43).

The novelty of our approach here is in the use of noninvasive ultrasound color M-mode tissue Doppler imaging (CMM-TDI) to accurately obtain instantaneous diameter of the pulmonary artery during cardiac catheterization (8, 13, 25); combining these two measurements, echo-assessed diameter and catheter-assessed pressure, the possibility of routinely obtaining modulus in a clinical setting emerges. From a structural mechanics viewpoint, the modulus of the arterial wall should be a useful parameter because it reflects basic material properties. Given the link between stiffness quantification (with pulmonary vascular input impedance) and improved prognostication (18), we speculate that this measure could provide similar prognostic benefits. In this study, we 1) develop two approximations to compute incremental elastic modulus from PD data; 2) apply these to obtain PD data in the chronically hypoxic neonatal calf model; and 3) compare the in vivo elastic modulus approximations thus obtained to direct (ex vivo) measures of tissue stiffness to assess the method's accuracy.


Equations to compute elastic modulus.

Equations for the approximation of elastic modulus from PD data are well known (10, 15); here, we briefly review our implementation specific for use with ultrasound data, which provides lumen diameter. The typical simplifying assumptions for use of these approximations are isotropy, indicating no variation in mechanical properties with direction; homogeneity, indicating similar material throughout the arterial wall; and local linearity in mechanical behavior, which assumes that the stress is a linear function of the strain for strain ranges of interest. In addition, for our echo-based acquisition of diameter, we must assume the thickness of the arterial wall (as a percentage of the arterial diameter); this will be addressed in greater detail below.

Thick-walled plane strain modulus.

The thick-walled plane-strain model is based on the two-dimensional plane strain approximation for a thick-walled elastic tube, which allows for variation in stress across the thickness and is appropriate when an artery's walls are greater than 10% of its diameter. The geometry of the tube cross section is described by thickness h, diastolic inner and outer diameters DID and DOD, respectively, and the instantaneous incremental transmural pressure PI. The instantaneous inner diameter at any time is referred to as DI. Simplification of this model's stress equilibrium equation provides a linear relationship between the incremental pressure and instantaneous diameter of the form DI = m PI + DID, in which m is a function (40) of the elastic modulus and the DID and DOD. Because the parameter m also corresponds to the slope of the incremental PD curve, we may equate this slope (obtained from experimental measurements) to the function and solve for the modulus, ETWPS=(1+ν)DIDDOD2+(1ν)DID3m(DOD2DID2) (1) where ETWPS is thick-walled plane strain modulus, and υ is Poisson's ratio (set to 0.45 to model the nearly incompressible nature of the artery wall). Because the diastolic diameters are constants in this expression for each PD curve, and the slope m may be extracted from linear regression of the PD data, ETWPS is readily obtained from Eq. 1.

Pressure-strain modulus.

The modulus Ep = (PS − PD)DID/(DISDID) has been used in prior studies of arterial mechanics (10, 15, 23, 31), where PS and PD are systolic and diastolic pressures (their difference being the pulse pressure), and DIS and DID are the artery inner diameter at systole and diastole, respectively. Applied to the outer diameter, it is known as the Peterson modulus (31). Ep is not a true elastic modulus (it lacks thickness information and is thus a ratio of pressure to strain rather than stress to strain) and is not investigated herein; however, an especially simple elastic modulus approximation based on Ep can be calculated as follows. The thin-walled model circumferential stress is defined as σθ = PID/2h; this yields an incremental change between systole and diastole of σθ = (PS − PD)DID/2h, while the same incremental change in the circumferential strain is εθ = (DISDID)/DID. We thus define the pressure-strain elastic modulus (EPS) as: EPS=(PSPD)DID/2h(DISDID)/DID=(PSPD)(DISDID)/DID×DID2h=Ep×DID2h (2) EPS is thus understood as Ep scaled by the wall-thickness-normalized diameter.

Ultrasound-based acquisition of diameter.

Institutional Animal Care and Use Committee approval was obtained for all animal studies. Sixteen newborn male dairy calves (Holstein), weighing between 35 and 50 kg, were used. Nine of the calves were put into a large hypobaric chamber at an atmospheric pressure of 445 mmHg, equivalent to an elevation of 4,300 m, between 1 and 3 days of age. These calves were maintained at this simulated altitude for 2 wk and developed hypoxic PH; the remaining animals remained at local ambient pressure (630 mmHg; Ft. Collins, CO) while awaiting measurement. Daily care for all animals proceeded as described previously (36) under the supervision of large-animal veterinarians. Both animal groups were at ambient pressure for at least 1 h before testing and underwent testing at ambient pressure to minimize the potential for residual hypoxic proximal vasoconstriction.

To begin each study, a calf was placed on its right side on an animal examination table. The animal's eyes were covered during the entire procedure to induce calm, and sedation was not required. A solid-state (Millar SPC-350) catheter was inserted percutaneously into the jugular vein to access the RV and MPA; placement was continuously assessed via evaluation of the pressure waveform. Pressure was recorded digitally by a laptop-computer-based data-acquisition system (LabView, National Instruments, Austin, TX) and by a commercial ultrasound machine through its PHYSIO port (Vivid 5, GE Medical Systems, Waukesha, WI). To allow echo measurements, a 15 × 15-cm patch was shaved on the right side of the chest, between the third and fifth intercostal spaces; a 3.5-MHz FPA probe accessed the animal through a hole in the underside of the examination table. A caudal, short-axis view at the fourth intercostal space, 3–6 cm dorsal to the elbow, was used to obtain CMM-TDI measurements of instantaneous diameter and velocity of both walls of the MPA; before acquisition, the ultrasound beam was swept through the long axis of the MPA to determine maximal diameter. Both pressure and ECG time histories were recorded with the echo images to ease offline processing.

Echo image processing.

Analysis of the combined CMM-TDI data pressure was performed offline using custom-written MATLAB V7.7 analysis software (The Mathworks, Natick, MA; code based on EchoMAT version 2.1, GE Medical Systems, Waukesha, WI). The program displays the data visually; an example screenshot is shown in Fig. 1 and displays the two-dimensional echocardiogram (top plot) and M-mode image from a normoxic calf (the TDI data is not shown), as well as the electrocardiogram trace and pressure time history (blue and yellow traces, respectively, bottom plot). The program uses an image decimation-based method [described in detail elsewhere (17, 18)] to initially find the upper and lower wall boundaries as seen in the M-mode image. After this first pass, the code then time integrates the velocity data (obtained from the TDI image) along the initial trace to obtain an improved, final wall definition. In this way, both the M-mode and TDI data are incorporated into calculation of the final trace, while the integration yields a smoother wall definition than would M-mode data alone. The typical results for the upper and lower boundary are shown in Fig. 1 as solid red lines, whereas the computed diameter is shown in blue floating between the two wall boundaries (not to scale).

Fig. 1.

Software interface for the processing of in vivo color M-mode tissue Doppler images of wall motion. The brightness-mode image appears as a triangular sector at the top of the image, M-mode data is at the center, while pressure (yellow) and ECG (blue) traces appear at the bottom. Also seen in the M-mode image are upper and lower main pulmonary artery wall definitions (red lines) and superimposed diameter (blue line).

After wall detection is complete, the PD data for each cardiac cycle are separated using the ECG trace and stored for analysis. PD extrema and linear regression slopes of the PD and stress-strain curves are computed for the end-diastole to end-systole portion of each cycle. During postprocessing, we found that using only this portion of the curve (corresponding to loading) was both consistent with our ex vivo protocol (see Ref. 22) and more reproducible cycle to cycle compared with using the whole cycle. Means of the PD extrema and the error-weighted averages (9) of the PD and stress-strain regression slopes are then found and used to compute the reported values of in vivo elastic modulus.

Calculation of in vivo wall thickness.

While the lumen of the proximal arteries is well defined in the CMM-TDI image, the limit of the adventitia is often more difficult to discern. As a result, we chose to estimate wall thickness for the elastic modulus approximations based on an average ex vivo measured value; this avoids the uncertainties associated with determining the outer wall in the echo image and simplifies postprocessing. We avoided using individual thicknesses measured from each artery, because such measurements will not be available in a clinical application of this method.

By processing high-resolution photographs (see Fig. 2) of the fresh MPA from the 16 animals and direct caliper measurements, we found diameter and thickness values and calculated wall thickness to diameter ratios (h/D) for each animal. While the hypertensive animal MPAs displayed a significant (P = 0.018) 13% increase in thickness with respect to the controls, their h/D values were not significantly different from the control h/D (P = 0.121), and thus we chose the mean of all animals, h/D = 0.178 (± 0.018), to represent both groups. An animal's estimated wall thickness was then obtained by multiplying its diameter at diastole (as determined from the echo image) by this group-average h/D.

Fig. 2.

Example cross-sectional image of hypertensive calf main pulmonary artery used for determination of diameter and wall thickness.

Ex vivo measurements.

To validate our elastic modulus approximations obtained from PD measurements, we compared them to ex vivo measurements of the tissue elastic modulus (22). Briefly, after death, the heart and lungs were removed from the animal in-block, and the right pulmonary artery, left pulmonary artery, and MPA were processed into circumferential strips and tested on a standard material testing system (MTS Insight II, Eden Prairie, MN) to obtain bench-top or “true” ex vivo elastic modulus (Eex) over a wide range of strains. As noted earlier, stiffness is the total measure of an elastic material's total resistance to mechanical deformation and may change due to both dimensional and inherent material factors, while elastic modulus quantifies only the inherent material property change. Here, Eex is the measured force normalized by the test sample area (22). We note that residual stress was not considered in the calculation of ex vivo parameters.

To compare these ex vivo measures to the in vivo approximations, values for Eex from the MPA must be found at relevant physiological strains (i.e., at the same relative operating conditions that existed in vivo). These physiological strains are found by relating an in vivo pressure measurement to an ex vivo strain state with Lame's equation, as described previously (22). We then find average “whole-cycle” values for Eex from the stress-strain linear regression slopes over the physiological strain range; our intent is to obtain an average modulus and stiffness representative of the physiological process of arterial distension. During this process, no regression had an r2 < 0.59, with the majority (88%) having r2 > 0.95, suggesting the use of this linear approximation to obtain slope was acceptable.

Statistical analysis.

The F-test and appropriate Student's t-test were used to detect differences between animal groups. Two tests were used to assess validation: 1) univariate linear regression analysis was performed using SAS 9.1 (SAS Institute, Cary, NC) to find relationships (and confidence limits) between the in vivo measured quantities and their corresponding ex vivo counterparts; and 2) Bland-Altman analysis (5) was used to assess limits of agreement (LoA). Analysis of covariance (32) was performed in MATLAB to examine differences in the correlations between the two approximations and ex vivo measures. All results are reported as statistically significant at the 95% confidence level, and P values refer to the two-tailed value. Finally, propagation of error analysis (9) was used to determine bias and random errors of the direct measurements and to estimate slope uncertainty (6); the details of this analysis are omitted. This propagated uncertainty is shown as error bars in each regression plot of the results.


Nine hypertensive (Hy1–Hy9) and seven control animals (Cn1–Cn7) were studied, and in vivo and ex vivo mechanical measures were obtained (see Table 1). As expected for this model of PH, the animal groups had significantly different (P < 0.001) mean pulmonary artery pressures, with the hypertensive group having a mean pressure (58 mmHg) more than double that of the control group (26 mmHg). Only one animal displayed small incremental strains (Hy2: 7.7%), although even this small deformation was over 22 times the M-mode resolution, and the TDI image data were readily used to obtain the wall definitions from echo for all animals. Interestingly, the incremental strain (assessed purely from the echo data) was also significantly different between groups (also see Table 1, P = 0.003).

View this table:
Table 1.

Calf hemodynamic, mechanical, and ultrasound image data grouped by condition (control, hypertensive)

Representative measured PD curves from control (Cn3) and hypertensive (Hy1) animals are shown in Fig. 3 during systole of multiple cardiac cycles, along with regression lines, obtained by averaging the individual linear regression slope and intercept parameters from each cardiac cycle. Notably, these individual correlations all had r2 > 0.98 (and for a majority of measured curves of all animals), indicating the PD curves may be reasonably represented by a linear approximation. The control curves are readily distinguished from the hypertensive curves by both their greater range of displacement (and thus greater incremental strain) and the lower operating pressure (and thus lower incremental stress).

Fig. 3.

Representative measured pressure-diameter (PD) curves obtained from control (○) and hypertensive (+) calves.

Error of prediction: validation of in vivo measurements.

As shown previously through postmortem bench-top measurements (22), hypoxia-induced PH induces statistically significant increases in the elastic modulus of the calf MPA (Fig. 4, ex vivo, P < 0.001). We see that our in vivo approximations allow these increases to be assessed in the intact animal: both show significant differences between groups (ETWPS, P = 0.001; EPS, P = 0.001). However, increases in in vivo modulus are not identical to increases in the Eex: from control to hypertensive condition, the means of ETWPS and EPS increase 125 and 132%, respectively, while Eex increases only 99%. The in vivo means from control and hypertensive animals are, on average, 46 and 38% smaller than the means of Eex for controls and hypertensive animals, respectively.

Fig. 4.

Group comparisons of two elastic modulus approximations [thick-walled plane strain (ETWPS), pressure-strain (EPS)] and the ex vivo elastic modulus (Eex). Values are means ± SD. *P < 0.005; #P < 0.001.

Table 2 provides Pearson correlation coefficients for regressions of the elastic modulus approximations against Eex. Regressions of the two approximations against Eex are shown in Fig. 5, A (ETWPS = 0.786 Eex − 9.69, r2 = 0.844, P < 0.001) and B (EPS = 0.516 Eex − 5.23, r2 = 0.811, P < 0.001). Also shown are dashed lines, showing the 95% confidence limits for the predicted value of Eex. Neither regression displays “exact validation” (e.g., slope = 1, intercept = 0), although both have significant correlation and goodness of fit is of good quality overall. Comparing Fig. 5, A and B, we see qualitatively that the thin-walled approximation (EPS) regression slope appears distinct from that of the ETWPS vs. Eex regression slope; this is confirmed by analysis of covariance, which reveals that the two slopes are significantly different (P = 0.023).

View this table:
Table 2.

Pearson correlation coefficients and P values for regressions comparing the in vivo elastic modulus approximations to the ex vivo elastic modulus (all for measurements on the main pulmonary artery)

Fig. 5.

The two elastic modulus approximations [ETWPS (A), EPS (B)] regressed against Eex of the calf MPA. Error bars refer to propagated measurement uncertainty. Open and shaded circles are indicative of control and hypertensive animals, respectively.

Bland-Altman plots are shown for ETWPS and EPS vs. Eex in Fig. 6, A and B, respectively. The thick-walled approximation has a roughly uniform bias of −42.0 kPa and LoA of 52.4 kPa, while the thin-walled approximation has both greater and nonuniform bias (mean of −78.2 kPa) that increases proportionally with average modulus and higher agreement limits (71.3 kPa). Each plot has a single outlier, indicated by a dashed circle: calf Hy4, which displayed substantially larger Eex compared with ETWPS (70.3% larger) and EPS (152.3% larger). When this outlier is removed, ETWPS bias and LoA improve to −37.8 and 42.6 kPa, respectively, while EPS bias and LoA improve to −73.2 and 62.1 kPa, respectively.

Fig. 6.

Bland-Altman agreement analysis between the elastic modulus approximations [ETWPS (A), EPS (B)] and the Eex of the calf MPA. Each solid middle line is the bias, and the bounding dashed lines indicate the limits of agreement (± 1.96 SD). Open and shaded circles are indicative of control and hypertensive animals, respectively. Each plot has a single outlier, indicated by a dashed circle (see results for details).


Current clinical evaluation of PAH relies primarily on measurement of PVR; however, incipient clinical analysis has begun to demonstrate the importance of vascular stiffening on the prediction of outcomes and mortality in this disease (14, 18, 29). In this work, we examined the utility of two equations that use PD data to estimate proximal pulmonary vascular elastic modulus in the intact animal; these estimates were validated against “true” elastic modulus of these tissues measured with an industry-standard uniaxial material tester. We obtained pressure and diameter data from invasive measurements of pressure and noninvasive CMM-TDI ultrasound, respectively. Our results suggest that both thick- and thin-walled approximations may be used to estimate Eex, although the thick-walled approximation is preferred. The implications and limitations of these topics are discussed below.

Estimation equations.

Two basic equations to obtain elastic modulus from PD data corresponding to thin-walled and thick-walled theory were derived to determine if a thick-walled approximation was required to accurately estimate elastic modulus in the intact animal. Both equations needed wall thickness, which was estimated using an average of data obtained from extracted (ex vivo) tissues. Our analyses indicate that the thick-walled equation results are more accurate, have tighter agreement bounds, and are significantly different from the thin-walled equation results. Both underestimate the true modulus; this could be expected, given this application of infinitesimal deformation theory, which assumes constant thickness, to our clearly finite-displacement in vivo conditions, in which the wall must undergo thinning. Furthermore, the underestimation of circumferential stress even at smaller strains by thin-walled theory accounts for the uniformly smaller values of elastic modulus values obtained from the EPS approximation and for its increasing agreement bias as the vessels become thicker due to PH.

In practice, EPS is the easier approximation for which to obtain measurements and to compute, only requiring PD extrema and simple arithmetic, respectively. A recommendation seems premature, however, given the limited number of animals examined in this study. Preliminary application of this method in pediatric PAH patients (19) suggests that the tighter agreement limits of ETWPS are likely more suitable for differentiation of pediatric patient groups. These tighter agreement limits, along with its better correlation coefficient, suggest that ETWPS is the more accurate measurement. Ongoing work by our group seeks to determine whether either approximation is a better predictor of patient outcomes compared with PVR alone.

Diameter acquisition.

Diameter data were relatively easy to obtain using a standard clinical digital ultrasound machine with color M-mode tissue-Doppler imaging. For the GE system we used, pressure data were digitized directly into the auxiliary input of the scanner, and our MATLAB postprocessing program allows direct extraction of both for simultaneous comparison and analysis. The use of transthoracic ultrasound is clearly advantageous over previous, more invasive methods to obtain diameter and has already allowed these techniques to be routinely used in a clinical setting when catheterization is performed.

Our group has previously explored the use of intravascular ultrasound to obtain diameter measurements (43). This previous work recommended the use of compliance measures (the inverse of stiffness) due to considerable measurement uncertainty that can result from a small strain appearing in the denominator of a computed elastic modulus or stiffness. The CMM-TDI method did not display this problem: in our studies, we found that the smallest deformations (≈7% strain for Hy2) were associated with uncertainties on the order of 30% (EPS), using only M-mode images to acquire the wall trace. While this is larger than desired, it is much smaller than what was predicted previously with the intravascular ultrasound modality (43).

In vivo measurement.

To the authors' knowledge, this work represents the first effort to compare in vivo and ex vivo mechanical measurements of the pulmonary vasculature in an animal model of PH. Given our group's recent findings of significant increases in proximal vascular elastic modulus and stiffness in hypoxic PH (22) and preliminary results from application of this method in pediatric PAH patients (19), we hypothesize that this measurement could have significant diagnostic and/or prognostic value. As site-specific proximal measures of elastic modulus, increases in the approximations could be indicative of proximal extracellular matrix remodeling, which is associated with vascular stiffening in the calf model. Stiffening of the proximal vessels should be most detrimental to the vascular windkessel (as opposed to distal vessels) (15) and would significantly change the pulse flow seen by the distal vasculature.

The potential presence of smooth muscle cell (SMC) tone in our in vivo measurements, which has been shown to alter arterial stiffness in muscular arteries (10, 15), could be a major source of the error seen in comparing these to ex vivo measurements, conducted on passive (dead) tissue. To the contrary, several very recent studies have suggested that conduit artery SMC tone does not appear elevated after hypoxia-induced PH and has little impact on proximal mechanics (38, 41). Regardless, such SMC-associated differences between in vivo and ex vivo tissue mechanics, if present, should be most apparent in the hypoxic group. However, groupwise Bland-Altman ETWPS and Eex agreement analyses suggest that any differences are small: the hypoxic group bias was slightly larger (48.9 kPa) than the control group bias (33.1 kPa), although the difference between these biases (15.8 kPa) is smaller than the standard deviation of either group mean. Such difference might be attributed to SMC tone, but could also be due to the infinitesimal displacement approximation inherent in the Lame equation. The impact of proximal SMC tone on proximal mechanics is an important area of future research.

Mechanical studies by our group have shown that proximal PVS increases in both rat (12) and calf (22) models of hypoxia-induced PH. The cellular and molecular mechanisms responsible for these increases are a subject of current research; however, many studies have previously examined histological and morphological features of the proximal pulmonary arteries in disease (37). Briefly, small mammals, such as the rat, primarily display immediate adventitial thickening due to accumulation of cells and collagen upon exposure to hypoxia, while the media responds more slowly. In contrast, large mammals, such as the calf, exhibit early and substantial medial thickening due to accumulation of both elastin and collagen. It is hypothesized that these differences are due to a more complex cellular composition of the conduit vessels in the larger mammals. Much more detail regarding proximal changes in hypoxic PH may be found in a recent review (37).

Clinical application.

Recent clinical studies have focused on the prognostic capabilities of invasive measurement of global stiffness [as assessed with pulmonary vascular input impedance (18) or pulmonary vascular compliance (29)] and MRI-based noninvasive measurement of pulmonary vascular relative area change (RAC) (14). These are all promising prognostics, but only the last is a site-specific measure, like the approximations explored herein. We used our echo-measured diameter alone to estimate RAC, which clearly correlates to our approximations (given the functional dependencies shown in Eqs. 1 and 2); however, RAC is a lesser correlate of Eex (r2 = 0.477, P = 0.003) in our animal model. Given the established link between proximal vascular remodeling and increasing vascular elastic modulus in animal models (7, 11, 12, 2022, 38, 41, 44), we speculate that our approximations, as better correlates to true elastic modulus, could be superior clinical prognostics.

Although coupled in vivo and in vitro mechanics work remains to be done to establish their correlative relationship in humans, preliminary evidence exists that such information may not be necessary for this method to have clinical utility (19). Indeed, the three stiffness prognostics above have been used without such correlation. The ease with which the ultrasound acquisition may be added to the typical clinical catheterization workflow (requiring only 2–3 additional min to obtain an image), the minimally invasive nature of this addition, and the simplicity of calculation (specifically EPS) all enable straightforward clinical adoption.


There are several limitations that must be acknowledged. First, our use of linear regressions to obtain averages of ex vivo modulus and stiffness may not best represent the typical in vivo situation; other possibilities for comparison include the incremental secant elastic modulus or standard or weighted averages of the elastic modulus over the physiological range. However, given our use of regression slopes with the PD data, ex vivo linear regressions were consistent, and the high correlation coefficients (r2) found suggests that an incremental, linear slope was representative across each physiological pressure range. We neglected finite-strain effects and longitudinal deformation in the derivation of our elastic modulus approximation equations; however, our goal here was to obtain an estimate of in vivo stiffness behavior that would represent heart load and could guide clinical diagnosis, not to obtain an exact measurement of elastic modulus, as can be obtained with bench-top measurement methods. Finally, we note that use of PD curves to estimate vascular stiffness is not new; our contribution is the use of ultrasound to assess diameter, which makes the method more practical for serial measurements in animals and minimally invasive in a clinical setting.


Through comparisons to true elastic modulus obtained from bench-top measurements, we have shown that PD data provide a reasonable means to estimate pulmonary vascular elastic modulus, a fundamental material property of the vasculature, in intact calves. Unlike previous work, we have used a noninvasive imaging modality, CMM-TDI ultrasound, to obtain arterial diameter, which should enable these measurements to be easily translated into a clinical environment. Our work provides evidence that simple mathematical models of the vessel are sufficient to extract an estimate of elastic modulus in these animals and that an estimate based on thick-walled theory provides the best accuracy. Finally, given the link between hypoxia-induced PH and proximal stiffening, we speculate that this measurement will be a useful prognostic tool in the treatment and management of human PH.


This project was supported in part by grants from the National Heart, Lung, and Blood Institute (R01-HL067393, T32-HL072738, K24-HL081506, and P50-HL084923) and the American Heart Association (09SDG2260194).


No conflicts of interest are declared by the author(s).


K. S. Hunter expresses gratitude to Maria Frid, Sandi Walchak, and the many others who assisted with various aspects of the calf experiments.


  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.
  42. 42.
  43. 43.
  44. 44.
  45. 45.
View Abstract