Coronary artery wall shear stress is associated with progression and transformation of atherosclerotic plaque and arterial remodeling in patients with coronary artery disease

Semi-Automatic Reconstruction of Patient-Specific Stented Coronaries based on Data Assimilation and Computer Aided Design

PURPOSE

The interplay between geometry and hemodynamics is a significant factor in the development of cardiovascular diseases. This is particularly true for stented coronary arteries. To elucidate this factor, an accurate patient-specific analysis requires the reconstruction of the geometry following the stent deployment for a computational fluid dynamics (CFD) investigation. The image-based reconstruction is troublesome for the different possible positions of the stent struts in the lumen and the coronary wall. However, the accurate inclusion of the stent footprint in the hemodynamic analysis is critical for detecting abnormal stress conditions and flow disturbances, particularly for thick struts like in bioresorbable scaffolds. Here, we present a novel reconstruction methodology that relies on Data Assimilation and Computer Aided Design.

METHODS

The combination of the geometrical model of the undeployed stent and image-based data assimilated by a variational approach allows the highly automated reconstruction of the skeleton of the stent. A novel approach based on computational mechanics defines the map between the intravascular frame of reference (called L-view) and the 3D geometry retrieved from angiographies. Finally, the volumetric expansion of the stent skeleton needs to be self-intersection free for the successive CFD studies; this is obtained by using implicit representations based on the definition of Nef-polyhedra.

RESULTS

We assessed our approach on a vessel phantom, with less than 10% difference (properly measured) vs. a customized manual (and longer) procedure previously published, yet with a significant higher level of automation and a shorter turnaround time. Computational hemodynamics results were even closer. We tested the approach on two patient-specific cases as well.

CONCLUSIONS

The method presented here has a high level of automation and excellent accuracy performances, so it can be used for larger studies involving patient-specific geometries.

Hemodynamics derived from computational fluid dynamics based on magnetic resonance angiography is associated with functional outcomes in atherosclerotic middle cerebral artery stenosis

BACKGROUND

To investigate the relationship between fluid-attenuated inversion recovery (FLAIR) vascular hyperintensity (FVH), hemodynamics, and functional outcome in atherosclerotic middle cerebral artery (MCA) stenosis using a computational fluid dynamics (CFD) model based on magnetic resonance angiography (MRA), according to a modified Rankin Scale (mRS) at 3 months.

METHODS

A total of 120 patients with 50-99% atherosclerotic MCA stenosis were included. The training and internal validation groups were composed of 99 participants and 21 participants, respectively. Demographic, imaging data, and functional outcome (mRS at 3 months) were collected. Hemodynamic parameters were obtained from the CFD model. The FVH score was based on the number of territories where FVH is positive, according to the spatial distribution in the Alberta Stroke Program Early Computed Tomography Score (ASPECTS). The prediction models were constructed according to clinical and hemodynamic parameters using multivariate logistic analysis. The DeLong test compared areas under the curves (AUCs) of the models.

RESULTS

The multivariable logistic regression analysis showed that the National Institute of Health Stroke Scale (NIHSS) at admission, hypertension, hyperlipidemia, the ratio of wall shear stress before treatment (WSSR_before), and difference in the ratio of wall shear stress (WSSR) were independently associated with functional outcome (all P<0.05). In the training group before treatment, the AUC of model 1a (only clinical variables) and 2a (clinical variables with addition of WSSR_before) were 0.750 and 0.802. After treatment, the AUC of model 1b (only clinical variables) and 2b (clinical variables with addition of difference in WSSR) were 0.815 and 0.883, respectively. The AUC of models with hemodynamic parameters was significantly higher than the models based on clinical variables only (all P<0.05, DeLong test). In the internal validation group before treatment, the AUC of the model (clinical variables) was 0.782, and that of the model (clinical variables and WSSR_before) was 0.800. After treatment, the AUC of the model (clinical variables) was 0.833, and that of the model (clinical variables and difference in WSSR) was 0.861. There were no significant differences between the good and the poor functional outcome group concerning FVH_before scores before treatment (0.30±0.81 vs. 0.26±0.97; P=0.321) and FVH_after scores after treatment (0.08±0.39 vs. 0.00±0.00; P=0.244).

CONCLUSIONS

Hemodynamics was associated with functional outcomes in patients with ischemic stroke attributed to atherosclerotic MCA stenosis, while FVH was not. Hemodynamic parameters were of great importance in the prediction models.

Smooth muscle cells affect differential nanoparticle accumulation in disturbed blood flow-induced murine atherosclerosis

Atherosclerosis is a lipid-driven chronic inflammatory disease that leads to the formation of plaques in the inner lining of arteries. Plaques form over a range of phenotypes, the most severe of which is vulnerable to rupture and causes most of the clinically significant events. In this study, we evaluated the efficacy of nanoparticles (NPs) to differentiate between two plaque phenotypes based on accumulation kinetics in a mouse model of atherosclerosis. This model uses a perivascular cuff to induce two regions of disturbed wall shear stress (WSS) on the inner lining of the instrumented artery, low (upstream) and multidirectional (downstream), which, in turn, cause the development of an unstable and stable plaque phenotype, respectively. To evaluate the influence of each WSS condition, in addition to the final plaque phenotype, in determining NP uptake, mice were injected with NPs at intermediate and fully developed stages of plaque growth. The kinetics of artery wall uptake were assessed in vivo using dynamic contrast-enhanced magnetic resonance imaging. At the intermediate stage, there was no difference in NP uptake between the two WSS conditions, although both were different from the control arteries. At the fully-developed stage, however, NP uptake was reduced in plaques induced by low WSS, but not multidirectional WSS. Histological evaluation of plaques induced by low WSS revealed a significant inverse correlation between the presence of smooth muscle cells and NP accumulation, particularly at the plaque-lumen interface, which did not exist with other constituents (lipid and collagen) and was not present in plaques induced by multidirectional WSS. These findings demonstrate that NP accumulation can be used to differentiate between unstable and stable murine atherosclerosis, but accumulation kinetics are not directly influenced by the WSS condition. This tool could be used as a diagnostic to evaluate the efficacy of experimental therapeutics for atherosclerosis.

Interactions Between Morphological Plaque Characteristics and Coronary Physiology: From Pathophysiological Basis to Clinical Implications

High-risk coronary plaque refers to a distinct set of plaque characteristics prone to future coronary events. Coronary physiology represents a group of indexes reflective of the local physiological environment and hemodynamic changes in the macrovascular and microvascular system. Although a large body of evidence has supported the clinical relevance of these 2 factors, currently, identifying plaque morphology cannot reliably capture the lesion subset that causes hard events. Also, the guideline-directed approach based on physiological indexes cannot fully predict and prevent clinical events. In parallel, there is accumulating evidence that these 2 aspects of coronary artery disease influence each other with significant clinical implications, despite traditionally being considered to have separate effects on significances, treatments, and outcomes. In this state-of-the-art review, we explore the clinical evidence of pathophysiological interplay of physiological indexes related to local hemodynamics, epicardial stenosis, and microvascular dysfunction with plaque morphological characteristics that provide a better understanding of the nature of coronary events. Furthermore, we examine the emerging data on the complementary role between plaque morphology and coronary physiology in prognostication and how to apply this concept to overcome the limitations of individual assessment alone. Finally, we propose the potential benefit of integrative assessment of coronary anatomy, plaque quantity and quality, and physiological aspects of a target lesion and vessels for personalized risk profiling and optimized treatment strategy.

The Association Between Time-Varying Wall Shear Stress and the Development of Plaque Ulcerations in Carotid Arteries From the Plaque at Risk Study

Background and Purpose: Shear stress (WSS) is involved in the pathophysiology of atherosclerotic disease and might affect plaque ulceration. In this case-control study, we compared carotid plaques that developed a new ulcer during follow-up and plaques that remained silent for their exposure to time-dependent oscillatory shear stress parameters at baseline.

Materials and Methods: Eighteen patients who underwent CTA and MRI of their carotid arteries at baseline and 2 years follow-up were included. These 18 patients consisted of six patients who demonstrated a new ulcer and 12 control patients selected from a larger cohort with similar MRI-based plaque characteristics as the ulcer group. (Oscillatory) WSS parameters [time average WSS, oscillatory shear index (OSI), and relative residence time (RRT)] were calculated using computational fluid dynamics applying the MRI-based geometry of the carotid arteries and compared among plaques (wall thickness>2 mm) with and without ulceration (Mann–Whitney U test) and ulcer-site vs. non-ulcer-site within the plaque (Wilcoxon signed rank test). More detailed analysis on ulcer cases was performed and the predictive value of oscillatory WSS parameters was calculated using linear and logistic mixed-effect regression models.

