In Vivo time harmonic elastography of the human heart

Explorative study using ultrasound time-harmonic elastography for stiffness-based quantification of skeletal muscle function

Time-harmonic elastography (THE) is an emerging ultrasound imaging technique that allows full-field mapping of the stiffness of deep biological tissues. THE's unique ability to rapidly capture stiffness in multiple tissues has never been applied for imaging skeletal muscle. Therefore, we addressed the lack of data on temporal changes in skeletal muscle stiffness while simultaneously covering stiffness of different muscles. Acquiring repeated THE scans every five seconds we quantified shear-wave speed (SWS) as a marker of stiffness of the long head (LHB) and short head (SHB) of biceps brachii and of the brachialis muscle (B) in ten healthy volunteers. SWS was continuously acquired during a 3-min isometric preloading phase, a 3-min loading phase with different weights (4, 8, and 12 kg), and a 9-min postloading phase. In addition, we analyzed temporal SWS standard deviation (SD) as a marker of muscle contraction regulation. Our results (median [min, max]) showed both SWS at preloading (LHB: 1.04 [0.94, 1.12] m/s, SHB: 0.86 [0.78, 0.94] m/s, B: 0.96 [0.87, 1.09] m/s, p < 0.001) and the increase in SWS with loading weight to be muscle-specific (LHB: 0.010 [0.002, 0.019] m/s/kg, SHB: 0.022 [0.017, 0.042] m/s/kg, B: 0.039 [0.019, 0.062] m/s/kg, p < 0.001). Additionally, SWS during loading increased continuously over time by 0.022 [0.004, 0.051] m/s/min (p < 0.01). Using an exponential decay model, we found an average relaxation time of 27 seconds during postloading. Analogously, SWS SD at preloading was also muscle-specific (LHB: 0.018 [0.011, 0.029] m/s, SHB: 0.021 [0.015, 0.027] m/s, B: 0.024 [0.018, 0.037] m/s, p < 0.05) and increased by 0.005 [0.003, 0.008] m/s/kg (p < 0.01) with loading. SWS SD did not change over loading time and decreased immediately in the postloading phase. Taken together, THE of skeletal muscle is a promising imaging technique for in vivo quantification of stiffness and stiffness changes in multiple muscle groups within seconds. Both the magnitude of stiffness changes and their temporal variation during isometric exercise may reflect the functional status of skeletal muscle and provide additional information to the morphological measures obtained by conventional imaging modalities.

Investigating the effects of microstructural changes induced by myocardial infarction on the elastic parameters of the heart

Within this work, we investigate how physiologically observed microstructural changes induced by myocardial infarction impact the elastic parameters of the heart. We use the LMRP model for poroelastic composites (Miller and Penta in Contin Mech Thermodyn 32:1533-1557, 2020) to describe the microstructure of the myocardium and investigate microstructural changes such as loss of myocyte volume and increased matrix fibrosis as well as increased myocyte volume fraction in the areas surrounding the infarct. We also consider a 3D framework to model the myocardium microstructure with the addition of the intercalated disks, which provide the connections between adjacent myocytes. The results of our simulations agree with the physiological observations that can be made post-infarction. That is, the infarcted heart is much stiffer than the healthy heart but with reperfusion of the tissue it begins to soften. We also observe that with the increase in myocyte volume of the non-damaged myocytes the myocardium also begins to soften. With a measurable stiffness parameter the results of our model simulations could predict the range of porosity (reperfusion) that could help return the heart to the healthy stiffness. It would also be possible to predict the volume of the myocytes in the area surrounding the infarct from the overall stiffness measurements.

Assessing cardiac stiffness using ultrasound shear wave elastography

Shear wave elastography offers a new dimension to echocardiography: it measures myocardial stiffness. Therefore, it could provide additional insights into the pathophysiology of cardiac diseases affecting myocardial stiffness and potentially improve diagnosis or guide patient treatment. The technique detects fast mechanical waves on the heart wall with high frame rate echography, and converts their propagation velocity into a stiffness value. A proper interpretation of shear wave data is required as the shear wave interacts with the intrinsic, yet dynamically changing geometrical and material characteristics of the heart under pressure. This dramatically alters the wave physics of the propagating wave, demanding adapted processing methods compared to other shear wave elastography applications as breast tumor and liver stiffness staging. Furthermore, several advanced analysis methods have been proposed to extract supplementary material features such as viscosity and anisotropy, potentially offering additional diagnostic value. This review explains the general mechanical concepts underlying cardiac shear wave elastography and provides an overview of the preclinical and clinical studies within the field. We also identify the mechanical and technical challenges ahead to make shear wave elastography a valuable tool for clinical practice.

