Wave intensity analysis of right ventricular and pulmonary vascular contributions to higher pulmonary than aortic blood pressure in fetal lambs

Mechanistic insights on age-related changes in heart-aorta-brain hemodynamic coupling using a pulse wave model of the entire circulatory system

Age-related changes in aortic biomechanics can impact the brain by reducing blood flow and increasing pulsatile energy transmission. Clinical studies have shown that impaired cardiac function in patients with heart failure is associated with cognitive impairment. Although previous studies have attempted to elucidate the complex relationship between age-associated aortic stiffening and pulsatility transmission to the cerebral network, they have not adequately addressed the effect of interactions between aortic stiffness and left ventricle (LV) contractility (neither on energy transmission nor on brain perfusion). In this study, we use a well-established and validated one-dimensional blood flow and pulse wave computational model of the circulatory system to address how age-related changes in cardiac function and vasculature affect the underlying mechanisms involved in the LV-aorta-brain hemodynamic coupling. Our results reveal how LV contractility affects pulsatile energy transmission to the brain, even with preserved cardiac output. Our model demonstrates the existence of an optimal heart rate (near the normal human heart rate) that minimizes pulsatile energy transmission to the brain at different contractility levels. Our findings further suggest that the reduction in cerebral blood flow at low levels of LV contractility is more prominent in the setting of age-related aortic stiffening. Maintaining optimal blood flow to the brain requires either an increase in contractility or an increase in heart rate. The former consistently leads to higher pulsatile power transmission, and the latter can either increase or decrease subsequent pulsatile power transmission to the brain.NEW & NOTEWORTHY We investigated the impact of major aging mechanisms of the arterial system and cardiac function on brain hemodynamics. Our findings suggest that aging has a significant impact on heart-aorta-brain coupling through changes in both arterial stiffening and left ventricle (LV) contractility. Understanding the underlying physical mechanisms involved here can potentially be a key step for developing more effective therapeutic strategies that can mitigate the contributions of abnormal LV-arterial coupling toward neurodegenerative diseases and dementia.

Fluid-structure interaction in a fully coupled three-dimensional mitral-atrium-pulmonary model

This paper aims to investigate detailed mechanical interactions between the pulmonary haemodynamics and left heart function in pathophysiological situations (e.g. atrial fibrillation and acute mitral regurgitation). This is achieved by developing a complex computational framework for a coupled pulmonary circulation, left atrium and mitral valve model. The left atrium and mitral valve are modelled with physiologically realistic three-dimensional geometries, fibre-reinforced hyperelastic materials and fluid–structure interaction, and the pulmonary vessels are modelled as one-dimensional network ended with structured trees, with specified vessel geometries and wall material properties. This new coupled model reveals some interesting results which could be of diagnostic values. For example, the wave propagation through the pulmonary vasculature can lead to different arrival times for the second systolic flow wave (S2 wave) among the pulmonary veins, forming vortex rings inside the left atrium. In the case of acute mitral regurgitation, the left atrium experiences an increased energy dissipation and pressure elevation. The pulmonary veins can experience increased wave intensities, reversal flow during systole and increased early-diastolic flow wave (D wave), which in turn causes an additional flow wave across the mitral valve (L wave), as well as a reversal flow at the left atrial appendage orifice. In the case of atrial fibrillation, we show that the loss of active contraction is associated with a slower flow inside the left atrial appendage and disappearances of the late-diastole atrial reversal wave (AR wave) and the first systolic wave (S1 wave) in pulmonary veins. The haemodynamic changes along the pulmonary vessel trees on different scales from microscopic vessels to the main pulmonary artery can all be captured in this model. The work promises a potential in quantifying disease progression and medical treatments of various pulmonary diseases such as the pulmonary hypertension due to a left heart dysfunction.

Minor impact of constraint from perivascular flow probes on wave intensity analysis

Perivascular flow probes are considered the gold-standard for measuring volumetric blood flow in animal studies. Although flow probes are generally placed non-constrictively around the vessel of interest, pressure-elevating interventions performed during an experiment may lead to vessel expansion and some probe-vessel impingement, particularly in highly compliant vessels such as adult sheep aorta or major pulmonary arteries in fetus lambs. This study assessed to what extent such mild flow probe constraint may impact on wave intensity analysis. We also investigated whether errors arising from flow probe constraint could explain apparent pressure reflection indices (R_p > 1) that have been observed in fetus lamb pulmonary arteries under some experimental conditions. These questions were investigated with one-dimensional models of an adult sheep aorta and fetus lamb pulmonary artery, with a virtual flow probe incorporated as a non-linear external constraint term in the vessel constitutive equation. Model-derived flow and pressure were subjected to standard analysis procedures that would be applied experimentally (correcting for apparent velocity lags and calculating wave speed via the PU-loop method). For the adult sheep model, simulations covering a wide range of haemodynamic conditions revealed a mostly minor effect (<10%) of probe constraint on the intensity and pressure effects of the three major waves (forward compression wave, forward decompression wave, backward compression wave). Moreover, flow probe constraint had essentially no impact on R_p in the fetus lamb model, suggesting that such constraint is unlikely to be responsible for an observed R_p > 1. Mild flow probe constraint is likely to have little impact on wave intensity analysis.

