In vivo study of myocardial elastography under graded ischemia conditions

Segmental and transmural motion of the rat myocardium estimated using quantitative ultrasound with new strategies for infarct detection

Introduction: The estimation of myocardial motion abnormalities has great potential for the early diagnosis of myocardial infarction (MI). This study aims to quantitatively analyze the segmental and transmural myocardial motion in MI rats by incorporating two novel strategies of algorithm parameter optimization and transmural motion index (TMI) calculation. Methods: Twenty-one rats were randomly divided into three groups (n = 7 per group): sham, MI, and ischemia-reperfusion (IR) groups. Ultrasound radio-frequency (RF) signals were acquired from each rat heart at 1 day and 28 days after animal model establishment; thus, a total of six datasets were represented as Sham1, Sham28, MI1, MI28, IR1, and IR28. The systolic cumulative displacement was calculated using our previously proposed vectorized normalized cross-correlation (VNCC) method. A semiautomatic regional and layer-specific myocardium segmentation framework was proposed for transmural and segmental myocardial motion estimation. Two novel strategies were proposed: the displacement-compensated cross-correlation coefficient (DCCCC) for algorithm parameter optimization and the transmural motion index (TMI) for quantitative estimation of the cross-wall transmural motion gradient. Results: The results showed that an overlap value of 80% used in VNCC guaranteed a more accurate displacement calculation. Compared to the Sham1 group, the systolic myocardial motion reductions were significantly detected (p < 0.05) in the middle anteroseptal (M-ANT-SEP), basal anteroseptal (B-ANT-SEP), apical lateral (A-LAT), middle inferolateral (M-INF-LAT), and basal inferolateral (B-INF-LAT) walls as well as a significant TMI drop (p < 0.05) in the M-ANT-SEP wall in the MI1 rats; significant motion reductions (p < 0.05) were also detected in the B-ANT-SEP and A-LAT walls in the IR1 group. The motion improvements (p < 0.05) were detected in the M-INF-LAT wall in the MI28 group and the apical septal (A-SEP) wall in the IR28 group compared to the MI1 and IR1 groups, respectively. Discussion: Our results show that the MI-induced reductions and reperfusion-induced recovery in systolic myocardial contractility could be successfully evaluated using our method, and most post-MI myocardial segments could recover systolic function to various extents in the remodeling phase. In conclusion, the ultrasound-based quantitative estimation framework for estimating segmental and transmural motion of the myocardium proposed in our study has great potential for non-invasive, novel, and early MI detection.

Preliminary Feasibility of Stress Myocardial Elastography for the Detection of Coronary Artery Disease

Myocardial elastography (ME) is a cardiac strain imaging technique that has been found capable of detecting a decrease in radial strain caused by ischemia or infarction in patients with coronary artery disease (CAD) as well as in a canine model. Prior studies have focused on rest imaging, but stress testing can reveal functional deficits caused by stenoses that are asymptomatic at rest. Therefore, it has been proposed that stress ME (S-ME) improves the detection of CAD. A novel strain difference (Δε) metric is presented and investigated in a canine model of induced ischemia, as well as in a study in human patients with CAD validated by myocardial perfusion imaging. In the canine model study, flow-limiting stenosis was induced by partial ligation in n = 2 canines, and stenosis was found to consistently reduce Δε in the affected myocardial regions compared with baseline, as well as compared to myocardial regions that are remote to the induced stenosis. In the clinical study, the median Δε was significantly lower (p < 0.05) in infarcted myocardial regions (-6.29%) than in those with normal perfusion (4.62%), with Δε in ischemic regions falling in between (-2.91%). The same trend was observed when considering radial strain during stress and, to a lesser degree, at rest alone. The results indicate that S-ME may be more sensitive to mild cases of CAD that are functionally asymptomatic at rest.

Research Progress of Myocardial Fibrosis and Atrial Fibrillation

With the aging population and the increasing incidence of basic illnesses such as hypertension and diabetes (DM), the incidence of atrial fibrillation (AF) has increased significantly. AF is the most common arrhythmia in clinical practice, which can cause heart failure (HF) and ischemic stroke (IS), increasing disability and mortality. Current studies point out that myocardial fibrosis (MF) is one of the most critical substrates for the occurrence and maintenance of AF. Although myocardial biopsy is the gold standard for evaluating MF, it is rarely used in clinical practice because it is an invasive procedure. In addition, serological indicators and imaging methods have also been used to evaluate MF. Nevertheless, the accuracy of serological markers in evaluating MF is controversial. This review focuses on the pathogenesis of MF, serological evaluation, imaging evaluation, and anti-fibrosis treatment to discuss the existing problems and provide new ideas for MF and AF evaluation and treatment.

Myocardial Strain Imaging With Electrocardiogram-Gated and Coherent Compounding for Early Diagnosis of Coronary Artery Disease

Myocardial elastography (ME) is an ultrasound-based technique that uses radiofrequency signals for 2-D cardiac motion tracking and strain imaging at a high frame rate. Early diagnosis of coronary artery disease (CAD) is critical for timely treatment and improvement of patient outcome. The objective of this study was to assess the performance of ME radial and circumferential strains in the detection and characterization of CAD in patients. In this study, 86 patients suspected of CAD were imaged with ME prior to invasive coronary angiography (ICA). End-systolic radial and circumferential left ventricular strains were estimated in all patients in each of their perfusion territories: left anterior descending (LAD), left circumflex (LCX) and right coronary artery (RCA). ME radial strains were capable of differentiating the obstructive CAD group (55.3 ± 29.8%) from the non-obstructive CAD (72.5 ± 46.8%, p < 0.05) and no CAD groups (73.4 ± 30.4%, p < 0.05) in the RCA territory. ME circumferential strains were capable of differentiating the obstructive CAD group (-3.1 ± 7.5%) from the non-obstructive CAD (-7.2 ± 6.8%, p < 0.05) and normal (-6.9 ± 8.0%, p < 0.05) groups in the LAD territory and to differentiate the normal group (-17.1 ± 8.2%) from the obstructive (-12.8 ± 7.2%, p < 0.05) and non-obstructive CAD (-13.6 ± 8.5%, p < 0.05) groups in the RCA territory. ME circumferential strain performed better than ME radial strain in differentiating normal, non-obstructive and obstructive perfusion territories. In the LCX territory, both ME radial and circumferential strains decreased when the level of stenosis was higher. However, it was not statistically significant. The findings presented herein indicate that ME radial and circumferential estimation obtained from ECG-gated and compounded acquisitions is a promising tool for early, non-invasive and radiation-free detection of CAD in patients.

Dynamic Myocardial Ultrasound Localization Angiography

