Effects of severity and location of stenosis on the hemodynamics in human aorta and its branches

Automatic Laplacian-based shape optimization for patient-specific vascular grafts

Cognitional heart disease is one of the leading causes of mortality among newborns. Tissue-engineered vascular grafts offer the potential to help treat cognitional heart disease through patient-specific vascular grafts. However, current methods often rely on non-personalized designs or involve significant human intervention. This paper presents a computational framework for the automatic shape optimization of patient-specific tissue-engineered vascular grafts for repairing the aortic arch, aimed at reducing the need for manual input and improving current treatment outcomes, which either use non-patient-specific geometry or require extensive human intervention to design the vascular graft. The paper's core innovation lies in an automatic shape optimization pipeline that combines Bayesian optimization techniques with the open-source finite volume solver, OpenFOAM, and a novel graft deformation algorithm. Specifically, our framework begins with Laplacian mode computation and the approximation of a computationally low-cost Gaussian process surrogate model to capture the minimum weighted combination of inlet-outlet pressure drop (PD) and maximum wall shear stress (WSS). Bayesian Optimization then performs a limited number of OpenFOAM simulations to identify the optimal patient-specific shape. We use imaging and flow data obtained from six patients diagnosed with cognitional heart disease to evaluate our approach. Our results showcase the potential of online training and hemodynamic surrogate model optimization for providing optimal graft shapes. These results show how our framework successfully reduces inlet-outlet PD and maximum WSS compared to pre-lofted models that include both the native geometry and human-designed grafts. Furthermore, we compare how the performance of each design optimized under steady-state simulation compares to that design's performance under transient simulation, and to what extent the optimal design remains similar under both conditions. Our findings underscore that the automated designs achieve at least a 16% reduction in blood flow pressure drop in comparison to geometries optimized by humans.

Distribution of Arteriosclerotic Vessels in Patients with Arteriosclerosis and the Differences of Serum Lipid Levels Classified by Different Sites

Objective

To investigate the distribution of arteriosclerotic vessels of arteriosclerosis, differential serum lipid profiles, and differences in the proportion of dyslipidaemia between patients with single-site arteriosclerosis and multi-site arteriosclerosis (significant hardening of ≥2 arteries).

Methods

The data of 6581 single-site arteriosclerosis patients and 5940 multi-site arteriosclerosis patients were extracted from the hospital medical record system. Serum total cholesterol (TC), triglycerides (TGs), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), apolipoprotein (Apo) A1, ApoB concentrations and C-reactive protein (CRP) between patients with single-site arteriosclerosis and multi-site arteriosclerosis were collected and analyzed.

Results

The most diseased arteries were coronary arteries (n=7099, 33.7%), limb arteries (n=6546, 31.1%), and carotid arteries (n=5279, 25.1%). TC, LDL-C, TC/HDL-C, and LDL-C/HDL-C levels were higher and CRP level was lower in multi-site arteriosclerosis patients than those in single-site arteriosclerosis patients. The TC, LDL-C levels in non-elderly (<65 years old) female patients were higher and TG/HDL-C, TC/HDL-C, LDL-C/HDL-C levels were lower than those in non-elderly male patients, while the TG, TC, LDL-C, and TG/HDL-C levels in elderly (≥65 years old) female patients were higher and LDL-C/HDL-C level was lower than those in elderly male patients. The proportion of dyslipidemia in descending order was as follows: low HDL-C (31.9%), elevated TG (16.9%), elevated TC (9.0%), and elevated LDL-C (4.2%). The levels of TC, LDL-C, TC/HDL-C, and LDL-C/HDL-C in patients with peripheral arteriosclerosis were higher than those in patients with cardio-cerebrovascular arteriosclerosis.

Conclusion

There were differences in serum lipid levels in patients with arteriosclerosis with different age, gender and distribution of arteriosclerotic vessels.

Blood flow effects in a patient with a thoracic aortic endovascular prosthesis