Reducing lung liquid volume increases biventricular outputs and systemic arterial blood flows despite decreased cardiac filling pressures in fetal lambs

As prior work has shown that reducing lung liquid volume 1) increases pulmonary arterial (PA) blood flow, 2) augments right ventricular (RV) output/power, and 3) decreases left atrial (LA) pressure, we tested the hypothesis that this perturbation has global cardiovascular effects. Ten anesthetized, open-chest fetal lambs (128 ± 2 days gestation, full term = 147 days) were acutely instrumented with 1) LA and right atrial (RA) catheters, 2) aortic and pulmonary trunk catheters, 3) brachiocephalic trunk, aortic isthmus, ductal, and left PA flow probes to obtain left ventricular (LV) and RV outputs and hydraulic power and flow in the descending thoracic aorta, and 4) an endotracheal tube to remove lung liquid. A 17 ± 7 ml/kg reduction of lung liquid volume 1) decreased LA and RA pressures similarly (1.5-1.6 mmHg, P < 0.001), 2) augmented LV and RV outputs (21-24%, P < 0.001) and total power (27-28%, P < 0.005), 3) increased systolic flows in the brachiocephalic trunk (18%, P < 0.001), aortic isthmus (29%, P < 0.005), ductus (12%, P < 0.005), and descending thoracic aorta (16%, P < 0.001), 4) increased mean PA flow via a higher systolic inflow (37%, P < 0.001) and lower diastolic backflow (-16%, P < 0.05), and 5) did not change systemic vascular conductance or arterial compliance but increased both pulmonary vascular conductance and arterial compliance (1.8-fold, P < 0.001). These data suggest that hemodynamic effects of lung liquid volume reduction are not confined to the lungs but extend to all cardiac chambers via rises in LV and RV outputs and power, despite falls in cardiac filling pressures, as well as the systemic circulation, via downstream increases in systolic flows of major central arteries.

Increased right ventricular power and ductal characteristic impedance underpin higher pulmonary arterial blood flow after betamethasone therapy in fetal lambs

BACKGROUND

The glucocorticosteroid betamethasone is routinely administered prior to anticipated preterm birth to enhance lung maturation. While betamethasone also increases fetal pulmonary blood flow and reduces pulmonary vascular resistance (PVR), we investigated whether alterations in right ventricular (RV) function and ductal characteristic impedance (Z_c) additionally contributed to rises in pulmonary flow.

METHODS

Anesthetized preterm fetal lambs with (n = 10) or without (n = 8) betamethasone pretreatment were instrumented with a pulmonary trunk micromanometer and ductus arteriosus and left pulmonary artery (PA) flow probes to calculate Z_c, and obtain RV output and hydraulic power.

RESULTS

Betamethasone (1) increased systolic and pulse arterial pressures (P ≤ 0.04), heart rate (P = 0.02), and lowered PVR (P = 0.04), (2) increased mean (P = 0.008) and systolic (P = 0.004), but not diastolic PA flow or PA Z_c, (3) increased ductal Z_c (P < 0.05), but not ductal flow, (4) increased RV output (P = 0.03) and the proportion of PT flow distributed to the lungs (P = 0.02), and (5) increased RV power (P ≤ 0.002).

CONCLUSION

An increased fetal PA blood flow after betamethasone therapy was confined to the systole and underpinned not only by decreased PVR, but also greater RV power and preferential distribution of an augmented RV systolic outflow to the lungs due to higher ductal Z_c.

Novel wave power analysis linking pressure-flow waves, wave potential, and the forward and backward components of hydraulic power

