Carotid plaque stress analysis by fluid structure interaction based on in-vivo MRI: Implications to plaque vulnerability assessment

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University, 2010.Stroke is one of the leading causes of death in the world, resulting mostly from the sudden rupture of atherosclerotic plaques. From a biomechanical view, plaque rupture can be considered as a mechanical failure caused by extremely high plaque stress. In this PhD project, we are aiming to predict 3D plaque stress based on in-vivo MRI by using fluid structure interaction (FSI) method, and provide information for plaque rupture risk assessment. Fluid structure interaction was implemented with ANSYS 11.0, followed by a parameter study on fibrous cap thickness and lipid core size with realistic carotid plaque geometry. Twenty patients with carotid plaques imaged by in-vivo MRI were provided in the project. A framework of reconstructing 3D plaque geometry from in-vivo multispectral MRI was designed. The followed reproducibility study on plaque geometry reconstruction procedure and its effect on plaque stress analysis filled the gap in the literature on imaging based plaque stress modeling. The results demonstrated that current MRI technology can provide sufficient information for plaque structure characterization; however stress analysis result is highly affected by MRI resolution and quality. The application of FSI stress analysis to 4 patients with different plaque burdens has showed that the whole procedure from plaque geometry reconstruction to FSI stress analysis was applicable. In the study, plaque geometries from three patients with recent transient ischemic attack were reconstructed by repairing ruptured fibrous cap. The well correlated relationship between local stress concentrations and plaque rupture sites indicated that extremely high plaque stress could be a factor responsible for plaque rupture. Based on the 20 reconstructed carotid plaques from two groups (symptomatic and asymptomatic), fully coupled fluid structure interaction was performed. It was found that there is a significant difference between symptomatic and asymptomatic patients in plaque stress levels, indicating plaque stress could be used as one of the factors for plaque vulnerability assessment. A corresponding plaque morphological feature study showed that plaque stress is significantly affected by fibrous cap thickness, lipid core size and fibrous cap surface irregularities (curvedness). A procedure was proposed for predicting plaque stress by using fibrous cap thickness and curvedness, which requires much less computational time, and has the potential for clinical routine application. The effects of residual stress on plaque stress analysis and arterial wall material property characterization by using in-vivo MRI data were also discussed for patient specific modeling. As the further development, histological study of plaque sample has been combined with conventional plaque stress analysis by assigning material properties to each computational element, based on the data from histological analysis. This method could bridge the gap between biochemistry and biomechanical study of atherosclerosis plaques. In conclusion, extreme stress distributions in the plaque region can be predicted by modern numerical methods, and used for plaque rupture risk assessment, which will be helpful in clinical practice. The combination of plaque MR imaging analysis, computational modelling, and clinical study/ validation would advance our understandings of plaque rupture, prediction of future rupture, and establish new procedures for patient diagnose, management, and treatment.Financial Support was obtained from British Heart Foundation, Brunel Institute for Bioengineering and Brunel Graduate School

Characterization of carotid artery plaques using noninvasive vascular ultrasound elastography

