3-D reconstruction of tissue components for atherosclerotic human arteries using ex vivo high-resolution MRI

Computing of Low Shear Stress-Driven Endothelial Gene Network Involved in Early Stages of Atherosclerotic Process

Background

In the pathogenesis of atherosclerosis, a central role is represented by endothelial inflammation with influx of chemokine-mediated leukocytes in the vascular wall. Aim of this study was to analyze the effect of different shear stresses on endothelial gene expression and compute gene network involved in atherosclerotic disease, in particular to homeostasis, inflammatory cell migration, and apoptotic processes.

Methods

HUVECs were subjected to shear stress of 1, 5, and 10 dyne/cm² in a Flow Bioreactor for 24 hours to compare gene expression modulation. Total RNA was analyzed by Affymetrix technology and the expression of two specific genes (CXCR4 and ICAM-1) was validated by RT-PCR. To highlight possible regulations between genes and as further validation, a bioinformatics analysis was performed.

Results

At low shear stress (1 dyne/cm²) we observed the following: (a) strong upregulation of CXCR4; (b) mild upregulation of Caspase-8; (c) mild downregulation of ICAM-1; (d) marked downexpression of TNFAIP3. Bioinformatics analysis showed the presence of network composed by 59 new interactors (14 transcription factors and 45 microRNAs) appearing strongly related to shear stress.

Conclusions

The significant modulation of these genes at low shear stress and their close relationships through transcription factors and microRNAs suggest that all may promote an initial inflamed endothelial cell phenotype, favoring the atherosclerotic disease.

Novel methodology for 3D reconstruction of carotid arteries and plaque characterization based upon magnetic resonance imaging carotid angiography data

In this study, we present a novel methodology that allows reliable segmentation of the magnetic resonance images (MRIs) for accurate fully automated three-dimensional (3D) reconstruction of the carotid arteries and semiautomated characterization of plaque type. Our approach uses active contours to detect the luminal borders in the time-of-flight images and the outer vessel wall borders in the T(1)-weighted images. The methodology incorporates the connecting components theory for the automated identification of the bifurcation region and a knowledge-based algorithm for the accurate characterization of the plaque components. The proposed segmentation method was validated in randomly selected MRI frames analyzed offline by two expert observers. The interobserver variability of the method for the lumen and outer vessel wall was -1.60%±6.70% and 0.56%±6.28%, respectively, while the Williams Index for all metrics was close to unity. The methodology implemented to identify the composition of the plaque was also validated in 591 images acquired from 24 patients. The obtained Cohen's k was 0.68 (0.60-0.76) for lipid plaques, while the time needed to process an MRI sequence for 3D reconstruction was only 30 s. The obtained results indicate that the proposed methodology allows reliable and automated detection of the luminal and vessel wall borders and fast and accurate characterization of plaque type in carotid MRI sequences. These features render the currently presented methodology a useful tool in the clinical and research arena.

Automated tracing of the adventitial contour of aortoiliac and peripheral arterial walls in CT angiography (CTA) to allow calculation of non-calcified plaque burden

Aortoiliac and lower extremity arterial atherosclerotic plaque burden is a risk factor for the development of visceral and peripheral ischemic and aneurismal vascular disease. While prior research allows automated quantification of calcified plaque in these body regions using CT angiograms, no automated method exists to quantify soft plaque. We developed an automatic algorithm that defines the outer wall contour and wall thickness of vessels to quantify non-calcified plaque in CT angiograms of the chest, abdomen, pelvis, and lower extremities. The algorithm encodes the search space as a constrained graph and calculates the outer wall contour by deriving a minimum cost path through the graph, following the visible outer wall contour while minimizing path tortuosity. Our algorithm was statistically equivalent to a reference standard made by two reviewers. Absolute error was 1.9 ± 2.3% compared to the inter-observer variability of 3.9 ± 3.6%. Wall thickness in vessels with atherosclerosis was 3.4 ± 1.6 mm compared to 1.2 ± 0.4 mm in normal vessels. The algorithm shows promise as a tool for quantification of non-calcified plaque in CT angiography. When combined with previous research, our method has the potential to quantify both non-calcified and calcified plaque in all clinically significant systemic arteries, from the thoracic aorta to the arteries of the calf, over a wide range of diameters. This algorithm has the potential to enable risk stratification of patients and facilitate investigations into the relationships between asymptomatic atherosclerosis and a variety of behavioral, physiologic, pathologic, and genotypic conditions.

Identification of heterogeneous elastic properties in stenosed arteries: a numerical plane strain study

Assessing the vulnerability of atherosclerotic plaques requires an accurate knowledge of the mechanical properties of the plaque constituents. It is possible to measure displacements in vivo inside a plaque using magnetic resonance imaging. An important issue is to solve the inverse problem that consists in estimating the elastic properties inside the plaque from measured displacements. This study focuses on the identifiability of elastic parameters, e.g. on the compromise between identification time and identification accuracy. An idealised plane strain finite element (FE) model is used. The effects of the FE mesh of the a priori assumptions about the constituents, of the measurement resolution and of the data noise are numerically investigated.

