Efficient sensitivity analysis for biomechanical models with correlated inputs.

In most variance-based sensitivity analysis (SA) approaches applied to biomechanical models, statistical independence of the model input is assumed. However, often the model inputs are correlated. This might alter the interpretation of the SA results, which may severely impact the guidance provided during model development and personalization. Potential reasons for the infrequent usage of SA techniques that account for input correlation are the associated high computational costs, especially for models with many parameters, and the fact that the input correlation structure is often unknown. The aim of this study was to propose an efficient correlated global sensitivity analysis method by applying a surrogate model-based approach. Furthermore, this article demonstrates how correlated SA should be interpreted and how the applied method can guide the modeler during model development and personalization, even when the correlation structure is not entirely known beforehand. The proposed methodology was applied to a typical example of a pulse wave propagation model and resulted in accurate SA results that could be obtained at a theoretically 27,000× lower computational cost compared to the correlated SA approach without employing a surrogate model. Furthermore, our results demonstrate that input correlations can significantly affect SA results, which emphasizes the need to thoroughly investigate the effect of input correlations during model development. We conclude that our proposed surrogate-based SA approach allows modelers to efficiently perform correlated SA to complex biomechanical models and allows modelers to focus on input prioritization, input fixing and model reduction, or assessing the dependency structure between parameters.

Estimation of central blood pressure waveform from femoral blood pressure waveform by blind sources separation.

Background: Central blood pressure (cBP) is a better indicator of cardiovascular morbidity and mortality than peripheral BP (pBP). However, direct cBP measurement requires invasive techniques and indirect cBP measurement is based on rigid and empirical transfer functions applied to pBP. Thus, development of a personalized and well-validated method for non-invasive derivation of cBP from pBP is necessary to facilitate the clinical routine. The purpose of the present study was to develop a novel blind source separation tool to separate a single recording of pBP into their pressure waveforms composing its dynamics, to identify the compounds that lead to pressure waveform distortion at the periphery, and to estimate the cBP. The approach is patient-specific and extracts the underlying blind pressure waveforms in pBP without additional brachial cuff calibration or any a priori assumption on the arterial model. Methods: The intra-arterial femoral BPfe and intra-aortic pressure BPao were anonymized digital recordings from previous routine cardiac catheterizations of eight patients at the German Heart Centre Berlin. The underlying pressure waveforms in BPfe were extracted by the single-channel independent component analysis (SCICA). The accuracy of the SCICA model to estimate the whole cBP waveform was evaluated by the mean absolute error (MAE), the root mean square error (RMSE), the relative RMSE (RRMSE), and the intraclass correlation coefficient (ICC). The agreement between the intra-aortic and estimated parameters including systolic (SBP), diastolic (DBP), mean arterial pressure (MAP), and pulse pressure (PP) was evaluated by the regression and Bland-Altman analyses. Results: The SCICA tool estimated the cBP waveform non-invasively from the intra-arterial BPfe with an MAE of 0.159 ± 1.629, an RMSE of 5.153 ± 0.957 mmHg, an RRMSE of 5.424 ± 1.304%, and an ICC of 0.94, as well as two waveforms contributing to morphological distortion at the femoral artery. The regression analysis showed a strong linear trend between the estimated and intra-aortic SBP, DBP, MAP, and PP with high coefficient of determination R2 of 0.98, 0.99, 0.99, and 0.97 respectively. The Bland-Altman plots demonstrated good agreement between estimated and intra-aortic parameters with a mean error and a standard deviation of difference of -0.54 ± 2.42 mmHg [95% confidence interval (CI): -5.28 to 4.20] for SBP, -1.97 ± 1.62 mmHg (95% CI: -5.14 to 1.20) for DBP, -1.49 ± 1.40 mmHg (95% CI: -4.25 to 1.26) for MAP, and 1.43 ± 2.79 mmHg (95% CI: -4.03 to 6.90) for PP. Conclusions: The SCICA approach is a powerful tool that identifies sources contributing to morphological distortion at peripheral arteries and estimates cBP.

