Estimation of nonlinear mechanical properties of vascular tissues via elastography

Mechanical and histological characteristics of aortic dissection tissues

AIMS

This study investigated the association between the macroscopic mechanical response of aortic dissection (AoD) flap, its fibre features, and patient physiological features and clinical presentations.

METHODS

Uniaxial test was performed with tissue strips in both circumferential and longitudinal directions from 35 patients with (AoD:CC) and without (AoD:w/oCC) cerebral/coronary complications, and 19 patients with rheumatic or valve-related heart diseases (RH). A Bayesian inference framework was used to estimate the expectation of material constants (C₁, D₁, and D₂) of the modified Mooney-Rivlin strain energy density function. Histological examination was used to visualise the elastin and collagen in the tissue strips and image processing was performed to quantify their area percentages, fibre misalignment and waviness.

RESULTS

The elastin area percentage was negatively associated with age (p = 0.008), while collagen increased about 6% from age 40 to 70 (p = 0.03). Elastin fibre was less dispersed and wavier in old patients and no significant association was found between patient age and collagen fibre dispersion or waviness. Features of fibrous microstructures, either elastin or collagen, were comparable between AoD:CC and AoD:w/oCC group. Elastin and collagen area percentages were positively correlated with C₁ and D₂, respectively, while the elastin and collagen waviness were negatively correlated with C₁ and D₂, respectively. Elastin dispersion was negatively correlated to D₂. Multivariate analysis showed that D₂ was an effective parameter which could differentiate patient groups with different age and clinical presentations, as well as the direction of tissue strip.

CONCLUSION

Fibre dispersion and waviness in the aortic dissection flap changed with patient age and clinical presentations, and these can be captured by the material constants in the strain energy density function.

STATEMENT OF SIGNIFICANCE

Aortic dissection (AoD) is a severe cardiovascular disease. Understanding the mechanical property of intimal flap is essential for its risk evaluation. In this study, mechanical testing and histology examination were combined to quantify the relationship between mechanical presentations and microstructure features. A Bayesian inference framework was employed to estimate the expectation of the material constants in the modified Mooney-Rivlin constitutive equation. It was found that fibre dispersion and waviness in the AoD flap changed with patient age and clinical presentations, and these could be captured by the material constants. This study firstly demonstrated that the expectation of material constants can be used to characterise tissue microstructures and differentiate patients with different clinical presentations.

An inverse method for mechanical characterization of heterogeneous diseased arteries using intravascular imaging

The increasing prevalence of finite element (FE) simulations in the study of atherosclerosis has spawned numerous inverse FE methods for the mechanical characterization of diseased tissue in vivo. Current approaches are however limited to either homogenized or simplified material representations. This paper presents a novel method to account for tissue heterogeneity and material nonlinearity in the recovery of constitutive behavior using imaging data acquired at differing intravascular pressures by incorporating interfaces between various intra-plaque tissue types into the objective function definition. Method verification was performed in silico by recovering assigned material parameters from a pair of vessel geometries: one derived from coronary optical coherence tomography (OCT); one generated from in silico-based simulation. In repeated tests, the method consistently recovered 4 linear elastic (0.1 ± 0.1% error) and 8 nonlinear hyperelastic (3.3 ± 3.0% error) material parameters. Method robustness was also highlighted in noise sensitivity analysis, where linear elastic parameters were recovered with average errors of 1.3 ± 1.6% and 8.3 ± 10.5%, at 5% and 20% noise, respectively. Reproducibility was substantiated through the recovery of 9 material parameters in two more models, with mean errors of 3.0 ± 4.7%. The results highlight the potential of this new approach, enabling high-fidelity material parameter recovery for use in complex cardiovascular computational studies.

Brillouin confocal microscopy to determine biomechanical properties of SULEEI-treated bovine pericardium for application in cardiac surgery

BACKGROUND