Wave intensity analysis provides detailed insights into factors influencing hemodynamics. However, wave intensity is not a conserved quantity, so it is sensitive to diameter variations and is not distributed among branches of a junction. Moreover, the fundamental relation between waves and hydraulic power is unclear. We, therefore, propose an alternative to wave intensity called “wave power,” calculated via incremental changes in pressure and flow (dPdQ) and a novel time-domain separation of hydraulic pressure power and kinetic power into forward and backward wave-related components (ΠP± and ΠQ±). Wave power has several useful properties: 1) it is obtained directly from flow measurements, without requiring further calculation of velocity; 2) it is a quasi-conserved quantity that may be used to study the relative distribution of waves at junctions; and 3) it has the units of power (Watts). We also uncover a simple relationship between wave power and changes in ΠP± and show that wave reflection reduces transmitted power. Absolute values of ΠP± represent wave potential, a recently introduced concept that unifies steady and pulsatile aspects of hemodynamics. We show that wave potential represents the hydraulic energy potential stored in a compliant pressurized vessel, with spatial gradients producing waves that transfer this energy. These techniques and principles are verified numerically and also experimentally with pressure/flow measurements in all branches of a central bifurcation in sheep, under a wide range of hemodynamic conditions. The proposed “wave power analysis,” encompassing wave power, wave potential, and wave separation of hydraulic power provides a potent time-domain approach for analyzing hemodynamics.

Retrograde lower body arterial reservoir discharge underlies rapid reversal of ductus arteriosus shunting after early cord clamping at birth in preterm lambs

Arterial reservoir (“windkessel”) function, whereby a part of left ventricular (LV) output is stored in elastic arteries during systole and discharged in diastole, is a well-established physiological phenomenon. However, its role in rapid reversal (to left-to-right) and a systolic-to-diastolic shift of shunting across the ductus arteriosus after birth is unknown. To address this question, ductal and aortic isthmus flows were measured with high-fidelity transit-time probes in six anesthetized preterm fetal lambs before and after cord clamping and subsequent early mechanical ventilation and for 30 min postbirth. Descending aortic flow was calculated as the sum of isthmus and ductal flows. Left-to-right ductal flow profiles were related to those of the isthmus and descending aorta, with upper body arterial reservoir discharge indicated by forward diastolic isthmus flow, and retrograde lower body arterial reservoir discharge by negative diastolic descending aortic flow. Left-to-right ductal shunting appeared immediately after cord clamping ( P < 0.001), due entirely to newly emergent retrograde lower body reservoir discharge, and rose with ventilation via increased lower body reservoir discharge ( P < 0.005), supplemented by upper body reservoir discharge after 45 s ( P < 0.025) and LV systolic flow after 3 min ( P = 0.025). The contribution of lower body reservoir discharge to left-to-right ductal shunting fell to 55 ± 8% at ≥15 min ( P < 0.001) but remained higher ( P < 0.002) than LV systolic flow (33 ± 8%) or upper body reservoir discharge (12 ± 5%). These results suggest that retrograde lower body arterial reservoir discharge plays a key role in rapid reversal and a systolic-to-diastolic shift of ductal shunting after cord clamping and early ventilation at birth.

A computational study of pressure wave reflections in the pulmonary arteries

Experiments using wave intensity analysis suggest that the pulmonary circulation in sheep and dogs is characterized by negative or open-end type wave reflections, that reduce the systolic pressure. Since the pulmonary physiology is similar in most mammals, including humans, we test and verify this hypothesis by using a subject specific one-dimensional model of the human pulmonary circulation and a conventional wave intensity analysis. Using the simulated pressure and velocity, we also analyse the performance of the P-U loop and sum of squares techniques for estimating the local pulse wave velocity in the pulmonary arteries, and then analyse the effects of these methods on linear wave separation in the main pulmonary artery. P-U loops are found to provide much better estimates than the sum of squares technique at proximal locations, but both techniques accumulate progressive error at distal locations away from heart, particularly near junctions. The pulse wave velocity estimated using the sum of squares method also gives rise to an artificial early systolic backward compression wave. Finally, we study the influence of three types of pulmonary hypertension viz. pulmonary arterial hypertension, chronic thromboembolic pulmonary hypertension and pulmonary hypertension associated with hypoxic lung disease. Simulating these conditions by changing the relevant parameters in the model and then applying the wave intensity analysis, we observe that for each group the early systolic backward decompression wave reflected from proximal junctions is maintained, whilst the initial forward compression and the late systolic backward compression waves amplify with increasing pathology and contribute significantly to increases in systolic pressure.

Noninvasive pulmonary artery wave intensity analysis in pulmonary hypertension