A coupling strategy for a first 3D-1D model of the cardiovascular system to study the effects of pulse wave propagation on cardiac function.

A key factor governing the mechanical performance of the heart is the bidirectional coupling with the vascular system, where alterations in vascular properties modulate the pulsatile load imposed on the heart. Current models of cardiac electromechanics (EM) use simplified 0D representations of the vascular system when coupling to anatomically accurate 3D EM models is considered. However, these ignore important effects related to pulse wave transmission. Accounting for these effects requires 1D models, but a 3D-1D coupling remains challenging. In this work, we propose a novel, stable strategy to couple a 3D cardiac EM model to a 1D model of blood flow in the largest systemic arteries. For the first time, a personalised coupled 3D-1D model of left ventricle and arterial system is built and used in numerical benchmarks to demonstrate robustness and accuracy of our scheme over a range of time steps. Validation of the coupled model is performed by investigating the coupled system's physiological response to variations in the arterial system affecting pulse wave propagation, comprising aortic stiffening, aortic stenosis or bifurcations causing wave reflections. Our first 3D-1D coupled model is shown to be efficient and robust, with negligible additional computational costs compared to 3D-0D models. We further demonstrate that the calibrated 3D-1D model produces simulated data that match with clinical data under baseline conditions, and that known physiological responses to alterations in vascular resistance and stiffness are correctly replicated. Thus, using our coupled 3D-1D model will be beneficial in modelling studies investigating wave propagation phenomena.

Numerical predictions of shear stress and cyclic stretch in pulmonary hypertension due to left heart failure.

Isolated post-capillary pulmonary hypertension (Ipc-PH) occurs due to left heart failure, which contributes to 1 out of every 9 deaths in the United States. In some patients, through unknown mechanisms, Ipc-PH transitions to combined pre-/post-capillary PH (Cpc-PH) and is associated with a dramatic increase in mortality. Altered mechanical forces and subsequent biological signaling in the pulmonary vascular bed likely contribute to the transition from Ipc-PH to Cpc-PH. However, even in a healthy pulmonary circulation, the mechanical forces in the smallest vessels (the arterioles, capillary bed, and venules) have not been quantitatively defined. This study is the first to examine this question via a computational fluid dynamics model of the human pulmonary arteries, arterioles, venules, and veins. Using this model, we predict temporal and spatial dynamics of cyclic stretch and wall shear stress with healthy and diseased hemodynamics. In the normotensive case for large vessels, numerical simulations show that large arteries have higher pressure and flow than large veins, as well as more pronounced changes in area throughout the cardiac cycle. In the microvasculature, shear stress increases and cyclic stretch decreases as vessel radius decreases. When we impose an increase in left atrial pressure to simulate Ipc-PH, shear stress decreases and cyclic stretch increases as compared to the healthy case. Overall, this model predicts pressure, flow, shear stress, and cyclic stretch that providing a way to analyze and investigate hypotheses related to disease progression in the pulmonary circulation.

Hemodynamic modeling of the circle of Willis reveals unanticipated functions during cardiovascular stress.