Heart valves are exposed to a highly dynamic environment and underlie high tensile and shear forces during opening and closing. Therefore, analysis of mechanical performance of novel heart valve bioprostheses materials, like SULEEI-treated bovine pericardium, is essential and usually carried out by uniaxial tensile tests. Nevertheless, major drawbacks are the unidirectional strain, which does not reflect the in vivo condition and the deformation of the sample material. An alternative approach for measurement of biomechanical properties is offered by Brillouin confocal microscopy (BCM), a novel, non-invasive and three-dimensional method based on the interaction of light with acoustic waves.

OBJECTIVE

BCM is a powerful tool to determine viscoelastic tissue properties and is, for the first time, applied to characterize novel biological graft materials, such as SULEEI-treated bovine pericardium. Therefore, the method has to be validated as a non-invasive alternative to conventional uniaxial tensile tests.

METHODS

Vibratome sections of SULEEI-treated bovine pericardium (decellularized, riboflavin/UV-cross-linked and low-energy electron irradiated) as well as native and GA-fixed controls (n = 3) were analyzed by BCM. In addition, uniaxial tensile tests were performed on equivalent tissue samples and Young's modulus as well as length of toe region were analyzed from stress-strain diagrams. The structure of the extracellular matrix (ECM), especially collagen and elastin, was investigated by multiphoton microscopy (MPM).

RESULTS

SULEEI-treated pericardium exhibited a significantly higher Brillouin shift and hence higher tissue stiffness in comparison to native and GA-fixed controls (native: 5.6±0.2 GHz; GA: 5.5±0.1 GHz; SULEEI: 6.3±0.1 GHz; n = 3, p < 0.0001). Similarly, a significantly higher Young's modulus was detected in SULEEI-treated pericardia in comparison to native tissue (native: 30.0±10.4 MPa; GA: 31.8±10.7 MPa; SULEEI: 42.1±7.0 MPa; n = 3, p = 0.027). Native pericardia showed wavy and non-directional collagen fibers as well as thin, linear elastin fibers generating a loose matrix. The fibers of GA-fixed and SULEEI-treated pericardium were aligned in one direction, whereat the SULEEI-sample exhibited a much denser matrix.

CONCLUSION

BCM is an innovative and non-invasive method to analyze elastic properties of novel pericardial graft materials with special mechanical requirements, like heart valve bioprostheses.

A multiphysics approach for modeling early atherosclerosis

This work is devoted to the development of a mathematical model of the early stages of atherosclerosis incorporating processes of all time scales of the disease and to show their interactions. The cardiovascular mechanics is modeled by a fluid-structure interaction approach coupling a non-Newtonian fluid to a hyperelastic solid undergoing anisotropic growth and a change of its constitutive equation. Additionally, the transport of low-density lipoproteins and its penetration through the endothelium is considered by a coupled set of advection-diffusion-reaction equations. Thereby, the permeability of the endothelium is wall-shear stress modulated resulting in a locally varying accumulation of foam cells triggering a novel growth and remodeling formulation. The model is calibrated and applied to an murine-specific case study, and a qualitative validation of the computational results is performed. The model is utilized to further investigate the influence of the pulsatile blood flow and the compliance of the artery wall to the atherosclerotic process. The computational results imply that the pulsatile blood flow is crucial, whereas the compliance of the aorta has only a minor influence on atherosclerosis. Further, it is shown that the novel model is capable to produce a narrowing of the vessel lumen inducing an adaption of the endothelial permeability pattern.

Optical coherence elastography - OCT at work in tissue biomechanics [Invited]

Optical coherence elastography (OCE), as the use of OCT to perform elastography has come to be known, began in 1998, around ten years after the rest of the field of elastography - the use of imaging to deduce mechanical properties of tissues. After a slow start, the maturation of OCT technology in the early to mid 2000s has underpinned a recent acceleration in the field. With more than 20 papers published in 2015, and more than 25 in 2016, OCE is growing fast, but still small compared to the companion fields of cell mechanics research methods, and medical elastography. In this review, we describe the early developments in OCE, and the factors that led to the current acceleration. Much of our attention is on the key recent advances, with a strong emphasis on future prospects, which are exceptionally bright.