Results: The ulcer group demonstrated no difference in maximum WSS [9.9 (6.6–18.5) vs. 13.6 (9.7–17.7) Pa, p = 0.349], a lower maximum OSI [0.04 (0.01–0.10) vs. 0.12 (0.06–0.20) p = 0.019] and lower maximum RRT [1.25 (0.78–2.03) Pa⁻¹ vs. 2.93 (2.03–5.28) Pa⁻¹, p = 0.011] compared to controls. The location of the ulcer (ulcer-site) within the plaque was not always at the maximal WSS, but demonstrated higher average WSS, lower average RRT and OSI at the ulcer-site compared to the non-ulcer-sites. High WSS (WSS>4.3 Pa) and low RRT (RRT < 0.25 Pa) were associated with ulceration with an odds ratio of 3.6 [CI 2.1–6.3] and 2.6 [CI 1.54–4.44] respectively, which remained significant after adjustment for wall thickness.

Conclusion: In this explorative study, ulcers were not exclusively located at plaque regions exposed to the highest WSS, OSI, or RRT, but high WSS and low RRT regions had a significantly higher odds to present ulceration within the plaque even after adjustment for wall thickness.

Forward-viewing estimation of 3D blood flow velocity fields by intravascular ultrasound: Influence of the catheter on velocity estimation in stenoses

Coronary artery disease is the most common type of cardiovascular disease, affecting > 18 million adults, and is responsible for > 365 k deaths per year in the U.S. alone. Wall shear stress (WSS) is an emerging indicator of likelihood of plaque rupture in coronary artery disease, however, non-invasive estimation of 3-D blood flow velocity and WSS is challenging due to the requirement for high spatial resolution at deep penetration depths in the presence of significant cardiac motion. Thus we propose minimally-invasive imaging with a catheter-based, 3-D intravascular forward-viewing ultrasound (FV US) transducer and present experiments to quantify the effect of the catheter on flow disturbance in stenotic vessel phantoms with realistic velocities and luminal diameters for both peripheral (6.33 mm) and coronary (4.74 mm) arteries. An external linear array ultrasound transducer was used to quantify 2-D velocity fields in vessel phantoms under various conditions of catheter geometry, luminal diameter, and position of the catheter relative to the stenosis at a frame rate of 5000 frames per second via a particle imaging velocimetry (PIV) approach. While a solid catheter introduced an underestimation of velocity measurement by > 20% relative to the case without a catheter, the hollow catheter introduced < 10% velocity overestimation, indicating that a hollow catheter design allowing internal blood flow reduces hemodynamic disturbance. In addition, for both peripheral and coronary arteries, the hollow catheter introduced < 3% deviation in flow velocity at the minimum luminal area compared to the control case. Finally, an initial comparison was made between velocity measurements acquired using a low frequency, catheter-based, 3-D intravascular FV US transducer and external linear array measurements, with relative error < 12% throughout the region of interest for a flow rate of 150 mL/min. While further system development is required, results suggest intravascular ultrasound characterization of blood flow velocity fields in stenotic vessels could be feasible with appropriate catheter design.

Multiscale Computational Modeling of Vascular Adaptation: A Systems Biology Approach Using Agent-Based Models

The widespread incidence of cardiovascular diseases and associated mortality and morbidity, along with the advent of powerful computational resources, have fostered an extensive research in computational modeling of vascular pathophysiology field and promoted in-silico models as a support for biomedical research. Given the multiscale nature of biological systems, the integration of phenomena at different spatial and temporal scales has emerged to be essential in capturing mechanobiological mechanisms underlying vascular adaptation processes. In this regard, agent-based models have demonstrated to successfully embed the systems biology principles and capture the emergent behavior of cellular systems under different pathophysiological conditions. Furthermore, through their modular structure, agent-based models are suitable to be integrated with continuum-based models within a multiscale framework that can link the molecular pathways to the cell and tissue levels. This can allow improving existing therapies and/or developing new therapeutic strategies. The present review examines the multiscale computational frameworks of vascular adaptation with an emphasis on the integration of agent-based approaches with continuum models to describe vascular pathophysiology in a systems biology perspective. The state-of-the-art highlights the current gaps and limitations in the field, thus shedding light on new areas to be explored that may become the future research focus. The inclusion of molecular intracellular pathways (e.g., genomics or proteomics) within the multiscale agent-based modeling frameworks will certainly provide a great contribution to the promising personalized medicine. Efforts will be also needed to address the challenges encountered for the verification, uncertainty quantification, calibration and validation of these multiscale frameworks.

Risk of myocardial infarction based on endothelial shear stress analysis using coronary angiography

BACKGROUND AND AIMS

Wall shear stress (WSS) has been associated with atherogenesis and plaque progression. The present study assessed the value of WSS analysis derived from conventional coronary angiography to detect lesions culprit for future myocardial infarction (MI).

METHODS AND RESULTS

Three-dimensional quantitative coronary angiography (3DQCA), was used to calculate WSS and pressure drop in 80 patients. WSS descriptors were compared between 80 lesions culprit of future MI and 108 non-culprit lesions (controls). Endothelium-blood flow interaction was assessed by computational fluid dynamics (10.8 ± 1.41 min per vessel). Median time between baseline angiography and MI was 25.9 (21.9-29.8) months. Mean patient age was 70.3 ± 12.7. Clinical presentation was STEMI in 35% and NSTEMI in 65%. Culprit lesions showed higher percent area stenosis (%AS), translesional vFFR difference (ΔvFFR), time-averaged WSS (TAWSS) and topological shear variation index (TSVI) compared to non-culprit lesions (p < 0.05 for all). TSVI was superior to TAWSS in predicting MI (AUC-TSVI = 0.77, 95%CI 0.71-0.84 vs. AUC-TAWSS = 0.61, 95%CI 0.53-0.69, p < 0.001). The addition of TSVI increased predictive and reclassification abilities compared to a model based on %AS and ΔvFFR (NRI = 1.04, p < 0.001, IDI = 0.22, p < 0.001).

CONCLUSIONS

A 3DQCA-based WSS analysis was feasible and can identify lesions culprit for future MI. The combination of area stenoses, pressure gradients and WSS predicted the occurrence of MI. TSVI, a novel WSS descriptor, showed strong predictive capacity to detect lesions prone to cause MI.

Quantifying Patient-Specific in vivo Coronary Plaque Material Properties for Accurate Stress/Strain Calculations: An IVUS-Based Multi-Patient Study

Introduction: Mechanical forces are closely associated with plaque progression and rupture. Precise quantifications of biomechanical conditions using in vivo image-based computational models depend heavily on the accurate estimation of patient-specific plaque mechanical properties. Currently, mechanical experiments are commonly performed on ex vivo cardiovascular tissues to determine plaque material properties. Patient-specific in vivo coronary material properties are scarce in the existing literature.

Methods:In vivo Cine intravascular ultrasound and virtual histology intravascular ultrasound (IVUS) slices were acquired at 20 plaque sites from 13 patients. A three-dimensional thin-slice structure-only model was constructed for each slice to obtain patient-specific in vivo material parameter values following an iterative scheme. Effective Young's modulus (YM) was calculated to indicate plaque stiffness for easy comparison purposes. IVUS-based 3D thin-slice models using in vivo and ex vivo material properties were constructed to investigate their impacts on plaque wall stress/strain (PWS/PWSn) calculations.

Results: The average YM values in the axial and circumferential directions for the 20 plaque slices were 599.5 and 1,042.8 kPa, respectively, 36.1% lower than those from published ex vivo data. The YM values in the circumferential direction of the softest and stiffest plaques were 103.4 and 2,317.3 kPa, respectively. The relative difference of mean PWSn on lumen using the in vivo and ex vivo material properties could be as high as 431%, while the relative difference of mean PWS was much lower, about 3.07% on average.

Conclusion: There is a large inter-patient and intra-patient variability in the in vivo plaque material properties. In vivo material properties have a great impact on plaque stress/strain calculations. In vivo plaque material properties have a greater impact on strain calculations. Large-scale-patient studies are needed to further verify our findings.

