Aortic hemodynamics assessment prior and after valve sparing reconstruction: A patient-specific 4D flow-based FSI model

Abnormal Cardiac Magnetic Resonance-Derived Ascending Aortic Area Strain Demonstrates Altered Ventriculo-Vascular Function in Marfan Syndrome

PURPOSE

There remains a need for improved imaging markers for risk stratification and treatment guidance in Marfan syndrome (MFS). After aortic root replacement (ARR), vascular remodeling and progressive aneurysm formation can occur due to alterations in up- and downstream wall biomechanics and hemodynamics. We aim to compare the ventriculo-vascular properties of patients with MFS with controls, and investigate the correlation between ascending aortic area strain and descending aortic area strain (DAAS) with other clinical variables.

PATIENTS AND METHODS

Nineteen patients with MFS (47% males), including 6 with ARR were studied. In 26 studies, aortic area strain was measured using cross-sectional cardiac magnetic resonance images at the ascending and proximal descending aortic levels. Left atrial, left ventricular longitudinal, and left ventricle circumferential strain (left atrial longitudinal strain, left ventricular longitudinal strain, and left ventricular circumferential strain, respectively) were measured using cardiac magnetic resonance-feature tracking.

RESULTS

Compared with healthy controls, patients with MFS had significantly impaired left ventricular longitudinal strain and left ventricular circumferential strain (-15.8 ± 4.7 vs -19.7 ± 4.8, P = 0.005, and -17.7 ± 4.0 vs -27.0 ± 4.1, P < 0.001). Left atrial longitudinal strain was comparable between patients with MFS and controls. AAAS was significantly reduced (19.0 [11.9, 23.7] vs 46.1 ± 11.3, P < 0.001), whereas DAAS was not significantly decreased. AAAS and DAAS were negatively correlated with age, whereas no significant associations were identified with left ventricle function indices. No significant differences were observed between the ventriculo-vascular properties of patients with MFS who underwent ARR and those who did not.

CONCLUSION

Patients with MFS demonstrated impaired ventricular and vascular function compared with healthy controls. Further investigations are warranted to determine clinical utility of aortic stiffness indices for predicting primary and repeat aortic events.

Patient-specific non-invasive estimation of the aortic blood pressure waveform by ultrasound and tonometry

BACKGROUND AND OBJECTIVE

Aortic blood pressure (ABP) is a more effective prognostic indicator of cardiovascular disease than peripheral blood pressure. A highly accurate algorithm for non-invasively deriving the ABP wave, based on ultrasonic measurement of aortic flow combined with peripheral pulse wave measurements, has been proposed elsewhere. However, it has remained at the proof-of-concept stage because it requires a priori knowledge of the ABP waveform to calculate aortic pulse wave velocity (PWV). The objective of this study is to transform this proof-of-concept algorithm into a clinically feasible technique.

METHODS

We used the Bramwell-Hill equation to non-invasively calculate aortic PWV which was then used to reconstruct the ABP waveform from non-invasively determined aortic blood flow velocity, aortic diameter, and radial pressure. The two aortic variables were acquired by an ultrasound system from 90 subjects, followed by recordings of radial pressure using a SphygmoCor device. The ABPs estimated by the new algorithm were compared with reference values obtained by cardiac catheterization (invasive validation, 8 subjects aged 62.3 ± 12.7 years) and a SphygmoCor device (non-invasive validation, 82 subjects aged 45.0 ± 17.8 years).

RESULTS

In the invasive comparison, there was good agreement between the estimated and directly measured pressures: the mean error in systolic blood pressure (SBP) was 1.4 ± 0.8 mmHg; diastolic blood pressure (DBP), 0.9 ± 0.8 mmHg; mean blood pressure (MBP), 1.8 ± 1.2 mmHg and pulse pressure (PP), 1.4 ± 1.1 mmHg. In the non-invasive comparison, the estimated and directly measured pressures also agreed well: the errors being: SBP, 2.0 ± 1.4 mmHg; DBP, 0.8 ± 0.1 mmHg; MBP, 0.1 ± 0.1 mmHg and PP, 2.3 ± 1.6 mmHg. The significance of the differences in mean errors between calculated and reference values for SBP, DBP, MBP and PP were assessed by paired t-tests. The agreement between the reference methods and those obtained by applying the new approach was also expressed by correlation and Bland-Altman plots.