Nonlinear characterization of elasticity using quantitative optical coherence elastography

Optical coherence elastography (OCE) has been used to perform mechanical characterization on biological tissue at the microscopic scale. In this work, we used quantitative optical coherence elastography (qOCE), a novel technology we recently developed, to study the nonlinear elastic behavior of biological tissue. The qOCE system had a fiber-optic probe to exert a compressive force to deform tissue under the tip of the probe. Using the space-division multiplexed optical coherence tomography (OCT) signal detected by a spectral domain OCT engine, we were able to simultaneously quantify the probe deformation that was proportional to the force applied, and to quantify the tissue deformation. In other words, our qOCE system allowed us to establish the relationship between mechanical stimulus and tissue response to characterize the stiffness of biological tissue. Most biological tissues have nonlinear elastic behavior, and the apparent stress-strain relationship characterized by our qOCE system was nonlinear an extended range of strain, for a tissue-mimicking phantom as well as biological tissues. Our experimental results suggested that the quantification of force in OCE was critical for accurate characterization of tissue mechanical properties and the qOCE technique was capable of differentiating biological tissues based on the elasticity of tissue that is generally nonlinear.

A Framework for Local Mechanical Characterization of Atherosclerotic Plaques: Combination of Ultrasound Displacement Imaging and Inverse Finite Element Analysis

Biomechanical models have the potential to predict plaque rupture. For reliable models, correct material properties of plaque components are a prerequisite. This study presents a new technique, where high resolution ultrasound displacement imaging and inverse finite element (FE) modeling is combined, to estimate material properties of plaque components. Iliac arteries with plaques were excised from 6 atherosclerotic pigs and subjected to an inflation test with pressures ranging from 10 to 120 mmHg. The arteries were imaged with high frequency 40 MHz ultrasound. Deformation maps of the plaques were reconstructed by cross correlation of the ultrasound radiofrequency data. Subsequently, the arteries were perfusion fixed for histology and structural components were identified. The histological data were registered to the ultrasound data to construct FE model of the plaques. Material properties of the arterial wall and the intima of the atherosclerotic plaques were estimated using a grid search method. The computed displacement fields showed good agreement with the measured displacement fields, implying that the FE models were able to capture local inhomogeneities within the plaque. On average, nonlinear stiffening of both the wall and the intima was observed, and the wall of the atheroslcerotic porcine iliac arteries was markedly stiffer than the intima (877 ± 459 vs. 100 ± 68 kPa at 100 mmHg). The large spread in the data further illustrates the wide variation of the material properties. We demonstrated the feasibility of a mixed experimental–numerical framework to determine the material properties of arterial wall and intima of atherosclerotic plaques from intact arteries, and concluded that, due to the observed variation, plaque specific properties are required for accurate stress simulations.

Influence of material property variability on the mechanical behaviour of carotid atherosclerotic plaques: a 3D fluid-structure interaction analysis

Mechanical analysis has been shown to be complementary to luminal stenosis in assessing atherosclerotic plaque vulnerability. However, patient-specific material properties are not available and the effect of material properties variability has not been fully quantified. Media and fibrous cap (FC) strips from carotid endarterectomy samples were classified into hard, intermediate and soft according to their incremental Young's modulus. Lipid and intraplaque haemorrhage/thrombus strips were classified as hard and soft. Idealised geometry-based 3D fluid-structure interaction analyses were performed to assess the impact of material property variability in predicting maximum principal stress (Stress-P1 ) and stretch (Stretch-P1 ). When FC was thick (1000 or 600 µm), Stress-P1 at the shoulder was insensitive to changes in material stiffness, whereas Stress-P1 at mid FC changed significantly. When FC was thin (200 or 65 µm), high stress concentrations shifted from the shoulder region to mid FC, and Stress-P1 became increasingly sensitive to changes in material properties, in particular at mid FC. Regardless of FC thickness, Stretch-P1 at these locations was sensitive to changes in material properties. Variability in tissue material properties influences both the location and overall stress/stretch value. This variability needs to be accounted for when interpreting the results of mechanical modelling.