The definition of low wall shear stress and its effect on plaque progression estimation in human coronary arteries

Wall shear stress (WSS), the frictional force of the blood on the vessel wall, plays a crucial role in atherosclerotic plaque development. Low WSS has been associated with plaque growth, however previous research used different approaches to define low WSS to investigate its effect on plaque progression. In this study, we used four methodologies to allocate low, mid and high WSS in one dataset of human coronary arteries and investigated the predictive power of low WSS for plaque progression. Coronary reconstructions were based on multimodality imaging, using intravascular ultrasound and CT-imaging. Vessel-specific flow was measured using Doppler wire and computational fluid dynamics was performed to calculate WSS. The absolute WSS range varied greatly between the coronary arteries. On the population level, the established pattern of most plaque progression at low WSS was apparent in all methodologies defining the WSS categories. However, for the individual patient, when using measured flow to determine WSS, the absolute WSS values range so widely, that the use of absolute thresholds to determine low WSS was not appropriate to identify regions at high risk for plaque progression.

Numerical modelling of blood rheology and platelet activation through a stenosed left coronary artery bifurcation

Coronary bifurcations are prone to atherosclerotic plaque growth, experiencing regions of reduced wall shear stress (WSS) and increased platelet adhesion. This study compares effects across different rheological approaches on hemodynamics, combined with a shear stress exposure history model of platelets within a stenosed porcine bifurcation. Simulations used both single/multiphase blood models to determine which approach best predicts phenomena associated with atherosclerosis and atherothrombosis. A novel Lagrangian platelet tracking model was used to evaluate residence time and shear history of platelets indicating likely regions of thrombus formation. Results show a decrease in area of regions with pathologically low time-averaged WSS with the use of multiphase models, particularly in a stenotic bifurcation. Significant non-Newtonian effects were observed due to low-shear and varying hematocrit levels found on the outer walls of the bifurcation and distal to the stenosis. Platelet residence time increased 11% in the stenosed artery, with exposure times to low-shear sufficient for red blood cell aggregation (>1.5 s). increasing the risk of thrombosis. This shows stenotic artery hemodynamics are inherently non-Newtonian and multiphase, with variations in hematocrit (0.163-0.617) and elevated vorticity distal to stenosis (+15%) impairing the function of the endothelium via reduced time-averaged WSS regions, rheological properties and platelet activation/adhesion.

Consistency in Geometry Among Coronary Atherosclerotic Plaques Extracted From Computed Tomography Angiography

Background: The three-dimensional (3D) geometry of coronary atherosclerotic plaques is associated with plaque growth and the occurrence of coronary artery disease. However, there is a lack of studies on the 3D geometric properties of coronary plaques. We aim to investigate if coronary plaques of different sizes are consistent in geometric properties.

Methods: Nineteen cases with symptomatic stenosis caused by atherosclerotic plaques in the left coronary artery were included. Based on attenuation values on computed tomography angiography images, coronary atherosclerotic plaques and calcifications were identified, 3D reconstructed, and manually revised. Multidimensional geometric parameters were measured on the 3D models of plaques and calcifications. Linear and non-linear (i.e., power function) fittings were used to investigate the relationship between multidimensional geometric parameters (length, surface area, volume, etc.). Pearson correlation coefficient (r), R-squared, and p-values were used to evaluate the significance of the relationship. The analysis was performed based on cases and plaques, respectively. Significant linear relationship was defined as R-squared > 0.25 and p < 0.05.

Results: In total, 49 atherosclerotic plaques and 56 calcifications were extracted. In the case-based analysis, significant linear relationships were found between number of plaques and number of calcifications (r = 0.650, p = 0.003) as well as total volume of plaques (r = 0.538, p = 0.018), between number of calcifications and total volume of plaques (r = 0.703, p = 0.001) as well as total volume of calcification (r = 0.646, p = 0.003), and between the total volumes of plaques and calcifications (r = 0.872, p < 0.001). In plaque-based analysis, the power function showed higher R-squared values than the linear function in fitting the relationships of multidimensional geometric parameters. Two presumptions of plaque geometry in different growth stages were proposed with simplified geometric models developed. In the proposed models, the exponents in the power functions of geometric parameters were in accordance with the fitted values.

Conclusion: In patients with coronary artery disease, coronary plaques and calcifications are positively related in number and volume. Different coronary plaques are consistent in the relationship between geometry parameters in different dimensions.

Coronary Computed Tomography Angiography Assessment of High-Risk Plaques in Predicting Acute Coronary Syndrome

Coronary computed tomography angiography (CCTA) is a comprehensive, non-invasive and cost-effective imaging assessment approach, which can provide the ability to identify the characteristics and morphology of high-risk atherosclerotic plaques associated with acute coronary syndrome (ACS). The development of CCTA and latest advances in emerging technologies, such as computational fluid dynamics (CFD), have made it possible not only to identify the morphological characteristics of high-risk plaques non-invasively, but also to assess the hemodynamic parameters, the environment surrounding coronaries and so on, which may help to predict the risk of ACS. In this review, we present how CCTA was used to characterize the composition and morphology of high-risk plaques prone to ACS and the current role of CCTA, including emerging CCTA technologies, advanced analysis, and characterization techniques in prognosticating the occurrence of ACS.

Association Among Local Hemodynamic Parameters Derived From CT Angiography and Their Comparable Implications in Development of Acute Coronary Syndrome

Background: Association among local hemodynamic parameters and their implications in development of acute coronary syndrome (ACS) have not been fully investigated.

Methods: A total of 216 lesions in ACS patients undergoing coronary CT angiography (CCTA) before 1–24 months from ACS event were analyzed. High-risk plaque on CCTA was defined as a plaque with ≥2 of low-attenuation plaque, positive remodeling, spotty calcification, and napkin-ring sign. With the use of computational fluid dynamics analysis, fractional flow reserve (FFR) derived from CCTA (FFR_CT) and local hemodynamic parameters including wall shear stress (WSS), axial plaque stress (APS), pressure gradient (PG) across the lesion, and delta FFR_CT across the lesion (ΔFFR_CT) were obtained. The association among local hemodynamics and their discrimination ability for culprit lesions from non-culprit lesions were compared.

Results: A total of 66 culprit lesions for later ACS and 150 non-culprit lesions were identified. WSS, APS, PG, and ΔFFR_CT were strongly correlated with each other (all p < 0.001). This association was persistent in all lesion subtypes according to a vessel, lesion location, anatomical severity, high-risk plaque, or FFR_CT ≤ 0.80. In discrimination of culprit lesions causing ACS from non-culprit lesions, WSS, PG, APS, and ΔFFR_CT were independent predictors after adjustment for lesion characteristics, high-risk plaque, and FFR_CT ≤ 0.80; and all local hemodynamic parameters significantly improved the predictive value for culprit lesions of high-risk plaque and FFR_CT ≤ 0.80 (all p < 0.05). The risk prediction model for culprit lesions with FFR_CT ≤ 0.80, high-risk plaque, and ΔFFR_CT had a similar or superior discrimination ability to that with FFR_CT ≤ 0.80, high-risk plaque, and WSS, APS, or PG; and the addition of WSS, APS, or PG into ΔFFR_CT did not improve the model performance.

Conclusions: Local hemodynamic indices were significantly intercorrelated, and all indices similarly provided additive and independent predictive values for ACS risk over high-risk plaque and impaired FFR_CT.

Comparison of Newtonian and Non-newtonian Fluid Models in Blood Flow Simulation in Patients With Intracranial Arterial Stenosis

Background

Newtonian fluid model has been commonly applied in simulating cerebral blood flow in intracranial atherosclerotic stenosis (ICAS) cases using computational fluid dynamics (CFD) modeling, while blood is a shear-thinning non-Newtonian fluid. We aimed to investigate the differences of cerebral hemodynamic metrics quantified in CFD models built with Newtonian and non-Newtonian fluid assumptions, in patients with ICAS.

Methods

We built a virtual artery model with an eccentric 75% stenosis and performed static CFD simulation. We also constructed CFD models in three patients with ICAS of different severities in the luminal stenosis. We performed static simulations on these models with Newtonian and two non-Newtonian (Casson and Carreau-Yasuda) fluid models. We also performed transient simulations on another patient-specific model. We measured translesional pressure ratio (PR) and wall shear stress (WSS) values in all CFD models, to reflect the changes in pressure and WSS across a stenotic lesion. In all the simulations, we compared the PR and WSS values in CFD models derived with Newtonian, Casson, and Carreau-Yasuda fluid assumptions.