A methodology to analyze changes in lipid core and calcification onto fibrous cap vulnerability: the human atherosclerotic carotid bifurcation as an illustratory example

A lipid core that occupies a high proportion of the plaque volume in addition to a thin fibrous cap is a predominant indicator of plaque vulnerability. Nowadays, noninvasive imaging modalities can identify such structural components, however, morphological criteria alone cannot reliably identify high-risk plaques. Information, such as stresses in the lesion's components, seems to be essential. This work presents a methodology able to analyze the effect of changes in the lipid core and calcification on the wall stresses, in particular, on the fibrous cap vulnerability. Using high-resolution magnetic resonance imaging and histology of an ex vivo human atherosclerotic carotid bifurcation, a patient-specific three-dimensional geometric model, consisting of four tissue components, is generated. The adopted constitutive model accounts for the nonlinear and anisotropic tissue behavior incorporating the collagen fiber orientation by means of a novel and robust algorithm. The material parameters are identified from experimental data. A novel stress-based computational cap vulnerability index is proposed to assess quantitatively the rupture-risk of fibrous caps. Nonlinear finite element analyses identify that the highest stress regions are located at the vicinity of the shoulders of the fibrous cap and in the stiff calcified tissue. A parametric analysis reveals a positive correlation between the increase in lipid core portion and the mechanical stress in the fibrous cap and, hence, the risk for cap rupture. The highest values of the vulnerability index, which correlate to more vulnerable caps, are obtained for morphologies for which the lipid cores were severe; heavily loaded fibrous caps were thus detected. The proposed multidisciplinary methodology is able to investigate quantitatively the mechanical behavior of atherosclerotic plaques in patient-specific stenoses. The introduced vulnerability index may serve as a more quantitative tool for diagnosis, treatment and prevention.

Reconstruction and finite element mesh generation of abdominal aortic aneurysms from computerized tomography angiography data with minimal user interactions

Evaluating rupture risk of abdominal aortic aneurysms is critically important in reducing related mortality without unnecessarily increasing the rate of elective repair. According to the current clinical practice aneurysm rupture risk is (mainly) estimated from its maximum diameter and/or expansion rate; an approach motivated from statistics but known to fail often in individuals. In contrast, recent research demonstrated that patient specific biomechanical simulations can provide more reliable diagnostic parameters, however current structural model development is cumbersome and time consuming. This paper used 2D and 3D deformable models to reconstruct aneurysms from computerized tomography angiography data with minimal user interactions. In particular, formulations of frames and shells, as known from structural mechanics, were used to define deformable modes, which in turn allowed a direct mechanical interpretation of the applied set of reconstruction parameters. Likewise, a parallel finite element implementation of the models allows the segmentation of clinical cases on standard personal computers within a few minutes. The particular topology of the applied 3D deformable models supports a fast and simple hexahedral-dominated meshing of the arising generally polyhedral domain. The variability of the derived segmentations (luminal: 0.50(SD 0.19) mm; exterior 0.89(SD 0.45) mm) with respect to large variations in elastic properties of the deformable models was in the range of the differences between manual segmentations as performed by experts (luminal: 0.57(SD 0.24) mm; exterior: 0.77(SD 0.58) mm), and was particularly independent from the algorithm's initialization. The proposed interaction of deformable models and mesh generation defines finite element meshes suitable to perform accurate and efficient structural analysis of the aneurysm using mixed finite element formulations.

A methodology to study the morphologic changes in lesions during in vitro angioplasty using MRI and image processing

The assessment of morphologic changes in atherosclerotic lesions during interventional procedures such as transluminal balloon angioplasty is an issue of highest clinical importance. We propose a methodology that allows realistic 3D morphomechanical modeling of the vessel, the plaque and the lumen at different stages of in vitro angioplasty. We elaborate on a novel device designed to guide angioplasty under controlled experimental conditions. The device allows to reproduce in vivo conditions as good as possible, i.e. axial in situ pre-stretch, 100mmHg intraluminal pressure, 37 degrees C Tyrode solution, balloon inflation without external constraints using a high-pressure syringe and contrast medium. With a standard 1.5T MR-system we accomplish multi-spectral images at different stages of the angioplasty experiment. After MR image acquisition the specimen is used for histopathological analysis and biomechanical tests. A segmentation process is used to generate NURBS-based 3D geometric models of the individual vessel and plaque components at different balloon pressures. Tissue components are segmented automatically using generalized gradient vector flow active contours. We investigated 10 human femoral arteries. The effects of balloon compression on the individual artery components is particularly described for two obstructed arteries with an intact collagenous cap, a pronounced lipid pool and with calcification. In both arteries we observe a significant increase in lumen area after angioplasty. Dissection between intima and media and reduction of the lipid pool are primary mechanisms of dilatation. This methodology provides a basis for studying plaque biomechanics under supra-physiological loading conditions. It has the potential to improve and validate finite element models of atherosclerotic plaques which may allow a better prediction of angioplasty procedures.

A numerical model to study the interaction of vascular stents with human atherosclerotic lesions