Carotid plaque elasticity estimation using ultrasound elastography, MRI, and inverse FEA - A numerical feasibility study

The material properties of atherosclerotic plaques govern the biomechanical environment, which is associated with rupture-risk. We investigated the feasibility of noninvasively estimating carotid plaque component material properties through simulating ultrasound (US) elastography and in vivo magnetic resonance imaging (MRI), and solving the inverse problem with finite element analysis. 2D plaque models were derived from endarterectomy specimens of nine patients. Nonlinear neo-Hookean models (tissue elasticity C1) were assigned to fibrous intima, wall (i.e., media/adventitia), and lipid-rich necrotic core. Finite element analysis was used to simulate clinical cross-sectional US strain imaging. Computer-simulated, single-slice in vivo MR images were segmented by two MR readers. We investigated multiple scenarios for plaque model elasticity, and consistently found clear separations between estimated tissue elasticity values. The intima C1 (160 kPa scenario) was estimated as 125.8 ± 19.4 kPa (reader 1) and 128.9 ± 24.8 kPa (reader 2). The lipid-rich necrotic core C1 (5 kPa) was estimated as 5.6 ± 2.0 kPa (reader 1) and 8.5 ± 4.5 kPa (reader 2). A scenario with a stiffer wall yielded similar results, while realistic US strain noise and rotating the models had little influence, thus demonstrating robustness of the procedure. The promising findings of this computer-simulation study stimulate applying the proposed methodology in a clinical setting.

Practical Challenges of Current Video Rate OCT Elastography: Accounting for Dynamic and Static Tissue Properties

Optical coherence tomography (OCT) elastography (OCTE) has the potential to be an important diagnostic tool for pathologies including coronary artery disease, osteoarthritis, malignancies, and even dental caries. Many groups have performed OCTE, including our own, using a wide range of approaches. However, we will demonstrate current OCTE approaches are not scalable to real-time, in vivo imaging. As will be discussed, among the most important reasons is current designs focus on the system and not the target. Specifically, tissue dynamic responses are not accounted, with examples being the tissue strain response time, preload variability, and conditioning variability. Tissue dynamic responses, and to a lesser degree static tissue properties, prevent accurate video rate modulus assessments for current embodiments. Accounting for them is the focus of this paper. A top-down approach will be presented to overcome these challenges to real time in vivo tissue characterization. Discussed first is an example clinical scenario where OTCE would be of substantial relevance, the prevention of acute myocardial infarction or heart attacks. Then the principles behind OCTE are examined. Next, constrains on in vivo application of current OCTE are evaluated, focusing on dynamic tissue responses. An example is the tissue strain response, where it takes about 20 msec after a stress is applied to reach plateau. This response delay is not an issue at slow acquisition rates, as most current OCTE approaches are preformed, but it is for video rate OCTE. Since at video rate each frame is only 30 msec, for essentially all current approaches this means the strain for a given stress is changing constantly during the B-scan. Therefore the modulus can't be accurately assessed. This serious issue is an even greater problem for pulsed techniques as it means the strain/modulus for a given stress (at a location) is unpredictably changing over a B-scan. The paper concludes by introducing a novel video rate approach to overcome these challenges.

Lipid crystals mechanically stimulate adjacent extracellular matrix in advanced atherosclerotic plaques

OBJECTIVE

Although lipid crystals (LCs) have received attention as a causative factor of plaque rupture, the mechanisms by which they increase plaque vulnerability are unknown. We examined whether solid-state LCs physically affect the adjacent extracellular matrix (ECM) using a combination of multimodal nonlinear optical (MNLO) imaging and finite element analysis (FEA).

METHODS

The changes of ECMs affected by lipids in atherosclerotic arteries in apolipoprotein E-deficient mice (n = 32) fed a high-fat diet for 20-30 weeks were micro-anatomically visualized by a 3D MNLO imaging platform including CARS for lipids, TPEF for elastin, and SHG for collagen.

RESULTS AND CONCLUSION