Dynamic Myocardial Ultrasound Localization Angiography (MULA) is an ultrasound-based imaging modality destined to enhance the diagnosis and treatment monitoring of coronary pathologies. Current diagnosis methods of coronary artery disease focus on the observation of vessel narrowing in the coronary vasculature to assess the organ's condition. However, we would strongly benefit from mapping and measuring flow from intramyocardial arterioles and capillaries as they are the direct vehicle of the myocardium blood income. With the advent of ultrafast ultrasound scanners, imaging modalities based on the localization and tracking of injected microbubbles allow for the subwavelength resolution imaging of an organ's vasculature. Yet, the application of these vascular imaging modalities relies on an accumulation of cine loops of a region of interest undergoing no or minimal tissue motion. This work introduces the MULA framework that combines 1) the mapping of the dynamics of the microvascular flow using an ultrasound sequence triggered by the electrocardiogram with a 2) novel Lagrangian beamformer based on non-rigid motion registration algorithm to form images directly in the myocardium's material coordinates and thus correcting for the large myocardial motion and deformation. Specifically, we show that this framework enables the non-invasive imaging of the angioarchitecture and dynamics of intramyocardial flow in vessels as small as a few tens of microns in the rat's beating heart in vivo.

Evaluation of strain averaging area and strain estimation errors in a spheroidal left ventricular model using synthetic image data and speckle tracking

BACKGROUND

In majority of studies on speckle tracking echocardiography (STE) the strain estimates are averaged over large areas of the left ventricle. This may impair the diagnostic capability of the STE in the case of e.g. local changes of the cardiac contractility. This work attempts to evaluate, how far one can reduce the averaging area, without sacrificing the estimation accuracy that could be important from the clinical point of view.

METHODS

Synthetic radio frequency (RF) data of a spheroidal left ventricular (LV) model were generated using FIELD II package and meshes obtained from finite element method (FEM) simulation. The apical two chamber (A2C) view and the mid parasternal short axis view (pSAXM) were simulated. The sector encompassed the entire cross-section (full view) of the LV model or its part (partial view). The wall segments obtained according to the American Heart Association (AHA17) were divided into subsegments of area decreasing down to 3 mm2. Longitudinal, circumferential and radial strain estimates, obtained using a hierarchical block-matching method, were averaged over these subsegments. Estimation accuracy was assessed using several error measures, making most use of the prediction of the maximal relative error of the strain estimate obtained using the FEM derived reference. Three limits of this predicted maximal error were studied, namely 16.7%, 33% and 66%. The smallest averaging area resulting in the strain estimation error below one of these limits was considered the smallest allowable averaging area (SAAA) of the strain estimation.

RESULTS

In all AHA17 segments, using the A2C projection, the SAAA ensuring maximal longitudinal strain estimates error below 33% was below 3 mm2, except for the segment no 17 where it was above 278 mm2. The SAAA ensuring maximal circumferential strain estimates error below 33% depended on the AHA17 segment position within the imaging sector and view type and ranged from below 3-287 mm2. The SAAA ensuring maximal radial strain estimates error below 33% obtained in the pSAXM projection was not less than 287 mm2. The SAAA values obtained using other maximal error limits differ from SAAA values observed for the 33% error limit only in limited number of cases. SAAA decreased when using maximal error limit equal to 66% in these cases. The use of the partial view (narrow sector) resulted in a decrease of the SAAA.

CONCLUSIONS

The SAAA varies strongly between strain components. In a vast part of the LV model wall in the A2C view the longitudinal strain could be estimated using SAAA below 3 mm2, which is smaller than the averaging area currently used in clinic, thus with a higher resolution. The SAAA of the circumferential strain estimation strongly depends on the position of the region of interest and the parameters of the acquisition. The SAAA of the radial strain estimation takes the highest values. The use of a narrow sector could increase diagnostic capabilities of 2D STE.

Speckle-Tracking Echocardiography with Novel Imaging Technique of Higher Frame Rate

Background: The accuracy of speckle-tracking echocardiography (STE) depends on temporal resolution. The goal of this study was to demonstrate the feasibility of relatively high frame rate (rHi-FR) (~200 fps) for STE. Methods: In this prospective study, echocardiographic images were acquired using clinical scanners on patients with normal left ventricular systolic function using rHi-FR and conventional frame rate (Reg-FR) (~50 FPS). GLS values were evaluated on apical 4-, 2- and 3-chamber images acquired in both rHi-FR and Reg-FR. Inter-observer and intra-observer variabilities were assessed in rHi-FR and Reg-FR. Results: There were 143 echocardiograms evaluated in this study. The frame rate of rHi-FR was 190 ± 25 and Reg-FR was 50 ± 3, and the heart rate was 71 ± 13. Absolute strain values measured in rHi-FR were significantly higher than those measured in Reg-FR (all p < 0.001). Inter-observer and intra-observer correlations were strong in both rHi-FR and Reg-FR. Conclusions: We demonstrated that absolute strain values were significantly higher using rHi-FR when compared with Reg-FR. It is plausible that higher temporal resolution enabled the measurement of myocardial strain at desired time point. Further investigations are necessary to evaluate the value of rHi-FR to assess myocardial strain in the setting of tachycardia.

Radial Motion Estimation of Myocardium in Rats with Myocardial Infarction: A Hybrid Method of FNCCGLAM and Polar Transformation

Ultrasound elastography is a novel approach of evaluating regional myocardial systolic function and detecting infarcted area. This study aims to evaluate the radial motion of myocardial infarction (MI) area in left ventricular parasternal short axis (PSAX) view using a hybrid method of fast normalized cross-correlation and global analytic minimization (FNCCGLAM) and polar transformation. Fifteen rats were randomly selected for sham group, MI group and ischemia-reperfusion (IR) group (N = 5 for each group). The ultrasound radiofrequency data of the PSAX view of rat heart were acquired. After polar transformation of the data, the infarcted myocardium with the change of mechanical property was tracked over one myocardial systolic phase by the proposed method in comparison with fast normalized cross-correlation (FNCC) and dynamic programming analytic minimization (DPAM). To obtain a clear visualization of the myocardium, the inverse polar transformation was performed. The results indicated that the use of FNCCGLAM refined the myocardial displacements to obtain high-quality myocardial elastographic map with a higher contrast-to-noise ratio and dynamically tracked the infarcted myocardial segment with a higher success rate in comparison with FNCC and DPAM. It was found that the radial systolic motion of the infarcted anterior segment in the MI group reduced significantly (p < 0.05) in comparison with the sham group, while the systolic function of that myocardial segment in the IR group recovered at some extent. The results in this study suggest that FNCCGLAM is superior to FNCC and DPAM with the improved accuracy and robustness of motion estimation and has potentials as displacement estimator in ultrasound elastography.

Monitoring Canine Myocardial Infarction Formation and Recovery via Transthoracic Cardiac Strain Imaging