2-D Ultrasound Shear Wave Elastography With Multi-Sphere-Source External Mechanical Vibration: Preliminary Phantom Results

Ultrasound shear wave elastography (SWE) imaging is emerging as a quantitative and non-invasive tissue characterization modality. Shear wave generation using external mechanical vibration (EMV) has received extensive research interest over acoustic radiation force impulse (ARFI) because of its low cost and potential for portability. In this paper, we propose an EMV concept with multiple spherical sources that can be easily reconfigured in three configurations to induce unique shear wave propagation patterns. We introduce two design embodiments of this concept bench test design for proof of concept and a clinically deployable design. The latter is designed to incorporate size, ergonomics, portability and power consumption considerations and constraints. Experimental validation on elasticity phantoms using both EMV designs demonstrates shear wave generation and elasticity reconstruction comparable in performance to ElastQ, a commercial ARFI-based shear elastography technology from Philips. In addition, the local displacement amplitude induced by EMV is 10 times greater than that induced by ARFI at the same given depth. Finally, the multiple configurations of the presented EMV design would allow exploration of advanced elastography methods such as tissue anisotropic elasticity.

Parasternal Versus Apical View in Cardiac Natural Mechanical Wave Speed Measurements

Shear wave speed measurements can potentially be used to noninvasively measure myocardial stiffness to assess the myocardial function. Several studies showed the feasibility of tracking natural mechanical waves induced by aortic valve closure in the interventricular septum, but different echocardiographic views have been used. This article systematically studied the wave propagation speeds measured in a parasternal long-axis and in an apical four-chamber view in ten healthy volunteers. The apical and parasternal views are predominantly sensitive to longitudinal or transversal tissue motion, respectively, and could, therefore, theoretically measure the speed of different wave modes. We found higher propagation speeds in apical than in the parasternal view (median of 5.1 m/s versus 3.8 m/s, , n = 9 ). The results in the different views were not correlated ( r = 0.26 , p = 0.49 ) and an unexpectedly large variability among healthy volunteers was found in apical view compared with the parasternal view (3.5-8.7 versus 3.2-4.3 m/s, respectively). Complementary finite element simulations of Lamb waves in an elastic plate showed that different propagation speeds can be measured for different particle motion components when different wave modes are induced simultaneously. The in vivo results cannot be fully explained with the theory of Lamb wave modes. Nonetheless, the results suggest that the parasternal long-axis view is a more suitable candidate for clinical diagnosis due to the lower variability in wave speeds.

Reproducibility of Natural Shear Wave Elastography Measurements

For the quantification of myocardial function, myocardial stiffness can potentially be measured non-invasively using shear wave elastography. Clinical diagnosis requires high precision. In 10 healthy volunteers, we studied the reproducibility of the measurement of propagation speeds of shear waves induced by aortic and mitral valve closure (AVC, MVC). Inter-scan was slightly higher but in similar ranges as intra-scan variability (AVC: 0.67 m/s (interquartile range [IQR]: 0.40-0.86 m/s) versus 0.38 m/s (IQR: 0.26-0.68 m/s), MVC: 0.61 m/s (IQR: 0.26-0.94 m/s) versus 0.26 m/s (IQR: 0.15-0.46 m/s)). For AVC, the propagation speeds obtained on different day were not statistically different (p = 0.13). We observed different propagation speeds between 2 systems (AVC: 3.23-4.25 m/s [Zonare ZS3] versus 1.82-4.76 m/s [Philips iE33]), p = 0.04). No statistical difference was observed between observers (AVC: p = 0.35). Our results suggest that measurement inaccuracies dominate the variabilities measured among healthy volunteers. Therefore, measurement precision can be improved by averaging over multiple heartbeats.

In vivo time-harmonic ultrasound elastography of the human brain detects acute cerebral stiffness changes induced by intracranial pressure variations

