Framework to co-register longitudinal virtual histology-intravascular ultrasound data in the circumferential direction

Automated finite element approach to generate anatomical patient-specific biomechanical models of atherosclerotic arteries from virtual histology-intravascular ultrasound

Despite advancements in early detection and treatment, atherosclerosis remains the leading cause of death across all cardiovascular diseases (CVD). Biomechanical analysis of atherosclerotic lesions has the potential to reveal biomechanically instable or rupture-prone regions. Treatment decisions rarely consider the biomechanics of the stenosed lesion due in-part to difficulties in obtaining this information in a clinical setting. Previous 3D FEA approaches have incompletely incorporated the complex curvature of arterial geometry, material heterogeneity, and use of patient-specific data. To address these limitations and clinical need, herein we present a user-friendly fully automated program to reconstruct and simulate the wall mechanics of patient-specific atherosclerotic coronary arteries. The program enables 3D reconstruction from patient-specific data with heterogenous tissue assignment and complex arterial curvature. Eleven arteries with coronary artery disease (CAD) underwent baseline and 6-month follow-up angiographic and virtual histology-intravascular ultrasound (VH-IVUS) imaging. VH-IVUS images were processed to remove background noise, extract VH plaque material data, and luminal and outer contours. Angiography data was used to orient the artery profiles along the 3D centerlines. The resulting surface mesh is then resampled for uniformity and tetrahedralized to generate the volumetric mesh using TetGen. A mesh convergence study revealed edge lengths between 0.04 mm and 0.2 mm produced constituent volumes that were largely unchanged, hence, to save computational resources, a value of 0.2 mm was used throughout. Materials are assigned and finite element analysis (FEA) is then performed to determine stresses and strains across the artery wall. In a representative artery, the highest average effective stress was in calcium elements with 235 kPa while necrotic elements had the lowest average stress, reaching as low as 0.79 kPa. After applying nodal smoothening, the maximum effective stress across 11 arteries remained below 288 kPa, implying biomechanically stable plaques. Indeed, all atherosclerotic plaques remained unruptured at the 6-month longitudinal follow up diagnosis. These results suggest our automated analysis may facilitate assessment of atherosclerotic plaque stability.

Framework for lumen-based nonrigid tomographic coregistration of intravascular images

Purpose: Modern medical imaging enables clinicians to effectively diagnose, monitor, and treat diseases. However, clinical decision-making often relies on combined evaluation of either longitudinal or disparate image sets, necessitating coregistration of multiple acquisitions. Promising coregistration techniques have been proposed; however, available methods predominantly rely on time-consuming manual alignments or nontrivial feature extraction with limited clinical applicability. Addressing these issues, we present a fully automated, robust, nonrigid registration method, allowing for coregistering of multimodal tomographic vascular image datasets using luminal annotation as the sole alignment feature. Approach: Registration is carried out by the use of the registration metrics defined exclusively for lumens shapes. The framework is primarily broken down into two sequential parts: longitudinal and rotational registration. Both techniques are inherently nonrigid in nature to compensate for motion and acquisition artifacts in tomographic images. Results: Performance was evaluated across multimodal intravascular datasets, as well as in longitudinal cases assessing pre-/postinterventional coronary images. Low registration error in both datasets highlights method utility, with longitudinal registration errors-evaluated throughout the paired tomographic sequences-of 0.29 ± 0.14    mm ( < 2 longitudinal image frames) and 0.18 ± 0.16    mm ( < 1 frame) for multimodal and interventional datasets, respectively. Angular registration for the interventional dataset rendered errors of 7.7 ° ± 6.7 ° , and 29.1 ° ± 23.2 ° for the multimodal set. Conclusions: Satisfactory results across datasets, along with additional attributes such as the ability to avoid longitudinal over-fitting and correct nonlinear catheter rotation during nonrigid rotational registration, highlight the potential wide-ranging applicability of our presented coregistration method.

Oscillatory wall shear stress is a dominant flow characteristic affecting lesion progression patterns and plaque vulnerability in patients with coronary artery disease

