Predicting Lipid-Rich Plaque Progression in Coronary Arteries Using Multimodal Imaging and Wall Shear Stress Signatures

BACKGROUND

Plaque composition and wall shear stress (WSS) magnitude act as well-established players in coronary plaque progression. However, WSS magnitude per se does not completely capture the mechanical stimulus to which the endothelium is subjected, since endothelial cells experience changes in the WSS spatiotemporal configuration on the luminal surface. This study explores WSS profile and lipid content signatures of plaque progression to identify novel biomarkers of coronary atherosclerosis.

METHODS

Thirty-seven patients with acute coronary syndrome underwent coronary computed tomography angiography, near-infrared spectroscopy intravascular ultrasound, and optical coherence tomography of at least 1 nonculprit vessel at baseline and 1-year follow-up. Baseline coronary artery geometries were reconstructed from intravascular ultrasound and coronary computed tomography angiography and combined with flow information to perform computational fluid dynamics simulations to assess the timeaveraged WSS magnitude (TAWSS) and the variability in the contraction/expansion action exerted by WSS on the endothelium, which can be assessed by the topological shear variation index (TSVI). Plaque progression was measured as intravascular ultrasound-derived percentage atheroma volume change at 1-year follow-up (Δplaque atheroma volume). Plaque composition information was extracted from near-infrared spectroscopy and optical coherence tomography.

RESULTS

Exposure to high TSVI and low TAWSS was associated with higher plaque progression (4.00±0.69% and 3.60±0.62%, respectively). Plaque composition acted synergistically with TSVI or TAWSS, resulting in the highest plaque progression (≥5.90%) at locations where lipid-rich content is exposed to high TSVI or low TAWSS.

CONCLUSIONS

Luminal exposure to high TSVI, solely or combined with a lipid-rich plaque phenotype, is associated with enhanced plaque progression at 1-year follow-up. Where plaque progression occurred, low TAWSS was also observed. These findings suggest TSVI, in addition to low TAWSS, as a potential biomechanical predictor for plaque progression, showing promise for clinical translation to improve patient prognosis.

Noninvasive Advanced Cardiovascular Magnetic Resonance-Derived Fontan Hemodynamics Are Associated With Reduced Kidney Function But Not Albuminuria

BACKGROUND

Kidney disease is the most important predictor of death in patients with a Fontan circulation, yet its clinical and hemodynamic correlates have not been well established.

METHODS AND RESULTS

A total of 53 ambulatory patients with a Fontan circulation (median age, 16.2 years, 52.8% male patients) underwent advanced cardiovascular magnetic resonance assessment, including 4-dimensional flow imaging and computational fluid dynamics. Estimated glomerular filtration rate (eGFR) <90 mL/min per 1.73 m2 was observed in 20.8% and albumin-to-creatinine ratio >3 mg/mmol in 39.6%. The average eGFR decline rate was -1.83 mL/min per 1.73 m2 per year (95% CI, -2.67 to -0.99; P<0.001). Lower eGFR was associated with older age, larger body surface area at examination, longer time since Fontan procedure, and lower systemic ventricular ejection fraction. Higher albumin-to-creatinine ratio was associated with absence of fenestration at the Fontan operation, and older age and lower systemic ventricular ejection fraction at the assessment. Lower cross-sectional area of the Fontan conduit indexed to flow (r=0.32, P=0.038), higher inferior vena cava-conduit velocity mismatch factor (r=-0.35, P=0.022), higher kinetic energy indexed to flow in the total cavopulmonary connection (r=-0.59, P=0.005), and higher total cavopulmonary connection resistance (r=-0.42, P=0.005 at rest; r=-0.43, P=0.004 during exercise) were all associated with lower eGFR but not with albuminuria.

CONCLUSIONS

Kidney dysfunction and albuminuria are common among clinically well adolescents and young adults with a Fontan circulation. Advanced cardiovascular magnetic resonance-derived metrics indicative of declining Fontan hemodynamics are associated with eGFR and might serve as targets to improve kidney health. Albuminuria might be driven by other factors that need further investigation.

Serial RV wall stress measurements: association with right ventricular function in repaired Tetralogy of Fallot patients

Background

Optimal timing of pulmonary valve replacement (PVR) in Tetralogy of Fallot (TOF) patients remains challenging. Ventricular wall stress is considered to be an early marker of right ventricular (RV) dysfunction.

Objectives

To investigate the association of RV wall stresses and their change over time with functional parameters in TOF patients.

Methods

Ten TOF patients after surgical repair with moderate/severe pulmonary regurgitation were included. At two timepoints (median follow-up time 7.2 years), patient-specific computational biventricular models for wall stress assessment were created using CMR short-axis cine images and echocardiography-based RV pressures. RV ejection fraction (RVEF), NT-proBNP and cardiopulmonary exercise tests were used as outcome measures reflecting RV function. Associations between regional RV diastolic wall stress and RV function were investigated using linear mixed models.

Results

Increased wall stress correlated with lower RV mass (rrm = -0.70, p = 0.017) and lower RV mass-to-volume (rrm = -0.80, p = 0.003) using repeated measures. Wall stress decreased significantly over time, especially in patients with a stable RVEF (p < 0.001). Higher wall stress was independently associated with lower RVEF, adjusted for left ventricular ejection fraction, RV end-diastolic volume and time since initial surgery (decrease of 1.27% RVEF per kPa increase in wall stress, p = 0.029) using repeated measurements. No association was found between wall stress, NT-proBNP, and exercise capacity.

Conclusions

Using a computational method to calculate wall stress locally in geometrically complex ventricles, we demonstrated that lower wall stress might be important to maintain ventricular function. RV wall stress assessment can be used in serial follow-up, and is potentially an early marker of impending RV dysfunction.

Plaque burden is associated with minimal intimal coverage following drug-eluting stent implantation in an adult familial hypercholesterolemia swine model

Safety and efficacy of coronary drug-eluting stents (DES) are often preclinically tested using healthy or minimally diseased swine. These generally show significant fibrotic neointima at follow-up, while in patients, incomplete healing is often observed. The aim of this study was to investigate neointima responses to DES in swine with significant coronary atherosclerosis. Adult familial hypercholesterolemic swine (n = 6) received a high fat diet to develop atherosclerosis. Serial OCT was performed before, directly after, and 28 days after DES implantation (n = 14 stents). Lumen, stent and plaque area, uncovered struts, neointima thickness and neointima type were analyzed for each frame and averaged per stent. Histology was performed to show differences in coronary atherosclerosis. A range of plaque size and severity was found, from healthy segments to lipid-rich plaques. Accordingly, neointima responses ranged from uncovered struts, to minimal neointima, to fibrotic neointima. Lower plaque burden resulted in a fibrotic neointima at follow-up, reminiscent of minimally diseased swine coronary models. In contrast, higher plaque burden resulted in minimal neointima and more uncovered struts at follow-up, similarly to patients' responses. The presence of lipid-rich plaques resulted in more uncovered struts, which underscores the importance of advanced disease when performing safety and efficacy testing of DES.

Wall shear stress-related plaque growth of lipid-rich plaques in human coronary arteries: an near-infrared spectroscopy and optical coherence tomography study

AIMS

Low wall shear stress (WSS) is acknowledged to play a role in plaque development through its influence on local endothelial function. Also, lipid-rich plaques (LRPs) are associated with endothelial dysfunction. However, little is known about the interplay between WSS and the presence of lipids with respect to plaque progression. Therefore, we aimed to study the differences in WSS-related plaque progression between LRPs, non-LRPs, or plaque-free regions in human coronary arteries.

METHODS AND RESULTS

In the present single-centre, prospective study, 40 patients who presented with an acute coronary syndrome successfully underwent near-infrared spectroscopy intravascular ultrasound (NIRS-IVUS) and optical coherence tomography (OCT) of at least one non-culprit vessel at baseline and completed a 1-year follow-up. WSS was computed applying computational fluid dynamics to a three-dimensional reconstruction of the coronary artery based on the fusion of the IVUS-segmented lumen with a CT-derived centreline, using invasive flow measurements as boundary conditions. For data analysis, each artery was divided into 1.5 mm/45° sectors. Plaque growth based on IVUS-derived percentage atheroma volume change was compared between LRPs, non-LRPs, and plaque-free wall segments, as assessed by both OCT and NIRS. Both NIRS- and OCT-detected lipid-rich sectors showed a significantly higher plaque progression than non-LRPs or plaque-free regions. Exposure to low WSS was associated with a higher plaque progression than exposure to mid or high WSS, even in the regions classified as a plaque-free wall. Furthermore, low WSS and the presence of lipids had a synergistic effect on plaque growth, resulting in the highest plaque progression in lipid-rich regions exposed to low shear stress.

CONCLUSION

This study demonstrates that NIRS- and OCT-detected lipid-rich regions exposed to low WSS are subject to enhanced plaque growth over a 1-year follow-up. The presence of lipids and low WSS proves to have a synergistic effect on plaque growth.