L'athérosclérose est une maladie vasculaire complexe qui affecte la paroi des artères (par l'épaississement) et les lumières (par la formation de plaques). La rupture d'une plaque de l'artère carotide peut également provoquer un accident vasculaire cérébral ischémique et des complications. Bien que plusieurs modalités d'imagerie médicale soient actuellement utilisées pour évaluer la stabilité d'une plaque, elles présentent des limitations telles que l'irradiation, les propriétés invasives, une faible disponibilité clinique et un coût élevé. L'échographie est une méthode d'imagerie sûre qui permet une analyse en temps réel pour l'évaluation des tissus biologiques. Il est intéressant et prometteur d’appliquer une échographie vasculaire pour le dépistage et le diagnostic précoces des plaques d’artère carotide. Cependant, les ultrasons vasculaires actuels identifient uniquement la morphologie d'une plaque en termes de luminosité d'écho ou l’impact de cette plaque sur les caractéristiques de l’écoulement sanguin, ce qui peut ne pas être suffisant pour diagnostiquer l’importance de la plaque. La technique d’élastographie vasculaire non-intrusive (« noninvasive vascular elastography (NIVE) ») a montré le potentiel de détermination de la stabilité d'une plaque. NIVE peut déterminer le champ de déformation de la paroi vasculaire en mouvement d’une artère carotide provoqué par la pulsation cardiaque naturelle. En raison des différences de module de Young entre les différents tissus des vaisseaux, différents composants d’une plaque devraient présenter différentes déformations, caractérisant ainsi la stabilité de la plaque. Actuellement, les performances et l’efficacité numérique sous-optimales limitent l’acceptation clinique de NIVE en tant que méthode rapide et efficace pour le diagnostic précoce des plaques vulnérables. Par conséquent, il est nécessaire de développer NIVE en tant qu’outil d’imagerie non invasif, rapide et économique afin de mieux caractériser la vulnérabilité liée à la plaque. La procédure à suivre pour effectuer l’analyse NIVE consiste en des étapes de formation et de post-traitement d’images. Cette thèse vise à améliorer systématiquement la précision de ces deux aspects de NIVE afin de faciliter la prédiction de la vulnérabilité de la plaque carotidienne. Le premier effort de cette thèse a été dédié à la formation d'images (Chapitre 5). L'imagerie par oscillations transversales a été introduite dans NIVE. Les performances de l’imagerie par oscillations transversales couplées à deux estimateurs de contrainte fondés sur un modèle de déformation fine, soit l’ « affine phase-based estimator (APBE) » et le « Lagrangian speckle model estimator (LSME) », ont été évaluées. Pour toutes les études de simulation et in vitro de ce travail, le LSME sans imagerie par oscillation transversale a surperformé par rapport à l'APBE avec imagerie par oscillations transversales. Néanmoins, des estimations de contrainte principales comparables ou meilleures pourraient être obtenues avec le LSME en utilisant une imagerie par oscillations transversales dans le cas de structures tissulaires complexes et hétérogènes. Lors de l'acquisition de signaux ultrasonores pour la formation d'images, des mouvements hors du plan perpendiculaire au plan de balayage bidimensionnel (2-D) existent. Le deuxième objectif de cette thèse était d'évaluer l'influence des mouvements hors plan sur les performances du NIVE 2-D (Chapitre 6). À cette fin, nous avons conçu un dispositif expérimental in vitro permettant de simuler des mouvements hors plan de 1 mm, 2 mm et 3 mm. Les résultats in vitro ont montré plus d'artefacts d'estimation de contrainte pour le LSME avec des amplitudes croissantes de mouvements hors du plan principal de l’image. Malgré tout, nous avons néanmoins obtenu des estimations de déformations robustes avec un mouvement hors plan de 2.0 mm (coefficients de corrélation supérieurs à 0.85). Pour un jeu de données cliniques de 18 participants présentant une sténose de l'artère carotide, nous avons proposé d'utiliser deux jeux de données d'analyses sur la même plaque carotidienne, soit des images transversales et longitudinales, afin de déduire les mouvements hors plan (qui se sont avérés de 0.25 mm à 1.04 mm). Les résultats cliniques ont montré que les estimations de déformations restaient reproductibles pour toutes les amplitudes de mouvement, puisque les coefficients de corrélation inter-images étaient supérieurs à 0.70 et que les corrélations croisées normalisées entre les images radiofréquences étaient supérieures à 0.93, ce qui a permis de démontrer une plus grande confiance lors de l'analyse de jeu de données cliniques de plaques carotides à l'aide du LSME. Enfin, en ce qui concerne le post-traitement des images, les algorithmes NIVE doivent estimer les déformations des parois des vaisseaux à partir d’images reconstituées dans le but d’identifier les tissus mous et durs. Ainsi, le dernier objectif de cette thèse était de développer un algorithme d'estimation de contrainte avec une résolution de la taille d’un pixel ainsi qu'une efficacité de calcul élevée pour l'amélioration de la précision de NIVE (Chapitre 7). Nous avons proposé un estimateur de déformation de modèle fragmenté (SMSE) avec lequel le champ de déformation dense est paramétré avec des descriptions de transformées en cosinus discret, générant ainsi des composantes de déformations affines (déformations axiales et latérales et en cisaillement) sans opération mathématique de dérivées. En comparant avec le LSME, le SMSE a réduit les erreurs d'estimation lors des tests de simulations, ainsi que pour les mesures in vitro et in vivo. De plus, la faible mise en oeuvre de la méthode SMSE réduit de 4 à 25 fois le temps de traitement par rapport à la méthode LSME pour les simulations, les études in vitro et in vivo, ce qui pourrait permettre une implémentation possible de NIVE en temps réel.Atherosclerosis is a complex vascular disease that affects artery walls (by thickening) and lumens (by plaque formation). The rupture of a carotid artery plaque may also induce ischemic stroke and complications. Despite the use of several medical imaging modalities to evaluate the stability of a plaque, they present limitations such as irradiation, invasive property, low clinical availability and high cost. Ultrasound is a safe imaging method with a real time capability for assessment of biological tissues. It is clinically used for early screening and diagnosis of carotid artery plaques. However, current vascular ultrasound technologies only identify the morphology of a plaque in terms of echo brightness or the impact of the vessel narrowing on flow properties, which may not be sufficient for optimum diagnosis. Noninvasive vascular elastography (NIVE) has been shown of interest for determining the stability of a plaque. Specifically, NIVE can determine the strain field of the moving vessel wall of a carotid artery caused by the natural cardiac pulsation. Due to Young’s modulus differences among different vessel tissues, different components of a plaque can be detected as they present different strains thereby potentially helping in characterizing the plaque stability. Currently, sub-optimum performance and computational efficiency limit the clinical acceptance of NIVE as a fast and efficient method for the early diagnosis of vulnerable plaques. Therefore, there is a need to further develop NIVE as a non-invasive, fast and low computational cost imaging tool to better characterize the plaque vulnerability. The procedure to perform NIVE analysis consists in image formation and image post-processing steps. This thesis aimed to systematically improve the accuracy of these two aspects of NIVE to facilitate predicting carotid plaque vulnerability. The first effort of this thesis has been targeted on improving the image formation (Chapter 5). Transverse oscillation beamforming was introduced into NIVE. The performance of transverse oscillation imaging coupled with two model-based strain estimators, the affine phase-based estimator (APBE) and the Lagrangian speckle model estimator (LSME), were evaluated. For all simulations and in vitro studies, the LSME without transverse oscillation imaging outperformed the APBE with transverse oscillation imaging. Nonetheless, comparable or better principal strain estimates could be obtained with the LSME using transverse oscillation imaging in the case of complex and heterogeneous tissue structures. During the acquisition of ultrasound signals for image formation, out-of-plane motions which are perpendicular to the two-dimensional (2-D) scan plane are existing. The second objective of this thesis was to evaluate the influence of out-of-plane motions on the performance of 2-D NIVE (Chapter 6). For this purpose, we designed an in vitro experimental setup to simulate out-of-plane motions of 1 mm, 2 mm and 3 mm. The in vitro results showed more strain estimation artifacts for the LSME with increasing magnitudes of out-of-plane motions. Even so, robust strain estimations were nevertheless obtained with 2.0 mm out-of-plane motion (correlation coefficients higher than 0.85). For a clinical dataset of 18 participants with carotid artery stenosis, we proposed to use two datasets of scans on the same carotid plaque, one cross-sectional and the other in a longitudinal view, to deduce the out-of-plane motions (estimated to be ranging from 0.25 mm to 1.04 mm). Clinical results showed that strain estimations remained reproducible for all motion magnitudes since inter-frame correlation coefficients were higher than 0.70, and normalized cross-correlations between radiofrequency images were above 0.93, which indicated that confident motion estimations can be obtained when analyzing clinical dataset of carotid plaques using the LSME. Finally, regarding the image post-processing component of NIVE algorithms to estimate strains of vessel walls from reconstructed images with the objective of identifying soft and hard tissues, we developed a strain estimation method with a pixel-wise resolution as well as a high computation efficiency for improving NIVE (Chapter 7). We proposed a sparse model strain estimator (SMSE) for which the dense strain field is parameterized with Discrete Cosine Transform descriptions, thereby deriving affine strain components (axial and lateral strains and shears) without mathematical derivative operations. Compared with the LSME, the SMSE reduced estimation errors in simulations, in vitro and in vivo tests. Moreover, the sparse implementation of the SMSE reduced the processing time by a factor of 4 to 25 compared with the LSME based on simulations, in vitro and in vivo results, which is suggesting a possible implementation of NIVE in real time