Cerebral stiffness (CS) reflects the biophysical environment in which neurons grow and function. While long-term CS changes can occur in the course of chronic neurological disorders and aging, little is known about acute variations of CS induced by intracranial pressure variations. Current gold standard methods for CS and intracranial pressure such as magnetic resonance elastography and direct pressure recordings are either expensive and slow or invasive. The study objective was to develop a real-time method for in vivo CS measurement and to demonstrate its sensitivity to physiological aging and intracranial pressure variations induced by the Valsalva maneuver in healthy volunteers. We used trans-temporal ultrasound time-harmonic elastography (THE) with external shear-wave stimulation by continuous and superimposed vibrations in the frequency range from 27 to 56 Hz. Multifrequency wave inversion generated maps of shear wave speed (SWS) as a surrogate maker of CS. On average, cerebral SWS was 1.56 ± 0.08 m/s with a tendency to reduce with age (R = -0.76, p < 0.0001) while Valsalva maneuver induced an immediate stiffening of the brain as reflected by a 10.8 ± 2.5% increase (p < 0.0001) in SWS. Our results suggest that CS is tightly linked to intracranial pressure and might be used in the future as non-invasive surrogate marker for intracranial pressure, which otherwise requires invasive measurements.

Time-Harmonic Ultrasound elastography of the Descending Abdominal Aorta: Initial Results

Stiffening of central large vessels is considered a key pathophysiologic factor within the cardiovascular system. Current diagnostic parameters such as pulse wave velocity (PWV) indirectly measure aortic stiffness, a hallmark of coronary diseases. The aim of the present study was to perform elastography of the proximal abdominal aorta based on externally induced time-harmonic shear waves. Experiments were performed in 30 healthy volunteers (25 young, 5 old, >50 y) and 5 patients with longstanding hypertension (PWV >10 m/s). B-Mode-guided sonographic time-harmonic elastography was used for measurement of externally induced shear waves at 30-Hz vibration frequency. Thirty-hertz shear wave amplitudes (SWAs) within the abdominal aorta were measured and displayed in real time and processed offline for differences in SWA between systole and diastole (ΔSWA). Data were analyzed using the Kruskal-Wallis test and receiver operating characteristic curve analysis. The change in SWA over the cardiac cycle was reduced significantly in all patients as assessed with ΔSWA (volunteers: mean = 10 ± 5 μm, patients: mean = 4 ± 1 μm; p < 0.001). The best separation of healthy volunteers from patients was obtained with a ΔSWA threshold of 4.7 μm, resulting in a sensitivity of 0.9 and a specificity of 1.0, with an overall area under the curve of 0.96. Time harmonic elastography of the abdominal aorta is feasible and shows promise for the exploitation of time-varying shear wave amplitudes as a diagnostic marker for aortic wall stiffening. Patients with elevated PWVs suggesting increased aortic wall stiffness were best identified by ΔSWA-a parameter that could be related to the ability of the vessel walls to distend on passages of the pulse wave.

Ultrasound palpation for fast in-situ quantification of articular cartilage stiffness, thickness and relaxation capacity

Most current cartilage testing devices require the preparation of excised samples and therefore do not allow intra-operative application for diagnostic purposes. The gold standard during open or arthroscopic surgery is still the subjective perception of manual palpation. This work presents a new diagnostic method of ultrasound palpation (USP) to acquire applied stress and strain data during manual palpation of articular cartilage. With the proposed method, we obtain cartilage thickness and stiffness. Moreover, repeated palpations allow the quantification of relaxation effects. USP measurements on elastomer phantoms demonstrated very good repeatability for both, stage-guided (97.2%) and handheld (96.0%) applications. The USP measurements were compared with conventional indentation experiments and revealed very good agreement on elastomer phantoms ([Formula: see text]) and good agreement on porcine cartilage samples ([Formula: see text]). Artificially degenerated cartilage samples showed reduced stiffness, weak capacity to relax after palpation and an increase of stiffness of approximately 50% with each single palpation. Intact cartilage was measured by USP directly at the patella (in situ) and after excision and removal of the subchondral bone (ex situ), leading to stiffness values of [Formula: see text] and [Formula: see text] ([Formula: see text]), respectively. The results demonstrate the potential of the USP system for cartilage testing, its sensitivity to degenerative changes and as a method for quantifying relaxation processes by means of repeated palpations. Furthermore, the differences in the results of in-situ and ex-situ measurements are of general interest, since such comparison has not been reported previously. We point out the limited comparability of ex-situ cartilage with its in-situ biomechanical behavior.

