3D fusion of intravascular ultrasound and coronary computed tomography for in-vivo wall shear stress analysis: a feasibility study

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 time-averaged WSS magnitude (TAWSS) and the variability in the contraction/expansion action exerted by WSS on the endothelium, quantifiable in terms of topological shear variation index (TSVI). Plaque progression was measured as intravascular ultrasound-derived percentage plaque atheroma volume change at 1-year follow-up. 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 plaque 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.

Advances in Intravascular Ultrasound

Since its inception, intravascular ultrasound (IVUS) and optical coherence tomography (OCT) have played a significant role in evaluating the pathophysiology of coronary artery disease (CAD) guiding the interventional and medical management of CAD improving outcomes in patients. Although the benefits of each of these modalities have been proven, due to some limitations, no single intravascular imaging technique has been proven to provide a detailed and complete evaluation of all CAD lesions. The use of different intravascular imaging modalities sequentially may lead to complications, which are cumbersome, consume time, and add financial burden to the patient. Recently, hybrid imaging catheters that combine OCT and IVUS benefits have been developed to limit these problems. Intravascular imaging techniques we are using currently have some drawbacks that hinder accurate assessment of plaque morphology and pathobiology as demonstrated in many histological studies, causing difficulty in identifying high-risk plaques. To overcome these limitations, great efforts have been put into developing hybrid, dual-probe catheters by combining imaging modalities to get an accurate analysis of plaque characteristics, and high-risk lesions. At present, many dual-probe catheters are available including combined IVUS-OCT, near-infrared spectroscopy-IVUS that is available commercially, the OCT-near infrared fluorescence (NIRF) molecular imaging, IVUS-NIRF, and combined fluorescence lifetime-IVUS imaging. Application of this combined multimodal imaging in clinical practice overcomes the limitations of standalone imaging and helps in providing a comprehensive and accurate visualization of plaque characteristics, composition, and plaque biology. The present article summarizes the advances in hybrid intravascular imaging, analyses the technical hindrances that should be known to have a use in the different clinical circumstances, and the till date shreds of evidence available from their first clinical application aiming to bring these modalities into the limelight and their potential role in the study of CAD.

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.

The role of cardiac computed tomography in predicting adverse coronary events

Cardiac computed tomography (CCT) is now considered a first-line diagnostic test for suspected coronary artery disease (CAD) providing a non-invasive, qualitative, and quantitative assessment of the coronary arteries and pericoronary regions. CCT assesses vascular calcification and coronary lumen narrowing, measures total plaque burden, identifies plaque composition and high-risk plaque features and can even assist with hemodynamic evaluation of coronary lesions. Recent research focuses on computing coronary endothelial shear stress, a potent modulator in the development and progression of atherosclerosis, as well as differentiating an inflammatory from a non-inflammatory pericoronary artery environment using the simple measurement of pericoronary fat attenuation index. In the present review, we discuss the role of the above in the diagnosis of coronary atherosclerosis and the prediction of adverse cardiovascular events. Additionally, we review the current limitations of cardiac computed tomography as an imaging modality and highlight how rapid technological advancements can boost its capacity in predicting cardiovascular risk and guiding clinical decision-making.

Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases

Hemodynamic factors, induced by pulsatile blood flow, play a crucial role in vascular health and diseases, such as the initiation and progression of atherosclerosis. Computational fluid dynamics, finite element analysis, and fluid-structure interaction simulations have been widely used to quantify detailed hemodynamic forces based on vascular images commonly obtained from computed tomography angiography, magnetic resonance imaging, ultrasound, and optical coherence tomography. In this review, we focus on methods for obtaining accurate hemodynamic factors that regulate the structure and function of vascular endothelial and smooth muscle cells. We describe the multiple steps and recent advances in a typical patient-specific simulation pipeline, including medical imaging, image processing, spatial discretization to generate computational mesh, setting up boundary conditions and solver parameters, visualization and extraction of hemodynamic factors, and statistical analysis. These steps have not been standardized and thus have unavoidable uncertainties that should be thoroughly evaluated. We also discuss the recent development of combining patient-specific models with machine-learning methods to obtain hemodynamic factors faster and cheaper than conventional methods. These critical advances widen the use of biomechanical simulation tools in the research and potential personalized care of vascular diseases.

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.

Study on the association of wall shear stress and vessel structural stress with atherosclerosis: An experimental animal study

BACKGROUND AND AIMS

Artery is subject to wall shear stress (WSS) and vessel structural stress (VSS) simultaneously. This study is designed to explore the role of VSS in development of atherosclerosis.

METHODS

Silastic collars were deployed on the carotid to create two constrictions on 13 rabbits for a distinct mechanical environment at the constriction. MRI was performed to visualize arteries' configuration. Animals with high fat (n = 9; Model-group) and normal diet (n = 4; Control-group) were sacrificed after 16 weeks. 3D fluid-structure interaction analysis was performed to quantify WSS and VSS simultaneously.

RESULTS