Haemodynamic performance of 16-20-mm extracardiac Goretex conduits in adolescent Fontan patients at rest and during simulated exercise

OBJECTIVES

To date, it is not known if 16-20-mm extracardiac conduits are outgrown during somatic growth from childhood to adolescence. This study aims to determine total cavopulmonary connection (TCPC) haemodynamics in adolescent Fontan patients at rest and during simulated exercise and to assess the relationship between conduit size and haemodynamics.

METHODS

Patient-specific, magnetic resonance imaging-based computational fluid dynamic models of the TCPC were performed in 51 extracardiac Fontan patients with 16-20-mm conduits. Power loss, pressure gradient and normalized resistance were quantified in rest and during simulated exercise. The cross-sectional area (CSA) (mean and minimum) of the vessels of the TCPC was determined and normalized for flow rate (mm2/l/min). Peak (predicted) oxygen uptake was assessed.

RESULTS

The median age was 16.2 years (Q1-Q3 14.0-18.2). The normalized mean conduit CSA was 35-73% smaller compared to the inferior and superior vena cava, hepatic veins and left/right pulmonary artery (all P < 0.001). The median TCPC pressure gradient was 0.7 mmHg (Q1-Q3 0.5-0.8) and 2.0 (Q1-Q3 1.4-2.6) during rest and simulated exercise, respectively. A moderate-strong inverse non-linear relationship was present between normalized mean conduit CSA and TCPC haemodynamics in rest and exercise. TCPC pressure gradients of ≥1.0 at rest and ≥3.0 mmHg during simulated exercise were observed in patients with a conduit CSA ≤ 45 mm2/l/min and favourable haemodynamics (<1 mmHg during both rest and exercise) in conduits ≥125 mm2/l/min. Normalized TCPC resistance correlated with (predicted) peak oxygen uptake.

CONCLUSIONS

Extracardiac conduits of 16-20 mm have become relatively undersized in most adolescent Fontan patients leading to suboptimal haemodynamics.

Divergence of the normalized wall shear stress as an effective computational template of low-density lipoprotein polarization at the arterial blood-vessel wall interface

BACKGROUND AND OBJECTIVE

Near-wall transport of low-density lipoproteins (LDL) in arteries plays a relevant role in the initiation of atherosclerosis. Although it can be modelled in silico by coupling the Navier-Stokes equations with the 3D advection-diffusion (AD) equation, the associated computational cost is high. As wall shear stress (WSS) represents a first-order approximation of the near-wall velocity in arteries, we aimed at identifying computationally convenient WSS-based quantities to infer LDL near-wall transport based on the underlying near-wall hemodynamics in five models of three human arterial districts (aorta, carotid bifurcations, coronary arteries). The simulated LDL transport and its WSS-based surrogates were qualitatively compared with in vivo longitudinal measurements of wall thickness growth on the coronary artery models.

METHODS

Numerical simulations of blood flow coupled with AD equations for LDL transport and blood-wall transfer were performed. The co-localization of the simulated LDL concentration polarization patterns with luminal surface areas characterized by low cycle-average WSS, near-wall flow stagnation and WSS attracting patterns was quantitatively assessed by the similarity index (SI). In detail, the latter two represent features of the WSS topological skeleton, obtained respectively through the Lagrangian tracking of surface-born particles, and the Eulerian analysis of the divergence of the normalized cycle-average WSS vector field.

RESULTS

Convergence of the solution of the AD problem required the simulation of 3 (coronary artery) to 10 (aorta) additional cardiac cycles with respect to the Navier-Stokes problem. Co-localization results underlined that WSS topological skeleton features indicating near-wall flow stagnation and WSS attracting patterns identified LDL concentration polarization profiles more effectively than low WSS, as indicated by higher SI values (SI range: 0.17-0.50 for low WSS; 0.24-0.57 for WSS topological skeleton features). Moreover, the correspondence between the simulated LDL uptake and WSS-based quantities profiles with the in vivo measured wall thickness growth in coronary arteries appears promising.

CONCLUSIONS

The recently introduced Eulerian approach for identifying WSS attracting patterns from the divergence of normalized WSS provides a computationally affordable template of the LDL polarization at the arterial blood-wall interface without simulating the AD problem. It thus candidates as an effective biomechanical tool for elucidating the mechanistic link amongst LDL transfer at the arterial blood-wall interface, WSS and atherogenesis.

Near-infrared spectroscopy to predict plaque progression in plaque-free artery regions

BACKGROUND

Positive near-infrared spectroscopy (NIRS) signals might be encountered in areas without evident artery wall thickening, being typically perceived as artefacts.

AIMS

We aimed to evaluate the utility of NIRS to identify artery wall regions associated with an increase in wall thickness (WT) as assessed by serial intravascular ultrasound (IVUS) and optical coherence tomography (OCT).

METHODS

In this prospective, single-centre study, patients presenting with acute coronary syndrome (ACS) underwent NIRS-IVUS and OCT assessment of a non-culprit artery at baseline and 12-month follow-up. For each vessel, 1.5 mm segments were identified, matched and divided into 45 sectors. The relationship between the change in IVUS-based WT (DWT) and the presence of NIRS-positive signals and OCT-detected lipid was evaluated using linear mixed models.

RESULTS

A total of 37 patients (38 vessels, 6,936 matched sectors) were analysed at baseline and 12 months. A total of 140/406 (34.5%) NIRS (+) sectors and 513/1,575 (32.6%) OCT-lipid (+) sectors were found to be located in thin (WT<0.5 mm) wall sectors. In the thin wall sectors, an increase in WT was significantly more pronounced in NIRS (+) vs NIRS (-) sectors (0.11 mm vs 0.06 mm, p<0.001). In the thick wall sectors, there was a decrease in WT observed that was less pronounced in the NIRS (+) versus NIRS (-) sectors (-0.08 mm vs -0.09 mm, p<0.001). Thin wall NIRS (+) OCT-lipid (+) sectors showed significant wall thickening (DWT=0.13 mm).

CONCLUSIONS

NIRS-positive signals in otherwise non-diseased arterial walls as assessed by IVUS could identify vessel wall regions prone to WT increase over 12-month follow-up. Our observations suggest that NIRS-positive signals in areas without evident wall thickening by IVUS should no longer be viewed as benign or imaging artefact.

Vascular biomechanics and molecular disease activity in the thoracic aorta: a novel imaging method

AIMS

The influence haemodynamics have on vessel wall pathobiology in aortic disease is incomplete. This aim of this study was to develop a repeatable method for assessing the relationship between aortic wall shear stress (WSS) and disease activity by fusing 4D flow cardiovascular magnetic resonance (CMR) with hybrid positron emission tomography (PET).

METHODS AND RESULTS

As part of an ongoing clinical trial, patients with bicuspid aortic valve (BAV) were prospectively imaged with both 18F-sodium fluoride (18F-NaF) PET, a marker of calcification activity, and 4D flow CMR. We developed novel software allowing accurate 3D co-registration and high-resolution comparison of aortic peak systolic WSS and 18F-NaF PET uptake (maximum tissue-to-background ratio). Intra-observer repeatability of both measurements was determined using Bland-Altman plots and intra-class correlation coefficients (ICCs). The relationship between localized WSS and 18F-NaF uptake was analysed using linear mixed-effect models. Twenty-three patients with BAV (median age 50 [44-55] years, 22% female) were included. Intra-observer repeatability for WSS (ICC = 0.92) and 18F-NaF (ICC = 0.91) measurements obtained within 1.4 ± 0.6 cm2 regions of interest was excellent. On multivariable analysis, 18F-NaF PET uptake was independently and negatively associated with WSS as well as diastolic blood pressure (both P < 0.05), adjusted for age.

CONCLUSION

Fused assessment of WSS and 18F-NaF PET uptake is feasible and repeatable, demonstrating a clear association between these two factors. This high spatial resolution approach has major potential to advance our understanding of the relationship between vascular haemodynamics and disease activity.

Wall shear stress angle is associated with aortic growth in bicuspid aortic valve patients

AIMS

Aortic wall shear stress (WSS) distributions in bicuspid aortic valve (BAV) patients have been associated with aortic dilatation, but prospective, longitudinal data are missing. This study assessed differences in aortic WSS distributions between BAV patients and healthy controls and determined the association of WSS with aortic growth in patients.

METHODS AND RESULTS

Sixty subjects underwent four-dimensional (4D) flow cardiovascular magnetic resonance of the thoracic aorta (32 BAV patients and 28 healthy controls). Peak velocity, pulse wave velocity, aortic distensibility, peak systolic WSS (magnitude, axial, and circumferential), and WSS angle were assessed. WSS angle is defined as the angle between the WSSmagnitude and WSSaxial component. In BAV patients, three-year computed tomography angiography-based aortic volumetric growth was determined in the proximal and entire ascending aorta. WSSaxial was significantly lower in BAV patients compared with controls (0.93 vs. 0.72 Pa, P = 0.047) and WSScircumferential and WSS angle were significantly higher (0.29 vs. 0.64 Pa and 18° vs. 40°, both P < 0.001). Significant volumetric growth of the proximal ascending aorta occurred in BAV patients (from 49.1 to 52.5 cm3, P = 0.003). In multivariable analysis corrected for baseline aortic volume and diastolic blood pressure, WSS angle was the only parameter independently associated with proximal aortic growth (P = 0.031). In the entire ascending aorta, besides the WSS angle, the WSSmagnitude was also independently associated with growth.

