AJP - Heart Information on EB 2010
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 292: H3006-H3018, 2007. First published February 16, 2007; doi:10.1152/ajpheart.01012.2006
0363-6135/07 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
292/6/H3006    most recent
01012.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ghorishi, Z.
Right arrow Articles by Eldridge, M. W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ghorishi, Z.
Right arrow Articles by Eldridge, M. W.

Shear stress paradigm for perinatal fractal arterial network remodeling in lambs with pulmonary hypertension and increased pulmonary blood flow

Zahra Ghorishi,1 Jay M. Milstein,1 Francis R. Poulain,1 Anita Moon-Grady,1,2 Theresa Tacy,2 Stephen H. Bennett,1 Jeffery R. Fineman,2 and Marlowe W. Eldridge3

1Neonatology Division, Department of Pediatrics, University of California, Davis and 2Department of Pediatrics, University of California, San Francisco, California; and 3Department of Pediatrics, Population Health Sciences and Biomedical Engineering, University of Wisconsin, Madison, Wisconsin

Submitted 14 September 2006 ; accepted in final form 12 February 2007


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX: ROLE OF ARTERIAL...
 REFERENCES
 
Congenital heart disease with increased blood flow commonly leads to the development of increased pulmonary vascular reactivity and pulmonary arterial hypertension by mechanisms that remain unclear. We hypothesized a shear stress paradigm of hemodynamic reactivity and network remodeling via the persistence and/or exacerbation of a fetal diameter bifurcation phenotype [parent diameter d0 and daughters d1 ≥ d2 with {alpha} < 2 in (d1/d0){alpha} + (d2/d0){alpha} and area ratio beta < 1 in beta = (d12+ d22)/ d02] that mechanically acts as a high resistance magnifier/shear stress amplifier to blood flow. Evidence of a hemodynamic influence on network remodeling was assessed with a lamb model of high-flow-induced secondary pulmonary hypertension in which an aortopulmonary graft was surgically placed in one twin in utero (Shunt twin) but not in the other (Control twin). Eight weeks after birth arterial casts were made of the left pulmonary arterial circulation. Bifurcation diameter measurements down to 0.010 mm in the Shunt and Control twins were then compared with those of an unoperated fetal cast. Network organization, cumulative resistance, and pressure/shear stress distributions were evaluated via a fractal model whose dimension D0 {approx} {alpha} delineates hemodynamic reactivity. Fetus and Control twin D0 differed: fetus D0 = 1.72, a high-resistance/shear stress amplifying condition; control twin D0 = 2.02, an area-preserving transport configuration. The Shunt twin (D0 = 1.72) maintained a fetal design but paradoxically remodeled diameter geometry to decrease cumulative resistance relative to the Control twin. Our results indicate that fetal/neonatal pulmonary hemodynamic reactivity remodels in response to shear stress, but the response to elevated blood flow and pulmonary hypertension involves the persistence and exacerbation of a fetal diameter bifurcation phenotype that facilitates endothelial dysfunction/injury.

pulmonary arterial morphometry; branching complexity; fractals


THE PULMONARY CIRCULATION undergoes a dramatic hemodynamic change at birth, from a state of high pressure/low flow to one of low pressure/high flow (50). This change occurs in conjunction with a decrease in pulmonary arterial resistance (8) and a transformation in the complexity of arterial diameter network organization (3). However, in cases of congenital heart disease with a systemic-to-pulmonary communication, an abnormally increased flow is superimposed on the usual decrease in resistance and commonly leads to development of pulmonary hypertension with functional and structural complications (25). While the mechanisms are not completely understood, the present consensus is that endothelial injury and dysfunction increase vascular tone and reactivity, leading to an increase in pulmonary artery pressure that is initially considered reversible (20). Pulmonary hypertension becomes problematic when vascular wall remodeling leads to the appearance of smooth muscle and alterations of medial wall thickness that encroach upon the lumen diameter (20). Since the forces of pressure and shear stress are known to be transduced and translated into molecular mechanisms that are instrumental in regulating diameter, vascular tone, and the remodeling of the vascular wall (11), further understanding of the factors aggravating these forces within the pulmonary arterial circulation should aid in understanding patterns of reactivity and the demarcation of reversibility in the disease process.

In general, mammalian vascular networks are thought to possess a universal fractal pattern of diameter network organization that appears independent of organ system, body size, and hemodynamic state (32, 33, 5557). However, the fetal pulmonary circulation is extraordinary. Fetal arterial bifurcations appear to possess a hemodynamic phenotype that is distinguished by an area ratio <1 (3), not >1 as typically found in the adult pulmonary circulation (28). This curious vessel diameter topology has a design that mechanically aggravates hemodynamic reactivity per unit blood flow, leading to an elevated driving pressure, a magnified resistance, and the concomitant amplification of shear stress in the course of branching (3). Although intended to be rapidly adapting to an adult pattern via perturbations in blood flow and shear stress in the fetal-neonatal transition after birth, the shear stress-amplifying characteristics of the fetal network design may be precarious to the pulmonary circulation in the presence of an elevated blood flow. Shear stress is a force known to influence the remodeling of vascular diameter and/or vascular wall properties in response to transient and chronic changes in blood flow (35). However, acute elevations beyond critical levels can induce endothelial dysfunction and/or injury (19, 35, 54). Thus the interaction of the fetal phenotype with amplified shear stress secondary to elevated blood flow, or any other source of fetal induced endothelial injury and/or pattern of dysfunction (7), may cause the fetal phenotype to persist or become exacerbated in the presence of elevated blood flow after birth (25).

While the changes in vascular diameter geometry and vascular wall morphology of the pulmonary circulation are well understood, the structural remodeling patterns of arterial network organization in response to the hemodynamic forces of pressure and shear stress are not. We hypothesized that hemodynamic reactivity and network remodeling via the persistence and/or exacerbation of a fetal diameter bifurcation phenotype might mechanically act as a high resistance magnifier/shear stress amplifier to blood flow. Therefore, one purpose of this study was to test the hypothesis that a fetal phenotype of arterial branching persists and/or is exacerbated after birth. We utilized a well-established lamb model of congenital heart disease with increased pulmonary blood flow, in which an aortopulmonary graft is surgically placed in one twin in utero (Shunt twin) but not in the other twin (Control twin) (40, 44). Diameter measurements in lung casts were made to characterize longitudinal branching patterns and design of bifurcations from the left pulmonary artery (LPA) in a Shunt twin and its Control twin, relative to an unoperated fetal lung (3). In addition, as the hemodynamic forces of pressure and shear stress are linked to multiple mechanisms of vascular diameter and wall remodeling, our second objective was to evaluate the fractal geometry, topology, and design of the fetal, Shunt, and Control pulmonary arterial networks as a means of predicting the distribution of pressure and shear stress, and to determine whether or not their respective fractal dimensions conform to a universal design independent of hemodynamic forces.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX: ROLE OF ARTERIAL...
 REFERENCES
 
Cast preparation. The respective studies were approved by animal use protocols issued by the University of California Davis (3) and the University of California San Francisco (40). Lung casts of 8-wk-old twin lambs and a term fetus (140 days gestation, 145 days term) were made from experimental animals originating from previous studies (3, 40). In the aortopulmonary shunt study, mixed-breed pregnant Western ewes with twin lambs (between 137 and 141 days gestation, term = 145 days) were operated on under sterile conditions (44). One twin had an aortopulmonary shunt surgically implanted as previously described (40) (Shunt twin). The other lamb was not operated on (Control twin). As previously described (40), after spontaneous birth the lambs were kept with their mothers and weighed daily, and respiratory rate and heart rate were measured. In Shunt lambs furosemide (1 mg/kg im) was administered daily and elemental iron (50 mg im) was given weekly. Eight weeks after delivery the lambs were euthanized with an intravenous injection of pentobarbital sodium (Euthanasia CII; Central City Medical, Union City, CA) and subjected to bilateral thoracotomy. In one set of twins, lung casts were made of the left pulmonary arterial circulation, while the right pulmonary arterial circulation was harvested for tissue. Lung casts were made as follows: a 4.5-Fr cannula was inserted into the proximal LPA. The airways were expanded with saline under a hydrostatic pressure gradient of 20 cmH2O via a tracheal cannula. The pulmonary arterial circulation was washed free of blood by saline perfusion. To remain consistent with the morphometric diameter measurement methods utilizing gelatin-barium-filled arteries, a high-viscosity methyl methacrylate plastic (Coe Tray Plastic, GC America, Chicago, IL) was injected into the pulmonary arterial system slowly with a syringe over a 2-min period under a maximum pressure of 75 mmHg. The plastic was allowed to polymerize overnight, whereupon the lung tissue was macerated in a 20% KOH bath for 3–5 days. The remaining tissue was washed away gently with distilled water, and the cast was allowed to dry before bifurcation diameter measurements.