Myocardial elastography (ME) is an ultrasound-based strain imaging method that aims to determine the degree of ischemia or infarction as a result of the change in the elastic properties of the myocardium. A survival canine model (n = 11) was employed to investigate the ability of ME to image myocardial infarction formation and recovery. Infarcts were generated by ligation of the left anterior descending coronary artery. Canines were survived and imaged for 4 days (n = 7) or 4 weeks (n = 4), allowing sufficient time for recovery via collateral perfusion. A radial strain-based metric, percentage of healthy myocardium by strain (PHMε), was developed as a marker for healthy myocardial tissue. PHMε was strongly linearly correlated with actual infarct size as determined by gross pathology (R2 = 0.80). Mean PHMε was reduced 1-3 days post-infarction (p < 0.05) at the papillary and apical short-axis levels; full infarct recovery was achieved by day 28, with mean PHMε returning to baseline levels. ME was capable of diagnosing individual myocardial segments as non-infarcted or infarcted with high sensitivity (82%), specificity (92%) and precision (85%) (area under the receiver operating characteristic curve = 0.90). The study therefore strengthens the ME premise that it can detect and assess myocardial infarction progression and recovery in vivo and could thus provide an important role in both disease diagnosis and treatment assesssment.

A comparison between unfocused and focused transmit strategies in cardiac strain imaging

Unfocused ultrasound imaging, particularly coherent compounding with diverging waves, is a commonly employed high-frame rate transmit strategy in cardiac strain imaging. However, the accuracy and precision of diverging wave imaging compared to focused-beam transmit approaches in human subjects is unknown. Three transmit strategies-coherent compounding imaging, composite focused imaging with ECG gating and narrow-beams, and focused imaging with wide-beams-were compared in simulation and in transthoracic imaging of healthy human subjects (n  =  7). The focused narrow-beam sequence estimated radial end-systolic cumulative strains of a simulated left ventricular deformation with 26%  ±  1.5% and 34%  ±  1.5% greater accuracy compared with compounding and wide-beam imaging, respectively. Strain estimation precision in transthoracic imaging was then assessed with the Strain Filter on cumulative end-systolic radial strains. Within the strain values where statistically significant differences in precision (E(SNR_e|ε)) were found between transmit strategies, the narrow-beam sequence estimated radial strain 13%  ±  0.71% and 34%  ±  8.9% more precisely on average compared to compounding or wide-beam imaging, respectively.

Optimization of Transmit Parameters in Cardiac Strain Imaging With Full and Partial Aperture Coherent Compounding

Coherent compounding methods using the full or partial transmit aperture have been investigated as a possible means of increasing strain measurement accuracy in cardiac strain imaging; however, the optimal transmit parameters in either compounding approach have yet to be determined. The relationship between strain estimation accuracy and transmit parameters-specifically the subaperture, angular aperture, tilt angle, number of virtual sources, and frame rate-in partial aperture (subaperture compounding) and full aperture (steered compounding) fundamental mode cardiac imaging was thus investigated and compared. Field II simulation of a 3-D cylindrical annulus undergoing deformation and twist was developed to evaluate accuracy of 2-D strain estimation in cross-sectional views. The tradeoff between frame rate and number of virtual sources was then investigated via transthoracic imaging in the parasternal short-axis view of five healthy human subjects, using the strain filter to quantify estimation precision. Finally, the optimized subaperture compounding sequence (25-element subperture, 90° angular aperture, 10 virtual sources, 300-Hz frame rate) was compared to the optimized steered compounding sequence (60° angular aperture, 15° tilt, 10 virtual sources, 300-Hz frame rate) via transthoracic imaging of five healthy subjects. Both approaches were determined to estimate cumulative radial strain with statistically equivalent precision (subaperture compounding E(SNRe %) = 3.56, and steered compounding E(SNRe %) = 4.26).

Simultaneous positron emission tomography and ultrafast ultrasound for hybrid molecular, anatomical and functional imaging

Positron emission tomography-computed tomography (PET-CT) is the most sensitive molecular imaging modality, but it does not easily allow for rapid temporal acquisition. Ultrafast ultrasound imaging (UUI)-a recently introduced technology based on ultrasonic holography-leverages frame rates of up to several thousand images per second to quantitatively map, at high resolution, haemodynamic, biomechanical, electrophysiological and structural parameters. Here, we describe a pre-clinical scanner that registers PET-CT and UUI volumes acquired simultaneously and offers multiple combinations for imaging. We demonstrate that PET-CT-UUI allows for simultaneous images of the vasculature and metabolism during tumour growth in mice and rats, as well as for synchronized multi-modal cardiac cine-loops. Combined anatomical, functional and molecular imaging with PET-CT-UUI represents a high-performance and clinically translatable technology for biomedical research.

Biomechanical assessment of myocardial infarction using optical coherence elastography

Myocardial infarction (MI) leads to cardiomyocyte loss, impaired cardiac function, and heart failure. Molecular genetic analyses of myocardium in mouse models of ischemic heart disease have provided great insight into the mechanisms of heart regeneration, which is promising for novel therapies after MI. Although biomechanical factors are considered an important aspect in cardiomyocyte proliferation, there are limited methods for mechanical assessment of the heart in the mouse MI model. This prevents further understanding the role of tissue biomechanics in cardiac regeneration. Here we report optical coherence elastography (OCE) of the mouse heart after MI. Surgical ligation of the left anterior descending coronary artery was performed to induce an infarction in the heart. Two OCE methods with assessment of the direction-dependent elastic wave propagation and the spatially resolved displacement damping provide complementary analyses of the left ventricle. In comparison with sham, the infarcted heart features a fibrotic scar region with reduced elastic wave velocity, decreased natural frequency, and less mechanical anisotropy at the tissue level at the sixth week post-MI, suggesting lower and more isotropic stiffness. Our results indicate that OCE can be utilized for nondestructive biomechanical characterization of MI in the mouse model, which could serve as a useful tool in the study of heart repair.

Cardiac Lesion Mapping In Vivo Using Intracardiac Myocardial Elastography