A methodology is proposed that identifies optimal stent devices for specific clinical criteria. It enables to predict the effect of stent designs on the mechanical environment of stenotic arteries. In particular, we present a numerical study which is based on the interaction of a vascular stent with a patient-specific, atherosclerotic human iliac lesion of type V. The stress evolution in four different tissue components during and after stenting is investigated. The geometric model of the artery is obtained through MRI, while anisotropic material models are applied to describe the behavior of tissues at finite strains. In order to model the observed fissuring and dissection of the plaque under dilation, the undeformed configuration of the arterial wall incorporates two initial tears. The 3D balloon-stent-artery interaction problem is modeled by means of a contact algorithm, which is based on a C(2)-continuous surface parametrization, hence avoiding numerical instabilities of standard facet-based techniques. In the simulations three different stent designs are studied. The performance of each stent is characterized by scalar quantities relating to stress changes in the artery, contact forces, and changes in lumen area after stenting. The study concludes by suggesting two optimal stent designs for two different clinically relevant parameters.

Modeling Plaque Fissuring and Dissection during Balloon Angioplasty Intervention

Balloon angioplasty intervention is traumatic to arterial tissue. Fracture mechanisms such as plaque fissuring and/or dissection occur and constitute major contributions to the lumen enlargement. However, these types of mechanically-based traumatization of arterial tissue are also contributing factors to both acute procedural complications and chronic restenosis of the treatment site. We propose physical and finite element models, which are generally useable to trace fissuring and/or dissection in atherosclerotic plaques during balloon angioplasty interventions. The arterial wall is described as an anisotropic, heterogeneous, highly deformable, nearly incompressible body, whereas tissue failure is captured by a strong discontinuity kinematics and a novel cohesive zone model. The numerical implementation is based on the partition of unity finite element method and the interface element method. The later is used to link together meshes of the different tissue components. The balloon angioplasty-based failure mechanisms are numerically studied in 3D by means of an atherosclerotic-prone human external iliac artery, with a type V lesion. Image-based 3D geometry is generated and tissue-specific material properties are considered. Numerical results show that in a primary phase the plaque fissures at both shoulders of the fibrous cap and stops at the lamina elastica interna. In a secondary phase, local dissections between the intima and the media develop at the fibrous cap location with the smallest thickness. The predicted results indicate that plaque fissuring and dissection cause localized mechanical trauma, but prevent the main portion of the stenosis from high stress, and hence from continuous tissue damage.

Computational stress-deformation analysis of arterial walls including high-pressure response

BACKGROUND

Changes in the mechanical behavior of arteries after balloon angioplasty cause cell reactions that may be responsible for restenosis. Hence, the study of the stress-deformation changes in arterial walls following supraphysiological tissue loading is an essential task.

METHODS

A normal LAD coronary artery was modeled and computationally analyzed as a two-layer, thick-walled, anisotropic and inelastic circular tube including residual strains. Each layer was treated as a fibre-matrix composite. The tube was subjected to an axial stretch of 1.1 and a transmural pressure of 750 mm Hg. Since overstretch of remnant non-diseased tissue in lesions is a primary mechanism of lumen enlargement this model approach represents a reasonable first step.

RESULTS

At physiological loading, the residual stresses led to a significant reduction of the high circumferential stress values at the inner wall, and the stress gradients. At low pressure level the media was the mechanically relevant layer, while at supraphysiological loading, the adventitia was the predominant load-carrying constituent providing a stiff support for 'redistribution' of soft plaque components by means of radial compression. After unloading to physiological loading conditions the stress state in the arterial wall differed significantly from that before inflation; the stress gradient in the media even changed its sign. Complete unloading indicated lumen enlargement, material softening and energy dissipation, which is in agreement with experimental studies.

CONCLUSIONS

This method may be useful to improve interventional protocols for reducing the dilatational trauma, and thereby the adverse biological reaction in arterial walls following balloon angioplasty.

3D geometric reconstruction of thoracic aortic aneurysms

Background

The thoracic aortic aneurysm (TAA) is a pathology that involves an expansion of the aortic diameter in the thoracic aorta, leading to risk of rupture. Recent studies have suggested that internal wall stress, which is affected by TAA geometry and the presence or absence of thrombus, is a more reliable predictor of rupture than the maximum diameter, the current clinical criterion. Accurate reconstruction of TAA geometry is a crucial step in patient-specific stress calculations.

Methods

In this work, a novel methodology was developed, which combines data from several sets of magnetic resonance (MR) images with different levels of detail and different resolutions. Two sets of images were employed to create the final model, which has the highest level of detail for each component of the aneurysm (lumen, thrombus, and wall). A reference model was built by using a single set of images for comparison. This approach was applied to two patient-specific TAAs in the descending thoracic aorta.

Results

The results of finite element simulations showed differences in stress pattern between the coarse and fine models: higher stress values were found with the coarse model and the differences in predicted maximum wall stress were 30% for patient A and 11% for patient B.

Conclusion

This paper presents a new approach to the reconstruction of an aneurysm model based on the use of several sets of MR images. This enables more accurate representation of not only the lumen but also the wall surface of a TAA taking account of intraluminal thrombus.