CONCLUSION

Increased WSScircumferential and especially WSS angle are typical in BAV patients. WSS angle was found to predict aortic growth. These findings highlight the potential role of WSS measurements in BAV patients to stratify patients at risk for aortic dilation.

Sex-related differences in plaque characteristics and endothelial shear stress related plaque-progression in human coronary arteries

BACKGROUND AND AIMS

Clinical atherosclerosis manifestations are different in women compared to men. Since endothelial shear stress (ESS) is known to play a critical role in coronary atherosclerosis development, we investigated differences in anatomical characteristics and endothelial shear stress (ESS)-related plaque growth in human coronary arteries in men compared to women.

METHODS

1183 coronary arteries (male/female: 944/239) from the PREDICTION study were studied for differences in artery/plaque and ESS characteristics, and ESS-related plaque progression (6-10 months follow-up) among men and women and after stratification for age. All characteristics were derived from IVUS-based vascular profiling and reported per 3 mm-segments (13,030 3-mm-segments (male/female: 10,465/2,565)).

RESULTS

Coronary arteries and plaques were significantly smaller in females compared to males; but no important differences were observed in plaque burden, ESS and rate of plaque progression. Change in plaque burden was inversely related to ESS (p<0.001) with no difference between women versus men (β: -0.62 ± 0.13 vs -0.68 ± 0.05, p=0.62). However, stratification for age demonstrated that ESS-related plaque growth was more marked in young women compared to men (<55 years, β: -2.02 ± 0.61 vs -0.33 ± 0.10, p=0.007), reducing in magnitude over the age-categories up till 75 years.

CONCLUSIONS

Coronary artery and plaque size are smaller in women compared to men, but ESS and ESS- related plaque progression were similar. Sex-related differences in ESS-related plaque growth were evident after stratification for age. These observations suggest that although the fundamental processes of atherosclerosis progression are similar in men versus women, plaque progression may be influenced by age within gender.

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.

Moving shear stress towards the clinic: preclinical comparison of optical coherence tomography-based versus angiography-based time-averaged wall shear stress estimations

Introduction

Identification of coronary atherosclerotic plaques at risk of causing future acute coronary syndromes remains a major unmet clinical challenge. The addition of vessel biomechanics to intracoronary imaging derived evaluation of plaque morphology, improves identification of plaques likely to develop high risk features. We and others have developed a framework for intracoronary imaging (optical coherence tomography [OCT]) based 3D reconstructions of coronary arteries for computational fluid dynamics (CFD) simulations of shear stress, which are considered the current gold standard approach for quantification of coronary arterial haemodynamics. However, these approaches are time consuming and computationally intensive, resulting in a barrier to clinical uptake.

Purpose

Determination of time averaged wall shear stress (TAWSS) based on 3D coronary geometries from non-invasive 3D Quantitative Coronary Angiography (3D-QCA) has recently been developed (Pie Medical Imaging, Netherlands), which enables results of shear stress simulations to be available within 30 minutes. We sought to compare TAWSS determined from 3D-QCA with gold standard OCT-based CFD simulations in both normal and stenotic arteries in minipigs.

Methods

15 normal and 5 stenotic minipig coronary arteries were studied. Anatomically matched 3D arterial geometries were reconstructed from 3D-QCA and OCT using common centrelines. Boundary conditions for simulations included directly measured inlet blood velocities; parabolic inlet flow profiles, zero pressure outlet; no-slip arterial walls; blood density: 1.05 g/ml; blood dynamic viscosity: 0.035 g/cm.s. Blood was modelled as Newtonian. 3D-QCA TAWSS was obtained with a Kratos Multi-Physics CFD solver. OCT-based simulations were performed using Abaqus/CFD v6.14. TAWSS was calculated for 80 axially matched segments for both methods (1200 and 400 paired comparisons for normal and stenotic arteries, respectively). Data were analysed using Bland-Altman and Wilcoxon matched-pairs signed ranked tests.

Results

Computation times for 3D-QCA and OCT-based CFD were approximately 30 minutes and 2 hours respectively. Axial profiles of TAWSS were similar between the two methods and there was agreement in TAWSS magnitudes and narrow 95% limits of agreement (Figure 1 and Figure 2). Using co-registered TAWSS maps generated by each method, we find similar spatial regional distributions of TAWSS in both normal and stenotic arteries.

Conclusions

Our data suggest that 3D-QCA based TAWSS is feasible in both normal and stenotic arteries. Spatial TAWSS distributions between the two methods are similar with agreement in matched TAWSS comparisons, though there are some small systematic differences in the absolute values of TAWSS, due to different resultant arterial geometries. These encouraging data suggest that further clinical evaluation of rapid TAWSS from 3D-QCA is warranted, which may facilitate clinical adoption of TAWSS assessment.

Funding Acknowledgement

Type of funding sources: Public grant(s) – EU funding. Main funding source(s): Medical Research Council

The Influence of Respiration on Blood Flow in the Fontan Circulation: Insights for Imaging-Based Clinical Evaluation of the Total Cavopulmonary Connection

Congenital heart disease is the most common birth defect and functionally univentricular heart defects represent the most severe end of this spectrum. The Fontan circulation provides an unique solution for single ventricle patients, by connecting both caval veins directly to the pulmonary arteries. As a result, the pulmonary circulation in Fontan palliated patients is characterized by a passive, low-energy circulation that depends on increased systemic venous pressure to drive blood toward the lungs. The absence of a subpulmonary ventricle led to the widely believed concept that respiration, by sucking blood to the pulmonary circulation during inspiration, is of great importance as a driving force for antegrade blood flow in Fontan patients. However, recent studies show that respiration influences pulsatility, but has a limited effect on net forward flow in the Fontan circulation. Importantly, since MRI examination is recommended every 2 years in Fontan patients, clinicians should be aware that most conventional MRI flow sequences do not capture the pulsatility of the blood flow as a result of the respiration. In this review, the unique flow dynamics influenced by the cardiac and respiratory cycle at multiple locations within the Fontan circulation is discussed. The impact of (not) incorporating respiration in different MRI flow sequences on the interpretation of clinical flow parameters will be covered. Finally, the influence of incorporating respiration in advanced computational fluid dynamic modeling will be outlined.

Comparison of Swine and Human Computational Hemodynamics Models for the Study of Coronary Atherosclerosis

Coronary atherosclerosis is a leading cause of illness and death in Western World and its mechanisms are still non completely understood. Several animal models have been used to 1) study coronary atherosclerosis natural history and 2) propose predictive tools for this disease, that is asymptomatic for a long time, aiming for a direct translation of their findings to human coronary arteries. Among them, swine models are largely used due to the observed anatomical and pathophysiological similarities to humans. However, a direct comparison between swine and human models in terms of coronary hemodynamics, known to influence atherosclerotic onset/development, is still lacking. In this context, we performed a detailed comparative analysis between swine- and human-specific computational hemodynamic models of coronary arteries. The analysis involved several near-wall and intravascular flow descriptors, previously emerged as markers of coronary atherosclerosis initiation/progression, as well as anatomical features. To do that, non-culprit coronary arteries (18 right–RCA, 18 left anterior descending–LAD, 13 left circumflex–LCX coronary artery) from patients presenting with acute coronary syndrome were imaged by intravascular ultrasound and coronary computed tomography angiography. Similarly, the three main coronary arteries of ten adult mini-pigs were also imaged (10 RCA, 10 LAD, 10 LCX). The geometries of the imaged coronary arteries were reconstructed (49 human, 30 swine), and computational fluid dynamic simulations were performed by imposing individualized boundary conditions. Overall, no relevant differences in 1) wall shear stress-based quantities, 2) intravascular hemodynamics (in terms of helical flow features), and 3) anatomical features emerged between human- and swine-specific models. The findings of this study strongly support the use of swine-specific computational models to study and characterize the hemodynamic features linked to coronary atherosclerosis, sustaining the reliability of their translation to human vascular disease.

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.

Wall shear stress angle determines aortic growth in patients with bicuspid aortic valves

Funding Acknowledgements

Type of funding sources: Private grant(s) and/or Sponsorship. Main funding source(s): The Dutch Heart Foundation

Background

Patients with bicuspid aortic valve (BAV) have altered flow velocity patterns with different wall shear stress (WSS) distributions in the ascending aorta compared to patients with tricuspid aortic valves. These WSS distributions are associated with aortic dilatation in cross sectional studies, however, longitudinal data demonstrating a potential causative role is missing.