Pulmonary wave reflections are a potential hemodynamic biomarker for pulmonary hypertension (PH) and can be analyzed using wave intensity analysis (WIA). In this study we used pulmonary vessel area and flow obtained using cardiac magnetic resonance (CMR) to implement WIA noninvasively. We hypothesized that this method could detect differences in reflections in PH patients compared with healthy controls and could also differentiate certain PH subtypes. Twenty patients with PH (35% CTEPH and 75% female) and 10 healthy controls (60% female) were recruited. Right and left pulmonary artery (LPA and RPA) flow and area curves were acquired using self-gated golden-angle, spiral, phase-contrast CMR with a 10.5-ms temporal resolution. These data were used to perform WIA on patients and controls. The presence of a proximal clot in CTEPH patients was determined from contemporaneous computed tomography/angiographic data. A backwards-traveling compression wave (BCW) was present in both LPA and RPA of all PH patients but was absent in all controls (P = 6e⁻⁸). The area under the BCW was associated with a sensitivity of 100% [95% confidence interval (CI) 63–100%] and specificity of 91% (95% CI 75–98%) for the presence of a clot in the proximal PAs of patients with CTEPH. In conclusion, WIA metrics were significantly different between patients and controls; in particular, the presence of an early BCW was specifically associated with PH. The magnitude of the area under the BCW showed discriminatory capacity for the presence of proximal PA clot in patients with CTEPH. We believe that these results demonstrate that WIA could be used in the noninvasive assessment of PH.

Enhanced ventricular pump function and decreased reservoir backflow sustain rise in pulmonary blood flow after reduction of lung liquid volume in fetal lambs

Although a reduction in lung liquid volume increases fetal pulmonary blood flow, the changes in central flow patterns that sustain this increased pulmonary perfusion are unknown. To address this issue, eight anesthetized late-gestation fetal sheep were instrumented with pulmonary trunk (PT), ductus arteriosus (DA), and left pulmonary artery (PA) micromanometer catheters and transit-time flow probes, with blood flow profile and wave intensity analyses performed at baseline and after withdrawal of lung liquid via an endotracheal tube. Reducing lung liquid volume by 19 ± 6 ml/kg (mean ± SD) augmented right ventricular power by 34% (P < 0.001), with distribution of an accompanying increase in mean PT blood flow (245 ± 63 ml/min, P < 0.001) to the lungs (169 ± 91 ml/min, P = 0.001) and across the DA (77 ± 92 ml/min, P = 0.04). However, although PT and DA flow increments were confined to systole and were related to an increased magnitude of flow-increasing, forward-running compression waves, the rise in PA flow spanned both systole (108 ± 66 ml/min) and diastole (61 ± 32 ml/min). Flow profile analysis showed that the step-up in PA diastolic flow was associated with diminished PA diastolic backflow and accompanied by a lesser degree of diastolic right-to-left DA shunting. These data suggest that an increased pulmonary blood flow after reduction of lung liquid volume is associated with substantial changes in PT-DA-PA interactions and underpinned by two main factors: 1) enhanced right ventricular pump function that increases PA systolic inflow and 2) decreased PA diastolic backflow that arises from a fundamental change in PA reservoir function, thereby resulting in greater passage of systolic inflow through the lungs.

Noninvasive determination of local pulse wave velocity and wave intensity: changes with age and gender in the carotid and femoral arteries of healthy human

We recently introduced noninvasive methods to assess local pulse wave velocity (PWV) and wave intensity (ndI) in arteries based on measurements of flow velocity (U) and diameter (D). Although the methods were validated in an experimental setting, clinical application remains lacking. The aim of this study was therefore to investigate the effect of age and gender on PWV and ndI in the carotid and femoral arteries of an existing population. We measured D and U in the carotid and femoral arteries of 1,774 healthy subjects aged 35-55 yr, a subgroup of the Asklepios population. With the use of the lnDU-loop method, we calculated local PWV, which was used to determine arterial distensibility (nDs). We then used the new algorithm to determine maximum forward and backward wave intensities (ndI+max and ndI−min, respectively) and the reflection index (nRI). On average, PWV was higher, and nDs was lower in the femoral than at the carotid arteries. At the carotid artery, PWV increased with age, but nDs, ndI+max, and ndI−min decreased; nRI did not change with age. At the femoral artery, PWV was higher, and nDs was lower in male, but all parameters did not change significantly with age in both women and men. We conclude that the carotid artery is more affected by the aging process than the femoral artery, even in healthy subjects. The new techniques provide mechanical and hemodynamic parameters, requiring only D and U measurements, both of which can be acquired using ultrasound equipment widely available today, hence their advantage for potential use in the clinical setting.

The reservoir-wave paradigm introduces error into arterial wave analysis: a computer modelling and in-vivo study

OBJECTIVES