Twenty plaques were found in Model-group and 3 in Control-group. In Model-group, 8 plaques located proximally to the first constriction (Region-1, close to the heart) and 7 distally to the second (Region-2, close to the head) and 5 plaques were found on the contralateral side of 3 rabbits. Plaques at Region-1 tended to be bigger than those at Region-2 and the macrophage density at these locations was comparable. Minimum time-averaged WSS (TAWSS) in Region-1 was significantly higher than that in Region-2, and both maximum oscillatory shear index (OSI) and particle relative residence time (RRT) were significantly lower. Peak and mean VSS in Region-1 were significantly higher than those in Region-2. Correlation analyses indicated that low TAWSS, high OSI and RRT were only associated with plaque in Region-2, while lesions in Region-1 were only associated with high VSS. Moreover, only VSS was associated with wall thickness of plaque-free regions in both regions.

CONCLUSIONS

VSS might contribute to the initialization and development of atherosclerosis solely or in combination with WSS.

3D reconstruction of coronary artery bifurcations from coronary angiography and optical coherence tomography: feasibility, validation, and reproducibility

The three-dimensional (3D) representation of the bifurcation anatomy and disease burden is essential for better understanding of the anatomical complexity of bifurcation disease and planning of stenting strategies. We propose a novel methodology for 3D reconstruction of coronary artery bifurcations based on the integration of angiography, which provides the backbone of the bifurcation, with optical coherence tomography (OCT), which provides the vessel shape. Our methodology introduces several technical novelties to tackle the OCT frame misalignment, correct positioning of the OCT frames at the carina, lumen surface reconstruction, and merging of bifurcation lumens. The accuracy and reproducibility of the methodology were tested in n = 5 patient-specific silicone bifurcations compared to contrast-enhanced micro-computed tomography (µCT), which was used as reference. The feasibility and time-efficiency of the method were explored in n = 7 diseased patient bifurcations of varying anatomical complexity. The OCT-based reconstructed bifurcation models were found to have remarkably high agreement compared to the µCT reference models, yielding r² values between 0.91 and 0.98 for the normalized lumen areas, and mean differences of 0.005 for lumen shape and 0.004 degrees for bifurcation angles. Likewise, the reproducibility of our methodology was remarkably high. Our methodology successfully reconstructed all the patient bifurcations yielding favorable processing times (average lumen reconstruction time < 60 min). Overall, our method is an easily applicable, time-efficient, and user-friendly tool that allows accurate and reproducible 3D reconstruction of coronary bifurcations. Our technique can be used in the clinical setting to provide information about the bifurcation anatomy and plaque burden, thereby enabling planning, education, and decision making on bifurcation stenting.

Validation of Wall Shear Stress Assessment in Non-invasive Coronary CTA versus Invasive Imaging: A Patient-Specific Computational Study

Endothelial shear stress (ESS) identifies coronary plaques at high risk for progression and/or rupture leading to a future acute coronary syndrome. In this study an optimized methodology was developed to derive ESS, pressure drop and oscillatory shear index using computational fluid dynamics (CFD) in 3D models of coronary arteries derived from non-invasive coronary computed tomography angiography (CTA). These CTA-based ESS calculations were compared to the ESS calculations using the gold standard with fusion of invasive imaging and CTA. In 14 patients paired patient-specific CFD models based on invasive and non-invasive imaging of the left anterior descending (LAD) coronary arteries were created. Ten patients were used to optimize the methodology, and four patients to test this methodology. Time-averaged ESS (TAESS) was calculated for both coronary models applying patient-specific physiological data available at the time of imaging. For data analysis, each 3D reconstructed coronary artery was divided into 2 mm segments and each segment was subdivided into 8 arcs (45°).TAESS and other hemodynamic parameters were averaged per segment as well as per arc. Furthermore, the paired segment- and arc-averaged TAESS were categorized into patient-specific tertiles (low, medium and high). In the ten LADs, used for optimization of the methodology, we found high correlations between invasively-derived and non-invasively-derived TAESS averaged over segments (n = 263, r = 0.86) as well as arcs (n = 2104, r = 0.85, p < 0.001). The correlation was also strong in the four testing-patients with r = 0.95 (n = 117 segments, p = 0.001) and r = 0.93 (n = 936 arcs, p = 0.001).There was an overall high concordance of 78% of the three TAESS categories comparing both methodologies using the segment- and 76% for the arc-averages in the first ten patients. This concordance was lower in the four testing patients (64 and 64% in segment- and arc-averaged TAESS). Although the correlation and concordance were high for both patient groups, the absolute TAESS values averaged per segment and arc were overestimated using non-invasive vs. invasive imaging [testing patients: TAESS segment: 30.1(17.1-83.8) vs. 15.8(8.8-63.4) and TAESS arc: 29.4(16.2-74.7) vs 15.0(8.9-57.4) p < 0.001]. We showed that our methodology can accurately assess the TAESS distribution non-invasively from CTA and demonstrated a good correlation with TAESS calculated using IVUS/OCT 3D reconstructed models.

The Evolution of Data Fusion Methodologies Developed to Reconstruct Coronary Artery Geometry From Intravascular Imaging and Coronary Angiography Data: A Comprehensive Review

Understanding the mechanisms that regulate atherosclerotic plaque formation and evolution is a crucial step for developing treatment strategies that will prevent plaque progression and reduce cardiovascular events. Advances in signal processing and the miniaturization of medical devices have enabled the design of multimodality intravascular imaging catheters that allow complete and detailed assessment of plaque morphology and biology. However, a significant limitation of these novel imaging catheters is that they provide two-dimensional (2D) visualization of the lumen and vessel wall and thus they cannot portray vessel geometry and 3D lesion architecture. To address this limitation computer-based methodologies and user-friendly software have been developed. These are able to off-line process and fuse intravascular imaging data with X-ray or computed tomography coronary angiography (CTCA) to reconstruct coronary artery anatomy. The aim of this review article is to summarize the evolution in the field of coronary artery modeling; we thus present the first methodologies that were developed to model vessel geometry, highlight the modifications introduced in revised methods to overcome the limitations of the first approaches and discuss the challenges that need to be addressed, so these techniques can have broad application in clinical practice and research.