Radio frequency (RF) ablation of the myocardium is used to treat various cardiac arrhythmias. The size, spacing, and transmurality of lesions have been shown to affect the success of the ablation procedure; however, there is currently no method to directly image the size and formation of ablation lesions in real time. Intracardiac myocardial elastography (ME) has been previously used to image the decrease in cardiac strain during systole in the ablated region as a result of the lesion formation. However, the feasibility of imaging multiple lesions and identifying the presence of gaps between lesions has not yet been investigated. In this paper, RF ablation lesions ( ) were generated in the left ventricular epicardium in three anesthetized canines. Two sets of two lesions each were created in close proximity to one another with small gaps (1.5 and 4 cm), while one set of two lesions was created directly next to each other with no gap. A clinical intracardiac echocardiography system was programmed to transmit a custom diverging beam sequence at 600 Hz and used to image the ablation site before and after the induction of ablation lesions. Cumulative strains were estimated over systole using a normalized cross-correlational displacement algorithm and a least-squares strain kernel. Afterward, lesions were excised and subjected to tetrazolium chloride staining. Results indicate that intracardiac ME was capable of imaging the reduction in systolic strain associated with the formation of an ablation lesion. Furthermore, lesion sets containing gaps were able to be distinguished from lesion sets created with no gaps. These results indicate that the end-systolic strain measured using intracardiac ME may be used to image the formation of lesions induced during an RF ablation procedure, in order to provide critical assessment of lesion viability during the interventional procedure.

Electromechanical wave imaging and electromechanical wave velocity estimation in a large animal model of myocardial infarction

Echocardiography is often used in the clinic for detection and characterization of myocardial infarction. Electromechanical wave imaging (EWI) is a non-invasive ultrasound-based imaging technique based on time-domain incremental motion and strain estimation that can evaluate changes in contractility in the heart. In this study, electromechanical activation is assessed in infarcted heart to determine whether EWI is capable of detecting and monitoring infarct formation. Additionally, methods for estimating electromechanical wave (EW) velocity are presented, and changes in the EW propagation velocity after infarct formation are studied. Five (n  =  5) adult mongrels were used in this study. Successful infarct formation was achieved in three animals by ligation of the left anterior descending (LAD) coronary artery. Dogs were survived for a few days after LAD ligation and monitored daily with EWI. At the end of the survival period, dogs were sacrificed and TTC (tetrazolium chloride) staining confirmed the formation and location of the infarct. In all three dogs, as soon as day 1 EWI was capable of detecting late-activated and non-activated regions, which grew over the next few days. On final day images, the extent of these regions corresponded to the location of infarct as confirmed by staining. EW velocities in border zones of infarct were significantly lower post-infarct formation when compared to baseline, whereas velocities in healthy tissues were not. These results indicate that EWI and EW velocity might help with the detection of infarcts and their border zones, which may be useful for characterizing arrhythmogenic substrate.

Ultrasound functional imaging in an ex vivo beating porcine heart platform

In recent years, novel ultrasound functional imaging (UFI) techniques have been introduced to assess cardiac function by measuring, e.g. cardiac output (CO) and/or myocardial strain. Verification and reproducibility assessment in a realistic setting remain major issues. Simulations and phantoms are often unrealistic, whereas in vivo measurements often lack crucial hemodynamic parameters or ground truth data, or suffer from the large physiological and clinical variation between patients when attempting clinical validation. Controlled validation in certain pathologies is cumbersome and often requires the use of lab animals. In this study, an isolated beating pig heart setup was adapted and used for performance assessment of UFI techniques such as volume assessment and ultrasound strain imaging. The potential of performing verification and reproducibility studies was demonstrated. For proof-of-principle, validation of UFI in pathological hearts was examined. Ex vivo porcine hearts (n  =  6, slaughterhouse waste) were resuscitated and attached to a mock circulatory system. Radio frequency ultrasound data of the left ventricle were acquired in five short axis views and one long axis view. Based on these slices, the CO was measured, where verification was performed using flow sensor measurements in the aorta. Strain imaging was performed providing radial, circumferential and longitudinal strain to assess reproducibility and inter-subject variability under steady conditions. Finally, strains in healthy hearts were compared to a heart with an implanted left ventricular assist device, simulating a failing, supported heart. Good agreement between ultrasound and flow sensor based CO measurements was found. Strains were highly reproducible (intraclass correlation coefficients  >0.8). Differences were found due to biological variation and condition of the hearts. Strain magnitude and patterns in the assisted heart were available for different pump action, revealing large changes compared to the normal condition. The setup provides a valuable benchmarking platform for UFI techniques. Future studies will include work on different pathologies and other means of measurement verification.

A Comparison of Finite Element-Based Inversion Algorithms, Local Frequency Estimation, and Direct Inversion Approach Used in MRE

In quantitative elastography, maps of the mechanical properties of soft tissue, or elastograms, are calculated from the measured displacement data by solving an inverse problem. The model assumptions have a significant effect on elastograms. Motivated by the high sensitivity of imaging results to the model assumptions for in vivo magnetic resonance elastography of the prostate, we compared elastograms obtained with four different methods. Two finite-element method (FEM)-based methods developed by our group were compared with two other commonly used methods, local frequency estimator (LFE) and curl-based direct inversion (c-DI). All the methods assume a linear isotropic elastic model, but the methods vary in their assumptions, such as local homogeneity or incompressibility, and in the specific approach used. We report results using simulations, phantom, and ex vivo and in vivo data. The simulation and phantom studies show, for regions with an inclusion, that the contrast to noise ratio (CNR) for the FEM methods is about three to five times higher than the CNR for the LFE and c-DI and the rms error is about half. The LFE method produces very smooth results (i.e., low CNR) and is fast. c-DI is faster than the FEM methods but it is only accurate in areas where elasticity variations are small. The artifacts resulting from the homogeneity assumption in c-DI is detrimental in regions with large variations. The ex vivo and in vivo results also show similar trends as the simulation and phantom studies. The c-FEM method is more sensitive to noise compared with the mixed-FEM due to higher orders derivatives. This is especially evident at lower frequencies, where the wave curvature is smaller and it is more prone to such error, causing a discrepancy in the absolute values between the mixed-FEM and c-FEM in our in vivo results. In general, the proposed FEMs use fewer simplifying assumptions and outperform the other methods but they are computationally more expensive.

A Systematic Investigation of Lateral Estimation Using Various Interpolation Approaches in Conventional Ultrasound Imaging

Accurate lateral displacement and strain estimation is critical for some applications of elasticity imaging. Typically, motion estimation in the lateral direction is challenging because of low sampling frequency and lack of phase information in conventional ultrasound imaging. Several approaches have been proposed to improve the performance of lateral estimation, such as lateral interpolation on the radio frequency (RF) signals (Interp_RF), lateral interpolation on the cross-correlation function (Interp_CCF), and lateral interpolation on both the RF signals and cross-correlation function (Interp_Both). In this paper, the estimation performances of the above-mentioned three approaches are compared systematically in simulations and phantom experiments. In the simulations, the root-mean-square error (RMSE) of axial/lateral displacement and strain is utilized to assess the accuracy of motion estimation. In the phantom experiments, the displacement quality metric (DQM), defined as the normalized cross-correlation between the motion-compensated reference frame and the comparison frame, and the contrast-to-noise ratio (CNR) of axial/lateral strain are used as the evaluation criteria. The results show that the three approaches have similar performance in axial estimation. For lateral estimation, if the line density of ultrasound imaging is relatively high (i.e., >4.2 lines/mm), Interp_CCF is comparable to Interp_Both, and Interp_RF performs the worst. However, if the line density is relatively low (i.e., <2.8 lines/mm), Interp_Both performs the best as indicated by the lowest RMSEs or highest DQMs and CNRs in lateral estimation. The trend is consistent at different window sizes, applied strains, and sonographic signal-to-noise ratios (>20 dB). Besides, Interp_Both with a small interpolation factor (e.g., 3-5) is found to obtain the best tradeoff between the estimation accuracy and the computational cost, and thus is suggested for lateral motion estimation in the case of a low line density (i.e., <2.8 lines/mm).