Although experimental studies suggest that low and oscillatory wall shear stress (WSS) promotes plaque transformation to a more vulnerable phenotype, this relationship has not been examined in human atherosclerosis progression. Thus, the aim of this investigation was to examine the association between oscillatory WSS, in combination with WSS magnitude, and coronary atherosclerosis progression. We hypothesized that regions of low and oscillatory WSS will demonstrate progression towards more vulnerable lesions, while regions exposed to low and non-oscillatory WSS will exhibit progression towards more stable lesions. Patients (n = 20) with non-flow-limiting coronary artery disease (CAD) underwent baseline and six-month follow-up angiography, Doppler velocity and radiofrequency intravascular ultrasound (VH-IVUS) acquisition. Computational fluid dynamics models were constructed to compute time-averaged WSS magnitude and oscillatory WSS. Changes in VH-IVUS-defined total plaque and constituent areas were quantified in focal regions (i.e. sectors; n = 14 235) and compared across haemodynamic categories. Compared with sectors exposed to low WSS magnitude, high WSS sectors demonstrated regression of total plaque area (p < 0.001) and fibrous tissue (p < 0.001), and similar progression of necrotic core. Sectors subjected to low and oscillatory WSS exhibited total plaque area regression, while low and non-oscillatory WSS sectors demonstrated total plaque progression (p < 0.001). Furthermore, compared with low and non-oscillatory WSS areas, sectors exposed to low and oscillatory WSS demonstrated regression of fibrous (p < 0.001) and fibrofatty (p < 0.001) tissue and similar progression of necrotic core (p = 0.82) and dense calcium (p = 0.40). Herein, we demonstrate that, in patients with non-obstructive CAD, sectors subjected to low and oscillatory WSS demonstrated regression of total plaque, fibrous and fibrofatty tissue, and progression of necrotic core and dense calcium, which suggest a transformation to a more vulnerable phenotype.

An agent-based model of leukocyte transendothelial migration during atherogenesis

A vast amount of work has been dedicated to the effects of hemodynamics and cytokines on leukocyte adhesion and trans-endothelial migration (TEM) and subsequent accumulation of leukocyte-derived foam cells in the artery wall. However, a comprehensive mechanobiological model to capture these spatiotemporal events and predict the growth and remodeling of an atherosclerotic artery is still lacking. Here, we present a multiscale model of leukocyte TEM and plaque evolution in the left anterior descending (LAD) coronary artery. The approach integrates cellular behaviors via agent-based modeling (ABM) and hemodynamic effects via computational fluid dynamics (CFD). In this computational framework, the ABM implements the diffusion kinetics of key biological proteins, namely Low Density Lipoprotein (LDL), Tissue Necrosis Factor alpha (TNF-α), Interlukin-10 (IL-10) and Interlukin-1 beta (IL-1β), to predict chemotactic driven leukocyte migration into and within the artery wall. The ABM also considers wall shear stress (WSS) dependent leukocyte TEM and compensatory arterial remodeling obeying Glagov's phenomenon. Interestingly, using fully developed steady blood flow does not result in a representative number of leukocyte TEM as compared to pulsatile flow, whereas passing WSS at peak systole of the pulsatile flow waveform does. Moreover, using the model, we have found leukocyte TEM increases monotonically with decreases in luminal volume. At critical plaque shapes the WSS changes rapidly resulting in sudden increases in leukocyte TEM suggesting lumen volumes that will give rise to rapid plaque growth rates if left untreated. Overall this multi-scale and multi-physics approach appropriately captures and integrates the spatiotemporal events occurring at the cellular level in order to predict leukocyte transmigration and plaque evolution.

Quantification of the focal progression of coronary atherosclerosis through automated co-registration of virtual histology-intravascular ultrasound imaging data