Diameter measurement. A large number of diameters within bifurcations of the left pulmonary arterial tree were sampled according to previous methods (3). In the Shunt and Control twins, the samples consisted of diameters of bifurcations along the LPA branch, plus an extensive sampling of bifurcations within the upper cranial lobe. Bifurcation diameters were measured with a methodology described previously (3). Briefly, a video microscope system (Infinivar Video Microscope with Zoom, Infinity Photo-optical, Boulder, CO) connected to a Cohu model 2600 solid-state camera (Cohu, San Diego, CA) and to a Macintosh computer (Cupertino, CA) was used to acquire and analyze images via a program (Object Image) developed by Norbert Vischer of the Faculty of Science, University of Amsterdam, Amsterdam, The Netherlands (http://simon.bio.uva.nl). Object Image is an extended version of the program NIH Image (http://rsb.info.nih.gov/nih-image), which allows diameter measurements and derived calculations to be recorded into an extensive database. To measure diameters, the three-dimensional branching aspect of each cast was broken into pieces to facilitate placement of bifurcations onto a plane. Assuming vessels were circular, vessel diameter was calculated as the average of two diameter measurements per vessel segment. The three averaged diameter measurements, consisting of a parent vessel, d0, and two daughter vessels, d1 and d2, where d0 > d1 ≥ d2, were stored in a database and used to make functional morphometric calculations within bifurcations.

Functional morphometric relationships. The measurements d0, d1, and d2 were utilized to calculate functional morphometric parameters that characterize local hemodynamic reactivity, via diameter geometry, topology, and design of bifurcations within the pulmonary arterial tree (APPENDIX; Ref. 3). Diameter branching geometry is defined by the diameter ratio

Formula 1(1)
along with the asymmetry ratio {gamma} = d2/d1 ≤ 1. The area ratio

Formula 2(2)
designates the topology between parent and daughter diameters in terms of the relative cross-sectional area available for perfusion. The design parameter {alpha} was iteratively solved (28) via the equation

Formula 3(3)
The values of {alpha} and/or beta have hemodynamic significance conferring reactivity to bifurcation network organization (Fig. 1; APPENDIX): a high-resistance/shear stress-amplifying fetal design (3) occurs when {alpha} < 2 and/or beta < 1; a lower resistance transport design (57) is evident when {alpha} ≥ 2 and/or beta ≥ 1.


Figure 1
View larger version (24K):
[in this window]
[in a new window]

 
Fig. 1. Impact of diameter network organization on magnifying cumulative resistance (Rcum) and amplifying shear stress ({tau}) to blood flow with branching order as seen in a symmetric dichotomous fractal arterial network via different area ratios (beta) and fractal dimensions (D0). A: hemodynamic reactivity ({kappa}), defined as the tendency to magnify and amplify properties with arterial branching, has interrelationships with the scaling of diameter geometry (via parent d0 and daughter d1 = d2 vessels), topology [via the area ratio beta = (d

12+ d22)/d02= 2d12/d02], and fractal design D0 = 2/(1 – lnbeta/ln2) according to {kappa} = 2R1/R0 = 2{tau}1/{tau}0 = (beta)3–D0/2–D0(APPENDIX). Hemodynamic reactivity {kappa} spans a spectrum of network designs and topologies that encompass different theoretical laws of network form and function: West-Brown-Enquist law D0 = 2.0 (57); Kurz-Sandau law D0 = 2.7 (32); Murray's law D0 = 3.0 (33, 51). Reactivity is neutral ({kappa} = 1) when the critical area ratio betac = 21/3 and D0 = 3 is achieved; increases ({kappa} > 1) when beta < betac; and decreases ({kappa} < 1) when beta < betac. B: influence of beta on the cumulative resistance/branching order Rcum(k) in successive orders of branching (k = 1, 2,..., 5, designated by different symbols: bullet, order 1; {square}, 2; {triangleup}, 3; {circ}, 4; {blacksquare}, 5), normalized by parent resistance R0: Rcum(k)/R0 = 1 + {kappa} + {kappa}2 + ... {kappa}d. C: influence of beta on amplifying shear stress in successive vessel orders k = 1, 2,..., 5 via {tau}(k)/{tau}0 = {kappa}k.

 
Diameter ordering/analysis of arterial tree. We categorized the parent diameters (d0) of bifurcations into orders within the arterial tree via a diameter ranking method that calculates an average diameter per order, Formula 3k, as previously described (3). The diameter data set for each bifurcation (d0, d1, d2, {alpha}, and beta) was first ranked from the largest parent diameter Formula 30 (designated order 0) down to the smallest diameter, by orders k = 0, 1... n categorized by the overall average diameter ratio Formula 3d according to Formula 3dk ≥ (di/Formula 30) > Formula 3d–(k+1). Here, di is a parent diameter of a bifurcation in the sorted list, d0 is the diameter of the largest (main or left) pulmonary artery, and Formula 3d represents the average diameter ratio, computed from the pool of Rd calculated from bifurcations. Once orders were assigned to each parent diameter, the average values of diameter Formula 3k, logFormula 3k, and logFormula 3k of each order in the arterial tree were calculated. The ranking procedure results in a nonoverlapping range of average parent diameters for each vessel diameter order, where each order is assumed to be statistically independent from another but where the log-transformed values of {alpha}k and betak are assumed to possess a statistically self-similar Gaussian distribution with bifurcation level independent of the average value of {alpha} or beta (3).

Fractal model. The average values of Formula 3k and Formula 3k for parent diameters within orders, from the largest diameter artery Formula 30 down to terminal arteries ~0.010 mm, were used to synthesize an equivalent homogeneous fractal branching tree to summarize and compare the hemodynamic reactive properties of arterial network organization (APPENDIX). The fractal properties of the arterial tree are summarized by global ratios: a network branching ratio Formula 3b, a network diameter ratio Formula 3d, and a fractal dimension D0, derived from the distribution of average Formula 3k and Formula 3k in the following way. To determine the number of vessels with each order of branching, a local branching ratio for each order was determined by

Formula 4(4)
where Rdk is the local diameter ratio between orders derived from

Formula 5(5)
and Formula 5k is the average value of {alpha} for the order k. The number of vessels per order was calculated according to Nd = N0Rb1Rb2... Rbd. The global network branching ratio Formula 5b was determined from the numbers of vessels per order via Nk = N0Formula 5bk, computed via the slope of the linear regression relationship between vessel number N and order

Formula 6(6)
where Formula 6b is the antilog of the slope and the intercept was comprised of the main parent N0 = 1. The global network diameter ratio Formula 6d, determining the decrement in average diameter with order Formula 6k = Formula 60Formula 6d–k and representative of an average of all local diameter ratios Rdk between orders (Eq. 5), was evaluated by linear regression of the logarithm of average diameter Formula 6k vs. order k,

Formula 7(7)
where Formula 7d is the negative of the antilog of the regression slope. The fractal dimension of the model, D0, designating the design of the arterial tree, was then evaluated from the power law relationship between vessel number Nk = N0Formula 7d–kD0 and average diameter Formula 7k = Formula 70Formula 7d–k via the slope of the linear regression of

Formula 8(8)
where the fractal dimension D0 is equal to the negative of the regression slope. The values of the fractal dimension D0, the network diameter ratio Formula 8d, and the network branching ratio Formula 8b in the fetal, Control twin, and Shunt twin lung casts were then used to calculate their respective cumulative resistance and shear stress distribution (APPENDIX; Fig. 1). The cumulative resistance (Rcum) for each order k (0–k+1) down to vessels ~0.020 mm was calculated according to

Formula 9(9)

The hydraulic resistance of the LPA to blood flow was computed via Poiseuille's law,

Formula 10(10)
where µ is blood viscosity (assumed to be 3.0 cP), l0 is vessel length, and Formula 100 is the LPA vessel diameter. Equation 9 assumes the presence of a length-diameter scaling relationship Formula 10k = aFormula 10k{alpha}2 where a and {alpha}2 are constants (30). We assumed that {alpha}2 = 1 (30) and that a = 1.7, where the constant a was evaluated on the basis of human pulmonary data (52) ordered according to the Horsfield method (27). The rectilinear pressure response to an increase in flow (13) was calculated according to

Formula 11(11)
where q is pulmonary arterial flow in the parent vessel normalized to body weight and pint is the intercept pressure, set to 25 mmHg in the fetal state (46) and 6 mmHg in the postnatal state. The associated pressure distribution in the fractal continuum equivalent tree at a steady-state flow qss was evaluated according to

Formula 12(12)
The distribution of shear stress with branching order k was calculated according to

Formula 13(13)
where

Formula 14(14)
and qss is the absolute steady-state blood flow (ml/min) with Formula 140 the diameter of the LPA. Steady-state values for flow q used in the calculations were taken from previously published studies in which LPA flow was measured: 139-day fetal lambs/145 days term (with average body wt 5.0 kg), where q = 9.5 ml·min–1·kg–1 or q = 38 ml/min (46); 8-wk-old Control twins (average body wt 22.6 kg), where q = 33 ml·min–1·kg–1 or 743 ml/min; and 8-wk-old Shunt twins (Ref. 40; average body wt 18.9 kg), where shunt open q =115 ml·min–1·kg–1 or 2,200 ml/min. The Reynolds number per order k associated with a given flow to the network was calculated according to Re(k) = Re(0)Rd(1–D0)k, where Re(0) = q{rho}d0, with the density of blood {rho} = 1.05 g/cm and blood viscosity µ = 0.03 g·cm–1·s–1 (APPENDIX).

Statistical analysis. The distributions of {alpha} and beta exhibit log-normal behavior; therefore, a two-way ANOVA was performed on log-transformed values of {alpha} and beta to evaluate differences between the fetal, Shunt twin, and Control twin data. One factor was the experimental group (Control twin, Shunt twin, fetus), and the other factor was branching order k. Our null hypothesis was that the mean values of {alpha} and beta, over all orders, are equal to a transport design ({alpha} = 2 and beta = 1), considered to be a universal theoretical pattern of arterial vascular network organization (57). Results of the analysis of the log-transformed values of {alpha} and beta are reported as means ± 95% confidence intervals (CI) (58), expressed in terms of their geometric mean adjusted values (37). The regression slopes D0 of Eq. 8 derived from the lung casts are also reported in terms of 95% CI. Significant differences comparing branching order design, topology, and group fractal dimensions were distinguished by nonoverlapping CI, indicative of a P > 0.05 level of significance. Statistical tests were performed with JMP 6.0.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX: ROLE OF ARTERIAL...
 REFERENCES
 
A total of n = 46,673 bifurcation diameter measurements for {alpha} and beta were evaluated in the three lung casts (Control twin n = 16,443 bifurcations, Shunt twin n = 19,816 bifurcations, and fetus n = 10,414 bifurcations), where 1 bifurcation = 1 parent diameter connected to 2 daughter diameters.

Distributions of arterial design ({alpha}) and topology (beta) with diameter and daughter asymmetry ({gamma}). ANOVA revealed significant differences in the averages between experimental groups and within orders for (Formula 14, Formula 14) at a P < 0.0001. In the fetal arterial cast Formula 14 =1.736 (CI 1.66–1.815) and Formula 14 = 0.846 (CI 0.816–0.878); in the Control twin cast Formula 14 =1.94 (CI 1.88–2.00) and Formula 14 = 0.934 (CI 0.911–0.957); and in the Shunt twin cast Formula 14 = 1.671 (CI 1.60–1.744) and Formula 14 = 0.850 (CI 0.821–0.881). These results indicate that, relative to the fetal cast, the control cast altered its bifurcation design and topology toward a West-Brown-Enquist design. However, the Shunt twin demonstrated the persistence of a fetal design and topology.

Figure 2 compares conventional histologically derived indices of vascular-diameter reactivity [percentage of smooth muscle and medial thickness-diameter ratio measured previously in this experimental model (40)] relative to indices of hemodynamic reactivity, average Formula 14 and Formula 14, relative to average diameter and branching order. The fetal-reactive regions (beta < 1, {alpha} < 2) are colored gray in Fig. 2, while the white area (beta ≥ 1, {alpha} ≥ 2) indicates a hemodynamically low-reactive adult transport criterion of bifurcation branching. Statistically significant differences of Formula 14 and Formula 14 between groups within orders according to ANOVA are earmarked within the Shunt twin markers as a white circle within a black circle in Fig. 2. Significant differences between any two groups within an order are characterized by nonoverlapping CI. Overlapping CI were observed for orders 0–6 in the Control twin and Shunt twin and orders 0–2 in the fetus, and they are not shown in Fig. 2. Figure 2 illustrates that the fetal pulmonary circulation exhibits a regional mechanical source of hemodynamic reactivity in intra-acinar arteries [0.020-mm to 0.100-mm parent inner diameter (ID)], a region we previously showed to be devoid of smooth muscle(40). The Shunt twin lamb demonstrated a similar fetuslike design profile of {alpha} < 2 and beta < 1 in preacinar arteries >0.150 mm. In addition, intra-acinar Shunt twin arteries (0.020- to 0.040-mm ID), which we previously showed to demonstrate the appearance of increased smooth muscle and medial wall thickness (40), were associated with significant differences in hemodynamic reactivity (smaller {alpha} and beta) relative to both the Control twin and the fetus.


Figure 2
View larger version (36K):
[in this window]
[in a new window]

 
Fig. 2. Comparison of classic histological indices of vascular reactivity (A) and hemodynamic reactivity (B and C) in 8-wk-old Control twin and Shunt twin lambs, both indices relative to the distribution of vascular diameter. A: vascular reactivity histologically is defined by appearance and/or elevation in % smooth muscle and/or increased % medial wall thickness relative to diameters ordered via accompanying airways classified as terminal bronchiole (TB), respiratory bronchioles (RB), alveolar ducts (AD), and the alveolar wall (AW) (40). B and C: hemodynamic reactivity expressed by degree that the average values of {alpha} (B) and beta (C) [±95% confidence intervals (CI)] deviate from a universal transport design state is demarcated by {alpha} = 2, beta = 1. ANOVA indicated that the average {alpha} and beta were significantly different between groups (fetus, Shunt twin, and Control twin) and between branching orders, indicating regional variations of hemodynamic reactivity. Significant differences within a branching order are indicated by a white dot superimposed upon the Shunt indicator (black dots). Hemodynamically reactive values ({alpha} < 2, beta < 1) are evident before the histological appearance of vascular reactivity in the fetal condition. Differences in {alpha} and beta between fetal and Control twin indicate that the pulmonary arterial circulation adapts to hemodynamic conditions of a high-flow/low-pressure state after birth by remodeling area ratios in a universal direction, toward area ratios equal to 1, which mitigates the hemodynamic reactivity of the fetal phenotype. However, the shunt exacerbates the fetal design in preacinar arteries and in small intra-acinar arteries, whose overall average indicates a persistence of the fetal configuration.

 
The fetal cast demonstrated significantly different average values of diameter asymmetry Formula 14 relative to the Control twin and the Shunt twin: Fetus 0.62 (SD 0.20), Control twin 0.71 (SD 0.18); Shunt twin 0.70 (SD 0.19). Figure 3 shows the interrelationship between area ratio beta, bifurcation design exponent {alpha}, and diameter asymmetry {gamma} as a bivariate frequency distribution. The overall range and spread of beta between the casts are similar and substantial. The fetal lung cast bifurcations show a predominant frequency toward acute asymmetry ({gamma} {approx} 0.5), with area ratios favoring 0.8 < beta < 1 rather than beta > 1. In the Control twin and Shunt twin casts, the greatest frequency of asymmetry occurred in a higher range, 0.6 < {gamma} < 0.8. The area ratios of the Control twin demonstrate their greatest frequencies about 0.8 < beta < 1.2, while the area ratios in the Shunt twin included relatively larger numbers of bifurcations 0.6 < beta < 0.8.


Figure 3
View larger version (23K):
[in this window]
[in a new window]

 
Fig. 3. Area ratio (beta) and asymmetry ratio ({gamma}) of fetus (left) and Control (center) and Shunt (right) twins expressed as a bivariate frequency distribution. Gray scale represents numbers of bifurcations within a particular range of {gamma}, beta, and {alpha}. The equation of the lines for different bifurcation scaling exponents {alpha} is given by Eq. 2 in text.

 
Fractal power law properties. The power law behavior synthesized via Eqs. 47 for the fetal, Control twin, and Shunt twin lung casts is shown in Fig. 4.


Figure 4
View larger version (9K):
[in this window]
[in a new window]

 
Fig. 4. Power law behavior of fractal models of arterial network organization derived from {alpha} and Rd and Rb = R

dD0 as calculated from observed values in Fig. 2. A: arterial diameter decreases with increasing order according to Formula A19k+1 = Rd–k+1Formula A19k. Rd (±CI): fetus = 1.30 ± 0.003; Control twin = 1.31 ± 0.001; Shunt twin = 1.41 ± 0.015, significantly different from fetus and Control twin. The differences in Rd slope between the Shunt twin and the Control twin indicate that the Shunt twin has remodeled and rescaled diameter geometry to larger diameters/order for the first 8 orders of vessels >1.0-mm internal diameter, and to smaller diameters for order < 11 or 0.5 mm. B: N increases with branching order k according to Nk = N0Rbk–1. Rb (±CI): fetus = 1.57 ± 0.03; Control twin = 1.73 ± 0.002; Shunt twin = 1.79 ± 0.01. The fetus is predicted to have a smaller number of arterial vessels than the 8-wk Control and Shunt twins. C: fractal dimension D0 as slope of regression line between Nk = (N0Formula A190D0)Formula A19kD0. Fractal dimensions D0 (±CI): fetus = 1.72 ± 0.033; Control twin = 2.02 ± 0.010; Shunt twin = 1.70 ± 0.018, with the fetus and the Shunt twin with statistically identical fractal dimensions of design. The relative values of Formula A19d, D0, and Formula A19b = Formula A19dD0 define the loci of hemodynamic reactivity inherent in network organization that impacts the distribution of cumulative resistance, volume, and shear stress.

 
In Fig. 4A, the diameter ratio (Formula 14d) of the Shunt twin was significantly different from that of the fetal and Control twin casts: Shunt twin Formula 14d = 1.41 ± 0.015 (CI), fetus Formula 14d = 1.30 ± 0.003 (CI), and Control twin Formula 14d = 1.31 ± 0.001 (CI). While the difference between the fetus and the Control twin was statistically significant, it was small. The relatively small difference between the fetus and Control twin cast Formula 14d values, plus the parallel lines in average diameter with order in Fig. 4A, indicate that growth and development scaled diameters upward with branching order. However, in the Shunt twin the average arterial diameters for orders 0–8 were larger relative to the Control twin or fetus cast, while orders 9–10 of the Shunt and Control twins were similar. Diameters of vessel orders 11 and larger were smaller in the Shunt twin compared with the Control twin. Alternatively, the fractal regression indicates that the intersection between Control twin and Shunt twin occurs at order 11, at a diameter of 0.5 mm. Thus the differences between the Shunt twin and other casts indicate that the Shunt twin remodeled and rescaled diameter geometry about a diameter pivot region of 0.5–1.0 mm, leading to smaller resistances in larger arteries and elevated resistances in smaller arteries.

In Fig. 4B, the branching ratios Formula 14b and their respective 95% CI were statistically different from one another: fetus Formula 14b = 1.572 ± 0.03 (CI), Shunt twin Formula 14b = 1.79 ± 0.03 (CI), and Control twin Formula 14b = 1.73 ± 0.002 (CI).

In Fig. 4C, the fractal dimension ± CI of the fetus was D0 = 1.72 ± 0.069, Shunt twin D0 =1.70 ± 0.038, and Control twin D0 = 2.02 ± 0.021. Thus the Shunt twin was not statistically different from the fetus, but both were both significantly different from the Control twin, indicating the persistence of a fetal fractal network design in the arterial tree at 8 wk of shunting.

Together these results indicate that the influence of a chronic shunt interrupts an apparent transition in arterial network organization at birth: the fetal-Control transition maintained the scaling of arterial diameter geometry uniformly increasing diameters but remodeled arterial design and topology to an adult transport design. Shunting remodels diameter geometry by rescaling large and small diameters to where large arterial vessels increase diameter to decrease resistance and small arterial vessels decrease diameter to increase resistance, while leaving the topology and fractal design of the pulmonary circulation in a hyperreactive fetuslike configuration.

Cumulative resistance distribution. Figure 5 illustrates the influence of fractal network organization in the fetal, Control twin, and Shunt twin casts on the cumulative resistance distribution versus diameter, as calculated from the regression values of Formula 14d and D0 according to Eqs. 9 and 10. In general, the smaller fractal dimension (D0 = 1.7) of the fetus and the Shunt twin demonstrates greater resistance reactivity with arterial branching, evidenced by a larger slope of cumulative resistance with decreasing diameter. However, the Shunt twin demonstrates an absolute resistance that is less than the fetus and the Control twin: the remodeling of arterial geometry (R0, Formula 14b) counteracts the increased hemodynamic reactivity of a smaller fractal dimension by decreasing the absolute resistance per order, leaving a smaller cumulative arterial resistance.


Figure 5
View larger version (13K):
[in this window]
[in a new window]

 
Fig. 5. Cumulative resistance of the fetus, Control twin, and Shunt twin fractal network in terms of diameter (d0) or order calculated from power law relationships of Fig. 4 according to Eq. 9. The Shunt and fetal lambs have similar slopes of Rcum change/order relative to the Control lamb because of smaller fractal dimensions of design. The values of R0 and Rd determine the vertical offsets of the Rcum graph: the fetus has the largest cumulative resistance with branching because of the size scale of the left pulmonary artery (LPA) and increased reactivity due to the smaller fractal dimension with arterial branching. The Shunt twin has a smaller cumulative resistance for all parent diameters of branching, which arises from the geometric scaling leading to a larger LPA diameter d0 = 2.12 cm and larger diameter ratio (Rd = 1.41) increasing diameter in orders 0–10 (Fig. 3A) relative to Control (LPA d0 = 1.15 cm and Rd = 1.3) and fetus (LPA d0 = 0.53 cm and Rd = 1.31). This graph indicates that by 8 wk hemodynamic forces of hypertension and increased blood flow of the shunt remodel vascular arterial diameter geometry (R0 and Rd) in a manner that paradoxically decreases the cumulative resistance relative to control (44), but maintains the hemodynamic resistance reactivity of the fetal state (slope), via the persistence of the fetal topology and design (fractal dimensions of fetus = Shunt).

 
Pressure-flow curves/pressure distribution. Figure 6 illustrates the influence of Rcum on the predicted pressure-flow curves (Fig. 6A) of the LPA circulation in the fetal, Control twin, and Shunt twin fractal configurations and (Fig. 6B) the corresponding distribution of pressure versus diameter at a given steady states of flow. The fetal state is the most pressure reactive, yielding hypertensive pressures in all arterial regions over a small flow range. Fractal network organization in the Shunt twin configuration is the least pressure reactive to flow of the three states, with hypertension arising secondary to the increased steady-state blood flow. As Fig. 6 indicates, the Control twin state also predicts pulmonary hypertension at a flow comparable to the Shunt twin, but with higher driving pressures and higher values in vessels >0.060 mm. Conversely, arterial network remodeling in the shunt condition predicts that shunt closure will result in nonhypertensive pressures with smaller pulmonary vascular resistances than the corresponding Control twin.


Figure 6
View larger version (19K):
[in this window]
[in a new window]

 
Fig. 6. A: influence of Rcum on pressure-flow relationships and pressure distribution vs. diameter. Rcum influences the slope of pressure-flow curves determining the pressure reactivity to a given increment in flow about its steady-state value [LPA steady-state flows Qlpa: fetus = 10 (46), Control twin = 33 (40), Shunt twin = 100 (40)]. The fetus has the greatest reactivity, followed by the Control twin and then the Shunt twin. Paradoxically, the remodeling of diameter geometry, topology, and design in the Shunt twin decreases the pressure response to all increments in flow. Increasing steady-state flow in the pulmonary arterial circulation to some threshold results in a clinical definition of pulmonary hypertension (Ppa = 25 mmHg). Structural remodeling in the Shunt twin elevates the steady-state flow threshold to 80 ml·min–1·kg–1, compared with the lower threshold of 50 ml·min–1·kg–1 in the Control twin. B: comparison of steady-state pressure-diameter profiles computed at low flow (33 ml·min–1·kg–1) and high flow (100 ml·min–1·kg–1) in the Control and Shunt twins. The structural remodeling in the Shunt twin lamb diminishes the hypertensive pressure stimulus compared with the Control twin at high flow.

 
Shear stress distribution. Figure 7 illustrates the influence of diameter network organization on the amplification of shear stress at steady-state blood flows in the fractal model. The calculated values for {tau}0 shear stress and Reynolds numbers in the LPA for the different states were as follows: fetus, LPA d0 = 5.4 mm at a flow of q = 38 ml/min yields {tau}0 = 1.25 dyn-cm–2, Re(0) = 82; Shunt twin with the shunt-open condition, LPA d0 = 21.2 mm and flow q0 = 2,200 ml/min yields {tau}0 = 1.16 dyn-cm–2 at Re(0) = 774; Control twin, LPA d0 = 11.5 mm and q0 =743 ml/min yields {tau}0 = 2.51 dyn-cm–2, Re(0) = 486. The common feature in all groups is that extremely large shear stress avalanches are developed with low Reynolds numbers (Re < 10) in arteries of 0.020 mm. For comparison purposes, the gray area in Fig. 7 represents the reported ranges of critical and/or erosion yield shear stresses (165–8,000 dyn-cm–2) of endothelial dysfunction/damage to an acute increase in shear stress as it occurs in large-diameter conduit vessels over a range of Reynolds numbers (0.21–2,300) (19, 35, 54). Interestingly, the fetal fractal model ({tau}0.020 = 1,200 dyn-cm–2 at Re0.020 = 1.4) and Control twin model ({tau}0.020 = 1,100 dyn-cm–2 at Re0.020 = 0.43) demonstrated similar magnitudes of shear stress experienced in arterial vessels of 0.020-mm ID as derived from their respective steady-state blood flows. However, the Shunt twin steady-state values of blood flow yielded a magnitude of shear stress at {tau}0.010 =7,100 dyn-cm–2 at Re0.020 = 2.6, which extends into the upper region of thresholds reported to induce endothelial injury.


Figure 7
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 7. Predicted arterial shear stress distribution in fetal, Control twin, and Shunt twin fractal models. Shear stress amplified with branching order k according to {tau}k = {tau}0Rd(3–D0)k using steady-state blood flows in the parent arteries to solve for shear stress in the parent {tau}0 = 32µq0/{pi}d03. Here, shear stress reactivity, the increase in shear stress with branching order, is determined primarily by the degree that D0 deviates from D0 = 3 (Murray's law). If the distribution with order is {alpha} = D0 = 3, then shear stress is expected to be constant, with branching and equal to {tau}0. Otherwise, if {alpha} = D0 < 3, shear stress is amplified in the course of branching, to a degree determined by how small {alpha} = D0 becomes. A: gray area demarcates reported ranges of critical values of acute shear stress within large-conduit arteries that induce endothelial dysfunction as a permeability change or endothelial mechanical injury (19, 35, 54). At steady-state flow rates, the predicted shear stress in intra-acinar arteries in the fetus and the Control twin at 8 wk are predicted to be similar. The shear stress in the Shunt twin is predicted to be greater for all arterial diameters. B: comparison of Shunt and Control twin shear stress distributions at corresponding steady-state flow rates, high-flow equivalent (100 ml·min–1·kg–1) vs. low flow (33 ml·min–1·kg–1). Increasing or decreasing flow increases/decreases amplitude of shear stress in small vessels, but amplification is always greater with smaller fractal dimensions (Shunt twin > Control twin) for any given flow rate.

 

    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX: ROLE OF ARTERIAL...
 REFERENCES
 
The present study is the first to demonstrate the existence of an unusual fetuslike pattern of arterial diameter network organization after birth in an experimental model of pulmonary hypertension with increased pulmonary blood flow. Persistence and/or exacerbation of the fetal pattern after birth suggests a mechanical basis of hemodynamic reactivity that may be linked to the origin of the disease and implicates a shear stress paradigm of structural remodeling in the pathology and progression of pulmonary hypertension in high-flow lesions. Although this study is limited to a morphometric analysis based on a sample of three lung casts only and not a population, our results indicate that a fetal phenotype may be a major factor that influences the state of structural remodeling in intra-acinar vessels when viewed in the context of classic grades of vascular remodeling determined by vessel geometry alone and is involved in two divergent patterns of adaptation in arterial network organization in a model of neonatal congenital heart disease with increased pulmonary blood flow.

Previous work regarding pulmonary vascular abnormalities in patients with pulmonary hypertension has focused primarily on histological appearance of the arterial wall. In 1958, Heath and Edwards (24a) described a classification sequence that ranged from reversible medial hypertrophy to terminal changes such as angiomatoid formation and fibrinoid necrosis. Rabinovitch et al. in 1978 (42a) described a modification that featured alterations in normal remodeling and growth in three progressive stages (A–C), characterized by progressive medial hypertrophy, abnormal muscle extension, and reduced artery size and concentration, that have been correlated with perioperative hemodynamics. The histological pattern of vascular remodeling in the Shunt twin at 8 wk of age occurring in this experimental model has been previously characterized as grade B alterations, with medial hypertrophy and abnormal appearance of muscle in the walls of intra-acinar arteries (40). A fundamental limitation of histologically derived indices of vascular reactivity is that it is not clear whether their impact on arterial diameter comprises a mechanical cause of pulmonary hypertension or represents an adaptation in response to it (47). Alternatively, our study was concerned with further characterization of the high-flow secondary hypertension experimental model with specific regard to the influence of elevated flow complicating the structural remodeling of a preexisting fetal phenotype of arterial diameter network organization to a postnatal adult pattern in the Shunt twin, but not in the Control twin. In this study we show that a fetal phenotypic source of hemodynamic reactivity ({alpha} < 2, beta < 1) is evident in arterial diameters of fetal and Shunt twin intra-acinar arteries. We note that fetal intra-acinar arteries are typically devoid of smooth muscle (2, 41), in which case the fetal phenotype represents a preexisting mechanical source and state of aggravated hemodynamic reactivity contributing to the high-resistance, hypertensive condition normally seen in the late-term fetal pulmonary circulation (50). Consequently, as Fig. 1 indicates, we would expect that the structural remodeling of the fetal phenotype to smaller area ratios would be of greater hemodynamic significance to the progression and timing of hypertension than that expected from the structural remodeling and encroachment of vessel diameter alone fixed at a given area ratio (47). Moreover, because the unusual fetal topological pattern between diameters represents a source of hemodynamic reactivity not considered in the classic vascular remodeling perspective, it is likely that changes in the topological form of hemodynamic reactivity in the fetal-neonatal state may precede the timing and appearance of more complex lesions after birth (43). Indeed, preferential patterns of structural remodeling of bifurcation topology and design in relationship to the appearance of smooth muscle and wall thickness are seen in our data in two ways. First, as seen in the 8-wk Control twin lung cast (Fig. 3), the appearance of a normal pattern of smooth muscle is related to the timing of significant increases in the area ratio. In this case, the migration of {alpha} and beta toward values of {alpha} {approx} 2 and beta {approx} 1, and to the fractal dimension of D0 {approx} 2, indicates that the degree of hemodynamic reactivity is significantly diminished relative to that of the fetus. Under these circumstances, the structural and functional impact of the residual fetal configuration in the Control twin on overall hemodynamic reactivity at this point may be comparatively minimal relative to the emerging influence of smooth muscle on vascular reactivity and tone. Second, Fig. 3 also illustrates that the abnormally increased muscle mass in Shunt twin intra-acinar vessels (0.020- to 0.040-mm ID) relative to the Control twin is associated with diameters of parent/daughter vessels with significantly smaller {alpha} and beta, relative to both the fetal and Control twin lungs. Under these circumstances, the increased muscle mass and the aggravated hemodynamic reactive design and topology indicate a preferential action of smooth muscle vascular remodeling that can exacerbate the fetal phenotype in the course of hypertension disease. Thus the regional pulmonary patterns of {alpha} and beta we observed in this study emphasize that there is a preexisting mechanical source of hemodynamic reactivity and remodeling in the fetal pulmonary circulation that extends beyond the apparently latent remodeling patterns of vascular wall properties at a vessel level of organization.

The conceptual advantage of the fractal model is that it simplifies the complexity of arterial network organization into a simple static parametric characterization of arterial geometry (Rd), topology (area ratio beta), and design (fractal dimension D0), whose impact on hemodynamics can be inferred from an equivalent symmetric bifurcation-branching arterial network (24, 30). The fractal dimension and/or area ratio have structural and functional significance in that their values define the global hemodynamic reactivity of the arterial tree, along with the pattern of transformation of pressure and shear stress with arterial branching (Fig. 1; APPENDIX) (4, 30). Our morphometric results demonstrated differences in the fractal dimension between the fetus (D0 = 1.72) and the Control twin (D0 = 2.02), suggesting a transformation in fractal arterial network organization at birth; furthermore, our results indicate that this transformation did not take place in the Shunt twin (D0 = 1.70). However, we emphasize that our study is limited in the sense that it is based on the results of three lung casts only, which are isolated samples of a term lamb fetus and a postnatal state at 8 wk of birth. Unfortunately, the limited number of lung casts we studied does not identify the exact timing of when the adult bifurcation transport design is achieved after birth. Also, the limited number of lung casts we utilized may not necessarily reflect the developmental dynamics of change in network properties representative of the population of human subjects or other species (39) experiencing the same hemodynamic conditions. Here, genotypic variation and heterogeneity in patterns of vascular wall remodeling may exist within a population of individuals and between species that also influences the distribution and reactivity of the fetal phenotype of arterial branching and its responsiveness to hemodynamic forces for arterial network structural remodeling (53). In addition, at a given developmental time period lung casts are static pictures of network organization that neglect the dynamic information about the influence of vasomotor activity on bifurcation constriction or dilation. Here it should be appreciated that the actions of vasomotion and vascular tone are now more complex beyond their influence on vascular diameter (12, 17), because they also modulate bifurcation design and topology, which can also influence the fractal dimension of design. Thus a spectrum of hemodynamic reactivity in network organization may be present at different times under different hemodynamic circumstances in the perinatal period than those observed in this study. Despite limitations, the pervasive nature of the unusual fetal design phenotype and its divergent directions of structural remodeling seen here after birth indicate that persistence and/or aggravation of the fetal phenotype after birth may be important in developing further understanding of the spectrum of abnormal patterns in hemodynamic reactivity and vascular wall remodeling seen in patients with pulmonary hypertension secondary to elevated blood flow (43).

The fractal dimensions observed in the fetus and the Shunt twin are curious because they are at odds with what is predicted to be universal for arterial network organization in mammalian vascular systems in general (32, 33, 5557). While hemodynamic forces, such as cyclic or mean pressure and shear stress secondary to increased flow, are stimuli known to influence endothelial dysfunction and vascular wall remodeling (10), various theories predict that arterial vascular network organization maintains a universal design, or a target configuration, with arterial branching ratios always beta ≥ 1, with the condition of beta < 1 improbable. An assertion made by West-Brown-Enquist (57) is that arterial network organization has a universal design based on a widely known metabolic scaling principle in adult species, which uniquely links the fractal dimensions of arterial network organization to the exponent of the 3/4 law metabolic scaling law [Kleiber's relationship (29)]. Their theory maintains that a fractal design of D0 = 2 and a diameter topology of beta = 1 exist for all arterial transport and capillary exchange networks, and that this design remains independent of animal size, ontogenic state, and hemodynamic condition (56, 57). Another classic design, called Murray's law, dictates that arterial network organization has a universal configuration of D0 = 3.0, a condition optimized for constant shear stress and resistance with arterial branching (33). Alternatively, Kurz and Sandau (32) predicted that the bifurcation scaling exponents of arterial network organization remodel during ontogenesis and converge to a target configuration of D0 = 2.7, a design indicative of a vascular system of constant wall tension. Whatever the target, the assumption of a fixed, universal design independent of hemodynamic conditions is relevant and crucial to the study of structural remodeling in pulmonary hypertension. Design invariance with respect to hemodynamic conditions is an operational premise that serves to justify the idea that the complexity of network organization is not a precursor to the generation of hemodynamic forces that are instrumental to the remodeling process of pulmonary hypertension disease. In addition, network design invariance provides a necessary theoretical justification that the complexity of the pulmonary arterial system can be exclusively reduced to a vessel level of organization to elucidate mechanisms controlling diameter restriction in pulmonary hypertension. What the theories of arterial network organization overlook is that the fetal pulmonary circulation does not function as a generalized transport/exchange organ system in utero. In addition, hemodynamic conditions in the fetal pulmonary circulation are also much different from those in the target adult pulmonary circulation: the fetal pulmonary bed is hypertensive in relation to the systemic bed and is subject to a relatively small blood flow (50). At birth the fetal configuration switches in response to an increase in blood flow, via the removal of the placenta and the closure of the ductus arteriosis and foramen ovale, in conjunction with a decrease in capillary resistance with airway expansion. Our morphometric results indicate that the hemodynamic transition at birth is associated with an important, necessary transformation in arterial network organization, such that the fetal configuration functions as a rapidly adapting high-impedance transformer, which is capable of dynamically altering its hemodynamic reactivity in response to shear stress perturbations via changes in blood flow.

Although our data are limited in not allowing us to generalize to the fractal state of arterial network organization in other fetal and neonatal states as a population or conditions of hypertension, we can address the general question of whether fractal network organization is independent of hemodynamic conditions in another way, via the influence of fractal network organization on the scaling of cumulative arterial resistance and slopes of pressure-flow curves relative to body weight. As Eq. 12 indicates, Rcum functionally dictates the slope of the pressure-flow curve (13), with increases in Rcum (Eq. 10) accomplished via smaller fractal dimensions and larger diameter ratios Rd, while a compensating decrease in Rcum can be accomplished via a smaller R0. However, arterial resistance Rcum is predicted to scale with body weight independent of species, state of development, and hemodynamic condition via a fixed fractal dimension according to different patterns: the West-Brown-Enquist law predicts resistance to universally scale according to Rcum {propto} M–3/4 with arterial D0 = 2 (57), while Murray's law predicts Rcum {propto} M–1 with arterial D0 = 3 (51). Consequently, if vascular network organization is independent of hemodynamic conditions and state of pulmonary development, then the slopes of pulmonary pressure-flow curves sampled from different species, ages, and conditions of pulmonary hypertension would be anticipated to be devoid of outliers and invariant about a given power law. Figure 8 summarizes the scaling of Rcum estimated via the slope of pulmonary pressure-flow curves drawn from the literature (1, 15, 16, 22, 31, 36, 38, 42, 45) and indicates that the resistance scaling law for the pulmonary circulation of the normotensive adult does not adhere to a Rcum {propto} M–3/4 as predicted by the West-Brown-Enquist theory, but instead migrates to an M–1 power law only in the adult case free of pulmonary hypertension. At the same time, under conditions of pulmonary hypertension in the fetal-neonatal state (8, 46) and the adult state with primary pulmonary hypertension (36), the deviations about the M–1 power law indicate that resistance does not adhere to any particular law at all. Under these circumstances, the degrees of freedom possible in a fractal pattern of network organization to influence resistance suggest that the elevated cumulative resistance seen in the fetus and adult with pulmonary hypertension may possibly share a common topology (beta < 1) and design (D0 < 2) in arterial network organization. At the same time, if arterial network organization is, in general, not independent of hemodynamic conditions, then there must be a common hemodynamic feedback signal to condition arterial design to target the fractal dimensions of arterial network organization from one resistance state of organization to another. As for targets and their feedback signals, LaBarbera (33) proposed constant shear stress for a target of D0 = 3 while Kurz and Sandau (32) theorized constant wall tension for a target of D0 = 2.7. Regardless of the signal and the target design, intermediate states of Rcum and its associated fractal dimension seen in Fig. 8 can be understood exclusively in terms of a network design that is dynamic, with a capacity to change its organization to act as a variable hemodynamic transformer (3, 4) that changes its reactivity under different hemodynamic conditions.


Figure 8
View larger version (22K):
[in this window]
[in a new window]

 
Fig. 8. A: Rcum power laws. Adult pulmonary cumulative resistance in absence of hypertension adheres to Rcum {propto} M–1 power law (51), not a universal Rcum {propto} M–3/4 law (57). Pressure-flow data for Rcum derived from Refs. 1, 15, 16, 22, 31, 36, 38, 42, 45. B: Rcum does not adhere to an invariant power law over developmental/physiological timescales. Rcum·M {propto} M–0 median trend line (0.08 mmHg·kg·ml–1·min–1) is value of impedance predicted by universal power law and allometric result of a mean pressure of 15 mmHg divided by allometric scaling of cardiac output with body size. Outliers to the law include a vertical allometric comparison between species at same body weight demonstrating significantly different resistances (39), not a singular impedance scaled to M–3/4 design. Also, high-resistance outliers to Rcum {propto} M–1 occur in fetal hypertension states (8, 46) and in the adult with primary pulmonary hypertension (PPHN) (36).

 
One major pattern of structural remodeling observed in this study was the transformation in fractal dimension between the fetal and postnatal states (fetus and Control twin), which was not seen in the Shunt twin. There is likely more than one contribution responsible for the structural remodeling of arterial network organization. As the hemodynamic transition at birth is associated with changes in vascular tone (17), we suspect that vascular tone is involved in remodeling of bifurcation design. Indeed, Griffith et al. (23) demonstrated in isolated bifurcations that vasodilation selectively increases the area ratio in bifurcations from a baseline state of {alpha} = 2 to {alpha} -> 3, a constant shear stress design, via a nitric oxide (NO) mechanism, while vasoconstriction from a baseline state to a fetal design of {alpha} -> 1.25 and beta < 1 is induced by endothelin (ET)-1. Another factor expected to influence the structural remodeling of arterial design and topology is the influence of shear stress. In general, the endothelium responds to shear stress by altering its functional phenotype and endothelial gene expression via several thousand genes over different timescales (21). The endothelium incorporates machinery for a rapid transduction-translation response in function and geometry coupled to shear stress (11), which includes potassium channels (9), vasodilation via NO and ET-1 (5, 18), and possibly lung Krüppel-like transcription factors to gradually remodel the state of tone (14). Over hours to days shear stress is capable of signaling a spectrum of remodeling changes within the vessel wall (endothelium, intima, smooth muscle, and adventitia) working as a syncytium that add or detract from vessel wall mass to either encroach or increase diameter (21, 34). Thus, insofar as the unusual configuration of the term fetal state is already preconditioned to hypertension, and small increments in flow lead to large transient avalanches of shear stress, causal reasoning implicates that it is the transduction of shear stress via the endothelium and its translation of endothelial/smooth muscle function to reorganize bifurcation topology, design, and hemodynamic reactivity that is a principal hemodynamic stimulus for remodeling to an adult configuration at birth, rather than the vasodilation of vascular diameter alone (17). We also suspect that shear stress is implicated in the abnormal structural remodeling of arterial diameter geometry, topology, and design in the Shunt twin. In postnatal experimental models of congenital heart disease, where the reduction in resistance has occurred in conjunction with an adult pulmonary arterial configuration, we would have anticipated that the chronically elevated pressure would have increased pulmonary vascular resistance (6, 48). However, while the reactive branching pattern persisted and was aggravated in smaller arterial vessels, the Shunt twin demonstrated a paradoxical decrease in absolute cumulative resistance relative to the Control twin, which is consistent with our previous findings (40, 44). In this example, the decrease in resistance can be attributed to geometric remodeling leading to a significant increase in larger diameter vessels. Figure 5 illustrates the regional influence of arterial diameter scaling in the Shunt twin relative to the fetus and the Control twin: about a pivot point of order 11, orders 0–8 possess larger diameters and orders 9–11 approximate Control twin diameters, while orders 12–22 are smaller than Control twin diameters. The impact of the increased diameters in the larger arteries is to increase arterial volume, relative to control, a common finding in this model and in patients with chronic high-flow hypertension (44). This geometric feature of structural remodeling may have arisen from the increased flow, rather than the pressure alone, via a negative feedback mechanism of diameter structural remodeling dependent on an endothelial shear stress mechanism, mediated by NO (34). However, the persistent pattern of fetal arterial network organization in the Shunt twin in combination with elevated blood flow suggests a potential initiating physical-mechanical cause and site of action of endothelial injury/dysfunction due to amplified shear stress (25). In Fig. 7 the flow in the Shunt twin may be elevated to a sufficient extent to lead to excessive chronic levels of shear stress known to transiently induce permeability changes and mechanical injury (19, 35, 54). Unfortunately, our study is limited in that it does not identify whether the amplitude of shear stress alone is an effective cause of the persistence of the fetal bifurcation design. Alternatively, other factors operating within the fetal state besides shear stress-induced endothelial dysfunction and/or injury, such as the disruption of endothelial/smooth muscle function and communication (10), may be instrumental effective causes interrupting the fetal-neonatal transformation to a low-impedance transport design in which elevated shear stress is a complication. Under these conditions, amplified shear stress itself may represent an additional complication factor that interferes with the structural remodeling of the fetal bifurcation design to an adult configuration. We submit that the nature and direction of structural remodeling in arterial networks with area ratios <1 under a shear stress paradigm may have implications for understanding the origins and course of vascular reactivity in pulmonary hypertension in the perinatal period.

In summary, our results indicate that the unusual hemodynamic reactivity of the perinatal pulmonary arterial circulation in response to elevated blood flow leading to pulmonary hypertension can be potentially explained via a shear stress paradigm of arterial diameter network remodeling. In an elevated pulmonary blood flow lamb model, the increased resistance and reactivity predicted in pulmonary hypertension has a preexisting mechanical basis evidenced by a fetal diameter branching phenotype that acts as a transformer/amplifier of hemodynamic forces that persists and/or becomes exacerbated after birth. This fetal phenotypic pattern appears before the emergence of abnormal vascular histological patterns of remodeling seen in the course of secondary pulmonary hypertension and serves to maintain and/or exacerbate the hemodynamic reactivity of the fetal phenotype. The levels of shear stress predicted in preacinar and intra-acinar arterial vessels in our model indicate that the critical levels of shear stress leading to endothelial dysfunction/injury may arise mechanically from the interaction of blood flow and the fetal phenotype of arterial network organization itself. While the actual mechanisms responsible for arterial network remodeling and adaptation in the perinatal period are not yet identified, our shear stress paradigm suggests possible remodeling pathways of arterial network diameter topology and design that extend beyond present understanding of a stationary pattern of arterial diameter geometry under a transmural pressure paradigm of vascular remodeling diameter and wall properties. Given the pressure-flow properties of the pulmonary circulation in patients with other forms of pulmonary hypertension, and the theoretical dependence of the area ratio on hemodynamic reactivity, we speculate that diverse forms of pulmonary arterial hypertension may share a common origin and formal cause embedded in replication of a high-impedance fetal pattern of branching complexity.


    APPENDIX: ROLE OF ARTERIAL FRACTAL NETWORK ORGANIZATION ON HEMODYNAMIC REACTIVITY
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX: ROLE OF ARTERIAL...
 REFERENCES
 
Our present concept of vascular reactivity is based primarily on factors that influence vascular diameter impacting vessel hemodynamic resistance via Poiseuille's law via Eq. 10. However, an additional source of hemodynamic reactivity arises from an influence of arterial network topology and design, which affects the distribution and scaling of hemodynamic forces via the manner in which diameters are interconnected. The purpose of this APPENDIX is to summarize the functional morphometric relationships of topology and design used to characterize hemodynamic reactivity in terms of dimensionless constants, such as the area ratio beta, the bifurcation scaling exponent {alpha}, the fractal dimension D0, and the resistance/shear stress reactivity constant {kappa} with arterial branching.

Bifurcation topology. The area ratio expresses the cross-sectional area of perfusion with arterial branching that can be approximated in terms of network topology and geometry

Formula A1(A1)
This relationship states that the change in cross-sectional area of perfusion to blood flow is determined by the product of the relative increase in the number of vessels with arterial branching (via the branching ratio Formula A1b = Nk+1/Nk) and the relative decrement in average diameter (Formula A1d = Formula A1k+1/Formula A1k), where the subscripts k and k+1 refer to the order of branching, with k the parent to daughters k+1. A bifurcation is the simplest network topology in which the area ratio can be explicitly expressed in terms of diameter geometry. In general, arterial networks branch asymmetrically: branching typically consists of a parent diameter d0 and daughter diameters d1 and d2 such that d0 > d1 ≥ d2. Here bifurcation summarizes the bifurcation's topology via a geometric relationship

Formula A2(A2)
where the local bifurcation diameter ratio is defined as Rd = d1/d0 and the asymmetry ratio is defined as {gamma} = d2/d1. Thus, in the context of the definition in Eq. A1, the diameter ratio and daughter asymmetry of the bifurcation are associated with hemodynamically equivalent average branching ratio and diameter ratio. In the case of a symmetric bifurcation {gamma} = 1 the area ratio has a simple topological-geometric interpretation, as

Formula A3(A3)
where Formula A3b = 2 and Formula A3d = Rd. Alone, the area ratio alone says nothing about its hemodynamic impact in response to flow.

Bifurcation design. A structure-function scaling relationship between flow and diameter within a vessel leads to the concept of a hemodynamic design for a given diameter network topology that can be used for hemodynamic analysis and prediction. Roy and Woldenberg (49) argued that the steady-state flow (q) in the vasculature has a general scaling relationship with diameter (d) according to

Formula A4(A4)
where k is a constant. The value of the scaling constant {alpha} has hemodynamic significance, imparting a design or plan for the steady components of flow and pressure that are optimal in some sense (3). In the case of a bifurcation, the conservation of flow dictates that q0 = q1 + q2. With the fraction of flow distributed in daughter vessels according to (p1 = q1/q0, p2 = q2/q0 = 1 – p1) the flow-diameter scaling relationship in the bifurcation is

Formula A5(A5)
In the case of an asymmetric bifurcation the area ratio is related to the diameter (49)

Formula A6(A6)
where the value for a symmetric bifurcation network is noted by {gamma} = 1.

The value of {alpha}, although dimensionless, has hemodynamic significance that reflects an underlying law or principle associated with the influence of topology of network organization on the distribution of blood flow (3): {alpha} = 4, beta > 1 confers constant resistance with branching; {alpha} = 3, beta > 1 constant shear stress (Murray's law); {alpha} = 2.7, beta > 1 constant wall tension (Kurz-Sandau law); and {alpha} = 2, beta = 1 constant velocity (West-Brown-Enquist law). The hemodynamic influence of a given value of {alpha} also impacts the scaling of hemodynamic properties associated with flow q, such as the Reynolds number (Re = q{rho}d), flow velocity ({nu} = q/{pi}d2), and/or shear stress ({tau} = 32q/{pi}d2). In essence, arterial network organization can influence a given hemodynamic property depending on the value of the design parameter {alpha}. In a symmetric bifurcation (Rb = 2 and Rd = 21/{alpha}):

Formula A7(A7)

Formula A8(A8)

Formula A9(A9)
These indices ignore the development of flow, the formation of fluid dynamic boundary layers, and flow-diameter asymmetry. However, what these scaling relationships have in common is that there is a critical value of {alpha}crit that influences whether network topology either attenuates or amplifies a given hemodynamic property.

Network topology and design. Horsfield and Woldenberg (28) argued that, akin to flow-diameter scaling in a bifurcation, there is a fractal flow-diameter scaling relationship in arterial trees. The number of branches in successive orders down a fractal tree, such as the pulmonary circulation, increases by a constant proportion (Formula A9b = Nk+1/Nk), while the mean flow in branches in successive orders will decrease by a proportion Formula A9 = Formula A9k+1/Formula A9k that is inversely proportional to the branching ratio, Formula A9b = Formula A9–1. The scaling relationship between the average decrease in diameter with branching, Formula A9d = Formula A9k+1/Formula A9k in conjunction with the vessel number or inverse fraction of flow distributed, leads to the seed of a fractal power law scaling relationship,

Formula A10(A10)
If Rb and Rd are constant with branching order k, then the arterial network can be formally considered as self-similar, with the scaling constant {alpha} approximately equal to a fractal dimension {alpha} {approx} D0 (30). Under these circumstances, the topology network organization via the area ratio is related to the fractal network scaling property according to

Formula A11(A11)
Consequently, network organization and design has an exponential scaling influence with branching order k on the attenuation/amplification properties according to

Formula A12(A12)

Formula A13(A13)

Formula A14(A14)
The idea that network organization influences hemodynamic properties in response to blood flow suggests that the diameter scaling properties can be summarized in terms of a reactivity parameter {kappa} that demonstrates its dependence on arterial topology and design with arterial branching. Thus, for shear stress reactivity, the relative increment of shear stress from the first order 0 to the first branching level 1 is given by

Formula A15(A15)
Figure 1 illustrates the relative influence of beta and D0 for {kappa} relative to one level of branching, as well as the amplification effects over several orders of branching.

Resistance reactivity. Under the assumptions that Poiseuille's law applies (Eq. 10), the cumulative resistance of a complex network depends on how resistance changes with arterial branching

Formula A16(A16)
where Formula A16k is the average resistance (Eq. 10) of the vessels of order k with average diameter Formula A16k and length Formula A16k, and Nk is the number of vessels within order k. Expressing Eq. A16 in terms of the presence of self-similar fractal scaling relationships, a constant diameter ratio Formula A16d = Formula A16k+1/Formula A16k and a constant length ratio Formula A16l = Formula A16k+1/Formula A16k, Rcum becomes

Formula A17(A17)
In asymmetric trees, the relative values of Formula A17d and Formula A17b depend on the methods of diameter ordering: Horsfield ordering yields Formula A17bH≤ 2, and a smaller Formula A17dH(26, 27) with more orders for the same hemodynamically equivalent tree than that ordered by a Strahler method with Formula A17bS≤ 2 (30). The branching ratios are related by 1 = (RdH)–1 + (RdS)–1 = (RdH)–1 + (RdH)–(1+{Delta}) via a symmetry factor {Delta} (26, 27). If an additional scaling relationship for vessel length and diameter is taken into account, assuming that arterial network organization has a common fractal length-diameter scaling relationship Formula A17k = aFormula A17k{alpha}2({alpha}2 = 1 and a is an arbitrary scaling constant that depends on the method of vessel ordering) and fractal scaling condition Formula A17b = Formula A17dD0, then Rcum can be expressed more simply as

Formula A18(A18)
In effect, the resistance of order k is R(k) = R0Formula A18d(3–{alpha})k, in which the fractal network organization imparts a contributing reactive factor {kappa}R = Formula A18d3–D0 within each order. Thus, under the assumptions of network organization in this fractal model, the resistance reactivity is equal to the shear stress reactivity {kappa} = {kappa}{tau} of Eq. A15.

Formula A19(A19)
Figure 1 illustrates the impact of hemodynamic reactivity with branching on normalized values of Rcum and {tau} for the simplest case of a symmetrically bifurcating network where Formula A19b = 2 and Formula A19d = 21/D0.


    FOOTNOTES
 

Address for reprint requests and other correspondence: S. H. Bennett, Neonatology Surge I Rm. 1121, Univ. of California, Davis, CA 95616 (e-mail: shbennett{at}ucdavis.edu)

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
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 APPENDIX: ROLE OF ARTERIAL...
 REFERENCES
 

  1. Barnard J, Wilson P, Moore T, Thompson W, Taylor A. Effect of nitric oxide and cyclooxygenase products on vascular resistance in dog and rat lungs. J Appl Physiol 74: 2940–2948, 1993.[Abstract/Free Full Text]
  2. Belik J, Keeley FW, Baldwin F, Rabinovitch M. Pulmonary hypertension and vascular remodeling in fetal sheep. Am J Physiol Heart Circ Physiol 266: H2303–H2309, 1994.[Abstract/Free Full Text]
  3. Bennett S, Eldridge M, Zaghi D, Zaghi S, Milstein J, Goetzman B. Form and function of fetal and neonatal pulmonary arterial bifurcations. Am J Physiol Heart Circ Physiol 279: H3047–H3057, 2000.[Abstract/Free Full Text]
  4. Bennett SH, Goetzman BW, Milstein JM, Pannu JS. Role of arterial design on pulse wave reflection in a fractal pulmonary network. J Appl Physiol 80: 1033–1056, 1996.[Abstract/Free Full Text]
  5. Black SM, Bekker JM, Johengen MJ, Parry AJ, Soifer SJ, Fineman JR. Altered regulation of the ET-1 cascade in lambs with increased pulmonary blood flow and pulmonary hypertension. Pediatr Res 47: 97–106, 2000.[Web of Science][Medline]
  6. Bousamra M, Rossi R, Jacobs E, Parviz M, Busch C, Nelin L, Haworth S, Dawson C. Systemic lobar shunting induces advanced pulmonary vasculopathy. J Thorac Cardiovasc Surg 120: 88–98, 2000.[Abstract/Free Full Text]
  7. Budhiraja R, Tuder R, Hassoun PM. Endothelial dysfunction in pulmonary hypertension. Circulation 109: 159–165, 2004.[Free Full Text]
  8. Cassin S, Dawes G, Mott J, Ross B, Strang L. The vascular resistance of the foetal and newly ventilated lung of the lamb. J Physiol 171: 61–79, 1963.[Web of Science]
  9. Cornfield D, Resnik E, Herron J, Reinhartz O, Fineman J. Pulmonary vascular K+ channel expression and vasoreactivity in a model of congenital heart disease. Am J Physiol Lung Cell Mol Physiol 283: L1210–L1219, 2002.[Abstract/Free Full Text]
  10. Dakshinamurti S. Pathophysiologic mechanisms of persistent pulmonary hypertension of the newborn. Pediatr Pulmonol 39: 492–503, 2005.[CrossRef][Web of Science][Medline]
  11. Davies PF. Flow-mediated endothelial mechanotransduction. Physiol Rev 75: 519–560, 1995.[Abstract/Free Full Text]
  12. Dawson C. Role of pulmonary vasomotion in physiology of the lung. Physiol Rev 64: 544–616, 1984.[Free Full Text]
  13. Dawson CA, Krenz GS, Karau KL, Haworth ST, Hanger CC, Linehan JH. Structure-function relationships in the pulmonary artery tree. J Appl Physiol 86: 569–583, 1999.[Abstract/Free Full Text]
  14. Dekker R, Thienen JV, Rohlena J, Jager SD, Elderkamp Y, Seppen J, deVries C, Blessen E, Berkel TV, Pannekeoek H, Horrevoets A. Endothelial KLF2 links local arterial shear stress levels to the expression of vascular tone-regulating genes. Am J Pathol 167: 609–618, 2005.[Abstract/Free Full Text]
  15. Emery C, Bee D, Barer G. Mechanical properties and reactivity of vessels in isolated perfused lungs of chronically hypoxic rats. Clin Sci (Lond) 61: 569–580, 1981.[Medline]
  16. Engelberg J, Dubois A. Mechanics of pulmonary circulation in isolated rabbit lungs. Am J Physiol 196: 401–414, 1959.[Abstract/Free Full Text]
  17. Fineman JR, Soifer SJ, Heyman MA. Regulation of pulmonary vascular tone in the perinatal period. Annu Rev Physiol 57: 115–134, 1995.[CrossRef][Web of Science][Medline]
  18. Fineman JR, Soifer SJ, Heymann MA. The role of pulmonary vascular endothelium in perinatal pulmonary circulatory regulation. Semin Perinatol 15: 58–62, 1991.[Web of Science][Medline]
  19. Fry D. Acute vascular endothelial changes associated with increased blood velocity gradients. Circ Res 22: 165–197, 1968.[Abstract/Free Full Text]
  20. Galie N, Torbicki A, Barst R, Dartevelle P, Haworth S, Higenbottam T, Olschewski H, Peacock A, Pietra G, Rubin LJ, Simonneau G, Priori SG, Garcia MA, Blanc JJ, Budaj A, Cowie M, Dean V, Deckers J, Burgos EF, Lekakis J, Lindahl B, Mazzotta G, McGregor K, Morais J, Oto A, Smiseth OA, Barbera JA, Gibbs S, Hoeper M, Humbert M, Naeije R, Pepke-Zaba J, Force T. Guidelines on diagnosis and treatment of pulmonary arterial hypertension. The Task Force on Diagnosis and Treatment of Pulmonary Arterial Hypertension of the European Society of Cardiology. Eur Heart J 25: 2243–2278, 2004.[Free Full Text]
  21. Garcia-Cardena G, Comander J, Anderson KR, Blackman BR, Gimbrone MA Jr. Biomechanical activation of vascular endothelium as a determinant of its functional phenotype. Proc Natl Acad Sci USA 98: 4478–4485, 2001.[Abstract/Free Full Text]
  22. Gregory T, Ellsworth M, Newell J. Hypoxic pulmonary vasoconstriction during steady and pulsatile flow in ferrets. J Appl Physiol 57: 205–212, 1984.[Abstract/Free Full Text]
  23. Griffith T, Edwards D, Randall M. Blood flow and optimal vascular topography: role of the endothelium. Basic Res Cardiol 86: 89–96, 1991.[Web of Science][Medline]
  24. Haworth S, Linehan J, Bronikowski T, Dawson C. A hemodynamic model representation of the dog lung. J Appl Physiol 79: 15–26, 1991.
  25. Heath D and Edwards J. The pathology of hypertensive pulmonary vascular disease. Circulation 18: 533–547, 1958.[Web of Science][Medline]
  26. Hoffman JIE, Rudolph AM, Heyman MA. Pulmonary vascular disease with congenital heart lesions: pathologic features and causes. Circulation 64: 873–877, 1981.[Abstract/Free Full Text]
  27. Horsfield K. Some mathematical properties of branching trees with application to the respiratory system. Bull Math Biol 38: 305–315, 1976.[Web of Science][Medline]
  28. Horsfield K, Cumming G. Morphometry of the bronchial tree in man. J Appl Physiol 24: 373–383, 1968.[Free Full Text]
  29. Horsfield K, Woldenberg MJ. Diameters and cross-sectional areas of branches in the human pulmonary arterial tree. Anat Rec 223: 245–251, 1989.[CrossRef][Medline]
  30. Kleiber M. The Fire of Life: an Introduction to Animal Energetics. New York: Wiley, 1961.
  31. Krenz GS, Linehan JH, Dawson CA. A fractal continuum model of the pulmonary arterial tree. J Appl Physiol 72: 2225–2237, 1992.[Abstract/Free Full Text]
  32. Krishnan A, Linehan JH, Rickaby DA, Dawson CA. Cat lung hemodynamics: comparison of experimental results and model predictions. J Appl Physiol 61: 2023–2034, 1986.[Abstract/Free Full Text]
  33. Kurz H, Sandau KS. Modeling of blood vessel development-bifurcation pattern and hemodynamics, optimality and allometry. Comments Theor Biol 4: 261–291, 1997.
  34. LaBarbera M. Principles of design of fluid transport systems in biology. Science 249: 992–1000, 1990.[Abstract/Free Full Text]
  35. Langille B, O'Donnell F. Reductions in arterial diameter produced by chronic decreases in blood flow are endothelium-dependent. Science 231: 405–407, 1986.[Abstract/Free Full Text]
  36. Langille L. Integrity of arterial endothelium following acute exposure to high shear stress. Biorheology 21: 333–346, 1984.[Web of Science][Medline]
  37. Laskey W, Ferrari VA, Palevsky HI, Kussmaul WG. Pulmonary artery hemodynamics in primary pulmonary hypertension. J Am Coll Cardiol 21: 406–412, 1993.[Abstract]
  38. Limpert E, Stahel W, Abbt M. Log-normal distributions across the sciences: keys and clues. Bioscience 51: 341–352, 2001.[CrossRef][Web of Science]
  39. Maggiorini M, Brimioulle S, De Canniere D, Delcroix M, Wauthy P, Naeije R. Pulmonary vascular impedance response to hypoxia in dogs and minipigs: effects of inhaled nitric oxide. J Appl Physiol 79: 1156–1162, 1995.[Abstract/Free Full Text]
  40. Maggiorini M, Melot C, Gilbert E, Vermeulen F, Naeije R. Pulmonary vascular resistance in dogs and minipigs—effects of hypoxia and inhaled nitric oxide. Respir Physiol 111: 213–222, 1998.[CrossRef][Web of Science][Medline]
  41. Mata-Greenwood E, Meyrick B, Steinhorn RH, Fineman JR, Black SM. Alterations in TGF-beta1 expression in lambs with increased pulmonary blood flow and pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol 285: L209–L221, 2003.[Abstract/Free Full Text]
  42. Murphy JD, Rabinovitch M, Goldstein JD, Reid LM. The structural basis of persistent pulmonary hypertension of the newborn infant. J Pediatr 98: 962–967, 1981.[CrossRef][Web of Science][Medline]
  43. Peake M, Harabin A, Brennan N, Sylvester J. Steady-state vascular responses to graded hypoxia in isolated lungs of five species. J Appl Physiol 51: 1214–1219, 1981.[Abstract/Free Full Text]
  44. Rabinovitch M, Haworth S, Castenada A, Nadas A, and Reid L. Lung biopsy in congenital heart disease: a morphometric approach to vascular disease. Circulation 58: 1107–1122, 1978.[Abstract/Free Full Text]
  45. Rabinovitch M, Keane J, Norwood W, Castenada A, Reid L. Vascular structure in lung biopsy tissue correlated with pulmonary hemodynamic findings after repair of congenital heart defects. Circulation 69: 655–667, 1984.[Abstract/Free Full Text]
  46. Reddy VM, Meyrick B, Wong J, Khoor A, Liddicoat JR, Hanley FL, Fineman JR. In utero placement of aortopulmonary shunts—a model of postnatal pulmonary hypertension with increased pulmonary blood flow in lambs. Circulation 92: 606–613, 1995.[Abstract/Free Full Text]
  47. Reeves J, Taylor A. Pulmonary hemodynamics and fluid exchange in the lungs during exercise. In: Handbook of Physiology. Exercise Regulation and Integration of Multiple Systems. Bethesda, MD: Am. Physiol. Soc., 1996, sect. 12, chapt. 13, p. 585–613.
  48. Reid DL, Thornberg KL. Pulmonary pressure-flow relationships in the fetal lamb during in utero ventilation. J Appl Physiol 69: 1630–1636, 1990.[Abstract/Free Full Text]
  49. Reid LM. Vascular remodeling. In: The Pulmonary Circulation, edited by Fishman A. Philadelphia: Univ. of Pennsylvania Press, 1990, p. 259–282.
  50. Rendas A, Lennox S, Reid L. Aortopulmonary shunts in growing pigs. J Thorac Cardiovasc Surg 77: 109–118, 1979.[Abstract]
  51. Roy AG, Woldenberg MJ. A generalization of the optimal models of arterial branching. Bull Math Biol 44: 349–360, 1982.[Web of Science][Medline]
  52. Rudolph AM. The changes in the circulation after birth: their importance in congenital heart disease. Circulation 41: 342–359, 1970.
  53. Sherman TF. On connecting large vessels to small: the meaning of Murray's law. J Gen Physiol 78: 431–453, 1981.[Abstract/Free Full Text]
  54. Singhal S, Henderson R, Horsfield K, Harding K, Cumming G. Morphometry of the human pulmonary arterial tree. Circ Res 33: 190–197, 1973.[Abstract/Free Full Text]
  55. Stenmark K, Mecham R. Cellular and molecular mechanisms of pulmonary vascular remodeling. Annu Rev Physiol 59: 89–144, 1997.[CrossRef][Web of Science][Medline]
  56. Vasinav S, Patel D, Atabek H, Deshpande M, Plowman F, Vossoughi J. Determination of the local erosion stress of the canine endothelium using a jet impingement method. J Biomech Eng 105: 77–83, 1983.[Web of Science][Medline]
  57. Weibel E. Design of biological organisms and fractal geometry. In: Fractals in Biology and Medicine, edited by Nonnenmacher T, Losa G, and Weibel E. Basel: Birkhauser, 1994, p. 68–85.
  58. West G, Brown J, Enquist B. A general model for ontogenic growth. Nature 413: 628–631, 2001.[CrossRef]
  59. West G, Brown J, Enquist B. A general model for the origin of allometric scaling laws in biology. Science 276: 122–126, 1997.[Abstract/Free Full Text]
  60. Zar H. Biostatistical Analysis. Englewood Cliffs, NJ: Prentice-Hall, 1984.




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
292/6/H3006    most recent
01012.2006v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ghorishi, Z.
Right arrow Articles by Eldridge, M. W.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ghorishi, Z.
Right arrow Articles by Eldridge, M. W.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2007 by the American Physiological Society.