Cardiac Strain Imaging With Coherent Compounding of Diverging Waves

Current methods of cardiac strain imaging at high frame rate suffer from motion matching artifacts or poor lateral resolution. Coherent compounding has been shown to improve echocardiographic image quality while maintaining a high frame rate, but has never been used to image cardiac strain. However, myocardial velocity can have an impact on coherent compounding due to displacements between frames. The objective of this paper was to investigate the feasibility and performance of coherent compounding for cardiac strain imaging at a low and a high myocardial velocity. Left-ventricular contraction in short-axis view was modeled as an annulus with radial thickening and circumferential rotation. Simulated radio-frequency channel data with a cardiac phased array were obtained using three different beamforming methods: single diverging wave, coherent compounding of diverging waves, and conventional focusing. Axial and lateral displacements and strains as well as radial strains were estimated and compared to their true value. In vivo feasibility of cardiac strain imaging with coherent compounding was performed and compared to single diverging wave imaging. At low myocardial velocities, the axial, lateral, and radial strain relative error for nine compounded waves (16.3%, 40.4%, and 18.9%) were significantly lower than those obtained with single diverging wave imaging (19.9%, 80.3%, and 30.6%) and closer to that obtained with conventional focusing (16.7%, 43.7%, and 16%). In vivo left-ventricular radial strains exhibited higher quality with nine compounded waves than with single diverging wave imaging. These results indicate that cardiac strain can be imaged using coherent compounding of diverging waves with a better performance than with single diverging wave imaging while maintaining a high frame rate, and therefore, has the potential to improve diagnosis of myocardial strain-based cardiac diseases.

Evaluation of Coronary Artery Disease Using Myocardial Elastography with Diverging Wave Imaging: Validation against Myocardial Perfusion Imaging and Coronary Angiography

Myocardial elastography (ME) is an ultrasound-based technique that can image 2-D myocardial strains. The objectives of this study were to illustrate that 2-D myocardial strains can be imaged with diverging wave imaging and differ, on average, between normal and coronary artery disease (CAD) patients. In this study, 66 patients with symptoms of CAD were imaged with myocardial elastography before a nuclear stress test or an invasive coronary angiography. Radial cumulative strains were estimated in all patients. The end-systolic radial strain in the total cross section of the myocardium was significantly higher in normal patients (17.9 ± 8.7%) than in patients with reversible perfusion defect (6.2 ± 9.3%, p < 0.001) and patients with significant (-0.9 ± 7.4%, p < 0.001) and non-significant (3.7 ± 5.7%, p < 0.01) lesions. End-systolic radial strain in the left anterior descending, left circumflex and right coronary artery territory was found to be significantly higher in normal patients than in CAD patients. These preliminary findings indicate that end-systolic radial strain measured with ME is higher on average in healthy persons than in CAD patients and that ME has the potential to be used for non-invasive, radiation-free early detection of CAD.

3D Myocardial Elastography In Vivo

Strain evaluation is of major interest in clinical cardiology as it can quantify the cardiac function. Myocardial elastography, a radio-frequency (RF)-based cross-correlation method, has been developed to evaluate the local strain distribution in the heart in vivo. However, inhomogeneities such as RF ablation lesions or infarction require a three-dimensional approach to be measured accurately. In addition, acquisitions at high volume rate are essential to evaluate the cardiac strain in three dimensions. Conventional focused transmit schemes using 2D matrix arrays, trade off sufficient volume rate for beam density or sector size to image rapid moving structure such as the heart, which lowers accuracy and precision in the strain estimation. In this study, we developed 3D myocardial elastography at high volume rates using diverging wave transmits to evaluate the local axial strain distribution in three dimensions in three open-chest canines before and after radio-frequency ablation. Acquisitions were performed with a 2.5 MHz 2D matrix array fully programmable used to emit 2000 diverging waves at 2000 volumes/s. Incremental displacements and strains enabled the visualization of rapid events during the QRS complex along with the different phases of the cardiac cycle in entire volumes. Cumulative displacement and strain volumes depict high contrast between non-ablated and ablated myocardium at the lesion location, mapping the tissue coagulation. 3D myocardial strain elastography could thus become an important technique to measure the regional strain distribution in three dimensions in humans.

Imaging the Propagation of the Electromechanical Wave in Heart Failure Patients with Cardiac Resynchronization Therapy

BACKGROUND

Current electrocardiographic and echocardiographic measurements in heart failure (HF) do not take into account the complex interplay between electrical activation and local wall motion. The utilization of novel technologies to better characterize cardiac electromechanical behavior may lead to improved response rates with cardiac resynchronization therapy (CRT). Electromechanical wave imaging (EWI) is a noninvasive ultrasound-based technique that uses the transient deformations of the myocardium to track the intrinsic EW that precedes myocardial contraction. In this paper, we investigate the performance and reproducibility of EWI in the assessment of HF patients and CRT.

METHODS

EWI acquisitions were obtained in five healthy controls and 16 HF patients with and without CRT pacing. Responders (n = 8) and nonresponders (n = 8) to CRT were identified retrospectively on the basis of left ventricular (LV) reverse remodeling. Electromechanical activation maps were obtained in all patients and used to compute a quantitative parameter describing the mean LV lateral wall activation time (LWAT).

RESULTS

Mean LWAT was increased by 52.1 ms in HF patients in native rhythm compared to controls (P < 0.01). For all HF patients, CRT pacing initiated a different electromechanical activation sequence. Responders exhibited a 56.4-ms ± 28.9-ms reduction in LWAT with CRT pacing (P < 0.01), while nonresponders showed no significant change.

CONCLUSION

In this initial feasibility study, EWI was capable of characterizing local cardiac electromechanical behavior as it pertains to HF and CRT response. Activation sequences obtained with EWI allow for quantification of LV lateral wall electromechanical activation, thus providing a novel method for CRT assessment.

Electromechanical wave imaging (EWI) validation in all four cardiac chambers with 3D electroanatomic mapping in canines in vivo