A comparison of three multimodality coronary 3D reconstruction methods

The assessment of the severity of arterial stenoses is of utmost importance in clinical practice. Several image modalities invasive and non-invasive are nowadays available and can be utilized for the 3-dimensional (3D) reconstruction of the arterial geometry. Following our previous study, the present study was conducted to further strengthen the evaluation of three reconstruction methodologies, namely: (i) the Quantitative Coronary Analysis (QCA), (ii) the Virtual Histology Intravascular Ultrasound VH-IVUS-Angiography hybrid method and (iii) the Coronary Computed Tomography Angiography (CCTA). Data from 13 patients were employed to perform a quantitative analysis using specific metrics, such as, the Mean Wall Shear Stress (mWSS), the Minimum Lumen diameter (MLD), the Reference Vessel Diameter (RVD), the Degree of stenosis (DS%), and the Lesion length (LL). A high correlation was observed for the mWSS metric between the three reconstruction methods, especially between the QCA and CCTA (r=0.974, P<; 0.001).

Biomechanical Forces and Atherosclerosis: From Mechanism to Diagnosis and Treatment

The article provides an overview of current views on the role of biomechanical forces in the pathogenesis of atherosclerosis. The importance of biomechanical forces in maintaining vascular homeostasis is considered. We provide descriptions of mechanosensing and mechanotransduction. The roles of wall shear stress and circumferential wall stress in the initiation, progression and destabilization of atherosclerotic plaque are described. The data on the possibilities of assessing biomechanical factors in clinical practice and the clinical significance of this approach are presented. The article concludes with a discussion on current therapeutic approaches based on the modulation of biomechanical forces.

Ex-vivo study in nephroureterectomy specimens defining the role of 3-D upper urinary tract visualization using optical coherence tomography and endoluminal ultrasound

Minimal invasive endoscopic treatment for upper urinary tract urothelial carcinoma (UUT-UC) is advocated in patients with low-risk disease and limited tumor volume. Diagnostic ureterorenoscopy combined with biopsy is the diagnostic standard. This study aims to evaluate two alternative diagnostic techniques for UUT-UC: optical coherence tomography (OCT) and endoluminal ultrasound (ELUS). Following nephroureterectomy, OCT, ELUS, and computed tomography (CT) were performed of the complete nephroureterectomy specimen. Visualization software (AMIRA^®) was used for reconstruction and coregistration of CT, OCT, and ELUS. Finally, CT was used to obtain exact probe localization. Coregistered OCT and ELUS datasets were compared with histology. Coregistration with three-dimensional CT makes exact data matching possible in this ex-vivo setting to compare histology with OCT and ELUS. In OCT images of normal-appearing renal pelvis and ureter, urothelium, lamina propria, and muscularis were visible. With ELUS, all anatomical layers of the ureter could be distinguished, besides the urothelial layer. ELUS identified suspect lesions, although exact staging and differentiation between noninvasive and invasive lesions were not possible. OCT provides high-resolution imaging of normal ureter and ureter lesions. ELUS, however, is of limited value as it cannot differentiate between noninvasive and invasive tumors.

Intracoronary Shear Stress and CT Characteristics of Vulnerable Coronary Plaques

Vulnerable coronary plaques are associated with a significant risk for rupture, and the ability to detect their characteristic features is of extreme importance, as timely detection of rupture-prone plaques could lead to the appropriate initiation of adequate therapeutic measures and prevent the evolution to an acute coronary event. The most common features of vulnerability in coronary plaques are represented by the presence of low density atheroma, a thin fibrous cap, spotty calcifications, and positive remodeling. However, there is still a huge amount of information to be learned about the role of local forces, represented by the shear stress, on the plaque vulnerability. This clinical update aims to present the most recent advances in the field of knowledge regarding the relation between shear stress and plaque vulnerability, starting from the hypothesis that shear stress significantly correlates with the CT features of plaque vulnerability and can represent a new marker of vulnerability in coronary artery plaques.

Effect of Electrical Stimulation on Blood Flow Velocity and Vessel Size

Interferential current electrical stimulation alters blood flow velocity and vessel size. We aimed to investigate the changes in the autonomic nervous system depending on electrical stimulation parameters. Forty-five healthy adult male and female subjects were studied. Bipolar adhesive pad electrodes were used to stimulate the autonomic nervous system at the thoracic vertebrae 1-4 levels for 20 min. Using Doppler ultrasonography, blood flow was measured to determine velocity and vessel size before, immediately after, and 30 min after electrical stimulation. Changes in blood flow velocity were significantly different immediately and 30 min after stimulation. The interaction between intervention periods and groups was significantly different between the exercise and pain stimulation groups immediately after stimulation (p<0.05). The vessel size was significantly different before and 30 min after stimulation (p<0.05). Imbalances in the sympathetic nervous system, which regulates balance throughout the body, may present with various symptoms. Therefore, in the clinical practice, the parameters of electrical stimulation should be selectively applied in accordance with various conditions and changes in form.