CONCLUSION

The new method proposed here can accurately estimate ABP, allowing this important variable to be obtained non-invasively, using standard, well validated measurement techniques. It thus has the potential to relocate ABP estimation from a research environment to more routine use in the cardiac clinic.

SHORT ABSTRACT

A highly accurate algorithm for non-invasively deriving the ABP wave has been proposed elsewhere. However, it has remained at the proof-of-concept stage because it requires a priori knowledge of the ABP waveform to calculate aortic pulse wave velocity (PWV). This study aims to transform this proof-of-concept algorithm into a clinically feasible technique. We used the Bramwell-Hill equation to non-invasively calculate aortic PWV which was then used to reconstruct the ABP waveform. The ABPs estimated by the new algorithm were compared with reference values obtained by cardiac catheterization or a SphygmoCor device. The results showed that there was good agreement between the estimated and directly measured pressures. The new method proposed can accurately estimate ABP, allowing this important variable to be obtained non-invasively, using standard, well validated measurement techniques. It thus has the potential to relocate ABP estimation from a research environment to more routine use in the cardiac clinic.

Patient-Specific Haemodynamic Analysis of Virtual Grafting Strategies in Type-B Aortic Dissection: Impact of Compliance Mismatch

INTRODUCTION

Compliance mismatch between the aortic wall and Dacron Grafts is a clinical problem concerning aortic haemodynamics and morphological degeneration. The aortic stiffness introduced by grafts can lead to an increased left ventricular (LV) afterload. This study quantifies the impact of compliance mismatch by virtually testing different Type-B aortic dissection (TBAD) surgical grafting strategies in patient-specific, compliant computational fluid dynamics (CFD) simulations.

MATERIALS AND METHODS

A post-operative case of TBAD was segmented from computed tomography angiography data. Three virtual surgeries were generated using different grafts; two additional cases with compliant grafts were assessed. Compliant CFD simulations were performed using a patient-specific inlet flow rate and three-element Windkessel outlet boundary conditions informed by 2D-Flow MRI data. The wall compliance was calibrated using Cine-MRI images. Pressure, wall shear stress (WSS) indices and energy loss (EL) were computed.

RESULTS

Increased aortic stiffness and longer grafts increased aortic pressure and EL. Implementing a compliant graft matching the aortic compliance of the patient reduced the pulse pressure by 11% and EL by 4%. The endothelial cell activation potential (ECAP) differed the most within the aneurysm, where the maximum percentage difference between the reference case and the mid (MDA) and complete (CDA) descending aorta replacements increased by 16% and 20%, respectively.

CONCLUSION

This study suggests that by minimising graft length and matching its compliance to the native aorta whilst aligning with surgical requirements, the risk of LV hypertrophy may be reduced. This provides evidence that compliance-matching grafts may enhance patient outcomes.

Early clinical outcomes and molecular smooth muscle cell phenotyping using a prophylactic aortic arch replacement strategy in Loeys-Dietz syndrome

OBJECTIVES

Patients with Loeys-Dietz syndrome demonstrate a heightened risk of distal thoracic aortic events after valve-sparing aortic root replacement. This study assesses the clinical risks and hemodynamic consequences of a prophylactic aortic arch replacement strategy in Loeys-Dietz syndrome and characterizes smooth muscle cell phenotype in Loeys-Dietz syndrome aneurysmal and normal-sized downstream aorta.

METHODS