Purpose

The aim of this study was to assess the differences in WSS distributions between BAV patients and healthy subjects and to determine the predictive value of WSS for aortic growth in patients with a BAV.

Methods

Forty patients with a BAV and 32 healthy matched subjects were prospectively studied by 4D-flow cardiovascular magnetic resonance (CMR). Peak velocity, pulse wave velocity, aortic distensibility, peak systolic WSS (magnitude), the different WSS components (axial and circumferential), and WSS angle were assessed in the proximal ascending aorta. WSS angle was defined as the angle between the WSSmagnitude and WSSaxial component. In the BAV patients, aortic volumetric growth over three years was determined in the proximal ascending aorta (first 5cm) based on CT angiography. Multivariate linear regression analysis was used to identify independent predictors of aortic volumetric growth.

Results

Of the BAV patients, 21 (53%) had a left-right fusion pattern and eight patients had Turner syndrome. WSSaxial was significantly lower in BAV patients compared to healthy subjects (p = 0.008) and WSScircumferential and WSS angle were significantly higher (both p &lt; 0.001, see Figure). WSSmagnitude, pulse wave velocity, and aorta distensibility were not statistically significant different. WSSmagnitude (0.69 N/m² [0.51-0.81] vs 1.08 N/m² [0.89-1.24], p = 0.005), WSSaxial (0.50 N/m² [0.39-0.61] vs 0.72 N/m²  [0.54-0.94], p = 0.015) and WSScircumferential (0.34 N/m² [0.32-0.46] vs 0.64 N/m² [0.47-0.81], p = 0.008) were significantly lower in BAV Turner patients compared to BAV non-Turner patients, while WSS angle (40° [34-41] vs 40° [32-48], p = 0.607) was not statistically significant different. During a follow-up of three years, there was a significant growth of the proximal ascending aorta in the BAV patients (1.2 cm3 [-0.2-2.5], p = 0.001). In multivariate analysis corrected for baseline aortic volume and diastolic blood pressure, WSS angle was the only independent predictor for proximal aortic volume growth (β=0.108, p = 0.030).

Conclusions

Increased WSScircumferential and especially WSS angle are present in patients with BAV. WSS angle was the only independent predictor of aortic growth. These findings highlight the potential role of  WSS measurements in patients with BAV to stratify patients at risk for aortic dilation.

Dolichoarteriopathies of the extracranial internal carotid artery: The Plaque At RISK study

BACKGROUND AND PURPOSE

Dolichoarteriopathies of the extracranial part of the internal carotid artery (ICA) are associated with cerebrovascular events, yet information on their prevalence and risk factors remains limited. The aim of the present study therefore was to study the prevalence and risk factors of dolichoarteriopathies in a sample of patients with cerebrovascular symptoms from the Plaque At RISK (PARISK) study.

METHODS

In a random sample of 100 patients from the PARISK study, multidetector computed tomography angiography (MDCTA) was performed as part of clinical workup. On MDCTA, we evaluated the extracranial trajectory of the ICA by measuring the length (in millimeters), the tortuosity index (TI; defined as the ICA length divided by the shortest possible distance from bifurcation to skull base), and dolichoarteriopathy type (tortuosity, coiling or kinking). Next, we investigated the association between cardiovascular risk factors and these measurements using linear and logistic regression analyses.

RESULTS

The mean (standard deviation) length of the ICA was 93 (± 14) mm, with a median (interquartile range) TI of 1.2 (1.1-1.3). The overall prevalence of dolichoarteriopathies was 69%, with tortuosity being the most common (72%), followed by coiling (20%), and kinking (8%). We found that age and obesity were associated with a higher TI: difference per 10-year increase in age: 0.05 (95% confidence interval [CI] 0.02-0.08) and 0.16 (95% CI 0.07-0.25) for obesity. Obesity and hypercholesterolemia were associated with a more severe type of dolichoarteriopathy (odds ratio [OR] 2.07 [95% CI 1.04-4.12] and OR 2.17 [95% CI 1.02-4.63], respectively).

CONCLUSION

Dolichoarteriopathies in the extracranial ICA are common in patients with cerebrovascular symptoms, and age, obesity and hypercholesterolemia may play an important role in the pathophysiology of these abnormalities.

Vulnerable plaques and patients: state-of-the-art

Despite advanced understanding of the biology of atherosclerosis, coronary heart disease remains the leading cause of death worldwide. Progress has been challenging as half of the individuals who suffer sudden cardiac death do not experience premonitory symptoms. Furthermore, it is well-recognized that also a plaque that does not cause a haemodynamically significant stenosis can trigger a sudden cardiac event, yet the majority of ruptured or eroded plaques remain clinically silent. In the past 30 years since the term 'vulnerable plaque' was introduced, there have been major advances in the understanding of plaque pathogenesis and pathophysiology, shifting from pursuing features of 'vulnerability' of a specific lesion to the more comprehensive goal of identifying patient 'cardiovascular vulnerability'. It has been also recognized that aside a thin-capped, lipid-rich plaque associated with plaque rupture, acute coronary syndromes (ACS) are also caused by plaque erosion underlying between 25% and 60% of ACS nowadays, by calcified nodule or by functional coronary alterations. While there have been advances in preventive strategies and in pharmacotherapy, with improved agents to reduce cholesterol, thrombosis, and inflammation, events continue to occur in patients receiving optimal medical treatment. Although at present the positive predictive value of imaging precursors of the culprit plaques remains too low for clinical relevance, improving coronary plaque imaging may be instrumental in guiding pharmacotherapy intensity and could facilitate optimal allocation of novel, more aggressive, and costly treatment strategies. Recent technical and diagnostic advances justify continuation of interdisciplinary research efforts to improve cardiovascular prognosis by both systemic and 'local' diagnostics and therapies. The present state-of-the-art document aims to present and critically appraise the latest evidence, developments, and future perspectives in detection, prevention, and treatment of 'high-risk' plaques occurring in 'vulnerable' patients.

Non-uniform mixing of hepatic venous flow and inferior vena cava flow in the Fontan conduit

Fontan patients require a balanced hepatic blood flow distribution (HFD) to prevent pulmonary arteriovenous malformations. Currently, HFD is quantified by tracking Fontan conduit flow, assuming hepatic venous (HV) flow to be uniformly distributed within the Fontan conduit. However, this assumption may be unvalid leading to inaccuracies in HFD quantification with potential clinical impact. The aim of this study was to (i) assess the mixing of HV flow and inferior vena caval (IVC) flow within the Fontan conduit and (ii) quantify HFD by directly tracking HV flow and quantitatively comparing results with the conventional approach. Patient-specific, time-resolved computational fluid dynamic models of 15 total cavopulmonary connections were generated, including the HV and subhepatic IVC. Mixing of HV and IVC flow, on a scale between 0 (no mixing) and 1 (perfect mixing), was assessed at the caudal and cranial Fontan conduit. HFD was quantified by tracking particles from the caudal (HFD_{caudal conduit}) and cranial (HFD_{cranial conduit}) conduit and from the hepatic veins (HFD_HV). HV flow was non-uniformly distributed at both the caudal (mean mixing 0.66 ± 0.13) and cranial (mean 0.79 ± 0.11) level within the Fontan conduit. On a cohort level, differences in HFD between methods were significant but small; HFD_HV (51.0 ± 20.6%) versus HFD_{caudal conduit} (48.2 ± 21.9%, p = 0.033) or HFD_{cranial conduit} (48.0 ± 21.9%, p = 0.044). However, individual absolute differences of 8.2–14.9% in HFD were observed in 4/15 patients. HV flow is non-uniformly distributed within the Fontan conduit. Substantial individual inaccuracies in HFD quantification were observed in a subset of patients with potential clinical impact.

Multidirectional wall shear stress promotes advanced coronary plaque development: comparing five shear stress metrics

Aims

Atherosclerotic plaque development has been associated with wall shear stress (WSS). However, the multidirectionality of blood flow, and thus of WSS, is rarely taken into account. The purpose of this study was to comprehensively compare five metrics that describe (multidirectional) WSS behaviour and assess how WSS multidirectionality affects coronary plaque initiation and progression.

Methods and results

Adult familial hypercholesterolaemic pigs (n = 10) that were fed a high-fat diet, underwent imaging of the three main coronary arteries at three-time points [3 (T₁), 9 (T₂), and 10–12 (T₃) months]. Three-dimensional geometry of the arterial lumen, in combination with local flow velocity measurements, was used to calculate WSS at T₁ and T₂. For analysis, arteries were divided into 3 mm/45° sectors (n = 3648). Changes in wall thickness and final plaque composition were assessed with near-infrared spectroscopy–intravascular ultrasound, optical coherence tomography imaging, and histology. Both in pigs with advanced and mild disease, the highest plaque progression rate was exclusively found at low time-averaged WSS (TAWSS) or high multidirectional WSS regions at both T₁ and T₂. However, the eventually largest plaque growth was located in regions with initial low TAWSS or high multidirectional WSS that, over time, became exposed to high TAWSS or low multidirectional WSS at T₂. Besides plaque size, also the presence of vulnerable plaque components at the last time point was related to low and multidirectional WSS. Almost all WSS metrics had good predictive values for the development of plaque (47–50%) and advanced fibrous cap atheroma (FCA) development (59–61%).