Results

In all the static and transient simulations, the Newtonian/non-Newtonian difference on PR value was negligible. As to WSS, in static models (virtual and patient-specific), the rheological difference was not obvious in areas with high WSS, but observable in low WSS areas. In the transient model, the rheological difference of WSS areas with low WSS was enhanced, especially during diastolic period.

Conclusion

Newtonian fluid model could be applicable for PR calculation, but caution needs to be taken when using the Newtonian assumption in simulating WSS especially in severe ICAS cases.

Optical Coherence Tomography of Coronary Plaque Progression and Destabilization: JACC Focus Seminar Part 3/3

The development of optical coherence tomography (OCT) has revolutionized our understanding of coronary artery disease. In vivo OCT research has paralleled with advances in computational fluid dynamics, providing additional insights in the various hemodynamic factors influencing plaque growth and stability. Recent OCT studies introduced a new concept of plaque healing in relation to clinical presentation. In addition to known mechanisms of acute coronary syndromes such as plaque rupture and plaque erosion, a new classification of calcified plaque was recently reported. This review will focus on important new insights that OCT has provided in recent years into coronary plaque development, progression, and destabilization, with a focus on the role of local hemodynamics and endothelial shear stress, the layered plaque (signature of previous subclinical plaque destabilization and healing), and the calcified culprit plaque.

Hemodynamic Responses in Carotid Bifurcation Induced by Enhanced External Counterpulsation Stimulation in Healthy Controls and Patients With Neurological Disorders

Enhanced external counterpulsation is a Food and Drug Administration–approved, non-invasive, assisted circulation therapy for ischemic cardiovascular and cerebrovascular diseases. Previous studies have confirmed that EECP stimulation induces largely different cerebral hemodynamic responses in patients with ischemic stroke and healthy controls. However, the underlying mechanisms remain uncertain. We hypothesize that different blood redistributions at the carotid bifurcation may play a key role. Ten subjects were enrolled in this study, namely, five patients with neurological disorders and five young healthy volunteers as controls. Magnetic resonance angiography (MRA) was performed on the carotid artery. All the subjects received a single session of EECP treatment, with external cuff pressures ranging from 20 to 40 kPa. Vascular ultrasound measurements were taken in the common carotid artery (CCA), external carotid artery (ECA) and internal carotid artery (ICA). Three-dimensional patient-specific numerical models were developed to calculate the WSS-derived hemodynamic factors. The results indicated that EECP increased CCA and ECA blood flow in both groups. The ICA blood flow in the patient group exhibited a mean increase of 6.67% during EECP treatment compared with the pre-EECP condition; a mean decrease of 9.2% was observed in the healthy control group. EECP increased the averaged wall shear stress (AWSS) throughout the carotid bifurcation in the patient group; the ICA AWSS of the healthy group decreased during EECP. In both groups, the oscillatory shear index (OSI) in the ICA increased proportionally with external cuff pressure. In addition, the relative resident time (RRT) was constant or slightly decreased in the CCA and ECA in both groups but increased in the ICA. We suggest that the benefits of EECP to patients with neurological disorders may partly result from blood flow promotion in the ICA and increase in WSS at the carotid bifurcation. In the healthy subjects, the ICA blood flow remained constant during EECP, although the CCA blood flow increased significantly. A relatively low external cuff pressure (20 kPa) is recommended as the optimal treatment pressure for better hemodynamic effects. This study may play an important role in the translation of physiological benefits of EECP treatment in populations with or without neurological disorders.

Effect of Subject-Specific, Spatially Reduced, and Idealized Boundary Conditions on the Predicted Hemodynamic Environment in the Murine Aorta

Mouse models of atherosclerosis have become effective resources to study atherogenesis, including the relationship between hemodynamics and lesion development. Computational methods aid the prediction of the in vivo hemodynamic environment in the mouse vasculature, but careful selection of inflow and outflow boundary conditions (BCs) is warranted to promote model accuracy. Herein, we investigated the impact of animal-specific versus reduced/idealized flow boundary conditions on predicted blood flow patterns in the mouse thoracic aorta. Blood velocities were measured in the aortic root, arch branch vessel, and descending aorta in ApoE^-/- mice using phase-contrast MRI. Computational geometries were derived from micro-CT imaging and combinations of high-fidelity or reduced/idealized MR-derived BCs were applied to predict the bulk flow field and hemodynamic metrics (e.g., wall shear stress, WSS; cross-flow index, CFI). Results demonstrate that pressure-free outlet BCs significantly overestimate outlet flow rates as compared to measured values. When compared to models that incorporate 3-component inlet velocity data [[Formula: see text]] and time-varying outlet mass flow splits [[Formula: see text]] (i.e., high-fidelity model), neglecting in-plane inlet velocity components (i.e., [Formula: see text])) leads to errors in WSS and CFI values ranging from 10 to 30% across the model domain whereas the application of a steady outlet mass flow splits results in negligible differences in these hemodynamics metrics. This investigation highlights that 3-component inlet velocity data and at least steady mass flow splits are required for accurate predictions of flow patterns in the mouse thoracic aorta.

A Review on the Effect of Temporal Geometric Variations of the Coronary Arteries on the Wall Shear Stress and Pressure Drop

Temporal variations of the coronary arteries during a cardiac cycle are defined as the superposition of the changes in the position, curvature, and torsion of the coronary artery axis markers and the variations in the lumen cross-sectional shape due to the distensible wall motion induced by the pulse pressure and contraction of the myocardium in a cardiac cycle. This review discusses whether modeling of the temporal variations of the coronary arteries is needed for the investigation of hemodynamics specifically in time-critical applications such as a clinical environment. The numerical modelings in the literature that model or disregard the temporal variations of the coronary arteries on the hemodynamic parameters are discussed. The results in the literature show that neglecting the effects of temporal geometric variations is expected to result in about 5% deviation of the time-averaged pressure drop and wall shear stress values and also about 20% deviation of the temporal variations of hemodynamic parameters, such as time-dependent wall shear stress and oscillatory shear index. This review study can be considered as a guide for future studies to outline the conditions in which temporal variations of the coronary arteries can be neglected while providing a reliable estimation of hemodynamic parameters.

Using Optical Coherence Tomography and Intravascular Ultrasound Imaging to Quantify Coronary Plaque Cap Stress/Strain and Progression: A Follow-Up Study Using 3D Thin-Layer Models

Accurate plaque cap thickness quantification and cap stress/strain calculations are of fundamental importance for vulnerable plaque research. To overcome uncertainties due to intravascular ultrasound (IVUS) resolution limitation, IVUS and optical coherence tomography (OCT) coronary plaque image data were combined together to obtain accurate and reliable cap thickness data, stress/strain calculations, and reliable plaque progression predictions. IVUS, OCT, and angiography baseline and follow-up data were collected from nine patients (mean age: 69; m: 5) at Cardiovascular Research Foundation with informed consent obtained. IVUS and OCT slices were coregistered and merged to form IVUS + OCT (IO) slices. A total of 114 matched slices (IVUS and OCT, baseline and follow-up) were obtained, and 3D thin-layer models were constructed to obtain stress and strain values. A generalized linear mixed model (GLMM) and least squares support vector machine (LSSVM) method were used to predict cap thickness change using nine morphological and mechanical risk factors. Prediction accuracies by all combinations (511) of those predictors with both IVUS and IO data were compared to identify optimal predictor(s) with their best accuracies. For the nine patients, the average of minimum cap thickness from IVUS was 0.17 mm, which was 26.08% lower than that from IO data (average = 0.23 mm). Patient variations of the individual errors ranged from ‒58.11 to 20.37%. For maximum cap stress between IO and IVUS, patient variations of the individual errors ranged from ‒30.40 to 46.17%. Patient variations of the individual errors of maximum cap strain values ranged from ‒19.90 to 17.65%. For the GLMM method, the optimal combination predictor using IO data had AUC (area under the ROC curve) = 0.926 and highest accuracy = 90.8%, vs. AUC = 0.783 and accuracy = 74.6% using IVUS data. For the LSSVM method, the best combination predictor using IO data had AUC = 0.838 and accuracy = 75.7%, vs. AUC = 0.780 and accuracy = 69.6% using IVUS data. This preliminary study demonstrated improved plaque cap progression prediction accuracy using accurate cap thickness data from IO slices and the differences in cap thickness, stress/strain values, and prediction results between IVUS and IO data. Large-scale studies are needed to verify our findings.