Patients with genetically confirmed Loeys-Dietz syndrome (n = 8) underwent prophylactic aortic arch replacement during valve-sparing aortic root replacement. Four-dimensional flow magnetic resonance imaging studies were performed in 4 patients with Loeys-Dietz syndrome (valve-sparing aortic root replacement + arch) and compared with patients with contemporary Marfan syndrome (valve-sparing aortic root replacement only, n = 5) and control patients (without aortopathy, n = 5). Aortic tissues from 4 patients with Loeys-Dietz syndrome and 2 organ donors were processed for anatomically segmented single-cell RNA sequencing and histologic assessment.

RESULTS

Patients with Loeys-Dietz syndrome valve-sparing aortic root replacement + arch had no deaths, major morbidity, or aortic events in a median of 2 years follow-up. Four-dimensional magnetic resonance imaging demonstrated altered flow parameters in patients with postoperative aortopathy relative to controls, but no clear deleterious changes due to arch replacement. Integrated analysis of aortic single-cell RNA sequencing data (>49,000 cells) identified a continuum of abnormal smooth muscle cell phenotypic modulation in Loeys-Dietz syndrome defined by reduced contractility and enriched extracellular matrix synthesis, adhesion receptors, and transforming growth factor-beta signaling. These modulated smooth muscle cells populated the Loeys-Dietz syndrome tunica media with gradually reduced density from the overtly aneurysmal root to the nondilated arch.

CONCLUSIONS

Patients with Loeys-Dietz syndrome demonstrated excellent surgical outcomes without overt downstream flow or shear stress disturbances after concomitant valve-sparing aortic root replacement + arch operations. Abnormal smooth muscle cell-mediated aortic remodeling occurs within the normal diameter, clinically at-risk Loeys-Dietz syndrome arch segment. These initial clinical and pathophysiologic findings support concomitant arch replacement in Loeys-Dietz syndrome.

Risk of Type B Dissection in Marfan Syndrome: The Cornell Aortic Aneurysm Registry

BACKGROUND

With preventive aortic grafting decreasing the incidence of type A dissections in Marfan syndrome (MFS), most dissections are now type B, for which risk factors remain largely uncertain.

OBJECTIVES

We explored the determinants of type B dissection risk in a large, single-center MFS registry.

METHODS

Demographic and anthropometric features, cardiovascular disease, and surgical history were compared in patients with MFS with and without type B dissection.

RESULTS

Of 336 patients with MFS, 47 (14%) experienced a type B dissection (vs type A in 9%). Patients with type B dissection were more likely to have undergone elective aortic root replacement (ARR) (79 vs 46%; P < 0.001). Of the patients, 55% had type B dissection a mean of 13.3 years after ARR, whereas 45% experienced type B dissection before or in the absence of ARR; 41 patients (87%) were aware of their MFS diagnosis before type B dissection. Among those with predissection imaging, the descending aorta was normal or minimally dilated (<4.0 cm) in 88%. In multivariable analyses, patients with type B dissection were more likely to have undergone ARR and independent mitral valve surgery, to have had a type II dissection, and to have lived longer.

CONCLUSIONS

In our contemporary cohort, type B dissections are more common than type A dissections and occur at traditional nonsurgical thresholds. The associations of type B dissection with ARR, independent mitral valve surgery, and type II dissection suggest a more severe phenotype in the setting of prolonged life expectancy.

Comparative study between 1-way and 2-way coupled fluid-structure interaction in numerical simulation of aortic arch aneurysms

Hemodynamic forces are related to pathological variations of the cardiovascular system, and numerical simulations for fluid-structure interaction have been systematically used to analyze the behavior of blood flow and the arterial wall in aortic aneurysms. This paper proposes a comparative analysis of 1-way and 2-way coupled fluid-structure interaction for aortic arch aneurysm. The coupling models of fluid-structure interaction were conducted using 3D geometry of the thoracic aorta from computed tomography. Hyperelastic anisotropic properties were estimated for the Holzapfel arterial wall model. The rheological behavior of the blood was modeled by the Carreau-Yasuda model. The results showed that the 1-way approach tends to underestimate von Mises stress, displacement, and strain over the entire cardiac cycle, compared to the 2-way approach. In contrast, the behavior of the variables of flow field, velocity, wall shear stress, and Reynolds number when coupled by the 1-way model was overestimated at the systolic moment and tends to be equal at the diastolic moment. The quantitative differences found, especially during the systole, suggest the use of 2-way coupling in numerical simulations of aortic arch aneurysms due to the hyperelastic nature of the arterial wall, which leads to a strong iteration between the fluid and the arterial wall.