Conclusion

This study demonstrates that low and multidirectional WSS promote both initiation and progression of coronary atherosclerotic plaques. The high-predictive values of the multidirectional WSS metrics for FCA development indicate their potential as an additional clinical marker for the vulnerable disease.

Serial MRI-based right ventricular mechanical wall stress measurements and their association with right ventricle function in patients with repaired Tetralogy of Fallot

Funding Acknowledgements

Type of funding sources: Foundation. Main funding source(s): Stichting Hartekind en Thorax Foundation

Background

Optimal timing of pulmonary valve replacement (PVR) in Tetralogy of Fallot (TOF) patients remains challenging. Wall stress is considered to be a possible early marker of right ventricular (RV) dysfunction. With patient-specific computational models, wall stress can be determined regionally and with high accuracy, especially in complex shaped ventricles such as in TOF patients. We aimed to 1) develop patient-specific computational models to assess RV diastolic wall stresses and 2) investigate the association of wall stresses and their change over time with functional parameters in TOF patients.

Methods

Repaired TOF patients with at least moderate pulmonary regurgitation (PR) and prior to PVR were included. MRI-based patient-specific computational ventricular models were created (figure). The ventricular geometry was created by stacking endo- and epicardial contours traced on short axis SSFP cine images. Pressure in the right ventricle was estimated from echocardiography. Mid-diastolic wall stress in the RV free wall was analysed globally and regionally (basal, mid, apical, anterior, lateral and posterior) at two time points. RV ejection fraction (RVEF), NT-proBNP and exercise tests (% maximum predicted workload) were used as outcomes for RV function. Associations between wall stresses and outcomes were investigated using linear mixed models adjusted for follow-up duration.

Results

Five males and five females were included with an age at baseline of 24 (IQR 16-28) years and RV end-diastolic volume of 140 (IQR 127-144) ml/m2. The period between the two time points was 7.0 (IQR 5.8-7.3) years. Global wall stress of the RV free wall combining both time points was 5.8 kPa (IQR 5.2-7.2). There was no statistical difference between baseline and follow-up global wall stress. The mean wall stresses in the mid region was 1.69 kPa (p &lt; 0.01) higher than in the basal region and was 1.05 kPa (p = 0.03) higher than in the apical region cross-sectionally. The wall stress also increased more in the mid region compared to basal and apical region, corrected for duration of follow-up. Patients with more severe PR at baseline demonstrated a higher increase of global wall stress over time (p = 0.02), especially in lateral free wall. Higher global free wall stresses were cross-sectionally independently associated with lower RVEF, adjusted for LVEF and RVEDV (β=-1.29 % RVEF per kPa increase in wall stress, p = 0.01). This association was most prominent in the anterior, basal and mid part. No statistically significant association was found between wall stress, NT-proBNP, and exercise capacity.

Conclusions

This study generated a novel MRI-based method to calculate wall stress in geometrically complex ventricles. Wall stress associated negatively with RVEF in patients with TOF and PR. This promising tool for RV wall stress analysis can be used in future larger studies to validate these preliminary findings and to assess the predictive value of wall stress in TOF.

Abstract Figure.

Plaque Composition as a Predictor of Plaque Ulceration in Carotid Artery Atherosclerosis: The Plaque At RISK Study

BACKGROUND AND PURPOSE

Plaque ulceration is a marker of previous plaque rupture. We studied the association between atherosclerotic plaque composition at baseline and plaque ulceration at baseline and follow-up.

MATERIALS AND METHODS

We included symptomatic patients with a carotid stenosis of <70% who underwent MDCTA and MR imaging at baseline (n = 180). MDCTA was repeated at 2 years (n = 73). We assessed the presence of ulceration using MDCTA. Baseline MR imaging was used to assess the vessel wall volume and the presence and volume of plaque components (intraplaque hemorrhage, lipid-rich necrotic core, and calcifications) and the fibrous cap status. Associations at baseline were evaluated with binary logistic regression and reported with an OR and its 95% CI. Simple statistical testing was performed in the follow-up analysis.

RESULTS

At baseline, the prevalence of plaque ulceration was 27% (49/180). Increased wall volume (OR  = 12.1; 95% CI, 3.5-42.0), higher relative lipid-rich necrotic core (OR = 1.7; 95% CI, 1.3-2.2), higher relative intraplaque hemorrhage volume (OR = 1.7; 95% CI, 1.3-2.2), and a thin-or-ruptured fibrous cap (OR = 3.4; 95% CI, 1.7-6.7) were associated with the presence of ulcerations at baseline. In 8% (6/73) of the patients, a new ulcer developed. Plaques with a new ulceration at follow-up had at baseline a larger wall volume (1.04 cm³ [IQR, 0.97-1.16 cm³] versus 0.86 cm³ [IQR, 0.73-1.00 cm³]; P = .029), a larger relative lipid-rich necrotic core volume (23% [IQR, 13-31%] versus 2% [IQR, 0-14%]; P = .002), and a larger relative intraplaque hemorrhage volume (14% [IQR, 8-24%] versus 0% [IQR, 0-5%]; P < .001).

CONCLUSIONS

Large atherosclerotic plaques and plaques with intraplaque hemorrhage and lipid-rich necrotic cores were associated with plaque ulcerations at baseline and follow-up.

The Story of Wall Shear Stress in Coronary Artery Atherosclerosis: Biochemical Transport and Mechanotransduction

Coronary artery atherosclerosis is a local, multifactorial, complex disease, and the leading cause of death in the US. Complex interactions between biochemical transport and biomechanical forces influence disease growth. Wall shear stress (WSS) affects coronary artery atherosclerosis by inducing endothelial cell mechanotransduction and by controlling the near-wall transport processes involved in atherosclerosis. Each of these processes is controlled by WSS differently and therefore has complicated the interpretation of WSS in atherosclerosis. In this paper, we present a comprehensive theory for WSS in atherosclerosis. First, a short review of shear stress-mediated mechanotransduction in atherosclerosis was presented. Next, subject-specific computational fluid dynamics (CFD) simulations were performed in ten coronary artery models of diseased and healthy subjects. Biochemical-specific mass transport models were developed to study low-density lipoprotein, nitric oxide, adenosine triphosphate, oxygen, monocyte chemoattractant protein-1, and monocyte transport. The transport results were compared with WSS vectors and WSS Lagrangian coherent structures (WSS LCS). High WSS magnitude protected against atherosclerosis by increasing the production or flux of atheroprotective biochemicals and decreasing the near-wall localization of atherogenic biochemicals. Low WSS magnitude promoted atherosclerosis by increasing atherogenic biochemical localization. Finally, the attracting WSS LCS's role was more complex where it promoted or prevented atherosclerosis based on different biochemicals. We present a summary of the different pathways by which WSS influences coronary artery atherosclerosis and compare different mechanotransduction and biotransport mechanisms.

Association between Intraplaque Hemorrhage and Vascular Remodeling in Carotid Arteries: The Plaque at RISK (PARISK) Study

INTRODUCTION

Vascular remodeling is a compensatory enlargement of the vessel wall in response to atherosclerotic plaque growth. We aimed to investigate the association between intraplaque hemorrhage (IPH), vascular remodeling, and luminal dimensions in recently symptomatic patients with mild to moderate carotid artery stenosis in which the differences in plaque size were taken into account.

MATERIALS AND METHODS

We assessed vessel dimensions on MRI of the symptomatic carotid artery in 164 patients from the Plaque At RISK study. This study included patients with recent ischemic neurological event and ipsilateral carotid artery stenosis <70%. The cross section with the largest wall area (WA) in the internal carotid artery (ICA) was selected for analysis. On this cross section, the following parameters were determined: WA, total vessel area (TVA), and lumen area (LA). Vascular remodeling was quantified as the remodeling ratio (RR) and was calculated as TVA at this position divided by the TVA in an unaffected distal portion of the ipsilateral ICA. Adjustment for WA was performed to correct for plaque size.

RESULTS

Plaques with IPH had a larger WA (0.56 vs. 0.46 cm2; p < 0.001), a smaller LA (0.17 vs. 0.22 cm2; p = 0.03), and a higher RR (2.0 vs. 1.9; p = 0.03) than plaques without IPH. After adjustment for WA, plaques containing IPH had a smaller LA (B = -0.052, p = 0.01) than plaques without IPH, but the RR was not different.

CONCLUSION

After correcting for plaque size, plaques containing IPH had a smaller LA than plaques without IPH. However, RR was not different.

Physiology and coronary artery disease: emerging insights from computed tomography imaging based computational modeling