Abnormal shear stress and residence time are associated with proximal coronary atheroma in the presence of myocardial bridging

BACKGROUND

Atheromatous plaques tend to form in the coronary segments proximal to a myocardial bridge (MB), but the mechanism of this occurrence remains unclear. This study evaluates the relationship between blood flow perturbations and plaque formation in patients with an MB.

METHODS AND RESULTS

A total of 92 patients with an MB in the mid left anterior descending artery (LAD) and 20 patients without an MB were included. Coronary angiography, intravascular ultrasound, and coronary physiology measurements were performed. A moving-boundary computational fluid dynamics algorithm was used to derive wall shear stress (WSS) and peak residence time (PRT). Patients with an MB had lower WSS (0.46 ± 0.21 vs. 0.96 ± 0.33 Pa, p < 0.001) and higher maximal plaque burden (33.6 ± 15.0 vs. 14.2 ± 5.8%, p < 0.001) within the proximal LAD compared to those without. Plaque burden in the proximal LAD correlated significantly with proximal WSS (r = -0.51, p < 0.001) and PRT (r = 0.60, p < 0.001). In patients with an MB, the site of maximal plaque burden occurred 23.4 ± 13.3 mm proximal to the entrance of the MB, corresponding to the site of PRT.

CONCLUSIONS

Regions of low WSS and high PRT occur in arterial segments proximal to an MB, and this is associated with the degree and location of coronary atheroma formation.

Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

High shear stress enhances endothelial permeability in the presence of the risk haplotype at 9p21.3

Single nucleotide polymorphisms (SNPs) are exceedingly common in non-coding loci, and while they are significantly associated with a myriad of diseases, their specific impact on cellular dysfunction remains unclear. Here, we show that when exposed to external stressors, the presence of risk SNPs in the 9p21.3 coronary artery disease (CAD) risk locus increases endothelial monolayer and microvessel dysfunction. Endothelial cells (ECs) derived from induced pluripotent stem cells of patients carrying the risk haplotype (R/R WT) differentiated similarly to their non-risk and isogenic knockout (R/R KO) counterparts. Monolayers exhibited greater permeability and reactive oxygen species signaling when the risk haplotype was present. Addition of the inflammatory cytokine TNFα further enhanced EC monolayer permeability but independent of risk haplotype; TNFα also did not substantially alter haplotype transcriptomes. Conversely, when wall shear stress was applied to ECs in a microfluidic vessel, R/R WT vessels were more permeable at lower shear stresses than R/R KO vessels. Transcriptomes of sheared cells clustered more by risk haplotype than by patient or clone, resulting in significant differential regulation of EC adhesion and extracellular matrix genes vs static conditions. A subset of previously identified CAD risk genes invert expression patterns in the presence of high shear concomitant with altered cell adhesion genes, vessel permeability, and endothelial erosion in the presence of the risk haplotype, suggesting that shear stress could be a regulator of non-coding loci with a key impact on CAD.

Hemodynamic Analysis of a Three-Point Suture During Tapering Technique for Microanastomosis Using Computational Fluid Dynamics

The tapering technique is one of the useful methods of anastomosing 2 vessels with large discrepancies during microanastomoses. When the tapering technique is used, a three-point suture is always present. The authors analyzed the most appropriate suture technique for this using computational fluid dynamics. This aspect has not previously been addressed. Three different suture techniques were simulated: (1)

Three single-knot sutures (Type I);

(2)

Two single-knot sutures forming an X-shape (Type II); and

(3)

A single continuous ligature through the vascular wall (Type III).

Vascular models of these 3 types were created. The streamline, wall shear stress, and oscillatory shear index at the anastomosis site were measured using a previously prepared venous model. Streamline disruption was most severe for Type II. In all 3 types, the highest wall shear stress was recorded at the suture peak protruding into the vessel. The maximum oscillatory shear index was highest in Type II, and lowest in Type III. The present results suggest that Type III is the best three-point suturing method for the tapering technique.

Relationship of Endothelial Shear Stress with Plaque Features with Coronary CT Angiography and Vasodilating Capability with PET

Background Advances in three-dimensional reconstruction techniques and computational fluid dynamics of coronary CT angiography (CCTA) data sets make feasible evaluation of endothelial shear stress (ESS) in the vessel wall. Purpose To investigate the relationship between CCTA-derived computational fluid dynamics metrics, anatomic and morphologic characteristics of coronary lesions, and their comparative performance in predicting impaired coronary vasodilating capability assessed by using PET myocardial perfusion imaging (MPI). Materials and Methods In this retrospective study, conducted between October 2019 and September 2020, coronary vessels in patients with stable chest pain and with intermediate probability of coronary artery disease who underwent both CCTA and PET MPI with oxygen 15-labeled water or nitrogen 13 ammonia and quantification of myocardial blood flow were analyzed. CCTA images were used in assessing stenosis severity, lesion-specific total plaque volume (PV), noncalcified PV, calcified PV, and plaque phenotype. PET MPI was used in assessing significant coronary stenosis. The predictive performance of the CCTA-derived parameters was evaluated by using area under the receiver operating characteristic curve (AUC) analysis. Results There were 92 coronary vessels evaluated in 53 patients (mean age, 65 years ± 7; 31 men). ESS was higher in lesions with greater than 50% stenosis versus those without significant stenosis (mean, 15.1 Pa ± 30 vs 4.6 Pa ± 4 vs 3.3 Pa ± 3; P = .004). ESS was higher in functionally significant versus nonsignificant lesions (median, 7 Pa [interquartile range, 5-23 Pa] vs 2.6 Pa [interquartile range, 1.8-5 Pa], respectively; P ≤ .001). Adding ESS to stenosis severity improved prediction (change in AUC, 0.10; 95% CI: 0.04, 0.17; P = .002) for functionally significant lesions. Conclusion The combination of endothelial shear stress with coronary CT angiography (CCTA) stenosis severity improved prediction of an abnormal PET myocardial perfusion imaging result versus CCTA stenosis severity alone. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Kusmirek and Wieben in this issue.

Patient-specific fluid-structure interaction model of bile flow: comparison between 1-way and 2-way algorithms

Gallbladder disease is one of the most spread pathologies in the world. Despite the number of operations dealing with biliary surgery increases, the number of postoperative complications is also high. The aim of this study is to show the influence of the biliary system pathology on bile flow character and to numerically assess the effect of surgical operation (cholecystectomy) on the fluid dynamics in the extrahepatic biliary tree, and also to reveal the difference between 1-way and 2-way FSI algorithms on the results. Moreover, the bile viscosity and biliary tree geometry influence on the choledynamics were evaluated. Bile velocity, pressure, wall shear stress (WSS), displacements and von Mises stress distributions in the extrahepatic biliary tree are presented, and comparison is made between a healthy and a lithogenic bile. The patient-specific biliary tree model is created using magnetic resonance imaging (MRI) and imported in a commercial finite element analysis software. It is found that in the case of lithogenic bile, velocities have lower magnitude while pressures are higher. Furthermore, stress analysis of the bile ducts shows that the WSS distribution is found mostly prevailing in the common hepatic duct and common bile duct areas. It is shown that when it is necessary to evaluate the bile flow dynamics in urgent medical situations, 1-way analysis is acceptable. Nevertheless, 2-way FSI provides more accurate data, if necessary to evaluate the stress-strain state of bile ducts. The proposed model can be applied to medical practice to reduce the number of post-operative complications.

Choroidal changes in lens-induced myopia in guinea pigs

INTRODUCTION

This study aimed to determine the role of the choroid in lens-induced myopia (LIM) in guinea pigs.

METHODS

Guinea pigs were randomly divided into two groups: a normal control (NC) group and a LIM group. Refraction and axial length (AL) were measured by streak retinoscopy and A-scan ultrasonography. The choroidal thickness (ChT), vessel density of the choriocapillaris (VDCC) and vessel density of the choroidal layer (VDCL) were assessed by Spectral-domain Optical Coherence Tomography Angiography (SD-OCT). In addition, the choroidal expression of nitric oxide synthase (NOS) enzymes at the mRNA and protein levels was analyzed by real-time fluorescence quantitative PCR, enzyme-linked immunosorbent assay (ELISA) and immunohistochemistry.

RESULTS