A Comparison of Acoustic Radiation Force-Derived Indices of Cardiac Function in the Langendorff Perfused Rabbit Heart

In the past decade, there has been an increased interest in characterizing cardiac tissue mechanics utilizing newly developed ultrasound-based elastography techniques. These methods excite the tissue mechanically and track the response. Two frequently used methods, acoustic radiation force impulse (ARFI) and shear-wave elasticity imaging (SWEI), have been considered qualitative and quantitative techniques providing relative and absolute measures of tissue stiffness, respectively. Depending on imaging conditions, it is desirable to identify indices of cardiac function that could be measured by ARFI and SWEI and to characterize the relationship between the measures. In this study, we have compared two indices (i.e., relaxation time constant used for diastolic dysfunction assessment and systolic/diastolic stiffness ratio) measured nearly simultaneously by M-mode ARFI and SWEI techniques. We additionally correlated ARFI-measured inverse displacements with SWEI-measured values of the shear modulus of stiffness. For the eight animals studied, the average relaxation time constant ( τ) measured by ARFI and SWEI were ([Formula: see text]) and ([Formula: see text]), respectively ([Formula: see text]). Average systolic/diastolic stiffness ratios for ARFI and SWEI measurements were 6.01±1.37 and 7.12±3.24, respectively ([Formula: see text]). Shear modulus of stiffness (SWEI) was linearly related to inverse displacement values (ARFI) with a 95% CI for the slope of 0.010-0.011 [Formula: see text] ( R(2)=0.73). In conclusion, the relaxation time constant and the systolic/diastolic stiffness ratio were calculated with good agreement between the ARFI- and SWEI-derived measurements. ARFI relative and SWEI absolute stiffness measurements were linearly related with varying slopes based on imaging conditions and subject tissue properties.

Time Harmonic Elastography Reveals Sensitivity of Liver Stiffness to Water Ingestion

The aim of the study was to test the sensitivity of liver stiffness (LS) measured by time harmonic elastography in large tissue windows to water uptake and post-prandial effects. Each subject gave written informed consent to participate in this institutional review board-approved prospective study. LS was measured by time harmonic elastography in 10 healthy volunteers pre- and post-prandially, as well as before, directly after and 2 h after drinking water. The LS-time function during water intake was measured in 14 scans over 3 h in five volunteers. LS increased by 10% (p = 0.0015) post-prandially and by 11% (p = 0.0024) after pure water ingestion, and decreased to normal values after 2 h. LS was lower after overnight fasting than after 2-h fasting (3%, p = 0.04). Over the time course, LS increased to post-water peak values 15 min after drinking 0.25 L water and remained unaffected by further ingestion of water. In conclusion, our study indicates that LS measured by time harmonic elastography represents an effective-medium property sensitive to physiologic changes in vascular load of the liver.

Investigation of the effects of myocardial anisotropy for shear wave elastography using impulsive force and harmonic vibration

The myocardium is known to be an anisotropic medium where the muscle fiber orientation changes through the thickness of the wall. Shear wave elastography methods use propagating waves which are measured by ultrasound or magnetic resonance imaging (MRI) techniques to characterize the mechanical properties of various tissues. Ultrasound- or MR-based methods have been used and the excitation frequency ranges for these various methods cover a large range from 24-500 Hz. Some of the ultrasound-based methods have been shown to be able to estimate the fiber direction. We constructed a model with layers of elastic, transversely isotropic materials that were oriented at different angles to simulate the heart wall in systole and diastole. We investigated the effect of frequency on the wave propagation and the estimation of fiber direction and wave speeds in the different layers of the assembled models. We found that waves propagating at low frequencies such as 30 or 50 Hz showed low sensitivity to the fiber direction but also had substantial bias in estimating the wave speeds in the layers. Using waves with higher frequency content (>200 Hz) allowed for more accurate fiber direction and wave speed estimation. These results have particular relevance for MR- and ultrasound-based elastography applications in the heart.

Multifrequency time-harmonic elastography for the measurement of liver viscoelasticity in large tissue windows