Multiscale model for blood flow after a bileaflet artificial aortic valve implantation

Cardiovascular diseases are the leading cause of mortality in the world, mainly due to atherosclerosis and its consequences. The article presents the numerical model of the blood flow through artificial aortic valve. The overset mesh approach was applied to simulate the valve leaflets motion and to realize the moving mesh, in the aortic arch and the main branches of cardiovascular system. To capture the cardiac system's response and the effect of vessel compliance on the outlet pressure, the lumped parameter model has been also included within the solution procedure. Three different turbulence modeling approaches were used and compared - the laminar, k-ϵ and k-ω model. The simulation results were also compared with the model excluding the moving valve geometry and the importance of the lumped parameter model for the outlet boundary condition was analyzed. Proposed numerical model and protocol was found as suitable for performing the virtual operations on the real patient vasculature geometry. The time-efficient turbulence model and overall solving procedure allows to support the clinicians in making decisions about the patient treatment and to predict the results of the future surgery.

Data-driven generation of 4D velocity profiles in the aneurysmal ascending aorta

BACKGROUND AND OBJECTIVE

Numerical simulations of blood flow are a valuable tool to investigate the pathophysiology of ascending thoratic aortic aneurysms (ATAA). To accurately reproduce in vivo hemodynamics, computational fluid dynamics (CFD) models must employ realistic inflow boundary conditions (BCs). However, the limited availability of in vivo velocity measurements, still makes researchers resort to idealized BCs. The aim of this study was to generate and thoroughly characterize a large dataset of synthetic 4D aortic velocity profiles sampled on a 2D cross-section along the ascending aorta with features similar to clinical cohorts of patients with ATAA.

METHODS

Time-resolved 3D phase contrast magnetic resonance (4D flow MRI) scans of 30 subjects with ATAA were processed through in-house code to extract anatomically consistent cross-sectional planes along the ascending aorta, ensuring spatial alignment among all planes and interpolating all velocity fields to a reference configuration. Velocity profiles of the clinical cohort were extensively characterized by computing flow morphology descriptors of both spatial and temporal features. By exploiting principal component analysis (PCA), a statistical shape model (SSM) of 4D aortic velocity profiles was built and a dataset of 437 synthetic cases with realistic properties was generated.

RESULTS

Comparison between clinical and synthetic datasets showed that the synthetic data presented similar characteristics as the clinical population in terms of key morphological parameters. The average velocity profile qualitatively resembled a parabolic-shaped profile, but was quantitatively characterized by more complex flow patterns which an idealized profile would not replicate. Statistically significant correlations were found between PCA principal modes of variation and flow descriptors.

CONCLUSIONS

We built a data-driven generative model of 4D aortic inlet velocity profiles, suitable to be used in computational studies of blood flow. The proposed software system also allows to map any of the generated velocity profiles to the inlet plane of any virtual subject given its coordinate set.

A systematic pan-cancer analysis reveals the clinical prognosis and immunotherapy value of C-X3-C motif ligand 1 (CX3CL1)