In the LIM group, refraction and AL were increased significantly compared with those in the NC group at 2 weeks (refraction: LIM vs. NC, -4.23 ± 0.43 D vs. 2.20 ± 0.48 D; AL: LIM vs. NC, 8.36 ± 0.05 mm vs. 8.22 ± 0.03 mm) and 4 weeks (refraction: LIM vs. NC, -5.88 ± 0.49 D vs. 1.63 ± 0.41 D; AL: 8.57 ± 0.06 mm vs. 8.40 ± 0.04 mm). The ChT and VDCC were decreased significantly compared with those in the NC group at 2 weeks (ChT: LIM vs. NC, 60.92 ± 8.15 μm vs. 79.11 ± 7.47 μm; VDCC: LIM vs. NC, 23.43 ± 3.85% vs. 28.74 ± 4.11%) and 4 weeks (ChT: LIM vs. NC, 48.43 ± 6.85 μm vs. 76.38 ± 7.84 μm; VDCC: LIM vs. NC, 21.29 ± 2.17% vs. 27.64 ± 2.91%). The VDCL was also decreased compared with that in the NC group at 2 weeks and 4 weeks (NC vs. LIM, 24.87 ± 5.16% vs. 22.45 ± 3.26%; 23.37 ± 5.85% vs. 21.39 ± 2.62%; all P > 0.05). Moreover, the ChT was positively correlated with the VDCC and VDCL. The mRNA and protein expression of NOS enzymes (eNOS and nNOS) was increased.

CONCLUSIONS

During the development of myopia, the ChT, VDCC and VDCL were decreased, while NOS expression in the choroid was increased. The expression of NOS was negatively correlated with the ChT, VDCC and VDCL. NO may play an important role in regulating the choroid during myopia development.

Functional assessment of coronary plaques using CT based hemodynamic simulations: Current status, technical principles and clinical value

In recent years, coronary computed tomography angiography (CCTA) has emerged as an accurate and safe non-invasive imaging modality in terms of detecting and excluding coronary artery disease (CAD). In the latest European Society of Cardiology Guidelines CCTA received Class I recommendation for the evaluation of patients with stable chest pain with low to intermediate clinical likelihood of CAD. Despite its high negative predictive value, the diagnostic performance of CCTA is limited by the relatively low specificity, especially in patients with heavily calcified lesions. The discrepancy between the degree of stenosis and ischemia is well established based on both invasive and non-invasive tests. The rapid evolution of computational flow dynamics has allowed the simulation of CCTA derived fractional flow reserve (FFR-CT), which improves specificity by combining anatomic and functional information regarding coronary atherosclerosis. FFR-CT has been extensively validated against invasively measured FFR as the reference standard. Due to recent technological advancements FFR-CT values can also be calculated locally, without offsite processing. Wall shear stress (WSS) and axial plaque stress (APS) are additional key hemodynamic elements of atherosclerotic plaque characteristics, which can also be measured using CCTA images. Current evidence suggests that WSS and APS are important hemodynamic features of adverse coronary plaques. CCTA based hemodynamic calculations could therefore improve prognostication and the management of patients with stable CAD.

A multi-objective optimization of stent geometries

Stents are scaffolding cardiovascular implants used to restore blood flow in narrowed arteries. However, the presence of the stent alters local blood flow and shear stresses on the surrounding arterial wall, which can cause adverse tissue responses and increase the risk of adverse outcomes. There is a need for optimization of stent designs for hemodynamic performance. We used multi-objective optimization to identify ideal combinations of design variables by assessing potential trade-offs based on common hemodynamic indices associated with clinical risk and mechanical performance of the stents. We studied seven design variables including strut cross-section, strut dimension, strut angle, cell alignment, cell height, connector type and connector arrangement. Optimization objectives were the percentage of vessel area exposed to adversely low time averaged WSS (TAWSS) and adversely high Wall Shear Stress (WSS) assessed using computational fluid dynamics modeling, as well as radial stiffness of the stent using FEA simulation. Two multi-objective optimization algorithms were used and compared to iteratively predict ideal designs. Out of 50 designs, three best designs with respect to each of the three objectives, and two designs in regard to overall performance were identified.

Ointment therapy and prevention of cannulation-induced superficial phlebitis

The primary aim of this prospective randomized trial is to determine the effect of clobetasol ointment with nitroglycerin ointment on the prevention of phlebitis caused by cannulation. The target sample is 144 patients admitted to surgical ward, randomized to 3 groups. the data collection tools included demographic information, information about intravenous treatment and phlebitis measurement scale. In the intervention groups, after venipuncture, 1.5 cm of ointments (clobetasol and nitroglycerin) was applied n three time periods of 24, 48 and 72 hours from the time of venipuncture. In the first 24 hours after cannulation, there was no significant difference between the three groups. But at 48 and 72 hours after placement, the difference between intervention and control groups was significant (P<0.0001). It is recommended to use clobetasol ointment and nitroglycerin ointment to prevent the occurrence of phlebitis in patients who need long-term use of cannula (more than 48 hours).

Assessment with clinical data of a coupled bio-hemodynamics numerical model to predict leukocyte adhesion in coronary arteries

Numerical simulations of coupled hemodynamics and leukocyte transport and adhesion inside coronary arteries have been performed. Realistic artery geometries have been obtained for a set of four patients from intravascular ultrasound and angiography images. The numerical model computes unsteady three-dimensional blood hemodynamics and leukocyte concentration in the blood. Wall-shear stress dependent leukocyte adhesion is also computed through agent-based modeling rules, fully coupled to the hemodynamics and leukocyte transport. Numerical results have a good correlation with clinical data. Regions where high adhesion is predicted by the simulations coincide to a good approximation with artery segments presenting plaque increase, as documented by clinical data from baseline and six-month follow-up exam of the same artery. In addition, it is observed that the artery geometry and, in particular, the tortuosity of the centerline are a primary factor in determining the spatial distribution of wall-shear stress, and of the resulting leukocyte adhesion patterns. Although further work is required to overcome the limitations of the present model and ultimately quantify plaque growth in the simulations, these results are encouraging towards establishing a predictive methodology for atherosclerosis progress.

Coronary Computer Tomography Angiography in 2021-Acquisition Protocols, Tips and Tricks and Heading beyond the Possible

Recent technological advances, together with an increasing body of evidence from randomized trials, have placed coronary computer tomography angiography (CCTA) in the center of the diagnostic workup of patients with coronary artery disease. The method was proven reliable in the diagnosis of relevant coronary artery stenosis. Furthermore, it can identify different stages of the atherosclerotic process, including early atherosclerotic changes of the coronary vessel wall, a quality not met by other non-invasive tests. In addition, newer computational software can measure the hemodynamic relevance (fractional flow reserve) of a certain stenosis. In addition, if required, information related to cardiac and valvular function can be provided with specific protocols. Importantly, recent trials have highlighted the prognostic relevance of CCTA in patients with coronary artery disease, which helped establishing CCTA as the first-line method for the diagnostic work-up of such patients in current guidelines. All this can be gathered in one relatively fast examination with minimal discomfort for the patient and, with newer machines, with very low radiation exposure. Herein, we provide an overview of the current technical aspects, indications, pitfalls, and new horizons with CCTA, providing examples from our own clinical practice.

Investigation of Wall Shear Stress in Cardiovascular Research and in Clinical Practice-From Bench to Bedside

In the 1900s, researchers established animal models experimentally to induce atherosclerosis by feeding them with a cholesterol-rich diet. It is now accepted that high circulating cholesterol is one of the main causes of atherosclerosis; however, plaque localization cannot be explained solely by hyperlipidemia. A tremendous amount of studies has demonstrated that hemodynamic forces modify endothelial athero-susceptibility phenotypes. Endothelial cells possess mechanosensors on the apical surface to detect a blood stream-induced force on the vessel wall, known as "wall shear stress (WSS)", and induce cellular and molecular responses. Investigations to elucidate the mechanisms of this process are on-going: on the one hand, hemodynamics in complex vessel systems have been described in detail, owing to the recent progress in imaging and computational techniques. On the other hand, investigations using unique in vitro chamber systems with various flow applications have enhanced the understanding of WSS-induced changes in endothelial cell function and the involvement of the glycocalyx, the apical surface layer of endothelial cells, in this process. In the clinical setting, attempts have been made to measure WSS and/or glycocalyx degradation non-invasively, for the purpose of their diagnostic utilization. An increasing body of evidence shows that WSS, as well as serum glycocalyx components, can serve as a predicting factor for atherosclerosis development and, most importantly, for the rupture of plaques in patients with high risk of coronary heart disease.