This work analyzes hemodynamic phenomena within the aorta of two elderly patients and their impact on blood flow behavior, particularly affected by an endovascular prosthesis in one of them (Patient II). Computational Fluid Dynamics (CFD) was utilized for this study, involving measurements of velocity, pressure, and wall shear stress (WSS) at various time points during the third cardiac cycle, at specific positions within two cross sections of the thoracic aorta. The first cross-section (Cross-Section 1, CS1) is located before the initial fluid bifurcation, just before the right subclavian artery. The second cross-section (Cross-Section 2, CS2) is situated immediately after the left subclavian artery. The results reveal that, under regular aortic geometries, velocity and pressure magnitudes follow the principles of fluid dynamics, displaying variations. However, in Patient II, an endoprosthesis near the CS2 and the proximal border of the endoprosthesis significantly disrupts fluid behavior owing to the pulsatile flow. The cross-sectional areas of Patient I are smaller than those of Patient II, leading to higher flow magnitudes. Although in CS1 of Patient I, there is considerable variability in velocity magnitudes, they exhibit a more uniform and predictable transition. In contrast, CS2 of Patient II, where magnitude variation is also high, displays irregular fluid behavior due to the endoprosthesis presence. This cross-section coincides with the border of the fluid bifurcation. Additionally, the irregular geometry caused by endovascular aneurysm repair contributes to flow disruption as the endoprosthesis adjusts to the endothelium, reshaping itself to conform with the vessel wall. In this context, significant alterations in velocity values, pressure differentials fluctuating by up to 10%, and low wall shear stress indicate the pronounced influence of the endovascular prosthesis on blood flow behavior. These flow disturbances, when compounded by the heart rate, can potentially lead to changes in vascular anatomy and displacement, resulting in a disruption of the prosthesis-endothelium continuity and thereby causing clinical complications in the patient.

Impact of disturbed flow and arterial stiffening on mechanotransduction in endothelial cells

Disturbed flow promotes progression of atherosclerosis at particular regions of arteries where the recent studies show the arterial wall becomes stiffer. Objective of this study is to show how mechanotransduction in subcellular organelles of endothelial cells (ECs) will alter with changes in blood flow profiles applied on ECs surface and mechanical properties of arterial wall where ECs are attached to. We will examine the exposure of ECs to atherogenic flow profiles (disturbed flow) and non-atherogenic flow profiles (purely forward flow), while stiffness and viscoelasticity of arterial wall will change. A multicomponent model of endothelial cell monolayer was applied to quantify the response of subcellular organelles to the changes in their microenvironment. Our results show that arterial stiffening alters mechanotransduction in intra/inter-cellular organelles of ECs by slight increase in the transmitted stresses, particularly over central stress fibers (SFs). We also observed that degradation of glycocalyx and exposure to non-atherogenic flow profiles result in significantly higher stresses in subcellular organelles, while degradation of glycocalyx and exposure to atherogenic flow profiles result in dramatically lower stresses in the organelles. Moreover, we show that increasing the arterial wall viscoelasticity leads to slight increase in the stresses transmitted to subcellular organelles. FAs are particularly influenced with the changes in the arterial wall properties and viscoelasticity. Our study suggests that changes in viscoelasticity of arterial wall and degradation state of glycocalyx have to be considered along with arterial stiffening in designing more efficient treatment strategies for atherosclerosis. Our study provides insight into significant role of mechanotransduction in the localization of atherosclerosis by quantifying the role of ECs mechanosensors and suggests that mechanotransduction may play a key role in design of more efficient and precision therapeutics to slow down or block the progression of atherosclerosis.

ALDH2 rs671 and MTHFR rs1801133 polymorphisms are risk factors for arteriosclerosis in multiple arteries

BACKGROUND

Arteriosclerosis in multiple arteries has long been associated with heightened cardiovascular risk. Acetaldehyde dehydrogenase 2 (ALDH2) and methylenetetrahydrofolate reductase (MTHFR) play an important role in the pathogenesis of arteriosclerosis by participating in the oxidation and reduction reactions in vascular endothelial cells. The purpose was to investigate the relationship of ALDH2 and MTHFR gene polymorphisms with arteriosclerosis in multiple arteries.

METHODS

410 patients with arteriosclerosis in single artery and 472 patients with arteriosclerosis in multiple arteries were included. The relationship between ALDH2 rs671 and MTHFR rs1801133 polymorphisms and arteriosclerosis in single artery and arteriosclerosis in multiple arteries was analyzed.

RESULTS

The proportion of ALDH2 rs671 A allele (35.6% vs. 30.9%, P = 0.038) and MTHFR rs1801133 T allele (32.6% vs. 27.1%, P = 0.012) in patients with arteriosclerosis in multiple arteries was significantly higher than that in arteriosclerosis in single artery, respectively. The proportion of history of alcohol consumption in patients with ALDH2 rs671 G/G genotype was higher than those in ALDH2 rs671 G/A genotype and A/A genotype (P < 0.001). The results of logistic regression analysis indicated that ALDH2 rs671 A/A genotype (A/A vs. G/G: OR 1.996, 95% CI: 1.258-3.166, P = 0.003) and MTHFR rs1801133 T/T genotype (T/T vs. C/C: OR 1.943, 95% CI: 1.179-3.203, P = 0.009) may be independent risk factors for arteriosclerosis in multiple arteries (adjusted for age, sex, smoking, drinking, hypertension, and diabetes).