The TPEF signal of elastin was increased at the peripheral regions of LCs (<10 μm) compared with foam cell regions. In order to confirm the increase of elastin, biochemical assay (western blot) was performed. The protein level of elastin was increased approximately 2.25-fold (p = 0.024) in LC-rich arteries. Under the hypothesis that the increase of elastin resulted from the mechanical stimulus from solid-state LCs, MNLO images were subjected to FEA to simulate the displacement according to the expanding magnitude of the vessel during cardiac cycles. We found that microscale focal stress was increased specifically around the LCs. These FEA results corresponded with the increase of elastin observed by TPEF. These data suggest that LCs mechanically stimulate the adjacent ECM to alter the composition of ECM and cause vessel remodeling. The combination of MNLO imaging and FEA has great potential to verify the mechanical predictions in cardiovascular diseases.

Vascular elastography: a validation study

Vascular elastography techniques are promising tools for mechanical characterization of diseased arteries. These techniques are usually validated with simulations or phantoms or by comparing results with histology or other imaging modalities. In the study described here, vascular elastography was applied to porcine aortas in vitro during inflation testing (n = 10) and results were compared with those of standard bi-axial tensile testing, a technique that directly measures the force applied to the tissue. A neo-Hookean model was fit to the stress-strain data, valid for large deformations. Results indicated good correspondence between the two techniques, with GUS = 110 ± 11 kPa and GTT = 108 ± 10 kPa for ultrasound and tensile testing, respectively. Bland-Altman analysis revealed little bias (GUS-GTT = 2 ± 20 kPa). The next step will be the application of a non-linear material model that is also adaptable for in vivo measurements.

Plaque hemorrhage in carotid artery disease: pathogenesis, clinical and biomechanical considerations

Stroke remains the most prevalent disabling illness today, with internal carotid artery luminal stenosis due to atheroma formation responsible for the majority of ischemic cerebrovascular events. Severity of luminal stenosis continues to dictate both patient risk stratification and the likelihood of surgical intervention. But there is growing evidence to suggest that plaque morphology may help improve pre-existing risk stratification criteria. Plaque components such a fibrous tissue, lipid rich necrotic core and calcium have been well investigated but plaque hemorrhage (PH) has been somewhat overlooked. In this review we discuss the pathogenesis of PH, its role in dictating plaque vulnerability, PH imaging techniques, marterial properties of atherosclerotic tissues, in particular, those obtained based on in vivo measurements and effect of PH in modulating local biomechanics.

Methods for robust in vivo strain estimation in the carotid artery

A hierarchical block-matching motion tracking algorithm for strain imaging is presented. Displacements are estimated with improved robustness and precision by utilizing a Bayesian regularization algorithm and an unbiased subsample interpolation technique. A modified least-squares strain estimator is proposed to estimate strain images from a noisy displacement input while addressing the motion discontinuity at the wall-lumen boundary. Methods to track deformation over the cardiac cycle incorporate a dynamic frame skip criterion to process data frames with sufficient deformation to produce high signal-to-noise displacement and strain images. Algorithms to accumulate displacement and/or strain on particles in a region of interest over the cardiac cycle are described. New methods to visualize and characterize the deformation measured with the full 2D strain tensor are presented. Initial results from patients imaged prior to carotid endarterectomy suggest that strain imaging detects conditions that are traditionally considered high risk including soft plaque composition, unstable morphology, abnormal hemodynamics and shear of plaque against tethering tissue can be exacerbated by neoangiogenesis. For example, a maximum absolute principal strain exceeding 0.2 is observed near calcified regions adjacent to turbulent flow, protrusion of the plaque into the arterial lumen and regions of low echogenicity associated with soft plaques. Non-invasive carotid strain imaging is therefore a potentially useful tool for detecting unstable carotid plaque.

A Quasi-Nonlinear Analysis of the Anisotropic Behaviour of Human Gallbladder Wall