Hybrid intravascular imaging: recent advances, technical considerations, and current applications in the study of plaque pathophysiology

Cumulative evidence from histology-based studies demonstrate that the currently available intravascular imaging techniques have fundamental limitations that do not allow complete and detailed evaluation of plaque morphology and pathobiology, limiting the ability to accurately identify high-risk plaques. To overcome these drawbacks, new efforts are developing for data fusion methodologies and the design of hybrid, dual-probe catheters to enable accurate assessment of plaque characteristics, and reliable identification of high-risk lesions. Today several dual-probe catheters have been introduced including combined near infrared spectroscopy-intravascular ultrasound (NIRS-IVUS), that is already commercially available, IVUS-optical coherence tomography (OCT), the OCT-NIRS, the OCT-near infrared fluorescence (NIRF) molecular imaging, IVUS-NIRF, IVUS intravascular photoacoustic imaging and combined fluorescence lifetime-IVUS imaging. These multimodal approaches appear able to overcome limitations of standalone imaging and provide comprehensive visualization of plaque composition and plaque biology. The aim of this review article is to summarize the advances in hybrid intravascular imaging, discuss the technical challenges that should be addressed in order to have a use in the clinical arena, and present the evidence from their first applications aiming to highlight their potential value in the study of atherosclerosis.

High shear stress induces atherosclerotic vulnerable plaque formation through angiogenesis

Rupture of atherosclerotic plaques causing thrombosis is the main cause of acute coronary syndrome and ischemic strokes. Inhibition of thrombosis is one of the important tasks developing biomedical materials such as intravascular stents and vascular grafts. Shear stress (SS) influences the formation and development of atherosclerosis. The current review focuses on the vulnerable plaques observed in the high shear stress (HSS) regions, which localizes at the proximal region of the plaque intruding into the lumen. The vascular outward remodelling occurs in the HSS region for vascular compensation and that angiogenesis is a critical factor for HSS which induces atherosclerotic vulnerable plaque formation. These results greatly challenge the established belief that low shear stress is important for expansive remodelling, which provides a new perspective for preventing the transition of stable plaques to high-risk atherosclerotic lesions.

Interaction between the Stress Phase Angle (SPA) and the Oscillatory Shear Index (OSI) Affects Endothelial Cell Gene Expression

Hemodynamic forces play an important role in the non-uniform distribution of atherosclerotic lesions. Endothelial cells are exposed simultaneously to fluid wall shear stress (WSS) and solid circumferential stress (CS). Due to variations in impedance (global factors) and geometric complexities (local factors) in the arterial circulation a time lag arises between these two forces that can be characterized by the temporal phase angle between CS and WSS (stress phase angle-SPA). Asynchronous flows (SPA close to -180°) that are most prominent in coronary arteries have been associated with localization of atherosclerosis. Reversing oscillatory flows characterized by an oscillatory shear index (OSI) that is great than zero are also associated with atherosclerosis localization. In this study we examined the relationship between asynchronous flows and reversing flows in altering the expression of 37 genes relevant to atherosclerosis development. In the case of reversing oscillatory flow, we observed that the asynchronous condition upregulated 8 genes compared to synchronous hemodynamics, most of them proatherogenic. Upregulation of the pro-inflammatory transcription factor NFκB p65 was confirmed by western blot, and nuclear translocation of NFκB p65 was confirmed by immunofluorescence staining. A comparative study between non-reversing flow and reversing flow found that in the case of synchronous hemodynamics, reversing flow altered the expression of 11 genes, while in the case of asynchronous hemodynamics, reversing flow altered the expression of 17 genes. Reversing flow significantly upregulated protein expression of NFκB p65 for both synchronous and asynchronous conditions. Nuclear translocation of NFκB p65 was confirmed for synchronous and asynchronous conditions in the presence of flow reversal. These data suggest that asynchronous hemodynamics and reversing flow can elicit proatherogenic responses in endothelial cells compared to synchronous hemodynamics without shear stress reversal, indicating that SPA as well as reversal flow (OSI) are important parameters characterizing arterial susceptibility to disease.

Automatic Coronary Artery Segmentation Using Active Search for Branches and Seemingly Disconnected Vessel Segments from Coronary CT Angiography

We propose a Bayesian tracking and segmentation method of coronary arteries on coronary computed tomographic angiography (CCTA). The geometry of coronary arteries including lumen boundary is estimated in Maximum A Posteriori (MAP) framework. Three consecutive sphere based filtering is combined with a stochastic process that is based on the similarity of the consecutive local neighborhood voxels and the geometric constraint of a vessel. It is also founded on the prior knowledge that an artery can be seen locally disconnected and consist of branches which may be seemingly disconnected due to plaque build up. For such problem, an active search method is proposed to find branches and seemingly disconnected but actually connected vessel segments. Several new measures have been developed for branch detection, disconnection check and planar vesselness measure. Using public domain Rotterdam CT dataset, the accuracy of extracted centerline is demonstrated and automatic reconstruction of coronary artery mesh is shown.

Coronary fractional flow reserve measurements of a stenosed side branch: a computational study investigating the influence of the bifurcation angle

Background