Characterization and mapping of arrhythmias is currently performed through invasive insertion and manipulation of cardiac catheters. Electromechanical wave imaging (EWI) is a non-invasive ultrasound-based imaging technique, which tracks the electromechanical activation that immediately follows electrical activation. Electrical and electromechanical activations were previously found to be linearly correlated in the left ventricle, but the relationship has not yet been investigated in the three other chambers of the heart. The objective of this study was to investigate the relationship between electrical and electromechanical activations and validate EWI in all four chambers of the heart with conventional 3D electroanatomical mapping. Six (n  =  6) normal adult canines were used in this study. The electrical activation sequence was mapped in all four chambers of the heart, both endocardially and epicardially using the St Jude's EnSite 3D mapping system (St. Jude Medical, Secaucus, NJ). EWI acquisitions were performed in all four chambers during normal sinus rhythm, and during pacing in the left ventricle. Isochrones of the electromechanical activation were generated from standard echocardiographic imaging views. Electrical and electromechanical activation maps were co-registered and compared, and electrical and electromechanical activation times were plotted against each other and linear regression was performed for each pair of activation maps. Electromechanical and electrical activations were found to be directly correlated with slopes of the correlation ranging from 0.77 to 1.83, electromechanical delays between 9 and 58 ms and R ² values from 0.71 to 0.92. The linear correlation between electrical and electromechanical activations and the agreement between the activation maps indicate that the electromechanical activation follows the pattern of propagation of the electrical activation. This suggests that EWI may be used as a novel non-invasive method to accurately characterize and localize sources of arrhythmias.

Dynamic estimation of myocardial deformation using ultrasound RF-data: A preliminary study

Analytic minimization (AM) is one of elastography algorithms. In this study, a healthy male volunteer lying in lateral position was examined by an ultrasound system. The left ventricular long-axis view of the heart was acquired and the ultrasound radio frequency (RF) data were stored for estimation of myocardial deformation. The methods of AM and Linear polynomial curve fitting were used to calculate the displacement and strain, respectively. The myocardial deformation was dynamically estimated. The strain maps at approximately 50ms, 140ms, and 240ms after the beginning of the ventricle systole indicate the contraction of different myocardial segments. This study provides a potential method to assess the survival myocardium in patients with myocardial infarction in future study.

Intracardiac myocardial elastography in canines and humans in vivo

Intracardiac echocardiography (ICE) is a useful imaging modality which is used during RF ablation procedures to identify anatomical structures. Utilizing ICE in conjunction with myocardial elastography (ME) can provide additional information on the mechanical properties of cardiac tissue and provide information on mechanical changes caused by ablation. The objective of this study was to demonstrate that ICE can be used at high frame rate using a diverging beam transmit sequence to image myocardial strain and differentiate myocardial tissue properties before, during, and after ablation for a clinical ablation procedure. In this feasibility study, three normal canines and eight patients with atrial fibrillation (AF) were studied in vivo. A 5.8-MHz ICE transducer was used to image the heart with a diverging beam transmit method achieving 1200 frames per second (fps). Cumulative axial displacement estimation was performed using 1-D cross-correlation with a window size of 2.7 mm and 95% overlap. Axial cumulative strains were estimated in the left atrium (LA) and right atrium (RA) using a least-squares estimator with a kernel of 2 mm on the axial displacements. In the canine case, radial thickening was detected in the lateral wall and in the interatrial septum during LA emptying. For AF patients, the mean absolute strain in the ablated region was lower (6.7 ± 3.1%) than before the ablation (17.4 ± 9.3%) in LA at the end of the LA emptying phase. In the cavotricuspid isthmus (CTI) region, mean absolute strain magnitude at the end of the RA emptying phase was found to be higher during ablation (43.0 ± 18.1%) compared with after ablation (33.7 ± 15.8%). Myocardial strains in the LA of an AF patient were approximately 2.6 times lower in the ablated region than before ablation. This initial feasibility indicates that ME can be used as a new imaging modality in conjunction with ICE in RF ablation guidance and lesion monitoring.

3D ultrafast ultrasound imaging in vivo

Very high frame rate ultrasound imaging has recently allowed for the extension of the applications of echography to new fields of study such as the functional imaging of the brain, cardiac electrophysiology, and the quantitative imaging of the intrinsic mechanical properties of tumors, to name a few, non-invasively and in real time. In this study, we present the first implementation of Ultrafast Ultrasound Imaging in 3D based on the use of either diverging or plane waves emanating from a sparse virtual array located behind the probe. It achieves high contrast and resolution while maintaining imaging rates of thousands of volumes per second. A customized portable ultrasound system was developed to sample 1024 independent channels and to drive a 32  ×  32 matrix-array probe. Its ability to track in 3D transient phenomena occurring in the millisecond range within a single ultrafast acquisition was demonstrated for 3D Shear-Wave Imaging, 3D Ultrafast Doppler Imaging, and, finally, 3D Ultrafast combined Tissue and Flow Doppler Imaging. The propagation of shear waves was tracked in a phantom and used to characterize its stiffness. 3D Ultrafast Doppler was used to obtain 3D maps of Pulsed Doppler, Color Doppler, and Power Doppler quantities in a single acquisition and revealed, at thousands of volumes per second, the complex 3D flow patterns occurring in the ventricles of the human heart during an entire cardiac cycle, as well as the 3D in vivo interaction of blood flow and wall motion during the pulse wave in the carotid at the bifurcation. This study demonstrates the potential of 3D Ultrafast Ultrasound Imaging for the 3D mapping of stiffness, tissue motion, and flow in humans in vivo and promises new clinical applications of ultrasound with reduced intra--and inter-observer variability.

New developments in paediatric cardiac functional ultrasound imaging

Ultrasound imaging can be used to estimate the morphology as well as the motion and deformation of tissues. If the interrogated tissue is actively deforming, this deformation is directly related to its function and quantification of this deformation is normally referred as 'strain imaging'. Tissue can also be deformed by applying an internal or external force and the resulting, induced deformation is a function of the mechanical tissue characteristics. In combination with the load applied, these strain maps can be used to estimate or reconstruct the mechanical properties of tissue. This technique was named 'elastography' by Ophir et al. in 1991. Elastography can be used for atherosclerotic plaque characterisation, while the contractility of the heart or skeletal muscles can be assessed with strain imaging. Rather than using the conventional video format (DICOM) image information, radio frequency (RF)-based ultrasound methods enable estimation of the deformation at higher resolution and with higher precision than commercial methods using Doppler (tissue Doppler imaging) or video image data (2D speckle tracking methods). However, the improvement in accuracy is mainly achieved when measuring strain along the ultrasound beam direction, so it has to be considered a 1D technique. Recently, this method has been extended to multiple directions and precision further improved by using spatial compounding of data acquired at multiple beam steered angles. Using similar techniques, the blood velocity and flow can be determined. RF-based techniques are also beneficial for automated segmentation of the ventricular cavities. In this paper, new developments in different techniques of quantifying cardiac function by strain imaging, automated segmentation, and methods of performing blood flow imaging are reviewed and their application in paediatric cardiology is discussed.