Differential immunological signature at the culprit site distinguishes acute coronary syndrome with intact from acute coronary syndrome with ruptured fibrous cap: results from the prospective translational OPTICO-ACS study

AIMS

Acute coronary syndromes with intact fibrous cap (IFC-ACS), i.e. caused by coronary plaque erosion, account for approximately one-third of ACS. However, the underlying pathophysiological mechanisms as compared with ACS caused by plaque rupture (RFC-ACS) remain largely undefined. The prospective translational OPTICO-ACS study programme investigates for the first time the microenvironment of ACS-causing culprit lesions (CL) with intact fibrous cap by molecular high-resolution intracoronary imaging and simultaneous local immunological phenotyping.

METHODS AND RESULTS

The CL of 170 consecutive ACS patients were investigated by optical coherence tomography (OCT) and simultaneous immunophenotyping by flow cytometric analysis as well as by effector molecule concentration measurements across the culprit lesion gradient (ratio local/systemic levels). Within the study cohort, IFC caused 24.6% of ACS while RFC-ACS caused 75.4% as determined and validated by two independent OCT core laboratories. The IFC-CL were characterized by lower lipid content, less calcification, a thicker overlying fibrous cap, and largely localized near a coronary bifurcation as compared with RFC-CL. The microenvironment of IFC-ACS lesions demonstrated selective enrichment in both CD4+ and CD8+ T-lymphocytes (+8.1% and +11.2%, respectively, both P < 0.05) as compared with RFC-ACS lesions. T-cell-associated extracellular circulating microvesicles (MV) were more pronounced in IFC-ACS lesions and a significantly higher amount of CD8+ T-lymphocytes was detectable in thrombi aspirated from IFC-culprit sites. Furthermore, IFC-ACS lesions showed increased levels of the T-cell effector molecules granzyme A (+22.4%), perforin (+58.8%), and granulysin (+75.4%) as compared with RFC plaques (P < 0.005). Endothelial cells subjected to culture in disturbed laminar flow conditions, i.e. to simulate coronary flow near a bifurcation, demonstrated an enhanced adhesion of CD8+T cells. Finally, both CD8+T cells and their cytotoxic effector molecules caused endothelial cell death, a key potential pathophysiological mechanism in IFC-ACS.

CONCLUSIONS

The OPTICO-ACS study emphasizes a novel mechanism in the pathogenesis of IFC-ACS, favouring participation of the adaptive immune system, particularly CD4+ and CD8+ T-cells and their effector molecules. The different immune signatures identified in this study advance the understanding of coronary plaque progression and may provide a basis for future development of personalized therapeutic approaches to ACS with IFC.

TRIAL REGISTRATION

The study was registered at clinicalTrials.gov (NCT03129503).

A Numerical Model for Simulating the Hemodynamic Effects of Enhanced External Counterpulsation on Coronary Arteries

Traditional enhanced external counterpulsation (EECP) used for the clinical treatment of patients with coronary heart disease only assesses diastolic/systolic blood pressure (Q = D/S > 1.2). However, improvement of the hemodynamic environment surrounding vascular endothelial cells of coronary arteries after long-term application of EECP is the basis of the treatment. Currently, the quantitative hemodynamic mechanism is not well understood. In this study, a standard 0D/3D geometric multi-scale model of the coronary artery was established to simulate the hemodynamic effects of different counterpulsation modes on the vascular endothelium. In this model, the neural regulation caused by counterpulsation was thoroughly considered. Two clinical trials were carried out to verify the numerical calculation model. The results demonstrated that the increase in counterpulsation pressure amplitude and pressurization duration increased coronary blood perfusion and wall shear stress (WSS) and reduced the oscillatory shear index (OSI) of the vascular wall. However, the impact of pressurization duration was the predominant factor. The results of the standard model and the two real individual models indicated that a long pressurization duration would cause more hemodynamic risk areas by resulting in excessive WSS, which could not be reflected by the change in the Q value. Therefore, long-term pressurization during each cardiac cycle therapy is not recommended for patients with coronary heart disease and clinical treatment should not just pay attention to the change in the Q value. Additional physiological indicators can be used to evaluate the effects of counterpulsation treatment.

Impact of Malapposed and Overlapping Stents on Hemodynamics: A 2D Parametric Computational Fluid Dynamics Study

Despite significant progress, malapposed or overlapped stents are a complication that affects daily percutaneous coronary intervention (PCI) procedures. These malapposed stents affect blood flow and create a micro re-circulatory environment. These disturbances are often associated with a change in Wall Shear Stress (WSS), Time-averaged WSS (TAWSS), relative residence time (RRT) and oscillatory character of WSS and disrupt the delicate balance of vascular biology, providing a possible source of thrombosis and restenosis. In this study, 2D axisymmetric parametric computational fluid dynamics (CFD) simulations were performed to systematically analyze the hemodynamic effects of malapposition and stent overlap for two types of stents (drug-eluting stent and a bioresorbable stent). The results of the modeling are mainly analyzed using streamlines, TAWSS, oscillatory shear index (OSI) and RRT. The risks of restenosis and thrombus are evaluated according to commonly accepted thresholds for TAWSS and OSI. The small malapposition distances (MD) cause both low TAWSS and high OSI, which are potential adverse outcomes. The region of low OSI decrease with MD. Overlap configurations produce areas with low WSS and high OSI. The affected lengths are relatively insensitive to the overlap distance. The effects of strut size are even more sensitive and adverse for overlap configurations compared to a well-applied stent.

Acute Hemodynamic Improvement by Thermal Vasodilation inside the Abdominal and Iliac Arterial Segments of Young Sedentary Individuals

OBJECTIVE

To study the hemodynamic response to lower leg heating intervention (LLHI) inside the abdominal and iliac arterial segments (AIAS) of young sedentary individuals.

METHODS

A Doppler measurement of blood flow was conducted for 5 young sedentary adults with LLHI. Heating durations of 0, 20, and 40 min were considered. A lumped parameter model (LPM) was used to ascertain the hemodynamic mechanism. The hemodynamics were determined via numerical approaches.

RESULTS

Ultrasonography revealed that the blood flow waveform shifted upwards under LLHI; in particular, the mean flow increased significantly (p < 0.05) with increasing heating duration. The LPM showed that its mechanism depends on the reduction in afterload resistance, not on the inertia of blood flow and arterial compliance. The time-averaged wall shear stress, time-averaged production rate of nitric oxide, and helicity in the external iliac arteries increased more significantly than in other segments as the heating duration increased, while the oscillation shear index (OSI) and relative residence time (RRT) in the AIAS declined with increasing heating duration. There was a more obvious helicity response in the bilateral external iliac arteries than the OSI and RRT responses.

CONCLUSION

LLHI can effectively induce a positive hemodynamic environment in the AIAS of young sedentary individuals.

Predicting plaque vulnerability change using intravascular ultrasound + optical coherence tomography image-based fluid-structure interaction models and machine learning methods with patient follow-up data: a feasibility study

BACKGROUND

Coronary plaque vulnerability prediction is difficult because plaque vulnerability is non-trivial to quantify, clinically available medical image modality is not enough to quantify thin cap thickness, prediction methods with high accuracies still need to be developed, and gold-standard data to validate vulnerability prediction are often not available. Patient follow-up intravascular ultrasound (IVUS), optical coherence tomography (OCT) and angiography data were acquired to construct 3D fluid-structure interaction (FSI) coronary models and four machine-learning methods were compared to identify optimal method to predict future plaque vulnerability.

METHODS

Baseline and 10-month follow-up in vivo IVUS and OCT coronary plaque data were acquired from two arteries of one patient using IRB approved protocols with informed consent obtained. IVUS and OCT-based FSI models were constructed to obtain plaque wall stress/strain and wall shear stress. Forty-five slices were selected as machine learning sample database for vulnerability prediction study. Thirteen key morphological factors from IVUS and OCT images and biomechanical factors from FSI model were extracted from 45 slices at baseline for analysis. Lipid percentage index (LPI), cap thickness index (CTI) and morphological plaque vulnerability index (MPVI) were quantified to measure plaque vulnerability. Four machine learning methods (least square support vector machine, discriminant analysis, random forest and ensemble learning) were employed to predict the changes of three indices using all combinations of 13 factors. A standard fivefold cross-validation procedure was used to evaluate prediction results.