Coronary hemodynamics and physiology specific for bifurcation lesions was not well understood. To investigate the influence of the bifurcation angle on the intracoronary hemodynamics of side branch (SB) lesions computational fluid dynamics simulations were performed.

Methods

A parametric model representing a left anterior descending—first diagonal coronary bifurcation lesion was created according to the literature. Diameters obeyed fractal branching laws. Proximal and distal main branch (DMB) stenoses were both set at 60 %. We varied the distal bifurcation angles (40°, 55°, and 70°), the flow splits to the DMB and SB (55 %:45 %, 65 %:35 %, and 75 %:25 %), and the SB stenoses (40, 60, and 80 %), resulting in 27 simulations. Fractional flow reserve, defined as the ratio between the mean distal stenosis and mean aortic pressure during maximal hyperemia, was calculated for the DMB and SB (FFR_SB) for all simulations.

Results

The largest differences in FFR_SB comparing the largest and smallest bifurcation angles were 0.02 (in cases with 40 % SB stenosis, irrespective of the assumed flow split) and 0.05 (in cases with 60 % SB stenosis, flow split 55 %:45 %). When the SB stenosis was 80 %, the difference in FFR_SB between the largest and smallest bifurcation angle was 0.33 (flow split 55 %:45 %). By describing the ΔP_SB−Q_SB relationship using a quadratic curve for cases with 80 % SB stenosis, we found that the curve was steeper (i.e. higher flow resistance) when bifurcation angle increases (ΔP = 0.451*Q + 0.010*Q² and ΔP = 0.687*Q + 0.017*Q² for 40° and 70° bifurcation angle, respectively). Our analyses revealed complex hemodynamics in all cases with evident counter-rotating helical flow structures. Larger bifurcation angles resulted in more pronounced helical flow structures (i.e. higher helicity intensity), when 60 or 80 % SB stenoses were present. A good correlation (R² = 0.80) between the SB pressure drop and helicity intensity was also found.

Conclusions

Our analyses showed that, in bifurcation lesions with 60 % MB stenosis and 80 % SB stenosis, SB pressure drop is higher for larger bifurcation angles suggesting higher flow resistance (i.e. curves describing the ΔP_SB−Q_SB relationship being steeper). When the SB stenosis is mild (40 %) or moderate (60 %), SB resistance is minimally influenced by the bifurcation angle, with differences not being clinically meaningful. Our findings also highlighted the complex interplay between anatomy, pressure drops, and blood flow helicity in bifurcations.

Performance of a thrombectomy device for aspiration of thrombus with various sizes based on a computational fluid dynamic modeling

It is important to thoroughly remove the thrombus within the course of aspiration thrombectomy; otherwise, it may lead to further embolization. The performance of the aspiration thrombectomy device with a generic geometry is studied through the computational approach. In order to model the thrombus aspiration, a real left coronary artery is chosen while thrombi with various sizes are located at the bifurcation area of the coronary artery and, depending on the size of the thrombus, it is stretched toward the side branches. The thrombus occupies the artery resembling the blood current obstruction in the coronary vessel similar to the situation that leads to heart attack. It is concluded that the aspiration ability of the thrombectomy device is not linked to the thrombus size; it is rather linked to the aspiration pressure and thrombus age (organized versus fresh thrombus). However, the aspiration time period correlates to the thrombus size. The minimum applicable aspiration pressure is also investigated in this study.

Fusion of fibrous cap thickness and wall shear stress to assess plaque vulnerability in coronary arteries: a pilot study

Purpose

Identification of rupture-prone plaques in coronary arteries is a major clinical challenge. Fibrous cap thickness and wall shear stress are two relevant image-based risk factors, but these two parameters are generally computed and analyzed separately. Accordingly, combining these two parameters can potentially improve the identification of at-risk regions. Therefore, the purpose of this study is to investigate the feasibility of the fusion of wall shear stress and fibrous cap thickness of coronary arteries in patient data.

Methods

Fourteen patients were included in this pilot study. Imaging of the coronary arteries was performed with optical coherence tomography and with angiography. Fibrous cap thickness was automatically quantified from optical coherence tomography pullbacks using a contour segmentation approach based on fast marching. Wall shear stress was computed by applying computational fluid dynamics on the 3D volume reconstructed from two angiograms. The two parameters then were co-registered using anatomical landmarks such as side branches.

Results

The two image modalities were successfully co-registered, with a mean (±SD) error corresponding to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$8.6\,\pm \,6.7\,\%$$\end{document}8.6±6.7% of the length of the analyzed region. For all the analyzed participants, the average thinnest portion of each fibrous cap was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$129\,\pm \,69\,\upmu \text {m}$$\end{document}129±69μm, and the average WSS value at the location of the fibrous cap was \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.46\,\pm \,1.16\,\text {Pa}$$\end{document}1.46±1.16Pa. A unique index was finally generated for each patient via the fusion of fibrous cap thickness and wall shear stress measurements, to translate all the measured parameters into a single risk map.

Conclusion

The introduced risk map integrates two complementary parameters and has potential to provide valuable information about plaque vulnerability.

Comparison of angiographic and IVUS derived coronary geometric reconstructions for evaluation of the association of hemodynamics with coronary artery disease progression