It is now widely known that C-X3-C motif ligand 1 (CX3CL1) plays an essential part in the process of regulating pro-inflammatory cells migration across a wide range of inflammatory disorders, including a number of malignancies. However, there has been no comprehensive study on the correlation between CX3CL1 and cancers on the basis of clinical features. In order to investigate the potential function of CX3CL1 in the clinical prognosis and immunotherapy, I evaluated the expression of CX3CL1 in numerous cancer types, methylation levels and genetic alterations. I found CX3CL1 was differentially expressed in numerous cancer types, which indicated CX3CL1 may plays a potential role in tumor progression. Furthermore, CX3CL1 was variably expressed in methylation levels and gene alterations in most cancers according to The Cancer Genome Atlas (TCGA). CX3CL1 was robustly associated with clinical characteristics and pathological stages, suggesting that it was related to the degree of tumor malignancy and the physical function of patients. As determined by the Kaplan-Meier method of estimating survival, high CX3CL1 expression was associated with either favorable or unfavorable outcomes depending on the different types of cancer. It suggests the correlation between CX3CL1 and tumor prognosis. Significant positive correlations of CX3CL1 expression with CD4⁺ T cells, M1 macrophage cells and activated mast cells have been established in the majority of TCGA malignancies. Which indicates CX3CL1 plays an important role in tumor immune microenvironment. Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis suggested that the chemokine signaling pathway may shed light on the pathway for CX3CL1 to exert function. In a conclusion, our study comprehensively summarizes the potential role of CX3CL1 in clinical prognosis and immunotherapy, suggesting that CX3CL1 may represent a promising pharmacological treatment target of tumors.

Synthesis of patient-specific multipoint 4D flow MRI data of turbulent aortic flow downstream of stenotic valves

We propose to synthesize patient-specific 4D flow MRI datasets of turbulent flow paired with ground truth flow data to support training of inference methods. Turbulent blood flow is computed based on the Navier-Stokes equations with moving domains using realistic boundary conditions for aortic shapes, wall displacements and inlet velocities obtained from patient data. From the simulated flow, synthetic multipoint 4D flow MRI data is generated with user-defined spatiotemporal resolutions and reconstructed with a Bayesian approach to compute time-varying velocity and turbulence maps. For MRI data synthesis, a fixed hypothetical scan time budget is assumed and accordingly, changes to spatial resolution and time averaging result in corresponding scaling of signal-to-noise ratios (SNR). In this work, we focused on aortic stenotic flow and quantification of turbulent kinetic energy (TKE). Our results show that for spatial resolutions of 1.5 and 2.5 mm and time averaging of 5 ms as encountered in 4D flow MRI in practice, peak total turbulent kinetic energy downstream of a 50, 75 and 90% stenosis is overestimated by as much as 23, 15 and 14% (1.5 mm) and 38, 24 and 23% (2.5 mm), demonstrating the importance of paired ground truth and 4D flow MRI data for assessing accuracy and precision of turbulent flow inference using 4D flow MRI exams.

The impact of 4D-Flow MRI spatial resolution on patient-specific CFD simulations of the thoracic aorta

Magnetic Resonance Imaging (MRI) is considered the gold standard of medical imaging technologies as it allows for accurate imaging of blood vessels. 4-Dimensional Flow Magnetic Resonance Imaging (4D-Flow MRI) is built on conventional MRI, and provides flow data in the three vector directions and a time resolved magnitude data set. As such it can be used to retrospectively calculate haemodynamic parameters of interest, such as Wall Shear Stress (WSS). However, multiple studies have indicated that a significant limitation of the imaging technique is the spatiotemporal resolution that is currently available. Recent advances have proposed and successfully integrated 4D-Flow MRI imaging techniques with Computational Fluid Dynamics (CFD) to produce patient-specific simulations that have the potential to aid in treatments,surgical decision making, and risk stratification. However, the consequences of using insufficient 4D-Flow MRI spatial resolutions on any patient-specific CFD simulations is currently unclear, despite being a recognised limitation. The research presented in this study aims to quantify the inaccuracies in patient-specific 4D-Flow MRI based CFD simulations that can be attributed to insufficient spatial resolutions when acquiring 4D-Flow MRI data. For this research, a patient has undergone four 4D-Flow MRI scans acquired at various isotropic spatial resolutions and patient-specific CFD simulations have subsequently been run using geometry and velocity data produced from each scan. It was found that compared to CFD simulations based on a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.5\,{\text {mm}} \times 1.5\,{\text {mm}} \times 1.5\,{\text {mm}}$$\end{document}1.5mm×1.5mm×1.5mm, using a spatial resolution of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4\,{\text {mm}} \times 4\,{\text {mm}} \times 4\,{\text {mm}}$$\end{document}4mm×4mm×4mm substantially underestimated the maximum velocity magnitude at peak systole by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$110.55\%$$\end{document}110.55%. The impacts of 4D-Flow MRI spatial resolution on WSS calculated from CFD simulations have been investigated and it has been shown that WSS is underestimated in CFD simulations that are based on a coarse 4D-Flow MRI spatial resolution. The authors have concluded that a minimum 4D-Flow MRI spatial resolution of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.5\,{\text {mm}} \times 1.5\,{\text {mm}} \times 1.5\,{\text {mm}}$$\end{document}1.5mm×1.5mm×1.5mm must be used when acquiring 4D-Flow MRI data to perform patient-specific CFD simulations. A coarser spatial resolution will produce substantial differences within the flow field and geometry.