The goal of this study was to evaluate the accuracy of a novel algorithm that circumferentially co-registers serial virtual histology-intravascular ultrasound (VH-IVUS) data for the focal assessment of coronary atherosclerosis progression. Thirty-three patients with an abnormal non-invasive cardiac stress test or stable angina underwent baseline and follow-up (6 or 12 months) invasive evaluation that included acquisition of VH-IVUS image data. Baseline and follow-up image pairs (n = 4194) were automatically co-registered in the circumferential direction via a multi-variate cross-correlation algorithm. Algorithm stability and accuracy were assessed by comparing results from multiple iterations of the algorithm (iteration 1 vs. iteration 2) and against values determined manually by two expert VH-IVUS readers (algorithm vs. two expert readers). Furthermore, focal plaque progression values were compared between the algorithm and expert readers following co-registration by the independently determined angles. Strong agreement in circumferential co-registration angles were observed across multiple iterations of the algorithm (stability) and between the algorithm and expert readers (accuracy; all concordance correlation coefficients >0.98). Furthermore, circumferential co-registration angles determined by the algorithm were not statistically when compared to values determined by two expert readers (p = 0. 99). Bland-Altman analysis indicated minimal bias when comparing focal VH-IVUS defined plaque progression in corresponding sectors following circumferential co-registration between the algorithm and expert readers. Finally, average differences in changes in total plaque and constituent areas between the algorithm and readers were within the average range of difference between readers (interobserver variability). We present a stable and validated algorithm to automatically circumferentially co-register serial VH-IVUS imaging data for the focal quantification of coronary atherosclerosis progression.

Evaluation of a framework for the co-registration of intravascular ultrasound and optical coherence tomography coronary artery pullbacks

A growing number of studies have used a combination of intravascular ultrasound (IVUS) and optical coherence tomography (OCT) for the assessment of atherosclerotic plaques. Given their respective strengths these imaging modalities highly complement each other. Correlations of hemodynamics and coronary artery disease (CAD) have been extensively investigated with both modalities separately, though not concurrently due to challenges in image registration. Manual co-registration of these modalities is a time expensive task subject to human error, and the development of an automatic method has not been previously addressed. We developed a framework that uses dynamic time warping for the longitudinal co-registration and dynamic programming for the circumferential co-registration of images and evaluated the methodology in a cohort (n = 12) of patients with moderate CAD. Excellent correlation was seen between the algorithm and two expert readers for longitudinal co-registration (CCC = 0.9964, CCC = 0.9959) and circumferential co-registration (CCC = 0.9688, CCC = 0.9598). The mean error of the circumferential co-registration angle was found to be within 10%. A framework for the co-registration of IVUS and OCT pullbacks has been developed which provides a foundation for comprehensive studies of CAD biomechanics.

Simultaneous Registration of Location and Orientation in Intravascular Ultrasound Pullbacks Pairs Via 3D Graph-Based Optimization

A novel method is reported for simultaneous registration of location (axial direction) and orientation (circumferential direction) of two intravascular ultrasound (IVUS) pullbacks of the same vessel taken at different times. Monitoring plaque progression or regression (e.g., during lipid treatment) is of high clinical relevance. Our method uses a 3D graph optimization approach, in which the cost function jointly reflects similarity of plaque morphology and plaque/perivascular image appearance. Graph arcs incorporate prior information about temporal correspondence of the two IVUS sequences and limited angular twisting between consecutive IVUS images. Additionally, our approach automatically identifies starting and ending frame pairs in the two IVUS pullbacks. Validation of our method was performed in 29 pairs of IVUS baseline/follow-up pullback sequences consisting of 8 622 IVUS image frames in total. In comparison to manual registration by three experts, the average location and orientation registration errors ranged from 0.72 mm to 0.79 mm and from 7.3(°) to 9.3(°), respectively, all close to the inter-observer variability with no difference being statistically significant (p = NS). Rotation angles determined by our automated approach and expert observers showed high correlation (r(2) of 0.97 to 0.98) and agreed closely (mutual bias between the automated method and expert observers ranged from -1.57(°) to 0.15(°)). Compared with state-of-the-art approaches, the new method offers lower errors in both location and orientation registration. Our method offers highly automated and accurate IVUS pullback registration and can be employed in IVUS-based studies of coronary disease progression, enabling more focal studies of coronary plaque development and transition of vulnerability.

An assessment of intra-patient variability on observed relationships between wall shear stress and plaque progression in coronary arteries

BACKGROUND

Wall shear stress (WSS) has been associated with sites of plaque localization and with changes in plaque composition in human coronary arteries. Different values have been suggested for categorizing WSS as low, physiologic or high; however, uncertainties in flow rates, both across subjects and within a given individual, can affect the classification of WSS and thus influence the observed relationships between local hemodynamics and plaque changes over time. This study examines the effects of uncertainties in flow rate boundary conditions upon WSS values and investigates the influence of this variability on the observed associations of WSS with changes in VH-IVUS derived plaque components.