Wall shear stress (WSS) has been investigated as a prognostic marker for the prospective identification of rapidly progressing coronary artery disease (CAD) and atherosclerotic lesions likely to gain high-risk (vulnerable) characteristics. The goal of this study was to compare biplane angiographic vs. intravascular ultrasound (IVUS) derived reconstructed coronary geometries to evaluate agreement in geometry, computed WSS, and association of WSS and CAD progression. Baseline and 6-month follow-up angiographic and IVUS imaging data were collected in patients with non-obstructive CAD (n = 5). Three-dimensional (3D) reconstructions of the coronary arteries were generated with each technique, and patient-specific computational fluid dynamics models were constructed to compute baseline WSS values. Geometric comparisons were evaluated in arterial segments (n = 9), and hemodynamic data were evaluated in circumferential sections (n = 468). CAD progression was quantified from serial IVUS imaging data (n = 277), and included virtual-histology IVUS (VH-IVUS) derived changes in plaque composition. There was no significant difference in reconstructed coronary segment lengths and cross-sectional areas (CSA), however, IVUS derived geometries exhibited a significantly larger left main CSA than the angiographic reconstructions. Computed absolute time-averaged WSS (TAWSS_ABS) values were significantly greater in the IVUS derived geometries, however, evaluations of relative TAWSS (TAWSS_REL) values revealed improved agreement and differences within defined zones of equivalence. Associations between VH-IVUS defined CAD progression and angiographic or IVUS derived WSS exhibited poor agreement when examining TAWSS_ABS data, but improved when evaluating the association with TAWSS_REL data. We present data from a small cohort of patients highlighting strong agreement between angiographic and IVUS derived coronary geometries, however, limited agreement is observed between computed WSS values and associations of WSS with CAD progression.

Validation study of a 3D-QCA coronary reconstruction method using a hybrid intravascular ultrasound and angiography reconstruction method and patient-specific Fractional Flow Reserve data

The estimation of the severity of coronary lesions is of utmost importance in today's clinical practice, since Cardiovascular diseases often have fatal consequences. The most efficient method to estimate the severity of a lesion is the calculation of the Fractional Flow Reserve. The necessary use of a pressure wire, however, makes this method invasive and strenuous for the patient. In this work, we present a novel 3-Dimensional Quantitative Coronary Analysis coronary reconstruction method and a framework for the computation of the virtual Functional Assessment Index (vFAI). In a dataset of 5 coronary arterial segments, we use the aforementioned method to reconstruct them in 3D, and compare them to the respective 3D models reconstructed from our already validated hybrid IVUS-angiography reconstruction method [2]. The obtained results indicate a high correlation between the two methods in terms of the calculated FFR values, presenting a difference of 3.19% in the worst case scenario. Furthermore, when compared to the actual FFR values that derive from a pressure wire, the differences were statistically insignificant.

Co-registration of angiography and intravascular ultrasound images through image-based device tracking

OBJECTIVES

To determine the feasibility of automated co-registration of angiography and intravascular ultrasound (IVUS) to facilitate integration of these two imaging modalities in a synchronous manner.

BACKGROUND

IVUS provides cross-sectional imaging of coronary arteries but lacks overview of the vascular territory provided by angiography. Co-registration of angiography and IVUS would increase utility of IVUS in the clinical setting.

METHODS

Forty-nine consecutive patients undergoing surveillance for cardiac allograft vasculopathy with angiography and IVUS of the left anterior descending artery (LAD) were enrolled. A pre-IVUS angiogram of the LAD was performed followed by an ECG-triggered fluoroscopy (ECGTF) during IVUS pullback at 0.5 mm/s using an automatic pullback device. ECGTF was used to track the IVUS catheter during pullback and establish a spatial relationship to the pre-IVUS angiogram. Angio-IVUS co-registration was performed with a research prototype (Siemens Healthcare, Germany) and accuracy was evaluated by distance mismatch between angiography and IVUS images at vessel bifurcations.

RESULTS

Median age was 54 (44.5, 67) years. The population was 82.6% male with minimal risk factors. The median (IQR) co-registration distance mismatch measured at 108 bifurcations in 42 (85%) patients was 0.35 (0.00-1.16) mm. Seven patients were excluded due to inappropriate data acquisition (n = 3) and failure of tracking (n = 4), e.g., due to overlapping sternal wires. Estimated effective radiation dose for ECGTF was 0.09 mSv.

CONCLUSION

This study demonstrates the feasibility of angio-IVUS co-registration which may be used as a clinical tool for localizing IVUS cross-sections along an angiographic roadmap. © 2015 Wiley Periodicals, Inc.

Computational assessment of the fractional flow reserve from intravascular ultrasound and coronary angiography data: a pilot study

Cardiovascular disease is one of the primary causes of morbidity and mortality around the globe. Thus, the diagnosis of critical lesions in coronary arteries is of utmost importance in clinical practice. One useful and efficient method to assess the functional severity of one or multiple lesions in a coronary artery is the calculation of the fractional flow reserve (FFR). In the current work, we present a method which allows the calculation of the FFR value computationally, without the use of a pressure wire and the induction of hyperemia, using intravascular ultrasound (IVUS) and biplane angiography images for three-dimensional (3D) coronary artery reconstruction and measurements of the volumetric flow rate derived from angiographic sequences. The simulated FFR values were compared to the invasively measured FFR values in 7 cases, presenting high correlation (r=0.85) and good agreement (mean difference=0.002). FFR assessment without employing a pressure wire and the induction of hyperemia is feasible using 3D reconstructed coronary artery models from angiographic and IVUS data coupled with computational fluid dynamics.