Estimation of biomechanical parameters of soft tissues from noninvasive measurements has clinical significance in patient-specific modeling and disease diagnosis. In this work, we present a quasi-nonlinear method that is used to estimate the elastic moduli of the human gallbladder wall. A forward approach based on a transversely isotropic membrane material model is used, and an inverse iteration is carried out to determine the elastic moduli in the circumferential and longitudinal directions between two successive ultrasound images of gallbladder. The results demonstrate that the human gallbladder behaves in an anisotropic manner, and constitutive models need to incorporate this. The estimated moduli are also nonlinear and patient dependent. Importantly, the peak stress predicted here differs from the earlier estimate from linear membrane theory. As the peak stress inside the gallbladder wall has been found to strongly correlate with acalculous gallbladder pain, reliable mechanical modeling for gallbladder tissue is crucial if this information is to be used in clinical diagnosis.

Investigating in vivo airway wall mechanics during tidal breathing with optical coherence tomography

Optical coherence tomography (OCT) is a nondestructive imaging technique offering high temporal and spatial resolution, which makes it a natural choice for assessing tissue mechanical properties. We have developed methods to mechanically analyze the compliance of the rabbit trachea in vivo using tissue deformations induced by tidal breathing, offering a unique tool to assess the behavior of the airways during their normal function. Four-hundred images were acquired during tidal breathing with a custom-built endoscopic OCT system. The surface of the tissue was extracted from a set of these images via image processing algorithms, filtered with a bandpass filter set at respiration frequency to remove cardiac and probe motion, and compared to ventilatory pressure to calculate wall compliance. These algorithms were tested on elastic phantoms to establish reliability and reproducibility. The mean tracheal wall compliance (in five animals) was 1.3±0.3×10(-5) (mm Pa)(-1). Unlike previous work evaluating airway mechanics, this new method is applicable in vivo, noncontact, and loads the trachea in a physiological manner. The technique may have applications in assessing airway mechanics in diseases such as asthma that are characterized by significant airway remodeling.

Experimental determination of circumferential properties of fresh carotid artery plaques

Carotid endarterectomy (CEA) is currently accepted as the gold standard for interventional revascularisation of diseased arteries belonging to the carotid bifurcation. Despite the proven efficacy of CEA, great interest has been generated in carotid angioplasty and stenting (CAS) as an alternative to open surgical therapy. CAS is less invasive compared with CEA, and has the potential to successfully treat lesions close to the aortic arch or distal internal carotid artery (ICA). Following promising results from two recent trials (CREST; Carotid revascularisation endarterectomy versus stenting trial, and ICSS; International carotid stenting study) it is envisaged that there will be a greater uptake in carotid stenting, especially amongst the group who do not qualify for open surgical repair, thus creating pressure to develop computational models that describe a multitude of plaque models in the carotid arteries and their reaction to the deployment of such interventional devices. Pertinent analyses will require fresh human atherosclerotic plaque material characteristics for different disease types. This study analysed atherosclerotic plaque characteristics from 18 patients tested on site, post-surgical revascularisation through endarterectomy, with 4 tissue samples being excluded from tensile testing based on large width-length ratios. According to their mechanical behaviour, atherosclerotic plaques were separated into 3 grades of stiffness. Individual and group material coefficients were then generated analytically using the Yeoh strain energy function. The ultimate tensile strength (UTS) of each sample was also recorded, showing large variation across the 14 atherosclerotic samples tested. Experimental Green strains at rupture varied from 0.299 to 0.588 and the Cauchy stress observed in the experiments was between 0.131 and 0.779 MPa. It is expected that this data may be used in future design optimisation of next generation interventional medical devices for the treatment and revascularisation of diseased arteries of the carotid bifurcation.

Optical coherence elastography: current status and future applications

Optical coherence tomography (OCT) has several advantages over other imaging modalities, such as angiography and ultrasound, due to its inherently high in vivo resolution, which allows for the identification of morphological tissue structures. Optical coherence elastography (OCE) benefits from the superior spatial resolution of OCT and has promising applications, including cancer diagnosis and the detailed characterization of arterial wall biomechanics, both of which are based on the elastic properties of the tissue under investigation. We present OCE principles based on techniques associated with static and dynamic tissue excitation, and their corresponding elastogram image-reconstruction algorithms are reviewed. OCE techniques, including the development of intravascular- or catheter-based OCE, are in their early stages of development but show great promise for surgical oncology or intravascular cardiology applications.