Arterial wave reflection has traditionally been quantified from pressure and flow measurements using wave separation and wave intensity (WI) analysis. In the recently proposed reservoir-wave paradigm, these analyses are performed after dividing pressure into 'reservoir' and 'excess' components, yielding a modified wave intensity (WI(RW)). This new approach has led to controversial conclusions about the nature and significance of arterial wave reflection. Our aim was to assess whether WI or WI(RW) more accurately represent wave phenomena.

METHODS

We studied two computer models (a simple network and a full model of the systemic arterial tree) in which all systolic forward waves and reflection properties were known a priori. Results of these models were compared with haemodynamic measurements in the ascending aorta of five adult sheep at baseline and after incremental arterial constriction.

RESULTS

The key findings of model studies were that the reservoir-wave approach markedly underestimated or eliminated reflected compression waves, overestimated or artefactually introduced forward and backward expansion waves, and displayed nonphysical interactions between distal reflection sites and early systolic waves. These errors arose because, contrary to a key assumption of the reservoir-wave approach, reservoir pressure was not spatially uniform during systole. In-vivo results were qualitatively similar to model results, with baseline WI and WI(RW) suggesting that the arterial network was dominated by positive and negative wave reflection, respectively, while under all conditions, reflected WI(RW) compression waves were substantially smaller than corresponding WI waves.

CONCLUSION

We conclude that the reservoir-wave paradigm introduces error into arterial wave analyses.

Enhanced central and conduit pulmonary arterial reservoir function offsets reduced ductal systolic outflow during constriction of the fetal ductus arteriosus

Constriction of the fetal ductus arteriosus (DA) has disparate effects on mean and phasic hemodynamics, as mean DA blood flow is preserved until constriction is severe, but DA systolic and diastolic blood velocities change with only mild constriction. To determine the basis of this disparity and its physiological significance, seven anesthetized late-gestation fetal sheep were instrumented with pulmonary trunk (PT), DA, and left pulmonary artery (PA) micromanometer catheters and transit-time flow probes. Blood flow profile and wave intensity analyses were performed at baseline and during mild, moderate, and severe DA constriction (defined as pulmonary-aortic mean pressure differences of 4, 8, and 14 mmHg, respectively), produced with an adjustable snare. With DA constriction, mean DA flow was initially maintained but decreased with severe constriction (P < 0.05) in conjunction with a reduction (P < 0.05) in PT flow (i.e., right ventricular output). By contrast, DA systolic flow fell progressively during DA constriction (P < 0.001), due to decreased transmission of both early and midsystolic proximal flow-enhancing forward-running compression waves into the DA. However, DA constriction was also accompanied by greater systolic storage of blood in the PT and main PA (P < 0.025), and increased retrograde diastolic flow from compliant major branch PA (P < 0.001). Transductal discharge of these central and conduit PA blood reservoirs in diastole offset systolic DA flow reductions. These data suggest that, during DA constriction in the fetus, enhanced central and conduit PA reservoir function constitutes an important compensatory mechanism that contributes to preservation of mean DA flow via a systolic-to-diastolic redistribution of phasic DA flow.

Pulmonary trunk, ductus arteriosus, and pulmonary arterial phasic blood flow interactions during systole and diastole in the fetus

Although the distribution of average fetal pulmonary trunk (PT) blood flow favors the ductus arteriosus (DA) over the lungs, the phasic aspects of this distribution during systole and diastole are not well understood. Accordingly, flow profile and wave intensity (WI) analyses were performed at baseline and during brief flow increases accompanying an extrasystole (ES) in 10 anesthetized late-gestation fetal sheep instrumented with PT, DA, and left pulmonary artery (PA) micromanometer catheters and transit-time flow probes. At baseline, 83% of mean PT flow crossed the DA and 17% entered the lungs. However, early systolic flow associated with a forward-running compression wave (FCW(is)) was higher in the PA and predominant DA flow only emerged in midsystole when a large PA backward-running compression wave (BCW(ms)), which reduced PA flow, was transmitted into the DA as a forward-running compression wave (FCW(ms)) that increased flow. Subsequent protodiastolic forward DA flow occurring during pulmonary valve closure was associated with substantial retrograde PA flow, but insignificant PT flow. Conversely, forward DA flow in the remainder of diastole occurred with forward PT but near-zero PA flow. These flow and WI patterns, in conjunction with the results of mathematical modeling, suggest that 1) fetal PT flow preferentially passes into the PA during early systole due to a lower PA-than-DA characteristic impedance, while DA flow predominates in mid- and late systole due to flow effects arising from the PA BCW(ms), and 2) forward DA flow is mainly sustained by reversal of PA flow in protodiastole but discharge of a more central reservoir in diastole.