The circle of Willis (CW) allows blood to be redistributed throughout the brain during local ischemia; however, it is unlikely that the anatomic persistence of the CW across mammalian species is driven by natural selection of individuals with resistance to cerebrovascular disease typically occurring in elderly humans. To determine the effects of communicating arteries (CoAs) in the CW on cerebral pulse wave propagation and blood flow velocity, we simulated young, active adult humans undergoing different states of cardiovascular stress (i.e., fear and aerobic exercise) using discrete transmission line segments with stress-adjusted cardiac output, peripheral resistance, and arterial compliance. Phase delays between vertebrobasilar and carotid pulses allowed bidirectional shunting through CoAs: both posteroanterior shunting before the peak of the pulse waveform and anteroposterior shunting after internal carotid pressure exceeded posterior cerebral pressure. Relative to an absent CW without intact CoAs, the complete CW blunted anterior pulse waveforms, although limited to 3% and 6% reductions in peak pressure and pulse pressure, respectively. Systolic rate of change in pressure (i.e., ∂P/∂t) was reduced 15%-24% in the anterior vasculature and increased 23%-41% in the posterior vasculature. Bidirectional shunting through posterior CoAs was amplified during cardiovascular stress and increased peak velocity by 25%, diastolic-to-systolic velocity range by 44%, and blood velocity acceleration by 134% in the vertebrobasilar arteries. This effect may facilitate stress-related increases in blood flow to the cerebellum (improving motor coordination) and reticular-activating system (enhancing attention and focus) via a nitric oxide-dependent mechanism, thereby improving survival in fight-or-flight situations.NEW & NOTEWORTHY Hemodynamic modeling reveals potential evolutionary benefits of the intact circle of Willis (CW) during fear and aerobic exercise. The CW equalizes pulse waveforms due to bidirectional shunting of blood flow through communicating arteries, which boosts vertebrobasilar blood flow velocity and acceleration. These phenomena may enhance perfusion of the brainstem and cerebellum via nitric oxide-mediated vasodilation, improving performance of the reticular-activating system and motor coordination in survival situations.

Arterial pulse wave propagation across stenoses and aneurysms: assessment of one-dimensional simulations against three-dimensional simulations and in vitro measurements.

One-dimensional (1-D) arterial blood flow modelling was tested in a series of idealized vascular geometries representing the abdominal aorta, common carotid and iliac arteries with different sizes of stenoses and/or aneurysms. Three-dimensional (3-D) modelling and in vitro measurements were used as ground truth to assess the accuracy of 1-D model pressure and flow waves. The 1-D and 3-D formulations shared identical boundary conditions and had equivalent vascular geometries and material properties. The parameters of an experimental set-up of the abdominal aorta for different aneurysm sizes were matched in corresponding 1-D models. Results show the ability of 1-D modelling to capture the main features of pressure and flow waves, pressure drop across the stenoses and energy dissipation across aneurysms observed in the 3-D and experimental models. Under physiological Reynolds numbers (Re), root mean square errors were smaller than 5.4% for pressure and 7.3% for the flow, for stenosis and aneurysm sizes of up to 85% and 400%, respectively. Relative errors increased with the increasing stenosis and aneurysm size, aneurysm length and Re, and decreasing stenosis length. All data generated in this study are freely available and provide a valuable resource for future research.

Improved biomechanical metrics of cerebral vasospasm identified via sensitivity analysis of a 1D cerebral circulation model.

Cerebral vasospasm (CVS) is a life-threatening condition that occurs in a large proportion of those affected by subarachnoid haemorrhage and stroke. CVS manifests itself as the progressive narrowing of intracranial arteries. It is usually diagnosed using Doppler ultrasound, which quantifies blood velocity changes in the affected vessels, but has low sensitivity when CVS affects the peripheral vasculature. The aim of this study was to identify alternative biomarkers that could be used to diagnose CVS. We used a 1D modelling approach to describe the properties of pulse waves that propagate through the cardiovascular system, which allowed the effects of different types of vasospasm on waveforms to be characterised at several locations within a simulated cerebral network. A sensitivity analysis empowered by the use of a Gaussian process statistical emulator was used to identify waveform features that may have strong correlations with vasospasm. We showed that the minimum rate of velocity change can be much more effective than blood velocity for stratifying typical manifestations of vasospasm and its progression. The results and methodology of this study have the potential not only to improve the diagnosis and monitoring of vasospasm, but also to be used in the diagnosis of many other cardiovascular diseases where cardiovascular waves can be decoded to provide disease characterisation.

A novel, FFT-based one-dimensional blood flow solution method for arterial network.