RESULTS

For LPI change prediction using support vector machine, wall thickness was the optimal single-factor predictor with area under curve (AUC) 0.883 and the AUC of optimal combinational-factor predictor achieved 0.963. For CTI change prediction using discriminant analysis, minimum cap thickness was the optimal single-factor predictor with AUC 0.818 while optimal combinational-factor predictor achieved an AUC 0.836. Using random forest for predicting MPVI change, minimum cap thickness was the optimal single-factor predictor with AUC 0.785 and the AUC of optimal combinational-factor predictor achieved 0.847.

CONCLUSION

This feasibility study demonstrated that machine learning methods could be used to accurately predict plaque vulnerability change based on morphological and biomechanical factors from multi-modality image-based FSI models. Large-scale studies are needed to verify our findings.

Analysis of Syk/PECAM-1 signaling pathway in low shear stress induced atherosclerosis based on ultrasound imaging

BACKGROUND AND OBJECTIVE

Low shear stress (LSS) has been demonstrated to be involved in function of vascular endothelial cells. Here we tested the hypothesis that activation of Syk played an important in LSS-induced atherosclerosis via PECAM-1 signaling pathway.

METHODS

In vitro, primary human umbilical vein endothelial cells (HUVECs) were stimulated with parallel plate flow chamber system for 12h under normal shear stress (NSS, 15dyne/cm²), LSS (5dyne/cm²) and high shear stress (HSS, 25dyne/cm²), respectively, followed by inflammatory response analysis. In vivo, animal models of rat fed atherogenic diet were treated with LSS stimulation by constricting abdominal aorta with a blunted needle (0.6mm in diameter). The spatial distribution of WSS of blood vessels was generated by WSS quantitative analysis software through color Doppler flow imaging with a high-frequency small animal ultrasound system. Small molecule R406, a well-demonstrated Syk inhibitor, was applied to animals as well as HUVEC cells.

RESULTS

In vivo, comparison with the control group was performed, the mean value of WSS distribution of blood vessels was lower in LSS model rat. LSS promoted expression of phosphorylated PECAM-1 (p-PECAM-1) and Syk in LSS model rats. Compared with control group, endothelial cells of the abdominal aorta become less elongated and more polygonal in LSS group, and had a slender shape in LSS with R406 group. In vitro, LSS increased the expression of p-PECAM-1, Syk and NF-κB in HUVECs. Inhibition of Syk attenuated LSS-induced inflammatory response.

CONCLUSIONS

Activation of Syk resulted in LSS-induced inflammatory response via PECAM-1 signaling pathway both in vitro and in vivo. Syk might be involved in morphological changes of ECs under the influence of LSS.

Association of Carotid Artery Endothelial Signal Intensity Gradient with Unilateral Large Artery Ischemic Stroke

BACKGROUND

Common carotid artery (CCA) and internal carotid artery (ICA) are aligned linearly, but their hemodynamic role in ischemic stroke has not been studied in depth.

OBJECTIVES

We aimed to investigate whether CCA and ICA endothelial shear stress (ESS) could be associated with the ischemic stroke of large artery atherosclerosis (LAA).

METHODS

We enrolled consecutive patients with unilateral ischemic stroke of LAA and healthy controls aged >60 years in the stroke center of Jeonbuk National University Hospital. All patients and controls were examined with carotid artery time-of-flight magnetic resonance angiography, and their endothelial signal intensity gradients (SIGs) were determined, as a measure of ESS. The effect of right or left unilateral stroke on the association between carotid artery endothelial SIG and ischemic stroke of LAA was assessed.

RESULTS

In total, the results from 132 patients with ischemic stroke of LAA and 121 controls were analyzed. ICA endothelial SIG showed significant and independent associations with the same-sided unilateral ischemic stroke of LAA, even after adjusting for the potential confounders including carotid stenosis, whereas CCA endothelial SIG showed a significant association with the presence of the ischemic stroke of LAA.

CONCLUSION

Although CCA and ICA are located with continuity, the hemodynamics and their roles in large artery ischemic stroke should be considered separately. Further studies are needed to delineate the pathophysiologic roles of ESS in CCA and ICA for large artery ischemic stroke.

Multi-patient study for coronary vulnerable plaque model comparisons: 2D/3D and fluid-structure interaction simulations

Several image-based computational models have been used to perform mechanical analysis for atherosclerotic plaque progression and vulnerability investigations. However, differences of computational predictions from those models have not been quantified at multi-patient level. In vivo intravascular ultrasound (IVUS) coronary plaque data were acquired from seven patients. Seven 2D/3D models with/without circumferential shrink, cyclic bending and fluid–structure interactions (FSI) were constructed for the seven patients to perform model comparisons and quantify impact of 2D simplification, circumferential shrink, FSI and cyclic bending plaque wall stress/strain (PWS/PWSn) and flow shear stress (FSS) calculations. PWS/PWSn and FSS averages from seven patients (388 slices for 2D and 3D thin-layer models) were used for comparison. Compared to 2D models with shrink process, 2D models without shrink process overestimated PWS by 17.26%. PWS change at location with greatest curvature change from 3D FSI models with/without cyclic bending varied from 15.07% to 49.52% for the seven patients (average = 30.13%). Mean Max-FSS, Min-FSS and Ave-FSS from the flow-only models under maximum pressure condition were 4.02%, 11.29% and 5.45% higher than those from full FSI models with cycle bending, respectively. Mean PWS and PWSn differences between FSI and structure-only models were only 4.38% and 1.78%. Model differences had noticeable patient variations. FSI and flow-only model differences were greater for minimum FSS predictions, notable since low FSS is known to be related to plaque progression. Structure-only models could provide PWS/PWSn calculations as good approximations to FSI models for simplicity and time savings in calculation.

Atherogenic index of plasma and the risk of rapid progression of coronary atherosclerosis beyond traditional risk factors

BACKGROUND AND AIMS

The atherogenic index of plasma (AIP) has been suggested as a marker of plasma atherogenicity. This study aimed to assess the association between AIP and the rapid progression of coronary atherosclerosis using serial coronary computed tomography angiography (CCTA).

METHODS

A total of 1488 adults (60.9 ± 9.2 years, 58.9% male) who underwent serial CCTA with a median inter-scan period of 3.4 years were included. AIP was defined as the base 10 logarithm of the ratio of the concentrations of triglyceride to high-density lipoprotein cholesterol. Rapid plaque progression (RPP) was defined as the change of percentage atheroma volume (PAV) ≥1.0%/year. All participants were divided into three groups based on AIP tertiles.

RESULTS

Baseline total PAV (median [interquartile range (IQR)]) (%) (group I [lowest]: 1.91 [0.00, 6.21] vs. group II: 2.82 [0.27, 8.83] vs. group III [highest]: 2.70 [0.41, 7.50]), the annual change of total PAV (median [IQR]) (%/year) (group I: 0.27 [0.00, 0.81] vs. group II: 0.37 [0.04, 1.11] vs. group III: 0.45 [0.06, 1.25]), and the incidence of RPP (group I: 19.7% vs. group II: 27.3% vs. group III: 31.4%) were significantly different among AIP tertiles (all p < 0.05). In multiple logistic regression analysis, the risk of RPP was increased in group III (odds ratio: 1.52, 95% confidence interval: 1.02-2.26; p = 0.042) compared to group I after adjusting for clinical factors and baseline total PAV.

CONCLUSIONS

Based on serial CCTA findings, AIP is an independent predictive marker for RPP beyond traditional risk factors.

Application of physics-based flow models in cardiovascular medicine: Current practices and challenges

Personalized physics-based flow models are becoming increasingly important in cardiovascular medicine. They are a powerful complement to traditional methods of clinical decision-making and offer a wealth of physiological information beyond conventional anatomic viewing using medical imaging data. These models have been used to identify key hemodynamic biomarkers, such as pressure gradient and wall shear stress, which are associated with determining the functional severity of cardiovascular diseases. Importantly, simulation-driven diagnostics can help researchers understand the complex interplay between geometric and fluid dynamic parameters, which can ultimately improve patient outcomes and treatment planning. The possibility to compute and predict diagnostic variables and hemodynamics biomarkers can therefore play a pivotal role in reducing adverse treatment outcomes and accelerate development of novel strategies for cardiovascular disease management.