Elastography of the liver for the non-invasive diagnosis of hepatic fibrosis is an established method. However, investigations of obese patients or patients with ascites are often limited by small and superficial elastographic windows. Therefore, multifrequency time-harmonic elastography (THE) based on time-resolved A-line ultrasound has recently been developed for measuring liver viscoelasticity in wide soft tissue windows and at greater depths. In this study, THE was integrated into a clinical B-mode scanner connected to a dedicated actuator bed driven by superimposed vibrations of 30- to 60-Hz frequencies. The resulting shear waves in the liver were captured along multiple profiles 7 to 14 cm in width and automatically processed for reconstruction of mean efficient shear wave speed and shear wave dispersion slope. This new modality was tested in healthy volunteers and 22 patients with clinically proven cirrhosis. Patients could be separated from controls by higher shear wave speeds (3.11 ± 0.64 m/s, 2.14-4.81 m/s, vs. 1.74 ± 0.10 m/s, 1.60-1.91 m/s) without significant degradation of data by high body mass index or ascites. Furthermore, the wave speed dispersion slope was significantly (p = 0.0025) lower in controls (5.2 ± 1.8 m/s/kHz) than in patients (39.1 ± 32.2 m/s/kHz). In conclusion, THE is useful for the diagnosis of cirrhosis in large tissue windows.

Shear-wave amplitudes measured with cardiac MR elastography for diagnosis of diastolic dysfunction

PURPOSE

To test whether shear-wave amplitudes (SWAs) in the myocardium measured with cardiac magnetic resonance (MR) elastography enable diagnosis of myocardial relaxation abnormalities in patients with diastolic dysfunction.

MATERIALS AND METHODS

Each subject gave written informed consent to participate in this institutional review board-approved prospective study. Electrocardiographically triggered SWA-based cardiac MR elastography with 24.13-Hz external vibration frequency was performed in 50 subjects grouped into asymptomatic young (n = 10, 18-39 years) and asymptomatic old (n = 10, 40-68 years) subjects and patients with echocardiographically proved mild, moderate, or severe diastolic dysfunction (n = 30, 44-73 years). SWA images were analyzed in the left ventricular (LV) region and were normalized against reference SWA of the thoracic wall. Analysis of variance with Bonferroni-corrected pairwise comparison and Pearson correlation were used for statistical evaluation.

RESULTS

Young and old control subjects had normalized mean LV SWA of 0.67 ± 0.04 (standard error of the mean) and 0.56 ± 0.04 (P = .18, F test), respectively. Compared with the control groups, patients with mild, moderate, and severe diastolic dysfunction displayed significantly reduced normalized mean LV SWA of 0.37 ± 0.04, 0.34 ± 0.04, and 0.29 ± 0.04 (P < .001, F test), respectively, which was inversely correlated to the severity of diastolic dysfunction (R = -0.61, P < .001). The best cutoff value to differentiate between asymptomatic volunteers and patients was 0.43, yielding an area under the receiver operating characteristic curve of 0.92, with 90% sensitivity and 89.7% specificity.

CONCLUSION

LV SWA measured with cardiac MR elastography provides image contrast sensitive to myocardial relaxation abnormalities and shows significantly lower values in patients with diastolic dysfunction.

Isovolumetric elasticity alteration in the human heart detected by in vivo time-harmonic elastography

Time harmonic elastography (THE) has recently been introduced for measurement of the periodic alteration in myocardial shear modulus based on externally induced low-frequency acoustic vibrations produced by a loudspeaker. In this study, we propose further developments of cardiac THE toward a clinical modality including integration of the vibration source into the patient bed and automated parameter extraction from harmonic shear wave amplitudes, wall motion profiles and synchronized electrocardiographic records. This method has enabled us to evaluate the delay between wall motion and wave amplitude alteration for the measurement of isovolumetric times of elasticity alteration during contraction (τ(C)) and relaxation (τ(R)) in a group of 32 healthy volunteers. On average, the wave amplitudes changed between systole and diastole by a factor of 1.7 ± 0.3, with a τ(C) of 137 ± 61 ms and a τ(R) of 68 ± 73 ms, which agrees with results obtained with the more time-consuming and expensive cardiac magnetic resonance elastography. Furthermore, because of the high sampling rate, elasto-morphometric parameters such as transition times and the area of wave amplitude-cardiac motion cycles can be processed in an automated way for the future clinical detection of myocardial relaxation abnormalities.