In the present work, we propose an FFT-based method for solving blood flow equations in an arterial network with variable properties and geometrical changes. An essential advantage of this approach is in correctly accounting for the vessel skin friction through the use of Womersley solution. To incorporate nonlinear effects, a novel approximation method is proposed to enable calculation of nonlinear corrections. Unlike similar methods available in the literature, the set of algebraic equations required for every harmonic is constructed automatically. The result is a generalized, robust and fast method to accurately capture the increasing pulse wave velocity downstream as well as steepening of the pulse front. The proposed method is shown to be appropriate for incorporating correct convection and diffusion coefficients. We show that the proposed method is fast and accurate and it can be an effective tool for 1D modelling of blood flow in human arterial networks.

Ultrasound simulation model incorporating incident and reflected wave propagations along a common carotid artery.

An ultrasound simulation model incorporating incident and reflected wave propagations is proposed to provide a realistic data source for validation of transit time (TT)-based local pulse wave velocity (PWV) estimation algorithms. First, the incident wave (IW) and reflected wave (RW) at a certain position over a common carotid artery (CCA) are estimated. Then, the propagating pulse waves (PWs) along the CCA are modelled with the synthesizations of the estimated IWs and RWs, whose occurrences are delayed in opposite sequences according to a preset PWV. In ultrasound simulation, a geometric model of a CCA is built, and the dynamic scatterer models are constructed by moving the scatterer positions according to the synthesized PWs. The RF signals are generated using Field II. To characterize the PW propagations of different arterial stiffnesses consistent with clinical ones in the model, 30 healthy subjects from young, middle-aged, and elderly groups are recruited for extractions of IWs and RWs. To quantitatively verify the effectiveness of the simulation model, the normalized root-mean-squared errors (NRMSEs) are used to compare the estimated and preset PWs, time delays (TDs), and PWVs. Results show that for the three age groups, the estimated PWs, TDs, and PWVs conform to the preset ones with the mean NRMSEs of 0.92%, 18.47%, and 8.55%, respectively. Moreover, the model can characterize the effect of the wave reflection on the local PW propagation as its clinical manifestation. Therefore, the proposed model can be effective as a data source for the validation of TT-based local PWV estimation algorithms, particularly the effects of RWs on the estimation performance.

Uncertainty quantification of inflow boundary condition and proximal arterial stiffness-coupled effect on pulse wave propagation in a vascular network.

This work aims at quantifying the effect of inherent uncertainties from cardiac output on the sensitivity of a human compliant arterial network response based on stochastic simulations of a reduced-order pulse wave propagation model. A simple pulsatile output form is used to reproduce the most relevant cardiac features with a minimum number of parameters associated with left ventricle dynamics. Another source of significant uncertainty is the spatial heterogeneity of the aortic compliance, which plays a key role in the propagation and damping of pulse waves generated at each cardiac cycle. A continuous representation of the aortic stiffness in the form of a generic random field of prescribed spatial correlation is then considered. Making use of a stochastic sparse pseudospectral method, we investigate the sensitivity of the pulse pressure and waves reflection magnitude over the arterial tree with respect to the different model uncertainties. Results indicate that uncertainties related to the shape and magnitude of the prescribed inlet flow in the proximal aorta can lead to potent variation of both the mean value and standard deviation of blood flow velocity and pressure dynamics due to the interaction of different wave propagation and reflection features. Lack of accurate knowledge in the stiffness properties of the aorta, resulting in uncertainty in the pulse wave velocity in that region, strongly modifies the statistical response, with a global increase in the variability of the quantities of interest and a spatial redistribution of the regions of higher sensitivity. These results will provide some guidance in clinical data acquisition and future coupling of arterial pulse wave propagation reduced-order model with more complex beating heart models.

On the impact of modelling assumptions in multi-scale, subject-specific models of aortic haemodynamics.