Hemodynamics in the Stenosed Carotid Bifurcation with Plaque Ulceration

The presence of irregular plaque surface morphology or ulceration of the atherosclerotic lesion has been identified as an independent risk factor for ischemic stroke. Doppler ultrasound (DUS) is the most commonly performed non-invasive technique used to assess patients suspected of having carotid artery disease, but currently does not incorporate the diagnosis of plaque ulceration. Advanced Doppler analyses incorporating quantitative estimates of flow disturbances may result in diagnostic indices that identify plaque ulcerative conditions. A technique for the fabrication of DUS-compatible flow phantoms was developed, using a direct-machining method that is amenable to comprehensive DUS investigations. In vitro flow studies in an ensemble of matched model vessel geometries determined that ulceration as small as 2 mm can generate significant disturbances in the downstream flow field in a moderately stenosed carotid artery, which are detectable using the DUS velocity-derived parameter turbulence intensity (TI) measured with a clinical system. Further experimental results showed that distal TI was significantly elevated (P \u3c 0.001) due to proximal plaque ulceration in the mild and moderately stenosed carotid bifurcation (30%, 50%, 60% diameter reduction), and also increased with stenosis severity. Pulsatile computational fluid dynamics (CFD) models, with simulated particle tracking, demonstrated enhanced flow disruption of the stenotic jet and slight elevations in path-dependent shear exposure parameters in a stenosed carotid bifurcation model with ulceration. In addition, CFD models were used to evaluate the DUS index TI using finite volume sampling