Radial basis functions for combining shape and speckle tracking in 4D echocardiography

Quantitative analysis of left ventricular deformation can provide valuable information about the extent of disease as well as the efficacy of treatment. In this work, we develop an adaptive multi-level compactly supported radial basis approach for deformation analysis in 3D+time echocardiography. Our method combines displacement information from shape tracking of myocardial boundaries (derived from B-mode data) with mid-wall displacements from radio-frequency-based ultrasound speckle tracking. We evaluate our methods on open-chest canines (N=8) and show that our combined approach is better correlated to magnetic resonance tagging-derived strains than either individual method. We also are able to identify regions of myocardial infarction (confirmed by postmortem analysis) using radial strain values obtained with our approach.

Review of MR elastography applications and recent developments

The technique of MR elastography (MRE) has emerged as a useful modality for quantitatively imaging the mechanical properties of soft tissues in vivo. Recently, MRE has been introduced as a clinical tool for evaluating chronic liver disease, but many other potential applications are being explored. These applications include measuring tissue changes associated with diseases of the liver, breast, brain, heart, and skeletal muscle including both focal lesions (e.g., hepatic, breast, and brain tumors) and diffuse diseases (e.g., fibrosis and multiple sclerosis). The purpose of this review article is to summarize some of the recent developments of MRE and to highlight some emerging applications.

Assessment of myocardial elastography performance in phantoms under combined physiologic motion configurations with preliminary in vivo feasibility

Myocardial elastography (ME) is a non-invasive, ultrasound-based strain imaging technique, which can detect and localize abnormalities in myocardial function. By acquiring radio-frequency (RF) frames at high frame rates, the deformation of the myocardium can be estimated, and used to identify regions of abnormal deformation indicative of cardiovascular disease. In this study, the primary objective is to evaluate the effect of torsion on the performance of ME, while the secondary objective is to image inclusions during different motion schemes. Finally, the phantom findings are validated with an in vivo human case. Phantoms of homogeneous stiffness, or containing harder inclusions, were fixed to a pump and motors, and imaged. Incremental displacements were estimated from the RF signals, and accumulated over a motion cycle, and rotation angle, radial strain and circumferential strain were estimated. Phantoms were subjected to four motion schemes: rotation, torsion, deformation, and a combination of torsion and deformation. Sonomicrometry was used as a gold standard during deformation and combined motion schemes. In the rotation scheme, the input and estimated rotation angle agree in both the homogeneous and inclusion phantoms. In the torsion scheme, the estimated rotation angle was found to be highest, closest to the source of torsion and lowest farthest from the source of torsion. In the deformation scheme, if an inclusion was not present, the estimated strain patterns accurately depicted homogeneity, while if an inclusion was present, abnormalities were observed which enabled detection of the inclusion. In addition, no significant rotation was detected. In the combined scheme, if an inclusion was not present, the estimated strain patterns accurately depicted homogeneity, while, if an inclusion was present, abnormalities were observed which enabled detection of the inclusion. Also, torsion was separated from the combined scheme and was found to be similar to the pure torsion findings. This study shows ME to be capable of accurately depicting and distinguishing between different types of motion schemes, and to be sensitive to stiffness changes in localized regions of tissue-mimicking phantoms under physiologic cardiac motion configurations, while strains estimated in the combined motion scheme were noisier than in individual motion schemes. Finally, ME was shown to be capable of distinguishing between deformation and rotation in a normal human heart in vivo.

Imaging feedback of histotripsy treatments using ultrasound shear wave elastography

Histotripsy is a cavitation-based ultrasound therapy that mechanically fractionates soft solid tissues into fluid-like homogenates. This paper investigates the feasibility of imaging the tissue elasticity change during the histotripsy process as a tool to provide feedback for the treatments. The treatments were performed on agar tissue phantoms and ex vivo kidneys using 3-cycle ultrasound pulses delivered by a 750-kHz therapeutic array at peak negative/positive pressure of 17/108 MPa and a repetition rate of 50 Hz. Lesions with different degrees of damage were created with increasing numbers of therapy pulses from 0 to 2000 pulses per treatment location. The elasticity of the lesions was measured with ultrasound shear wave elastography, in which a quasi-planar shear wave was induced by acoustic radiation force generated by the therapeutic array, and tracked with ultrasound imaging at 3000 frames per second. Based on the shear wave velocity calculated from the sequentially captured frames, the Young's modulus was reconstructed. Results showed that the lesions were more easily identified on the shear wave velocity images than on B-mode images. As the number of therapy pulses increased from 0 to 2000 pulses/location, the Young's modulus decreased exponentially from 22.1 ± 2.7 to 2.1 ± 1.1 kPa in the tissue phantoms (R2 = 0.99, N = 9 each), and from 33.0 ± 7.1 to 4.0 ± 2.5 kPa in the ex vivo kidneys (R2 = 0.99, N = 8 each). Correspondingly, the tissues transformed from completely intact to completely fractionated as examined via histology. A good correlation existed between the lesions' Young's modulus and the degree of tissue fractionation as examined with the percentage of remaining structurally intact cell nuclei (R2 = 0.91, N = 8 each). These results indicate that lesions produced by histotripsy can be detected with high sensitivity using shear wave elastography. Because the decrease in the tissue elasticity corresponded well with the morphological and histological change, this study provides a basis for predicting the local treatment outcomes from tissue elasticity change.

Physiologic cardiovascular strain and intrinsic wave imaging

Cardiovascular disease remains the primary killer worldwide. The heart, essentially an electrically driven mechanical pump, alters its mechanical and electrical properties to compensate for loss of normal mechanical and electrical function. The same adjustment also is performed in the vessels, which constantly adapt their properties to accommodate mechanical and geometrical changes related to aging or disease. Real-time, quantitative assessment of cardiac contractility, conduction, and vascular function before the specialist can visually detect it could be feasible. This new physiologic data could open up interactive therapy regimens that are currently not considered. The eventual goal of this technology is to provide a specific method for estimating the position and severity of contraction defects in cardiac infarcts or angina. This would improve care and outcomes as well as detect stiffness changes and overcome the current global measurement limitations in the progression of vascular disease, at little more cost or risk than that of a clinical ultrasound.

In Vivo time harmonic elastography of the human heart