Simulation of haemodynamics has become increasingly popular within the research community. Irrespective of the modelling approach (zero-dimensional (0D), one-dimensional (1D) or three-dimensional (3D)), in vivo measurements are required to personalize the arterial geometry, material properties and boundary conditions of the computational model. Limitations in in vivo data acquisition often result in insufficient information to determine all model parameters and, hence, arbitrary modelling assumptions. Our goal was to minimize and understand the impact of modelling assumptions on the simulated blood pressure, flow and luminal area waveforms by studying a small region of the systemic vasculature-the upper aorta-and acquiring a rich array of non-invasive magnetic resonance imaging and tonometry data from a young healthy volunteer. We first investigated the effect of different modelling assumptions for boundary conditions and material parameters in a 1D/0D simulation framework. Strategies were implemented to mitigate the impact of inconsistencies in the in vivo data. Average relative errors smaller than 7% were achieved between simulated and in vivo waveforms. Similar results were obtained in a 3D/0D simulation framework using the same inflow and outflow boundary conditions and consistent geometrical and mechanical properties. We demonstrated that accurate subject-specific 1D/0D and 3D/0D models of aortic haemodynamics can be obtained using non-invasive clinical data while minimizing the number of arbitrary modelling decisions.

A Novel Analytical Approach to Pulsatile Blood Flow in the Arterial Network.

Haemodynamic simulations using one-dimensional (1-D) computational models exhibit many of the features of the systemic circulation under normal and diseased conditions. We propose a novel linear 1-D dynamical theory of blood flow in networks of flexible vessels that is based on a generalized Darcy's model and for which a full analytical solution exists in frequency domain. We assess the accuracy of this formulation in a series of benchmark test cases for which computational 1-D and 3-D solutions are available. Accordingly, we calculate blood flow and pressure waves, and velocity profiles in the human common carotid artery, upper thoracic aorta, aortic bifurcation, and a 20-artery model of the aorta and its larger branches. Our analytical solution is in good agreement with the available solutions and reproduces the main features of pulse waveforms in networks of large arteries under normal physiological conditions. Our model reduces computational time and provides a new approach for studying arterial pulse wave mechanics; e.g.,  the analyticity of our model allows for a direct identification of the role played by physical properties of the cardiovascular system on the pressure waves.

Fluid friction and wall viscosity of the 1D blood flow model.

We study the behavior of the pulse waves of water into a flexible tube for application to blood flow simulations. In pulse waves both fluid friction and wall viscosity are damping factors, and difficult to evaluate separately. In this paper, the coefficients of fluid friction and wall viscosity are estimated by fitting a nonlinear 1D flow model to experimental data. In the experimental setup, a distensible tube is connected to a piston pump at one end and closed at another end. The pressure and wall displacements are measured simultaneously. A good agreement between model predictions and experiments was achieved. For amplitude decrease, the effect of wall viscosity on the pulse wave has been shown as important as that of fluid viscosity.

A benchmark study of numerical schemes for one-dimensional arterial blood flow modelling.

Haemodynamical simulations using one-dimensional (1D) computational models exhibit many of the features of the systemic circulation under normal and diseased conditions. Recent interest in verifying 1D numerical schemes has led to the development of alternative experimental setups and the use of three-dimensional numerical models to acquire data not easily measured in vivo. In most studies to date, only one particular 1D scheme is tested. In this paper, we present a systematic comparison of six commonly used numerical schemes for 1D blood flow modelling: discontinuous Galerkin, locally conservative Galerkin, Galerkin least-squares finite element method, finite volume method, finite difference MacCormack method and a simplified trapezium rule method. Comparisons are made in a series of six benchmark test cases with an increasing degree of complexity. The accuracy of the numerical schemes is assessed by comparison with theoretical results, three-dimensional numerical data in compatible domains with distensible walls or experimental data in a network of silicone tubes. Results show a good agreement among all numerical schemes and their ability to capture the main features of pressure, flow and area waveforms in large arteries. All the information used in this study, including the input data for all benchmark cases, experimental data where available and numerical solutions for each scheme, is made publicly available online, providing a comprehensive reference data set to support the development of 1D models and numerical schemes.