Evolution and rupture of vulnerable plaques: a review of mechanical effects

Atherosclerosis occurs as a result of the buildup and infiltration of lipid streaks in artery walls, leading to plaques. Understanding the development of atherosclerosis and plaque vulnerability is of critical importance, since plaque rupture can result in heart attack or stroke. Plaques can be divided into two distinct types: those that rupture (vulnerable) and those that are less likely to rupture (stable). In the last few decades, researchers have been interested in studying the influence of the mechanical effects (blood shear stress, pressure forces, and structural stress) on the plaque formation and rupture processes. In the literature, physiological experimental studies are limited by the complexity of in vivo experiments to study such effects, whereas the numerical approach often uses simplified models compared with realistic conditions, so that no general agreement of the mechanisms responsible for plaque formation has yet been reached. In addition, in a large number of cases, the presence of plaques in arteries is asymptomatic. The prediction of plaque rupture remains a complex question to elucidate, not only because of the interaction of numerous phenomena involved in this process (biological, chemical, and mechanical) but also because of the large time scale on which plaques develop. The purpose of the present article is to review the current mechanical models used to describe the blood flow in arteries in the presence of plaques, as well as reviewing the literature treating the influence of mechanical effects on plaque formation, development, and rupture. Finally, some directions of research, including those being undertaken by the authors, are described

Carotid Atheroma Rupture Observed In Vivo and FSI-Predicted Stress Distribution Based on Pre-rupture Imaging

Atherosclerosis at the carotid bifurcation is a major risk factor for stroke. As mechanical forces may impact lesion stability, finite element studies have been conducted on models of diseased vessels to elucidate the effects of lesion characteristics on the stresses within plaque materials. It is hoped that patient-specific biomechanical analyses may serve clinically to assess the rupture potential for any particular lesion, allowing better stratification of patients into the most appropriate treatments. Due to a sparsity of in vivo plaque rupture data, the relationship between various mechanical descriptors such as stresses or strains and rupture vulnerability is incompletely known, and the patient-specific utility of biomechanical analyses is unclear. In this article, we present a comparison between carotid atheroma rupture observed in vivo and the plaque stress distribution from fluid–structure interaction analysis based on pre-rupture medical imaging. The effects of image resolution are explored and the calculated stress fields are shown to vary by as much as 50% with sub-pixel geometric uncertainty. Within these bounds, we find a region of pronounced elevation in stress within the fibrous plaque layer of the lesion with a location and extent corresponding to that of the observed site of plaque rupture