Improvements in spatial and temporal resolution now permit robust high quality characterization of presence, morphology and composition of coronary atherosclerosis in computed tomography (CT). These characteristics include high risk features such as large plaque volume, low CT attenuation, napkin-ring sign, spotty calcification and positive remodeling. Because of the high image quality, principles of patient-specific computational fluid dynamics modeling of blood flow through the coronary arteries can now be applied to CT and allow the calculation of local lesion-specific hemodynamics such as endothelial shear stress, fractional flow reserve and axial plaque stress. This review examines recent advances in coronary CT image-based computational modeling and discusses the opportunity to identify lesions at risk for rupture much earlier than today through the combination of anatomic and hemodynamic information.

Contemporary rationale for non-invasive imaging of adverse coronary plaque features to identify the vulnerable patient: a Position Paper from the European Society of Cardiology Working Group on Atherosclerosis and Vascular Biology and the European Association of Cardiovascular Imaging

Atherosclerotic plaques prone to rupture may cause acute myocardial infarction (MI) but can also heal without causing an event. Certain common histopathological features, including inflammation, a thin fibrous cap, positive remodelling, a large necrotic core, microcalcification, and plaque haemorrhage are commonly found in plaques causing an acute event. Recent advances in imaging techniques have made it possible to detect not only luminal stenosis and overall coronary atherosclerosis burden but also to identify such adverse plaque characteristics. However, the predictive value of identifying individual adverse atherosclerotic plaques for future events has remained poor. In this Position Paper, the relationship between vulnerable plaque imaging and MI is addressed, mainly for non-invasive assessments but also for invasive imaging of adverse plaques in patients undergoing invasive coronary angiography. Dynamic changes in atherosclerotic plaque development and composition may indicate that an adverse plaque phenotype should be considered at the patient level rather than for individual plaques. Imaging of adverse plaque burden throughout the coronary vascular tree, in combination with biomarkers and biomechanical parameters, therefore holds promise for identifying subjects at increased risk of MI and for guiding medical and invasive treatment.

Lipid-rich Plaques Detected by Near-infrared Spectroscopy Are More Frequently Exposed to High Shear Stress

High wall shear stress (WSS) and near-infrared spectroscopy (NIRS) detected lipid-rich plaque (LRP) are both known to be associated with plaque destabilization and future adverse cardiovascular events. However, knowledge of spatial co-localization of LRP and high WSS is lacking. This study investigated the co-localization of LRP based on NIRS and high WSS. Fifty-three patients presenting acute coronary syndrome underwent NIRS-intravascular-ultrasound (NIRS-IVUS) imaging of a non-culprit coronary artery. WSS was obtained using WSS profiling in 3D-reconstructions of the coronary arteries based on fusion of IVUS-segmented lumen and CT-derived 3D-centerline. Thirty-eight vessels were available for final analysis and divided into 0.5 mm/45° sectors. LRP sectors, as identified by NIRS, were more often colocalized with high WSS than sectors without LRP. Moreover, there was a dose-dependent relationship between lipid content and high WSS exposure. This study is a first step in understanding the evolution of LRPs to vulnerable plaques.

Automated Quantitative Assessment of Coronary Calcification Using Intravascular Ultrasound

Coronary calcification represents a challenge in the treatment of coronary artery disease by stent placement. It negatively affects stent expansion and has been related to future adverse cardiac events. Intravascular ultrasound (IVUS) is known for its high sensitivity in detecting coronary calcification. At present, automated quantification of calcium as detected by IVUS is not available. For this reason, we developed and validated an optimized framework for accurate automated detection and quantification of calcified plaque in coronary atherosclerosis as seen by IVUS. Calcified lesions were detected by training a supported vector classifier per IVUS A-line on manually annotated IVUS images, followed by post-processing using regional information. We applied our framework to 35 IVUS pullbacks from each of the three commonly used IVUS systems. Cross-validation accuracy for each system was >0.9, and the testing accuracy was 0.87, 0.89 and 0.89 for the three systems. Using the detection result, we propose an IVUS calcium score, based on the fraction of calcium-positive A-lines in a pullback segment, to quantify the extent of calcified plaque. The high accuracy of the proposed classifier suggests that it may provide a robust and accurate tool to assess the presence and amount of coronary calcification and, thus, may play a role in image-guided coronary interventions.

The clinical impact of phase offset errors and different correction methods in cardiovascular magnetic resonance phase contrast imaging: a multi-scanner study

BACKGROUND

Cardiovascular magnetic resonance (CMR) phase contrast (PC) flow measurements suffer from phase offset errors. Background subtraction based on stationary phantom measurements can most reliably be used to overcome this inaccuracy. Stationary tissue correction is an alternative and does not require additional phantom scanning. The aim of this study was 1) to compare measurements with and without stationary tissue correction to phantom corrected measurements on different GE Healthcare CMR scanners using different software packages and 2) to evaluate the clinical implications of these methods.

METHODS

CMR PC imaging of both the aortic and pulmonary artery flow was performed in patients on three different 1.5 T CMR scanners (GE Healthcare) using identical scan parameters. Uncorrected, first, second and third order stationary tissue corrected flow measurement were compared to phantom corrected flow measurements, our reference method, using Medis QFlow, Circle cvi42 and MASS software. The optimal (optimized) stationary tissue order was determined per scanner and software program. Velocity offsets, net flow, clinically significant difference (deviation > 10% net flow), and regurgitation severity were assessed.

RESULTS

Data from 175 patients (28 (17-38) years) were included, of which 84% had congenital heart disease. First, second and third order and optimized stationary tissue correction did not improve the velocity offsets and net flow measurements. Uncorrected measurements resulted in the least clinically significant differences in net flow compared to phantom corrected data. Optimized stationary tissue correction per scanner and software program resulted in net flow differences (> 10%) in 19% (MASS) and 30% (Circle cvi42) of all measurements compared to 18% (MASS) and 23% (Circle cvi42) with no correction. Compared to phantom correction, regurgitation reclassification was the least common using uncorrected data. One CMR scanner performed worse and significant net flow differences of > 10% were present both with and without stationary tissue correction in more than 30% of all measurements.

CONCLUSION

Phase offset errors had a significant impact on net flow quantification, regurgitation assessment and varied greatly between CMR scanners. Background phase correction using stationary tissue correction worsened accuracy compared to no correction on three GE Healthcare CMR scanners. Therefore, careful assessment of phase offset errors at each individual scanner is essential to determine whether routine use of phantom correction is necessary.

TRIAL REGISTRATION

Observational Study.

Exploring wall shear stress spatiotemporal heterogeneity in coronary arteries combining correlation-based analysis and complex networks with computational hemodynamics

Atherosclerosis at the early stage in coronary arteries has been associated with low cycle-average wall shear stress magnitude. However, parallel to the identification of an established active role for low wall shear stress in the onset/progression of the atherosclerotic disease, a weak association between lesions localization and low/oscillatory wall shear stress has been observed. In the attempt to fully identify the wall shear stress phenotype triggering early atherosclerosis in coronary arteries, this exploratory study aims at enriching the characterization of wall shear stress emerging features combining correlation-based analysis and complex networks theory with computational hemodynamics. The final goal is the characterization of the spatiotemporal and topological heterogeneity of wall shear stress waveforms along the cardiac cycle. In detail, here time-histories of wall shear stress magnitude and wall shear stress projection along the main flow direction and orthogonal to it (a measure of wall shear stress multidirectionality) are analyzed in a representative dataset of 10 left anterior descending pig coronary artery computational hemodynamics models. Among the main findings, we report that the proposed analysis quantitatively demonstrates that the model-specific inlet flow-rate shapes wall shear stress time-histories. Moreover, it emerges that a combined effect of low wall shear stress magnitude and of the shape of the wall shear stress–based descriptors time-histories could trigger atherosclerosis at its earliest stage. The findings of this work suggest for new experiments to provide a clearer determination of the wall shear stress phenotype which is at the basis of the so-called arterial hemodynamic risk hypothesis in coronary arteries.

In-vitro and in-vivo imaging of coronary artery stents with Heartbeat OCT

To quantify the impact of cardiac motion on stent length measurements with Optical Coherence Tomography (OCT) and to demonstrate in vivo OCT imaging of implanted stents, without motion artefacts. The study consists of: clinical data evaluation, simulations and in vivo tests. A comparison between OCT-measured and nominal stent lengths in 101 clinically acquired pullbacks was carried out, followed by a simulation of the effect of cardiac motion on stent length measurements, experimentally and computationally. Both a commercial system and a custom OCT, capable of completing a pullback between two consecutive ventricular contractions, were employed. A 13 mm long stent was implanted in the left anterior descending branch of two atherosclerotic swine and imaged with both OCT systems. The analysis of the clinical OCT images yielded an average difference of 1.1 ± 1.6 mm, with a maximum difference of 7.8 mm and the simulations replicated the statistics observed in clinical data. Imaging with the custom OCT, yielded an RMS error of 0.14 mm at 60 BPM with the start of the acquisition synchronized to the cardiac cycle. In vivo imaging with conventional OCT yielded a deviation of 1.2 mm, relative to the length measured on ex-vivo micro-CT, while the length measured in the pullback acquired by the custom OCT differed by 0.20 mm. We demonstrated motion artefact-free OCT-imaging of implanted stents, using ECG triggering and a rapid pullback.