A 1D pulse wave propagation model of the hemodynamics of calf muscle pump function.

The calf muscle pump is a mechanism which increases venous return and thereby compensates for the fluid shift towards the lower body during standing. During a muscle contraction, the embedded deep veins collapse and venous return increases. In the subsequent relaxation phase, muscle perfusion increases due to increased perfusion pressure, as the proximal venous valves temporarily reduce the distal venous pressure (shielding). The superficial and deep veins are connected via perforators, which contain valves allowing flow in the superficial-to-deep direction. The aim of this study is to investigate and quantify the physiological mechanisms of the calf muscle pump, including the effect of venous valves, hydrostatic pressure, and the superficial venous system. Using a one-dimensional pulse wave propagation model, a muscle contraction is simulated by increasing the extravascular pressure in the deep venous segments. The hemodynamics are studied in three different configurations: a single artery-vein configuration with and without valves and a more detailed configuration including a superficial vein. Proximal venous valves increase effective venous return by 53% by preventing reflux. Furthermore, the proximal valves shielding function increases perfusion following contraction. Finally, the superficial system aids in maintaining the perfusion during the contraction phase and reduces the refilling time by 37%.

A systematic comparison between 1-D and 3-D hemodynamics in compliant arterial models.

We present a systematic comparison of computational hemodynamics in arteries between a one-dimensional (1-D) and a three-dimensional (3-D) formulation with deformable vessel walls. The simulations were performed using a series of idealized compliant arterial models representing the common carotid artery, thoracic aorta, aortic bifurcation, and full aorta from the arch to the iliac bifurcation. The formulations share identical inflow and outflow boundary conditions and have compatible material laws. We also present an iterative algorithm to select the parameters for the outflow boundary conditions by using the 1-D theory to achieve a desired systolic and diastolic pressure at a particular vessel. This 1-D/3-D framework can be used to efficiently determine material and boundary condition parameters for 3-D subject-specific arterial models with deformable vessel walls. Finally, we explore the impact of different anatomical features and hemodynamic conditions on the numerical predictions. The results show good agreement between the two formulations, especially during the diastolic phase of the cycle.

Novel wave intensity analysis of arterial pulse wave propagation accounting for peripheral reflections.

We present a novel analysis of arterial pulse wave propagation that combines traditional wave intensity analysis with identification of Windkessel pressures to account for the effect on the pressure waveform of peripheral wave reflections. Using haemodynamic data measured in vivo in the rabbit or generated numerically in models of human compliant vessels, we show that traditional wave intensity analysis identifies the timing, direction and magnitude of the predominant waves that shape aortic pressure and flow waveforms in systole, but fails to identify the effect of peripheral reflections. These reflections persist for several cardiac cycles and make up most of the pressure waveform, especially in diastole and early systole. Ignoring peripheral reflections leads to an erroneous indication of a reflection-free period in early systole and additional error in the estimates of (i) pulse wave velocity at the ascending aorta given by the PU-loop method (9.5% error) and (ii) transit time to a dominant reflection site calculated from the wave intensity profile (27% error). These errors decreased to 1.3% and 10%, respectively, when accounting for peripheral reflections. Using our new analysis, we investigate the effect of vessel compliance and peripheral resistance on wave intensity, peripheral reflections and reflections originating in previous cardiac cycles.

Descending aorta subject-specific one-dimensional model validated against in vivo data.