2-D and 3-D high frame-rate Pulse Wave Imaging for the characterization of focal vascular disease

Cardiovascular diseases are major causes of morbidity and mortality in Western-style populations. Atherosclerosis and Abdominal Aortic Aneurysms (AAAs) are two prevalent vascular diseases that may progress without symptoms and contribute to acute cardiovascular events such as stroke and AAA rupture, which are consistently among the leading causes of death worldwide. The imaging methods used in the diagnosis of these diseases, have been reported to present several limitations. Given that both are associated with mechanical changes in the arterial wall, imaging of the arterial mechanical properties may improve early disease detection and patient care. Pulse wave velocity (PWV) refers to the velocity at which arterial waves generated by ventricular ejection travel along the arterial tree. PWV is a surrogate marker of arterial stiffness linked to cardiovascular mortality. The foot-to-foot method that is typically used to calculate PWV suffers from errors of distance measurements and time-delay measurements. Additionally, a single PWV estimate is provided over a relatively long distance, thus inherently lacking the capability to provide regional arterial stiffness measurements. Pulse Wave Imaging (PWI) is a noninvasive, ultrasound-based technique for imaging the propagation of pulse waves along the wall of major arteries and providing a regional PWV value for the imaged artery. The aim of this work was to enable PWI to provide more localized PWV and stiffness measurements within the imaged arterial segment and to further extend it into a 2-D and 3-D technique for the detection and monitoring of focal vascular disease at high temporal and spatial resolution. The improved modality was integrated with blood flow imaging modalities aiming to render PWI a comprehensive methodology for the study of arterial biomechanics in vivo. Spatial information was increased with the introduction of piecewise PWI. This novel technique was used to measure PWV within small sub-regions of the imaged vessel in murine aneurysmal (n = 8) and atherosclerotic aortas (n = 11) in vivo. It provided PWV and stiffness maps while capturing the progressive arterial stiffening caused by atherosclerosis. PWI was further augmented with a sophisticated adaptive algorithm, enabling it to optimally partition the imaged artery into relatively homogeneous segments, automatically isolating arterial stiffness inhomogeneities. Adaptive PWI was validated in silicone phantoms consisting of segments of varying stiffness and then tested in murine aortas in vivo. Subsequently, the conventional tradeoff between spatial and temporal resolution was addressed with a plane wave compounding implementation of PWI, allowing the acquisition of full field of view frames at over 2000 Hz. A GPU-accelerated PWI post-processing framework was developed for the processing of the big bulk of generated data. The parameters of coherent compounding were optimized in vivo. The optimized sequences were then used in the clinic to assess the mechanical properties of atherosclerotic carotids (n=10) and carotids of patients after endarterectomy (n=7), a procedure to remove the plaque and restore blood flow to the brain. In the case of atherosclerotic patients undergoing carotid endarterectomy, the results were compared against the histology of the excised plaques. Investigation of the mechanical properties of plaques was also conducted for the first time with a high-frequency transducer (18.5 MHz). Additionally, 4-D PWI was introduced, utilizing high frame rate 3-D plane wave acquisitions with a 2-D matrix array transducer (16x16 elements, 2.5 MHz). A novel methodology for PWV estimation along the direction of pulse wave propagation was implemented and validated in silicone phantoms. 4-D PWI provided comprehensive views of the pulse wave propagation in a plaque phantom and the results were compared against conventional PWI. Finally, its feasibility was tested in the carotid arteries of healthy human subjects (n=6). PWVs derived in 3-D were within the physiological range and showed good agreement with the results of conventional PWI. Finally, PWI was integrated with flow imaging modalities (Color and Vector Doppler). Thus, full field-of-view, high frame-rate, simultaneous and co-localized imaging of the arterial wall dynamics and color flow as well as 2-D vector flow was implemented. The feasibility of both techniques was tested in healthy subjects (n=6) in vivo. The relationship between the timings of the flow and wall velocities was investigated at multiple locations of the imaged artery. Vector flow velocities were found to be aligned with the vessel’s centerline during peak systole in the common carotid artery and interesting flow patterns were revealed in the case of the carotid bifurcation Consequently, with the aforementioned improvements and the inclusion of 3-D imaging, PWI is expected to provide comprehensive information on the mechanical properties of pathological arteries, providing clinicians with a powerful tool for the early detection of vascular abnormalities undetectable on the B-mode, while also enabling the monitoring of fully developed vascular pathology and of the recovery of post-operated vessels