Coronary Atherosclerosis and Interventional Cardiology

The atherosclerotic process in coronary arteries begins with endothelial dysfunction and may provoke thrombotic total occlusion and myocardial infarction. In this state-of-the-art review, we discuss recent evidence of atheroslerosis, vulnerable plaque, and hemodynamic changes in the coronary tree, as well as the current techniques we implement in the catheterization lab to evaluate coronary stenosis. It is clear that atherosclerosis is a chronic inflammatory condition with several consequences in the coronary tree, however, we are able now to characterize the plaque and to select the appropriate treatment for many patients.

A new methodology for accurate 3-dimensional coronary artery reconstruction using routine intravascular ultrasound and angiographic data: implications for widespread assessment of endothelial shear stress in humans

AIMS

To develop and validate a new methodology that allows accurate 3-dimensional (3-D) coronary artery reconstruction using standard, simple angiographic and intravascular ultrasound (IVUS) data acquired during routine catheterisation enabling reliable assessment of the endothelial shear stress (ESS) distribution.

METHODS AND RESULTS

Twenty-two patients (22 arteries: 7 LAD; 7 LCx; 8 RCA) who underwent angiography and IVUS examination were included. The acquired data were used for 3-D reconstruction using a conventional method and a new methodology that utilised the luminal 3-D centreline to place the detected IVUS borders and anatomical landmarks to estimate their orientation. The local ESS distribution was assessed by computational fluid dynamics. In corresponding consecutive 3 mm segments, lumen, plaque and ESS measurements in the 3-D models derived by the centreline approach were highly correlated to those derived from the conventional method (r>0.98 for all). The centreline methodology had a 99.5% diagnostic accuracy for identifying segments exposed to low ESS and provided similar estimations to the conventional method for the association between the change in plaque burden and ESS (centreline method: slope= -1.65%/Pa, p=0.078; conventional method: slope= -1.64%/Pa, p=0.084; p =0.69 for difference between the two methodologies).

CONCLUSIONS

The centreline methodology provides geometrically correct models and permits reliable ESS computation. The ability to utilise data acquired during routine coronary angiography and IVUS examination will facilitate clinical investigation of the role of local ESS patterns in the natural history of coronary atherosclerosis.

Assessing the hemodynamic influence between multiple lesions in a realistic right coronary artery segment: A computational study

Coronary artery disease is the primary cause of morbidity and mortality worldwide. Therefore, detailed assessment of lesions in the coronary vasculature is critical in current clinical practice. Fractional flow reserve (FFR) has been proven as an efficient method for assessing the hemodynamic severity of a coronary stenosis. However, functional assessment of a coronary segment with multiple stenoses (≥ 2) remains complex for guiding the strategy of percutaneous coronary intervention due to the hemodynamic interplay between adjacent stenoses. In this work, we created four 3-dimensional (3D) arterial models that derive from a healthy patient-specific right coronary artery segment. The initial healthy model was reconstructed using fusion of intravascular ultrasound (IVUS) and biplane angiographic patient data. The healthy 3D model presented a measured FFR value of 0.96 (pressure-wire) and a simulated FFR value of 0.98. We then created diseased models with two artificial sequential stenoses of 90% lumen area reduction or with the proximal and distal stenosis separately. We calculated the FFR value for each case: 0.65 for the case with the two stenoses, 0.73 for the case with the distal stenosis and 0.90 for the case with the proximal stenosis. This leads to the conclusion that although both stenoses had the same degree of lumen area stenosis, there was a large difference in hemodynamic severity, thereby indicating that angiographic lumen assessment by itself is often not adequate for accurate assessment of coronary lesions.

A fast seed detection using local geometrical feature for automatic tracking of coronary arteries in CTA

We propose a fast seed detection for automatic tracking of coronary arteries in coronary computed tomographic angiography (CCTA). To detect vessel regions, Hessian-based filtering is combined with a new local geometric feature that is based on the similarity of the consecutive cross-sections perpendicular to the vessel direction. It is in turn founded on the prior knowledge that a vessel segment is shaped like a cylinder in axial slices. To improve computational efficiency, an axial slice, which contains part of three main coronary arteries, is selected and regions of interest (ROIs) are extracted in the slice. Only for the voxels belonging to the ROIs, the proposed geometric feature is calculated. With the seed points, which are the centroids of the detected vessel regions, and their vessel directions, vessel tracking method can be used for artery extraction. Here a particle filtering-based tracking algorithm is tested. Using 19 clinical CCTA datasets, it is demonstrated that the proposed method detects seed points and can be used for full automatic coronary artery extraction. ROC (receiver operating characteristic) curve analysis shows the advantages of the proposed method.

Computational fluid dynamics in coronary artery disease

Computational fluid dynamics (CFD) is a widely used method in mechanical engineering to solve complex problems by analysing fluid flow, heat transfer, and associated phenomena by using computer simulations. In recent years, CFD has been increasingly used in biomedical research of coronary artery disease because of its high performance hardware and software. CFD techniques have been applied to study cardiovascular haemodynamics through simulation tools to predict the behaviour of circulatory blood flow in the human body. CFD simulation based on 3D luminal reconstructions can be used to analyse the local flow fields and flow profiling due to changes of coronary artery geometry, thus, identifying risk factors for development and progression of coronary artery disease. This review aims to provide an overview of the CFD applications in coronary artery disease, including biomechanics of atherosclerotic plaques, plaque progression and rupture; regional haemodynamics relative to plaque location and composition. A critical appraisal is given to a more recently developed application, fractional flow reserve based on CFD computation with regard to its diagnostic accuracy in the detection of haemodynamically significant coronary artery disease.