CONCLUSIONS

ALDH2 rs671 A/A and MTHFR rs1801133 T/T genotypes may be independent risk factors for arteriosclerosis in multiple arteries.

Non-invasive diagnostics of blockage growth in the descending aorta-computational approach

Atherosclerosis causes blockages to the main arteries such as the aorta preventing blood flow from delivering oxygen to the organs. Non-invasive diagnosis of these blockages is difficult, particularly in primary healthcare. In this paper, the effect of arterial blockage development and growth is investigated at the descending aorta on some possible non-invasive assessment parameters including the blood pressure waveform, wall shear stress (WSS), time-average WSS (TAWSS) and the oscillation shear index (OSI). Blockage severity growth is introduced in a simulation model as 25%, 35%, 50% and 65% stenosis at the descending aorta based on specific healthy control aorta data clinically obtained. A 3D aorta model with invasive pulsatile waveforms (blood flow and pressure) is used in the CFD simulation. Blockage severity is assessed by using blood pressure measurements at the left subclavian artery. An arterial blockage growth more than 35% of the lumen diameter significantly affects the pressure. A strong correlation is also observed between the ascending aorta pressure values, pressure at the left subclavian artery and the relative residence time (RRT). An increase of RRT downstream from the stenosis indicates a 35% stenosis at the descending aorta which results in high systolic and diastolic pressure readings. The findings of this study could be further extended by transferring the waveform reading from the left subclavian artery to the brachial artery.

Graphical abstract

Numerical investigation of patient-specific thoracic aortic aneurysms and comparison with normal subject via computational fluid dynamics (CFD)

Vascular hemodynamics play an important role in cardiovascular diseases. This work aimed to investigate the effects of an increase in ascending aortic diameter (AAD) on hemodynamics throughout a cardiac cycle for real patients. In this study, two scans of thoracic aortic aneurysm (TAA) subject with different AADs (42.94 mm and 48.01 mm) and a scan of a normal subject (19.81 mm) were analyzed to assess the effects of hemodynamics on the progression of TAA with the same flow rate. Real-patient aortic geometries were scanned by computed tomography angiography (CTA), and steady and pulsatile flow conditions were used to simulate real patient aortic geometries. Aortic arches were obtained from routine clinical scans. Computational fluid dynamics (CFD) simulations were performed with in vivo boundary conditions, and 3D Navier-Stokes equations were solved by a UDF (user-defined function) code defining a real cardiac cycle of one patient using Fourier series (FS). Wall shear stress (WSS) and pressure distributions were presented from normal subject to TAA cases. The results show that during the peak systolic phase pressure load increased by 18.56% from normal subject to TAA case 1 and by 23.8% from normal subject to TAA case 2 in the aneurysm region. It is concluded that although overall WSS increased in aneurysm cases but was low in dilatation areas. As a result, abnormal changes in WSS and higher pressure load may lead to rupture and risk of further dilatation. CFD simulations were highly effective to guide clinical predictions and assess the progress of aneurysm regions in case of early surgical intervention. Graphical abstract.

Fluid Structure Interaction on Paravalvular Leakage of Transcatheter Aortic Valve Implantation Related to Aortic Stenosis: A Patient-Specific Case

This study investigated the impact of paravalvular leakage (PVL) in relation to the different valve openings of the transcatheter aortic valve implantation (TAVI) valve using the fluid structure interaction (FSI) approach. Limited studies were found on the subject of FSI with regards to TAVI-PVL condition, which involves both fluid and structural responses in coupling interaction. Hence, further FSI simulation with the two-way coupling method is implemented to investigate the effects of hemodynamics blood flow along the patient-specific aorta model subjected to the interrelationship between PVL and the different valve openings using the established FSI software ANSYS 16.1. A 3D patient-specific aorta model is constructed using MIMICS software. The TAVI valve identical to Edward SAPIEN XT 26 (Edwards Lifesciences, Irvine, California), at different Geometrical Orifice Areas (GOAs), is implanted into the patient's aortic annulus. The leaflet opening of the TAVI valve is drawn according to severity of GOA opening represented in terms of 100%, 80%, 60%, and 40% opening, respectively. The result proved that the smallest percentage of GOA opening produced the highest possibility of PVL, increased the recirculatory flow proximally to the inner wall of the ascending aorta, and produced lower backflow velocity streamlines through the side area of PVL region. Overall, 40% GOA produced 89.17% increment of maximum velocity magnitude, 19.97% of pressure drop, 65.70% of maximum WSS magnitude, and a decrement of 33.62% total displacement magnitude with respect to the 100% GOA.