METHODS

Three patients with coronary artery disease underwent baseline and 12 month follow-up angiography and virtual histology-intravascular ultrasound (VH-IVUS) measurements. Coronary artery models were reconstructed from the data and models with and without side-branches were created. Patient-specific Doppler ultrasound (DUS) data were employed as inflow boundary conditions and computational fluid dynamics was used to calculate the WSS in each model. Further, the influence of representative coronary artery flow waveforms upon WSS values was investigated and the concept of treating WSS using relative, rather than actual, values was explored.

RESULTS

Models that included side-branch outflows and subject-specific DUS velocities were considered to be the reference cases. Hemodynamic differences were caused by the exclusion of side-branches and by imposing alternative velocity waveforms. One patient with fewer side-branches and a scaled generic waveform had little deviation from the reference case, while another patient with several side-branches excluded showed much larger departures from the reference situation. Differences between models and the respective reference cases were reduced when data were analyzed using relative, rather than actual, WSS.

CONCLUSIONS

When considering individual subjects, large variations in patient-specific flow rates and exclusion of multiple side-branches in computational models can cause significant differences in observed associations between plaque evolution and ranges of computed WSS. These differences may contribute to the large variability typically found among subjects in pooled populations. Relative WSS may be more useful than actual WSS as a correlative variable when there is a large degree of uncertainty in flow rate data.

Biomechanics and inflammation in atherosclerotic plaque erosion and plaque rupture: implications for cardiovascular events in women

Objective

Although plaque erosion causes approximately 40% of all coronary thrombi and disproportionally affects women more than men, its mechanism is not well understood. The role of tissue mechanics in plaque rupture and regulation of mechanosensitive inflammatory proteins is well established, but their role in plaque erosion is unknown. Given obvious differences in morphology between plaque erosion and rupture, we hypothesized that inflammation in general as well as the association between local mechanical strain and inflammation known to exist in plaque rupture may not occur in plaque erosion. Therefore, our objective was to determine if similar mechanisms underlie plaque rupture and plaque erosion.

Methods and Results

We studied a total of 74 human coronary plaque specimens obtained at autopsy. Using lesion-specific computer modeling of solid mechanics, we calculated the stress and strain distribution for each plaque and determined if there were any relationships with markers of inflammation. Consistent with previous studies, inflammatory markers were positively associated with increasing strain in specimens with rupture and thin-cap fibroatheromas. Conversely, overall staining for inflammatory markers and apoptosis were significantly lower in erosion, and there was no relationship with mechanical strain. Samples with plaque erosion most closely resembled those with the stable phenotype of thick-cap fibroatheromas.

Conclusions

In contrast to classic plaque rupture, plaque erosion was not associated with markers of inflammation and mechanical strain. These data suggest that plaque erosion is a distinct pathophysiological process with a different etiology and therefore raises the possibility that a different therapeutic approach may be required to prevent plaque erosion.

Focal association between wall shear stress and clinical coronary artery disease progression

Wall shear stress (WSS) has been investigated as a potential prospective marker to identify rapidly progressing coronary artery disease (CAD) and potential for lesions to acquire vulnerable characteristics. Previous investigations, however, are limited by a lack of understanding of the focal association between WSS and CAD progression (i.e., data are notably spatially averaged). Thus, the aim of this investigation was to examine the focal association between WSS and coronary atherosclerosis progression, and compare these results to those determined by spatial averaging. Five patients with CAD underwent baseline and 6-month follow-up angiographic and virtual histology-intravascular ultrasound imaging to quantify CAD progression. Patient-specific computational fluid dynamics models were constructed to compute baseline WSS values, which were either averaged around the entire artery circumference or examined in focal regions (sectors). Analysis of data within each sector (n = 3871) indicated that circumferentially averaged and sector WSS values were statistically different (p < 0.05) and exhibited poor agreement (concordance correlation coefficient = 0.69). Furthermore, differences were observed between the analysis techniques when examining the association of WSS and CAD progression. This investigation highlights the importance of examining spatially heterogeneous variables at a focal level to reduce the affect of data reduction and warrants implementation in a larger clinical study to determine the predictive power in prospectively identifying rapidly progressing and/or vulnerable coronary plaques.