Time harmonic elastography is introduced as a modality for assessing myocardial elasticity changes during the cardiac cycle. It is based on external stimulation and real-time analysis of 30-Hz harmonic shear waves in axial direction of a parasternal line of sight through the lateral heart wall. In 20 healthy volunteers, the externally induced waves showed smaller amplitudes during systole (76.0 ± 30.8 μm) and higher amplitudes during diastole (126.7 ± 52.1 μm). This periodic wave amplitude alteration preceded ventricular contraction and dilation by about 100 ms. The amplitude ratio of 1.75 ± 0.49 indicates a relative change in myocardial shear elasticity on the order of 14 ± 11. These results well agree with observations made by cardiac magnetic resonance elastography for a similar displacement component and region of the heart. The proposed method provides reproducible elastodynamic information on the heart in real-time and may help in diagnosing myocardial relaxation abnormalities in the future.

Elevated Strain and Structural Disarray Occur at the Right Ventricular Apex

The right ventricular apex (RVA) is a potential hot spot for development of cardiac rhythm anomalies. Many conditions, including arrhythmogenic right ventricular cardiomyopathy and Brugada's syndrome affect the RVA, and further, the RVA remains an incompletely characterized pacing region. Whether there are structural reasons underlying these conduction properties remains unsettled. In the current study, we characterize the mechanical strains and structural attributes of the right ventricular wall, and test the hypothesis that the right ventricular apex experiences heterogeneous strain distributions and altered fiber organization, and is thus susceptible to conduction alterations. Electromechanical wave imaging (EWI), or elastography, of hearts was used to quantify mechanical strains occurring through a cardiac cycle. Histological and immunofluorescence imaging techniques were used to examine cardiac wall structure and arrangement of junctional proteins. Right ventricular mechanical strains were elevated and sustained throughout systole, compared to the left ventricle and septum. Heterogeneous strain distributions, myocardial fiber disarray, and altered junctional protein localization occured at the RVA. Disarray and altered strain distributions suggest decreased structural strength at the right ventricular apex in particular and increased mechanical impositions in the right ventricle, respectively. Thus, these data demonstrate why the right ventricular apex may be particularly vulnerable to conduction abnormalities.

Imaging of wall motion coupled with blood flow velocity in the heart and vessels in vivo: a feasibility study

The mechanical property and geometry changes as a result of cardiovascular disease affect both the wall motion and blood flow in the heart and vessels, whereas the latter two are also coupled and therefore continuously influence one another. Simultaneous and registered imaging of both cardiovascular wall motion and blood velocity may thus contribute to more complete computational models of cardiovascular mechanical and fluid dynamics as well as provide additional diagnostic information. The objective of this paper was to determine the feasibility of imaging cardiovascular wall motion coupled with blood flow in vivo. Normal (n = 6) and infarcted (n = 5) murine left ventricles, and normal (n = 5) and aneurysmal (n = 4) murine abdominal aortas, were imaged in longitudinal views with a 30-MHz ultrasound probe. Using electrocardiogram (ECG) gating, 2-D radio-frequency (RF) data were acquired at a frame rate of 8 kHz. The axial wall velocity and blood velocity were estimated using a speckle-tracking technique. Spatially and temporally registered imaging of both cardiovascular wall motion and blood flow was shown to be feasible. Reduced wall motion was detected in the infarcted region, whereas vortex flow patterns were imaged in diastolic phases of both normal and infarcted left ventricles. The myocardial wall motion and blood flow were found to be more synchronous in the normal heart, where the blood moves toward the anteroseptal wall after the mitral valve opens (i.e., rapid filling phase), and the anteroseptal wall simultaneously undergoes outward motion. In the infarcted heart, however, in the rapid filling phase, the basal anteroseptal wall starts moving about 20 ms before the mitral valve opens and the blood enters the left ventricle. In the normal aorta, the wall motion and blood velocity were uniform and synchronous. In the aneurysmal aorta, reduced and spatially varied wall motion and vortex flow patterns in the aneurysmal sac were found. The wall motion and blood velocity were thus less synchronous in the aneurysmal aorta. Cardiovascular wall motion and blood flow were both imaged in mice in vivo. This dual information may provide important insights for the diagnosis of cardiovascular disease as well as essential parameters for its biomechanical modeling.

Imaging the electromechanical activity of the heart in vivo

Cardiac conduction abnormalities remain a major cause of death and disability worldwide. However, as of today, there is no standard clinical imaging modality that can noninvasively provide maps of the electrical activation. In this paper, electromechanical wave imaging (EWI), a novel ultrasound-based imaging method, is shown to be capable of mapping the electromechanics of all four cardiac chambers at high temporal and spatial resolutions and a precision previously unobtainable in a full cardiac view in both animals and humans. The transient deformations resulting from the electrical activation of the myocardium were mapped in 2D and combined in 3D biplane ventricular views. EWI maps were acquired during five distinct conduction configurations and were found to be closely correlated to the electrical activation sequences. EWI in humans was shown to be feasible and capable of depicting the normal electromechanical activation sequence of both atria and ventricles. This validation of EWI as a direct, noninvasive, and highly translational approach underlines its potential to serve as a unique imaging tool for the early detection, diagnosis, and treatment monitoring of arrhythmias through ultrasound-based mapping of the transmural electromechanical activation sequence reliably at the point of care, and in real time.

Mapping of cardiac electrical activation with electromechanical wave imaging: an in silico-in vivo reciprocity study

BACKGROUND

Electromechanical wave imaging (EWI) is an entirely noninvasive, ultrasound-based imaging method capable of mapping the electromechanical activation sequence of the ventricles in vivo. Given the broad accessibility of ultrasound scanners in the clinic, the application of EWI could constitute a flexible surrogate for the 3-dimensional electrical activation.

OBJECTIVE

The purpose of this report is to reproduce the electromechanical wave (EW) using an anatomically realistic electromechanical model, and establish the capability of EWI to map the electrical activation sequence in vivo when pacing from different locations.

METHODS

EWI was performed in 1 canine during pacing from 3 different sites. A high-resolution dynamic model of coupled cardiac electromechanics of the canine heart was used to predict the experimentally recorded electromechanical wave. The simulated 3-dimensional electrical activation sequence was then compared with the experimental EW.

RESULTS

The electrical activation sequence and the EW were highly correlated for all pacing sites. The relationship between the electrical activation and the EW onset was found to be linear, with a slope of 1.01 to 1.17 for different pacing schemes and imaging angles.

CONCLUSION

The accurate reproduction of the EW in simulations indicates that the model framework is capable of accurately representing the cardiac electromechanics and thus testing new hypotheses. The one-to-one correspondence between the electrical activation and the EW sequences indicates that EWI could be used to map the cardiac electrical activity. This opens the door for further exploration of the technique in assisting in the early detection, diagnosis, and treatment monitoring of rhythm dysfunction.