The aorta plays a major role in the cardiovascular system and its function and structure are primarily affected by aging, eating habits, life style and other cardiovascular risk factors, inducing increased stiffness which is associated with cardiovascular and cerebral morbi-mortality. Our objective was to develop and validate a robust subject-specific one-dimensional wave propagation numerical model of the descending aorta. This model with a cross-sectional area, velocity and pressure formulation is built using geometric and hemodynamic data measured on a specific person and is validated against in vivo data acquired on the same subject at three distinct anatomical locations along the thoracic aorta. We studied seven healthy volunteers, who underwent carotid applanation tonometry and aortic cardiovascular magnetic resonance (CMR). Responses of our model in terms of changes in central pressure waveform with arterial alterations were consistent with previously described physiological knowledge. Quantitative validation averaged over the three descending aortic locations and the seven subjects provided low rms errors (given in percentage of the maximal clinical value) between simulated and CMR data, i.e. area: 10±6%, velocity: 11±3%, flow rate: 9±3%. Finally, we also found low rms (5±2%) when comparing simulated pressure in the proximal aortic location against tonometric carotid pressure curves. In conclusion, this simple model performs similar to more complex models of the entire systemic arterial tree at a fraction of the cost, and could be of major usefulness in the non-invasive and local estimation of proximal biomechanical and hemodynamic indices.

Age-dependent changes in elastic properties of thoracic aorta evaluated by magnetic resonance in normal subjects.

OBJECTIVES: Aortic stiffness is an independent cardiovascular risk factor. Cardiac magnetic resonance (CMR) allows evaluation of aortic elastic properties by different indexes such as distensibility, the maximum rate of systolic distension (MRSD) and pulse wave velocity (PWV). We sought to define age-dependent changes of indexes of elastic properties of the thoracic aorta in healthy subjects. METHODS: We enrolled 85 healthy subjects (53 males) free of overt cardiovascular disease subdivided into 6 classes of age (from 15 to >60 years). Distensibility, MRSD and PWV were measured by the analysis of CMR images acquired using a 1.5 T clinical scanner. RESULTS: MRSD and distensibility decreased progressively through the classes of age (P < 0.001) after an initial plateau between 20 and 30 years in males and 15 and 20 years in females. Pulse wave velocity increased progressively with the age (P < 0.001). Distensibility was related to body mass index (P = 0.002), surface area (P < 0.005), weight (P = 0.005) and to left ventricular parameters such as mass index (P < 0.001) and end-diastolic volume index (P = 0.002). MRSD was related to end-diastolic volume index (P < 0.001) but not to body parameters. PWV was not related to body and ventricular parameters. CONCLUSIONS: This study confirmed that physiological ageing is associated with a progressive impairment of the elastic properties of the aortic wall. Results of this study may be useful for the early identification of subjects with impaired aortic wall properties providing referral values of elasticity indexes assessed by CMR in different classes of age.

Alterations in pulse wave propagation reflect the degree of outflow tract banding in HH18 chicken embryos.

Hemodynamic conditions play a critical role in embryonic cardiovascular development, and altered blood flow leads to congenital heart defects. Chicken embryos are frequently used as models of cardiac development, with abnormal blood flow achieved through surgical interventions such as outflow tract (OFT) banding, in which a suture is tightened around the heart OFT to restrict blood flow. Banding in embryos increases blood pressure and alters blood flow dynamics, leading to cardiac malformations similar to those seen in human congenital heart disease. In studying these hemodynamic changes, synchronization of data to the cardiac cycle is challenging, and alterations in the timing of cardiovascular events after interventions are frequently lost. To overcome this difficulty, we used ECG signals from chicken embryos (Hamburger-Hamilton stage 18, ∼3 days of incubation) to synchronize blood pressure measurements and optical coherence tomography images. Our results revealed that, after 2 h of banding, blood pressure and pulse wave propagation strongly depend on band tightness. In particular, while pulse transit time in the heart OFT of control embryos is ∼10% of the cardiac cycle, after banding (35% to 50% band tightness) it becomes negligible, indicating a faster OFT pulse wave velocity. Pulse wave propagation in the circulation is likewise affected; however, pulse transit time between the ventricle and dorsal aorta (at the level of the heart) is unchanged, suggesting an overall preservation of cardiovascular function. Changes in cardiac pressure wave propagation are likely contributing to the extent of cardiac malformations observed in banded hearts.