DYNAMIC OPTICAL COHERENCE ELASTOGRAPHY: A REVIEW

With the development of optical coherence tomography, the application optical coherence elastography (OCE) has gained more and more attention in biomechanics for its unique features including micron-scale resolution, real-time processing, and non-invasive imaging. In this review, one group of OCE techniques, namely dynamic OCE, are introduced and discussed including external dynamic OCE mapping and imaging of ex vivo breast tumor, external dynamic OCE measurement of in vivo human skin, and internal dynamic OCE including acoustomotive OCE and magnetomotive OCE. These techniques overcame some of the major drawbacks of traditional static OCE, and broadened the OCE application fields. Driven by scientific needs to engineer new quantitative methods that utilize the high micron-scale resolution achievable with optics, results of biomechanical properties were obtained from biological tissues. The results suggest potential diagnostic and therapeutic clinical applications. Results from these studies also help our understanding of the relationship between biomechanical variations and functional tissue changes in biological systems.

On the potential of a new IVUS elasticity modulus imaging approach for detecting vulnerable atherosclerotic coronary plaques: in vitro vessel phantom study

Peak cap stress amplitude is recognized as a good indicator of vulnerable plaque (VP) rupture. However, such stress evaluation strongly relies on a precise, but still lacking, knowledge of the mechanical properties exhibited by the plaque components. As a first response to this limitation, our group recently developed, in a previous theoretical study, an original approach, called iMOD (imaging modulography), which reconstructs elasticity maps (or modulograms) of atheroma plaques from the estimation of strain fields. In the present in vitro experimental study, conducted on polyvinyl alcohol cryogel arterial phantoms, we investigate the benefit of coupling the iMOD procedure with the acquisition of intravascular ultrasound (IVUS) measurements for detection of VP. Our results show that the combined iMOD-IVUS strategy: (1) successfully detected and quantified soft inclusion contours with high positive predictive and sensitivity values of 89.7 ± 3.9% and 81.5 ± 8.8%, respectively, (2) estimated reasonably cap thicknesses larger than ∼300 µm, but underestimated thinner caps, and (3) quantified satisfactorily Young's modulus of hard medium (mean value of 109.7 ± 23.7 kPa instead of 145.4 ± 31.8 kPa), but overestimated the stiffness of soft inclusions (mean Young`s moduli of 31.4 ± 9.7 kPa instead of 17.6 ± 3.4 kPa). All together, these results demonstrate a promising benefit of the new iMOD-IVUS clinical imaging method for in vivo VP detection.

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.

Magnetic resonance elastography: Inversions in bounded media

Magnetic resonance elastography is a noninvasive imaging technique capable of quantifying and spatially resolving the shear stiffness of soft tissues by visualization of synchronized mechanical wave displacement fields. However, magnetic resonance elastography inversions generally assume that the measured tissue motion consists primarily of shear waves propagating in a uniform, infinite medium. This assumption is not valid in organs such as the heart, eye, bladder, skin, fascia, bone and spinal cord, in which the shear wavelength approaches the geometric dimensions of the object. The aim of this study was to develop and test mathematical inversion algorithms capable of resolving shear stiffness from displacement maps of flexural waves propagating in bounded media such as beams, plates, and spherical shells, using geometry-specific equations of motion. Magnetic resonance elastography and finite element modeling of beam, plate, and spherical shell phantoms of various geometries were performed. Mechanical testing of the phantoms agreed with the stiffness values obtained from finite element modeling and magnetic resonance elastography data, and a linear correlation of r(2) >or= 0.99 was observed between the stiffness values obtained using magnetic resonance elastography and finite element modeling data. In conclusion, we have demonstrated new inversion methods for calculating shear stiffness that may be more appropriate for waves propagating in bounded media.