Geometric determinants of local hemodynamics in severe carotid artery stenosis

In cases of severe carotid artery stenosis (CAS), carotid endarterectomy (CEA) is performed to recover lumen patency and alleviate stroke risk. Under current guidelines, the decision to surgically intervene relies primarily on the percent loss of native arterial lumen diameter within the stenotic region (i.e. the degree of stenosis). An underlying premise is that the degree of stenosis modulates flow-induced wall shear stress elevations at the lesion site, and thus indicates plaque rupture potential and stroke risk. Here, we conduct a retrospective study on pre-CEA computed tomography angiography (CTA) images from 50 patients with severe internal CAS (>60% stenosis) to better understand the influence of plaque and local vessel geometry on local hemodynamics, with geometrical descriptors that extend beyond the degree of stenosis. We first processed CTA images to define a set of multipoint geometric metrics characterizing the stenosed region, and next performed computational fluid dynamics simulations to quantify local wall shear stress and associated hemodynamic metrics. Correlation and regression analyses were used to relate obtained geometric and hemodynamic metrics, with inclusion of patient sub-classification based on the degree of stenosis. Our results suggest that in the context of severe CAS, prediction of shear stress-based metrics can be enhanced by consideration of readily available, multipoint geometric metrics in addition to the degree of stenosis.

Mechanotransmission in endothelial cells subjected to oscillatory and multi-directional shear flow

Local haemodynamics are linked to the non-uniform distribution of atherosclerosic lesions in arteries. Low and oscillatory (reversing in the axial flow direction) wall shear stress (WSS) induce inflammatory responses in endothelial cells (ECs) mediating disease localization. The objective of this study is to investigate computationally how the flow direction (reflected in WSS variation on the EC surface over time) influences the forces experienced by structural components of ECs that are believed to play important roles in mechanotransduction. A three-dimensional, multi-scale, multi-component, viscoelastic model of focally adhered ECs is developed, in which oscillatory WSS (reversing or non-reversing) parallel to the principal flow direction, or multi-directional oscillatory WSS with reversing axial and transverse components are applied over the EC surface. The computational model includes the glycocalyx layer, actin cortical layer, nucleus, cytoskeleton, focal adhesions (FAs), stress fibres and adherens junctions (ADJs). We show the distinct effects of atherogenic flow profiles (reversing unidirectional flow and reversing multi-directional flow) on subcellular structures relative to non-atherogenic flow (non-reversing flow). Reversing flow lowers stresses and strains due to viscoelastic effects, and multi-directional flow alters stress on the ADJs perpendicular to the axial flow direction. The simulations predict forces on integrins, ADJ filaments and other substructures in the range that activate mechanotransduction.

Microparticles in patients undergoing transcatheter aortic valve implantation (TAVI)

Degenerative aortic stenosis (AS) is the most frequent form of acquired valvular heart disease. AS is known to entail endothelial dysfunction caused by increased mechanical shear stress leading to elevated circulatory levels of microparticles. Endothelial and platelet microparticles (EMP and PMP) are small vesicles that originate from activated cells and thrombocytes. We sought to evaluate whether transcatheter aortic valve implantation (TAVI) procedure would elicit effects on circulating EMP and PMP. 92 patients undergoing TAVI procedure for severe AS were included in this study. Samples were obtained at each visit before TAVI, 1 week post-procedure and at 1, 3 and after 6 months after TAVI and were evaluated using flow cytometry. A 12 month clinical follow-up was also performed. CD62E+ EMP concentration before TAVI was 21.11 % (±6.6 % SD) and declined to 20.99 % (±6.8 % SD) after 1 week, to 16.63 % (±5.4 % SD, p < 0.0001) after 1 month, to 17.08 % (±4.6 % SD, p < 0.0001) after 3 months and to 15.94 % (±5.4 % SD, p < 0.0001) after 6 months. CD31+/CD42b-, CD31+/Annexin+/- EMP remained unchanged. CD31+/CD41b+ PMP evidenced a slight, but statistically significant increase after TAVI and remained elevated during the entire follow-up. Apart from a procedure-related improvement in echocardiographic parameters, TAVI procedure led also to a decline in CD62E+ EMP. The reduction in pressure gradients with less hemodynamic shear stress seems also to have beneficially affected endothelial homeostasis.