3D reconstruction techniques of human coronary bifurcations for shear stress computations

BACKGROUND

Heterogeneity in plaque composition in human coronary artery bifurcations is associated with blood flow induced shear stress. Shear stress is generally determined by combing 3D lumen data and computational fluid dynamics (CFD). We investigated two new procedures to generate 3D lumen reconstructions of coronary artery bifurcations for shear stress computations.

METHODS

We imaged 10 patients with multislice computer tomography (MSCT) and intravascular ultrasound (IVUS). The 3D reconstruction of the main branch was based on the fusion of MSCT and IVUS. The proximal part of side branch was reconstructed using IVUS data or MSCT data, resulting in two different reconstructions of the bifurcation region. The distal part of the side branch was based on MSCT data alone. The reconstructed lumen was combined with CFD to determine the shear stress. Low and high shear stress regions were defined and shear stress patterns in the bifurcation regions were investigated.

RESULTS

The 3D coronary bifurcations were successfully generated with both reconstruction procedures. The geometrical features of the bifurcation region for the two reconstruction procedures did not reveal appreciable differences. The shear stress maps showed a qualitative agreement, and the low and high shear stress regions were similar in size and average shear stress values were identical. The low and high shear stress regions showed an overlap of approximately 75%.

CONCLUSION

Reconstruction of the side branch with MSCT data alone is an adequate technique to study shear stress and wall thickness in the bifurcation region. The reconstruction procedure can be applied to further investigate the effect of shear stress on atherosclerosis in coronary bifurcations.

Does low and oscillatory wall shear stress correlate spatially with early atherosclerosis? A systematic review

Low and oscillatory wall shear stress is widely assumed to play a key role in the initiation and development of atherosclerosis. Indeed, some studies have relied on the low shear theory when developing diagnostic and treatment strategies for cardiovascular disease. We wished to ascertain if this consensus is justified by published data. We performed a systematic review of papers that compare the localization of atherosclerotic lesions with the distribution of haemodynamic indicators calculated using computational fluid dynamics. The review showed that although many articles claim their results conform to the theory, it has been interpreted in different ways: a range of metrics has been used to characterize the distribution of disease, and they have been compared with a range of haemodynamic factors. Several studies, including all of those making systematic point-by-point comparisons of shear and disease, failed to find the expected relation. The various pre- and post-processing techniques used by different groups have reduced the range of shears over which correlations were sought, and in some cases are mutually incompatible. Finally, only a subset of the known patterns of disease has been investigated. The evidence for the low/oscillatory shear theory is less robust than commonly assumed. Longitudinal studies starting from the healthy state, or the collection of average flow metrics derived from large numbers of healthy vessels, both in conjunction with point-by-point comparisons using appropriate statistical techniques, will be necessary to improve our understanding of the relation between blood flow and atherogenesis.

Coronary pressure-derived fractional flow reserve in the assessment of coronary artery stenoses

OBJECTIVES

Catheter-based angiography is the reference-standard to establish coronary anatomy. While routinely employed clinically, lumen assessment correlates poorly with physiological measures of ischaemia. Moreover, functional studies to identify and localise ischaemia before elective angiography are often not available. This article reviews fractional flow reserve (FFR) and its role in guiding patient management for patients with a potentially haemodynamic significant coronary lesion.

METHODS

This review discusses the theory, evidence, indications, and limitations of FFR. Also included are emerging non-invasive imaging FFR surrogates currently under evaluation for accuracy with respect to standard FFR.

RESULTS

Coronary pressure-derived fractional flow reserve (FFR) rapidly assesses the haemodynamic significance of individual coronary artery lesions and can readily be performed in the catheterisation laboratory. The use of FFR has been shown to effectively guide coronary revascularization procedures leading to improved patient outcomes.

CONCLUSIONS

FFR is an invaluable modality in guiding coronary disease treatment decisions. It is safe, cost-effective and leads to improved patient outcomes. Non-invasive imaging modalities to assess the physiologic significance of CAD are currently being  developed and evaluated.

Modeling stented coronary arteries: where we are, where to go

In the last two decades, numerical models have become well-recognized and widely adopted tools to investigate stenting procedures. Due to limited computational resources and modeling capabilities, early numerical studies only involved simplified cases and idealized stented arteries. Nowadays, increased computational power allows for numerical models to meet clinical needs and include more complex cases such as the implantation of multiple stents in bifurcations or curved vessels. Interesting progresses have been made in the numerical modeling of stenting procedures both from a structural and a fluid dynamics points of view. Moreover, in the drug eluting stents era, new insights on drug elution capabilities are becoming essential in the stent development. Lastly, image-based methods able to reconstruct realistic geometries from medical images have been proposed in the recent literature aiming to better describe the peculiar anatomical features of coronary vessels and increase the accuracy of the numerical models. In this light, this review provides a comprehensive analysis of the current state-of-the-art in this research area, discussing the main methodological advances and remarkable results drawn from a number of significant studies.