Impact of extra-anatomical bypass on coarctation fluid dynamics using patient-specific lumped parameter and Lattice Boltzmann modeling

Accurate hemodynamic analysis is not only crucial for successful diagnosis of coarctation of the aorta (COA), but intervention decisions also rely on the hemodynamics assessment in both pre and post intervention states to minimize patient risks. Despite ongoing advances in surgical techniques for COA treatments, the impacts of extra-anatomic bypass grafting, a surgical technique to treat COA, on the aorta are not always benign. Our objective was to investigate the impact of bypass grafting on aortic hemodynamics. We investigated the impact of bypass grafting on aortic hemodynamics using a patient-specific computational-mechanics framework in three patients with COA who underwent bypass grafting. Our results describe that bypass grafting improved some hemodynamic metrics while worsened the others: (1) Doppler pressure gradient improved (decreased) in all patients; (2) Bypass graft did not reduce the flow rate substantially through the COA; (3) Systemic arterial compliance increased in patients #1 and 3 and didn't change (improve) in patient 3; (4) Hypertension got worse in all patients; (5) The flow velocity magnitude improved (reduced) in patient 2 and 3 but did not improve significantly in patient 1; (6) There were elevated velocity magnitude, persistence of vortical flow structure, elevated turbulence characteristics, and elevated wall shear stress at the bypass graft junctions in all patients. We concluded that bypass graft may lead to pseudoaneurysm formation and potential aortic rupture as well as intimal hyperplasia due to the persistent abnormal and irregular aortic hemodynamics in some patients. Moreover, post-intervention, exposures of endothelial cells to high shear stress may lead to arterial remodeling, aneurysm, and rupture.

Ultrasound-based method for individualized estimation of central aortic blood pressure from flow velocity and diameter

Central aortic blood pressure (CABP) is a better predictor for cardiovascular events than brachial blood pressure. However, direct CABP measurement is invasive. The objective of this paper is to develop an ultrasound-based method using individualized Windkessel (WK) models for non-invasive estimation of CABP. Three WK models (with two-, three- and four-element WK, named, WK2, WK3 and WK4, respectively) were created and the model parameters were individualized based on aortic flow velocity and diameter waveforms measured by ultrasound (US). Experimental data were acquired in 42 subjects aged 21-67 years. The CABP estimated by WK models was compared with the reference CABP obtained using a commercial system. The results showed that the overall performance of the WK3 and WK4 models was similar, outperforming the WK2 model. The estimated CABP based on WK3/WK4 model showed good agreement with the reference CABP: the absolute errors of systolic blood pressure (SBP), 2.4 ± 2.1/2.4 ± 2.0 mmHg; diastolic blood pressure (DBP), 1.4 ± 1.1/1.7 ± 1.5 mmHg; mean blood pressure (MBP), 1.3 ± 0.8/1.3 ± 0.8 mmHg; pulse pressure (PP), 3.0 ± 2.3/3.2 ± 2.6 mmHg; the root mean square error (RMSE) of the waveforms, 2.5 ± 1.0/2.6 ± 1.1 mmHg. Therefore, the proposed method can provide a non-invasive CABP estimation during routine cardiac US examination.