P5280The impact of background phase offset errors in cardiovascular magnetic resonance phase contrast imaging: a multi-scanner study

Background

Phase contrast (PC) CMR flow measurements (FM) are widely used for blood flow assessment, but they suffer from phase offset errors (POE). Stationary phantom correction limits these inaccuracies, however, this adds scan time. Stationary tissue (ST) correction is an alternative method that does not require additional scanning. The aim of this study was to evaluate the impact of POE, to assess interscanner variation, and to evaluate the ST correction usage.

Methods

We included 166 patients in which both aorta and main pulmonary artery FM were acquired including static gelatin phantom data. Subjects were scanned on three types of 1.5T scanners from the one vendor. Uncorrected and ST corrected FM were compared with phantom corrected FM, our reference value, and corrected for BSA. A difference of >10% in net flow was defined as clinically relevant. Regurgitation fraction was calculated and POE influences were assessed. Regurgitation severity was graded and POE influence on severity grading was assessed.

Results

Of the 166 cases included, the median age was 27 (5–74) years. Overall, the median difference between no corrected and phantom corrected FM was ≤6%, however, with a wide range of over- and underestimation (−155%–78% change) (figure). ST correction resulted in larger differences compared to no correction (p<0.01). Clinically significant differences were seen in 19% of all FM with no correction and in 30% of with ST correction (p<0.01). Furthermore, there were significant differences between scanners (no correction 10%, p<0.01; ST correction, p<0.01). Regurgitation severity indexing changed in 38 (11%) cases with no correction and in 48 (48%) with ST correction.

Magnitude of flow change with and without offset corrections (n=332) Flow (ml/m2) Δ no correction and phantom correction (%) Δ ST correction and phantom correction (%) Clinically significant difference (>10%) Mean ± SD Median IQR Range Median IQR Range No correction, N (%) ST correction , N (%) MRI 1 (n=126) 50±12 3 0 to 6 −8 to 30 5 −3 to 9 −26 to 28 13 (10%) 34 (27%) MRI 2 (n=102) 48±13 −2 −15 to 6 −155 to 78 5 −3 to 11 −74 to 52 50 (49%) 50 (49%) MRI 3 (n=104) 48±12 −1 −1 to 0 −7 to 14 2 −2 to 5 −39 to 29 1 (1%) 16 (15%) Total (n=332) 49±12 0 −2 to 4 −155 to 78 3 −2 to 8 −74 to 52 64 (19%) 100 (30%)

Conclusion

Background POE have a significant impact on flow quantification and regurgitation severity. Unexpectedly, background correction using ST correction worsens accuracy compared to no correction. POE vary greatly between scanners. Therefore, careful assessment of FM at each scanner is essential to determine if routine phantom scanning is necessary.

5222A distinct LDL profile to predict the risk of cardiovascular disease in familial hypercholesterolemia subjects: initial pre-clinical results

Aim

Cardiovascular disease development is highly variable between patients with familial hypercholesterolemia (FH). Since current risk prediction methods fail to estimate the risk of individual patients, all patients are treated with high-intensity lipid-lowering medication like statins and PCSK9-inhibitors. To prevent overtreatment (with its associated costs and side effects) of patients with a low risk, reliable biomarkers are urgently needed. While studying atherosclerosis development in an FH porcine model, we discovered a specific LDL profile that was directly associated with the severity of atherosclerosis development.

Methods

10 adult pigs with an LDLR mutation received a high-fat diet for 12 months. Atherosclerosis development in the three main coronary arteries was monitored with intravascular ultrasound (IVUS) and optical coherence tomography at three time points. After sacrifice, tissue was harvested for detailed histological analysis of the coronary plaque composition. Size-exclusion chromatography (SEC, n=10) and LC-MSMS (n=4) were used to assess the plasma lipoprotein profile, and the sphingolipid content of LDL, relative to cholesterol.

Results

Imaging and histology revealed a marked difference in pigs that developed large, lumen intruding plaques (IVUS-based plaque burden 13–77%) within 9 months (n=5) and pigs (n=5) that only developed early lesions (IVUS-based plaque burden 8–34%), even after 12 months of follow-up. The plaques seen in the fast responding pigs possessed distinct characteristics of advanced disease (i.e. heavy macrophage infiltration, large lipid-rich necrotic cores, neovascularisation, micro- and macrocalcifications and intraplaque haemorrhage, see figure). In these fast responders, fibrous cap atheroma occupied 34% of the total artery.

SEC revealed two distinct LDL subclasses: regular and “larger” LDL particles. Fast responding pigs with advanced atherosclerosis displayed a significantly higher ratio in cholesterol concentration of regular/'larger' LDL than slow responding pigs (1.7 (1.3–1.9) vs. 0.8 (0.6–1.2); (p=0.004), see figure). Compared to regular LDL, “larger” LDL contained relatively more sphingolipids in the fast responding than in the slow responding animals (regular LDL/'larger' LDL: S1P 0.5 (0.5–0.5) vs. 1.0 (0.8–1.2); Cer16:0 0.70 (0.67–0.73) vs. 1.04 (0.95–1.13); Cer18:0 0.60 (0.58–0.61) vs. 1.15 (1.13–1.16); Cer20:0 0.73 (0.73–0.74) vs. 0.94 (0.94–0.94)). “Larger” LDL particles and comparable sphingolipid ratios were also observed in FH patients. Cardiovascular data from our FH patient cohort, coupled to the LDL subclass distribution, will provide more insight into the potential of this novel biomarker.

Conclusion

A distinct difference in LDL subclasses, including a new “larger” LDL particle, was found in fast versus slow responding FH animals. This finding can potentially be used to identify FH patients at the highest risk of CVD to avoid overtreatement of low risk patients.

Acknowledgement/Funding

ERC- starter grant (grant agreement 310457)

P3588The synergistic effect of NIRS-detected lipid-rich plaque and 5 different multidirectional wall shear stress metrics on human coronary plaque growth

Introduction

Local wall shear stress (WSS) metrics, high local lipid levels (as detected by near-infrared spectroscopy (NIRS)), as well as systemic lipid levels, have been individually associated with atherosclerotic disease progression. However, a possible synergistic effect remains to be elucidated. This study is the first study to combine WSS metrics with NIRS-detected local lipid content to investigate a potential synergistic effect on plaque progression in human coronary arteries.

Methods

The IMPACT study is a prospective, single centre study investigating the relation between atherosclerotic plaque progression and WSS in human coronary arteries. Patients with ACS treated with PCI were included. At baseline and after 1-year follow-up, patients underwent near-infrared spectroscopy intravascular ultrasound (NIRS-IVUS) imaging and intravascular doppler flow measurements of at least one non-culprit coronary artery. After one month, a CT angiography was made. CT derived centreline combined with IVUS lumen contours resulted in a 3D reconstruction of the vessel. The following WSS metrics were computed using computational fluid dynamics applying the vessel specific invasive flow measurements: time-average wall shear stress (TAWSS), relative residence time (RRT), cross-flow index, oscillatory shear index and transverse wall shear stress. Low TAWSS is known as pro atherogenic, in contrast to all the other shear stress metrics, at which a high magnitude is pro-atherogenic. The arteries were divided into 1.5mm/45° sectors. Based on NIRS-IVUS, wall thickness change over time was determined and NIRS positive sectors detected. Furthermore, per vessel the shear stress was divided into tertiles (low, intermediate, high). To investigate the synergistic effect of local lipids on shear stress related plaque growth, wall thickness change over time was related to the different shear stress metrics comparing the NIRS-positive with the NIRS-negative sectors.

Results

15 non-culprit coronary arteries from the first 14 patients were analyzed (age 62±10 years old and 92.9% male). A total of 2219 sectors were studied (5.2%, N=130, NIRS-positive) for wall thickness changes. After studying all five shear stress metrics, we found for TAWSS and RRT that presence of lipids, as detected by NIRS, amplified the effect of shear stress on plaque progression (see figure). Sectors presenting with lipid-rich plaque, compared to NIRS-negative sectors, showed more progression when they were exposed to low TAWSS (p=0.07) or high RRT (p=0.012) and more regression in sectors exposed to high TAWSS (p=0.10) or low RRT (p=0.06).

Delta wall thickness vs shear stress

Conclusion

We presented the first preliminary results of the IMPACT study, showing the synergistic effect of lipid rich plaque and shear stress on plaque progression. Therefore, intravascular lipid-rich plaque (NIRS) assessment has added value to shear stress profiling for the prediction of plaque growth, leading to improved risk stratification.