Novel ultrasound features for the identification of the vulnerable carotid plaque

Background: The identification of the vulnerable carotid plaque is of paramount importance in order to prevent the significant stroke-related mortality and morbidity. Currently the clinical decision-making around this condition is based on the traditional ultrasound evaluation of the degree of stenosis. However, there is emerging evidence supporting that this is not sufficient for all patients. Aim of this thesis: The evaluation of novel carotid plaque features for the characterisation of plaque composition, volume and motion using 2 and 3 dimensional ultrasound technology. The ultimate goal is to identify novel sensitive imaging markers for carotid plaque characterisation and stroke-risk stratification. Methods: The Asymptomatic Carotid Stenosis and Risk of Stroke (ACSRS) Study was a large prospective multicentre trial that was recently completed. A post-hoc analysis of the sonographic and clinical data from this study was performed in order to evaluate the effectiveness of novel ultrasound texture features, such as second order statics, on stroke-risk prediction. In addition, the change of specific texture features and degree of stenosis during the ACSRS follow-up time (8 years) and their importance for stroke prediction was evaluated. In order to assess the potential of 3D ultrasound carotid imaging we also developed a special methodology using a 3D broadband, linear array probe and the Q-lab software. This methodology was then applied in a clinical, cross-sectional study of patients with symptomatic and asymptomatic carotid disease. Finally we developed a carotid plaque motion analysis methodology that we tested on a feasibility study. Results: The post-hoc analysis of more than 1, 000 patients from the ACSRS database showed that there are novel ultrasound features of plaque homogeneity that can contribute to plaque characterisation and improve stroke-risk prediction. Similarly our results suggest that the change of degree of stenosis or plaque’s composition through time might have significant predictive value when combined with the above novel features. The study in 3D ultrasound prospectively assessed more than 80 people with symptomatic and asymptomatic carotid disease with both 2 and 3D carotid ultrasound without, though, revealing any significant benefit from the use of 3D imaging in terms of stroke-risk prediction. Finally, our feasibility study on plaque motion analysis showed that it is possible to objectively characterise plaque motion, using ultrasound and dedicated software without complicated reconstructions. Conclusion: The use of novel 2D ultrasound texture features in combination with traditional ones can improve the stroke-risk stratification. 3D ultrasound is a promising new approach, however, the current technology does not appear to offer a significant benefit in comparison to cheaper traditional 2D ultrasound for carotid plaque evaluation. Further research is warranted on this issue.Open Acces

An Efficient Hemodynamic Workflow in Computational Surgery

For few decades, it has been shown that atherosclerosis is the cause of the majority of clinical cardiovascular diseases including peripheral arterial diseases. The diagnosis and treatment for vascular disease has evolved significantly over the past years considering the rapid advances in imaging technologies. In recent years, computational fluid dynamics has been increasingly used as a simulation tool for blood flows. Numerous researches connect wall shear stress quantities to endovascular diseases such as stenosis, aneurism, and atherosclerosis. A thorough knowledge of vascular anatomy and hemodynamic would be beneficial for understanding the development and progression of the disease, the therapeutic decision process and follow up. The objective of this dissertation is to propose a computational fluid dynamic framework that includes: Understanding how streamline efficiently hemodynamic simulation for main arteries to produce database for clinical study/Providing some confidence estimate on numerical results/Extending the state of the art of clinical study by including motion and particles analysis.Computer Science, Department o