Acknowledgement/Funding

ERC starting grant 310457

Variation in Coronary Atherosclerosis Severity Related to a Distinct LDL (Low-Density Lipoprotein) Profile: Findings From a Familial Hypercholesterolemia Pig Model

OBJECTIVE

In an adult porcine model of familial hypercholesterolemia (FH), coronary plaque development was characterized. To elucidate the underlying mechanisms of the observed inter-individual variation in disease severity, detailed lipoprotein profiles were determined. Approach and Results: FH pigs (3 years old, homozygous LDLR R84C mutation) received an atherogenic diet for 12 months. Coronary atherosclerosis development was monitored using serial invasive imaging and histology. A pronounced difference was observed between mildly diseased pigs which exclusively developed early lesions (maximal plaque burden, 25% [23%-34%]; n=5) and advanced-diseased pigs (n=5) which developed human-like, lumen intruding plaques (maximal plaque burden, 69% [57%-77%]) with large necrotic cores, intraplaque hemorrhage, and calcifications. Advanced-diseased pigs and mildly diseased pigs displayed no differences in conventional risk factors. Additional plasma lipoprotein profiling by size-exclusion chromatography revealed 2 different LDL (low-density lipoprotein) subtypes: regular and larger LDL. Cholesterol, sphingosine-1-phosphate, ceramide, and sphingomyelin levels were determined in these LDL-subfractions using standard laboratory techniques and high-pressure liquid chromatography mass-spectrometry analyses, respectively. At 3 months of diet, regular LDL of advanced-diseased pigs contained relatively more cholesterol (LDL-C; regular/larger LDL-C ratio 1.7 [1.3-1.9] versus 0.8 [0.6-0.9]; P=0.008) than mildly diseased pigs, while larger LDL contained more sphingosine-1-phosphate, ceramides, and sphingomyelins. Larger and regular LDL was also found in plasma of 3 patients with homozygous FH with varying LDL-C ratios.

CONCLUSIONS

In our adult FH pig model, inter-individual differences in atherosclerotic disease severity were directly related to the distribution of cholesterol and sphingolipids over a distinct LDL profile with regular and larger LDL shortly after the diet start. A similar LDL profile was detected in patients with homozygous FH.

In vivo intravascular photoacoustic imaging of plaque lipid in coronary atherosclerosis

Prospective identification of lipid-rich vulnerable plaque has remained an elusive goal. Intravascular photoacoustics, a hybrid optical and ultrasonic technology, was developed as a tool for lipid-rich plaque imaging. Here, we present the first in vivo images of lipid-rich coronary atherosclerosis acquired with this new technology in a large animal model, and relate them to independent catheter-based imaging and histology.

Flow Patterns in Carotid Webs: A Patient-Based Computational Fluid Dynamics Study

BACKGROUND AND PURPOSE

Carotid webs are increasingly recognized as an important cause of (recurrent) ischemic stroke in patients without other cardiovascular risk factors. Hemodynamic flow patterns induced by these lesions might be associated with thrombus formation. The aim of our study was to evaluate flow patterns of carotid webs using computational fluid dynamics.

MATERIALS AND METHODS

Patients with a carotid web in the Multicenter Randomized Clinical Trial of Endovascular Treatment of Acute Ischemic Stroke in the Netherlands (MR CLEAN) were selected for hemodynamic evaluation with computational fluid dynamics models based on lumen segmentations obtained from CT angiography scans. Hemodynamic parameters, including the area of recirculation zone, time-averaged wall shear stress, transverse wall shear stress, and the oscillatory shear index, were assessed and compared with the contralateral carotid bifurcation.

RESULTS

In our study, 9 patients were evaluated. Distal to the carotid webs, recirculation zones were significantly larger compared with the contralateral bifurcation (63 versus 43 mm², P = .02). In the recirculation zones of the carotid webs and the contralateral carotid bifurcation, time-averaged wall shear stress values were comparable (both: median, 0.27 Pa; P = .30), while transverse wall shear stress and oscillatory shear index values were significantly higher in the recirculation zone of carotid webs (median, 0.25 versus 0.21 Pa; P = .02 and 0.39 versus 0.30 Pa; P = .04). At the minimal lumen area, simulations showed a significantly higher time-averaged wall shear stress in the web compared with the contralateral bifurcation (median, 0.58 versus 0.45 Pa; P = .01).

CONCLUSIONS

Carotid webs are associated with increased recirculation zones and regional increased wall shear stress metrics that are associated with disturbed flow. These findings suggest that a carotid web might stimulate thrombus formation, which increases the risk of acute ischemic stroke.

Spatial correlations between MRI-derived wall shear stress and vessel wall thickness in the carotid bifurcation

Background

To explore the possibility of creating three-dimensional (3D) estimation models for patient-specific wall thickness (WT) maps using patient-specific and cohort-averaged WT, wall shear stress (WSS), and vessel diameter maps in asymptomatic atherosclerotic carotid bifurcations.

Methods

Twenty subjects (aged 75 ± 6 years [mean ± standard deviation], eight women) underwent a 1.5-T MRI examination. Non-gated 3D phase-contrast gradient-echo images and proton density-weighted echo-planar images were retrospectively assessed for WSS, diameter estimation, and WT measurements. Spearman’s ρ and scatter plots were used to determine correlations between individual WT, WSS, and diameter maps. A bootstrapping technique was used to determine correlations between 3D cohort-averaged WT, WSS, and diameter maps. Linear regression between the cohort-averaged WT, WSS, and diameter maps was used to predict individual 3D WT.

Results

Spearman’s ρ averaged over the subjects was − 0.24 ± 0.18 (p < 0.001) and 0.07 ± 0.28 (p = 0.413) for WT versus WSS and for WT versus diameter relations, respectively. Cohort-averaged ρ, averaged over 1000 bootstraps, was − 0.56 (95% confidence interval [− 0.74,− 0.38]) for WT versus WSS and 0.23 (95% confidence interval [− 0.06, 0.52]) for WT versus diameter. Scatter plots did not reveal relationships between individual WT and WSS or between WT and diameter data. Linear relationships between these parameters became apparent after averaging over the cohort. Spearman’s ρ between the original and predicted WT maps was 0.21 ± 0.22 (p < 0.001).

Conclusions

With a combination of bootstrapping and cohort-averaging methods, 3D WT maps can be predicted from the individual 3D WSS and diameter maps. The methodology may help to elucidate pathological processes involving WSS in carotid atherosclerosis.

Electronic supplementary material

The online version of this article (10.1186/s41747-018-0058-1) contains supplementary material, which is available to authorized users.

Future directions for therapeutic strategies in post-ischaemic vascularization: a position paper from European Society of Cardiology Working Group on Atherosclerosis and Vascular Biology

Modulation of vessel growth holds great promise for treatment of cardiovascular disease. Strategies to promote vascularization can potentially restore function in ischaemic tissues. On the other hand, plaque neovascularization has been shown to associate with vulnerable plaque phenotypes and adverse events. The current lack of clinical success in regulating vascularization illustrates the complexity of the vascularization process, which involves a delicate balance between pro- and anti-angiogenic regulators and effectors. This is compounded by limitations in the models used to study vascularization that do not reflect the eventual clinical target population. Nevertheless, there is a large body of evidence that validate the importance of angiogenesis as a therapeutic concept. The overall aim of this Position Paper of the ESC Working Group of Atherosclerosis and Vascular biology is to provide guidance for the next steps to be taken from pre-clinical studies on vascularization towards clinical application. To this end, the current state of knowledge in terms of therapeutic strategies for targeting vascularization in post-ischaemic disease is reviewed and discussed. A consensus statement is provided on how to optimize vascularization studies for the identification of suitable targets, the use of animal models of disease, and the analysis of novel delivery methods.

Atherosclerosis at arterial bifurcations: evidence for the role of haemodynamics and geometry

Atherosclerotic plaques are found at distinct locations in the arterial system, despite the exposure to systemic risk factors of the entire vascular tree. From the study of arterial bifurcation regions, emerges ample evidence that haemodynamics are involved in the local onset and progression of the atherosclerotic disease. This observed co-localisation of disturbed flow regions and lesion prevalence at geometrically predisposed districts such as arterial bifurcations has led to the formulation of a ‘haemodynamic hypothesis’, that in this review is grounded to the most current research concerning localising factors of vascular disease. In particular, this review focuses on carotid and coronary bifurcations because of their primary relevance to stroke and heart attack. We highlight reported relationships between atherosclerotic plaque location, progression and composition, and fluid forces at vessel’s wall, in particular shear stress and its ‘easier-tomeasure’ surrogates, i.e. vascular geometric attributes (because geometry shapes the flow) and intravascular flow features (because they mediate disturbed shear stress), in order to give more insight in plaque initiation and destabilisation. Analogous to Virchow’s triad for thrombosis, atherosclerosis must be thought of as subject to a triad of, and especially interactions among, haemodynamic forces, systemic risk factors, and the biological response of the wall.