Search

Objective. In proton therapy, range uncertainties prevent optimal benefit from the superior depth-dose characteristics of proton beams over conventional photon-based radiotherapy. To reduce these uncertainties we recently proposed the use of phase-change ultrasound contrast agents as an affordable and effective range verification tool. In particular, superheated nanodroplets can convert into echogenic microbubbles upon proton irradiation, whereby the resulting ultrasound contrast relates to the proton range with high reproducibility. Here, we provide a first in vivo proof-of-concept of this technology. Approach. First, the in vitro biocompatibility of radiation-sensitive poly(vinyl alcohol) perfluorobutane nanodroplets was investigated using several colorimetric assays. Then, in vivo ultrasound contrast was characterized using acoustic droplet vaporization (ADV) and later using proton beam irradiations at varying energies (49.7 MeV and 62 MeV) in healthy Sprague Dawley rats. A preliminary evaluation of the in vivo biocompatibility was performed using ADV and a combination of physiology monitoring and histology. Main results. Nanodroplets were non-toxic over a wide concentration range (<1 mM). In healthy rats, intravenously injected nanodroplets primarily accumulated in the organs of the reticuloendothelial system, where the lifetime of the generated ultrasound contrast (<30 min) was compatible with a typical radiotherapy fraction (<5 min). Spontaneous droplet vaporization did not result in significant background signals. Online ultrasound imaging of the liver of droplet-injected rats demonstrated an energy-dependent proton response, which can be tuned by varying the nanodroplet concentration. However, caution is warranted when deciding on the exact nanodroplet dose regimen as a mild physiological response (drop in cardiac rate, granuloma formation) was observed after ADV. Significance. These findings underline the potential of phase-change ultrasound contrast agents for in vivo proton range verification and provide the next step towards eventual clinical applications., (© 2024 Institute of Physics and Engineering in Medicine. All rights, including for text and data mining, AI training, and similar technologies, are reserved.)

Aims: By combining temporal changes in left ventricular (LV) global longitudinal strain (GLS) with LV volume, LV strain-volume loops can assess cardiac function across the cardiac cycle. This study compared LV strain-volume loops between bicuspid aortic valve (BAV) patients and controls, and investigated the loop's prognostic value for clinical events., Methods and Results: From a prospective cohort of congenital heart disease patients, BAV patients were selected and compared with healthy volunteers, who were matched for age and sex at group level. GLS analysis from apical views was used to construct strain-volume loops. Associations with clinical events, i.e. a composite of all-cause mortality, heart failure, arrhythmias, and aortic valve replacement, were assessed by Cox regression. A total of 113 BAV patients were included (median age 32 years, 40% female). BAV patients demonstrated lower Sslope (0.21%/mL, [Q1-Q3: 0.17-0.28] vs. 0.27%/mL [0.24-0.34], P < 0.001) and ESslope (0.19%/mL [0.12-0.25] vs. 0.29%/mL [0.21-0.43], P < 0.001) compared with controls, but also greater uncoupling during early (0.48 ± 1.29 vs. 0.05 ± 1.21, P = 0.04) and late diastole (0.66 ± 1.02 vs. -0.07 ± 1.07, P < 0.001). Median follow-up duration was 9.9 [9.3-10.4] years. Peak aortic jet velocity (HR 1.22, P = 0.03), enlarged left atrium (HR 3.16, P = 0.003), E/e' ratio (HR 1.17, P = 0.002), GLS (HR 1.16, P = 0.008), and ESslope (HR 0.66, P = 0.04) were associated with the occurrence of clinical events., Conclusion: Greater uncoupling and lower systolic and diastolic slopes were observed in BAV patients compared with healthy controls, suggesting presence of altered LV cardiomechanics. Moreover, lower ESslope was associated with clinical events, highlighting the strain-volume loop's potential as prognostic marker., Competing Interests: Conflict of interest: None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Objective . An elevated interstitial fluid pressure (IFP) can lead to strain-induced stiffening of poroelastic biological tissues. As shear wave elastography (SWE) measures functional tissue stiffness based on the propagation speed of acoustically induced shear waves, the shear wave velocity (SWV) can be used as an indirect measurement of the IFP. The underlying biomechanical principle for this stiffening behavior with pressurization is however not well understood, and we therefore studied how IFP affects SWV through SWE experiments and numerical modeling. Approach . For model set-up and verification, SWE experiments were performed while dynamically modulating IFP in a chicken breast. To identify the confounding factors of the SWV-IFP relationship, we manipulated the material model (linear poroelastic versus porohyperelastic), deformation assumptions (geometric linearity versus nonlinearity), and boundary conditions (constrained versus unconstrained) in a finite element model mimicking the SWE experiments. Main results . The experiments demonstrated a statistically significant positive correlation between the SWV and IFP. The model was able to reproduce a similar SWV-IFP relationship by considering an unconstrained porohyperelastic tissue. Material nonlinearity was identified as the primary factor contributing to this relationship, whereas geometric nonlinearity played a smaller role. The experiments also highlighted the importance of the dynamic nature of the pressurization procedure, as indicated by a different observed SWV-IFP for pressure buildup and relaxation, but its clinical relevance needs to be further investigated. Significance . The developed model provides an adaptable framework for SWE of poroelastic tissues and paves the way towards non-invasive measurements of IFP., (© 2024 Institute of Physics and Engineering in Medicine.)

Objective: Post-operative brain injury in neonates may result from disturbed cerebral perfusion, but accurate peri-operative monitoring is lacking. High-frame-rate (HFR) cerebral ultrasound could visualize and quantify flow in all detectable vessels using spectral Doppler; however, automated quantification in small vessels is challenging because of low signal amplitude. We have developed an automatic envelope detection algorithm for HFR pulsed wave spectral Doppler signals, enabling neonatal brain quantitative parameter maps during and after surgery., Methods: HFR ultrasound data from high-risk neonatal surgeries were recorded with a custom HFR mode (frame rate = 1000 Hz) on a Zonare ZS3 system. A pulsed wave Doppler spectrogram was calculated for each pixel containing blood flow in the image, and spectral peak velocity was tracked using a max-likelihood estimation algorithm of signal and noise regions in the spectrogram, where the most likely cross-over point marks the blood flow velocity. The resulting peak systolic velocity (PSV), end-diastolic velocity (EDV) and resistivity index (RI) were compared with other detection schemes, manual tracking and RIs from regular pulsed wave Doppler measurements in 10 neonates., Results: Envelope detection was successful in both high- and low-quality arterial and venous flow spectrograms. Our technique had the lowest root mean square error for EDV, PSV and RI (0.46 cm/s, 0.53 cm/s and 0.15, respectively) when compared with manual tracking. There was good agreement between the clinical pulsed wave Doppler RI and HFR measurement with a mean difference of 0.07., Conclusion: The max-likelihood algorithm is a promising approach to accurate, automated cerebral blood flow monitoring with HFR imaging in neonates., Competing Interests: Conflict of interest The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

The advent of high-frame rate imaging in ultrasound allowed the development of shear wave elastography as a noninvasive alternative for myocardial stiffness assessment. It measures mechanical waves propagating along the cardiac wall with speeds that are related to stiffness. The use of cardiac shear wave elastography in clinical studies is increasing, but a proper understanding of the different factors that affect wave propagation is required to correctly interpret results because of the heart's thin-walled geometry and intricate material properties. The aims of this review are to give an overview of the general concepts in cardiac shear wave elastography and to discuss in depth the effects of age, hemodynamic loading, cardiac morphology, fiber architecture, contractility, viscoelasticity, and system-dependent factors on the measurements, with a focus on clinical application. It also describes how these factors should be considered during acquisition, analysis, and reporting to ensure an accurate, robust, and reproducible measurement of the shear wave., Competing Interests: Funding Support and Author Disclosures This work was supported by the Research Foundation Flanders under grants 1211620N and 12B3124N to Dr Caenen and grant G092318N to Dr Bézy. Dr Voigt holds a personal research mandate of the Research Foundation Flanders (1832922N). Dr Nightingale has intellectual property related to radiation force-based imaging technologies that has been licensed to Siemens, Samsung, and MicroElastic Ultrasound Systems. Dr Pernot has intellectual property related to SWE that has been licensed to eMyoSound. Dr Vos has an ongoing research collaboration with Mindray Ultrasound. Drs Voigt and Jan D’hooge have an ongoing research collaboration with GE Vingmed. All other authors have reported that they have no relationships relevant to the contents of this paper to disclose., (Copyright © 2024. Published by Elsevier Inc.)

Objective: Described here is the development of an ultrasound matrix transducer prototype for high-frame-rate 3-D intra-cardiac echocardiography., Methods: The matrix array consists of 16 × 18 lead zirconate titanate elements with a pitch of 160 µm × 160 µm built on top of an application-specific integrated circuit that generates transmission signals and digitizes the received signals. To reduce the number of cables in the catheter to a feasible number, we implement subarray beamforming and digitization in receive and use a combination of time-division multiplexing and pulse amplitude modulation data transmission, achieving an 18-fold reduction. The proposed imaging scheme employs seven fan-shaped diverging transmit beams operating at a pulse repetition frequency of 7.7 kHz to obtain a high frame rate. The performance of the prototype is characterized, and its functionality is fully verified., Results: The transducer exhibits a transmit efficiency of 28 Pa/V at 5 cm per element and a bandwidth of 60% in transmission. In receive, a dynamic range of 80 dB is measured with a minimum detectable pressure of 10 Pa per element. The element yield of the prototype is 98%, indicating the efficacy of the manufacturing process. The transducer is capable of imaging at a frame rate of up to 1000 volumes/s and is intended to cover a volume of 70° × 70° × 10 cm., Conclusion: These advanced imaging capabilities have the potential to support complex interventional procedures and enable full-volumetric flow, tissue, and electromechanical wave tracking in the heart., Competing Interests: Conflict of interest The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Ultrafast cerebral Doppler ultrasound enables simultaneous quantification and visualization of cerebral blood flow velocity. The aim of this study is to compare the use of conventional and ultrafast spectral Doppler during anesthesia and their potential to show the effect of anesthesiologic procedures on cerebral blood flow velocities, in relation to blood pressure and cerebral oxygenation in infants undergoing inguinal hernia repair., Methods: A single-center prospective observational cohort study in infants up to six months of age. We evaluated conventional and ultrafast spectral Doppler cerebral ultrasound measurements in terms of number of successful measurements during the induction of anesthesia, after sevoflurane induction, administration of caudal analgesia, a fluid bolus and emergence of anesthesia. Cerebral blood flow velocity was quantified in pial arteries using conventional spectral Doppler and in the cerebral cortex using ultrafast Doppler by peak systolic velocity, end diastolic velocity and resistivity index., Results: Twenty infants were included with useable conventional spectral Doppler images in 72/100 measurements and ultrafast Doppler images in 51/100 measurements. Intraoperatively, the success rates were 53/60 (88.3%) and 41/60 (68.3%), respectively. Cerebral blood flow velocity increased after emergence for both conventional (end diastolic velocity, from 2.01 to 2.75 cm/s, p < 0.001) and ultrafast spectral Doppler (end diastolic velocity, from 0.59 to 0.94 cm/s), whereas cerebral oxygenation showed a reverse pattern with a decrease after the emergence of the infant (85% to 68%, p < 0.001)., Conclusion: It is possible to quantify cortical blood flow velocity during general anesthesia using conventional and ultrafast spectral Doppler cerebral ultrasound. Cerebral blood flow velocity and blood pressure decreased, while regional cerebral oxygenation increased during general anesthesia. Ultrafast spectral Doppler ultrasound offers novel insights into perfusion within the cerebral cortex, unattainable through conventional spectral ultrasound. Yet, ultrafast Doppler is curtailed by a lower success rate and a more rigorous learning curve compared to conventional method., Competing Interests: Declaration of Competing Interest The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (Copyright © 2023. Published by Elsevier Inc.)

Objective: The aim of this study was to assess the feasibility and imaging options of contrast-enhanced volumetric ultrasound kidney vasculature imaging in a porcine model using a prototype sparse spiral array., Methods: Transcutaneous freehand in vivo imaging of two healthy porcine kidneys was performed according to three protocols with different microbubble concentrations and transmission sequences. Combining high-frame-rate transmission sequences with our previously described spatial coherence beamformer, we determined the ability to produce detailed volumetric images of the vasculature. We also determined power, color and spectral Doppler, as well as super-resolved microvasculature in a volume. The results were compared against a clinical 2-D ultrasound machine., Results: Three-dimensional visualization of the kidney vasculature structure and blood flow was possible with our method. Good structural agreement was found between the visualized vasculature structure and the 2-D reference. Microvasculature patterns in the kidney cortex were visible with super-resolution processing. Blood flow velocity estimations were within a physiological range and pattern, also in agreement with the 2-D reference results., Conclusion: Volumetric imaging of the kidney vasculature was possible using a prototype sparse spiral array. Reliable structural and temporal information could be extracted from these imaging results., Competing Interests: Conflict of interest The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Assessing the coronary circulation with contrast-enhanced echocardiography has high clinical relevance. However, it is not being routinely performed in clinical practice because the current clinical tools generally cannot provide adequate image quality. The contrast agent's visibility in the myocardium is generally poor, impaired by motion and nonlinear propagation artifacts. The established multipulse contrast schemes (MPCSs) and the more experimental singular value decomposition (SVD) filter also fall short to solve these issues. Here, we propose a scheme to process amplitude modulation/amplitude-modulated pulse inversion (AM/AMPI) echoes with higher order SVD (HOSVD) instead of conventionally summing the complementary pulses. The echoes from the complementary pulses form a separate dimension in the HOSVD algorithm. Then, removing the ranks in that dimension with dominant coherent signals coming from tissue scattering would provide the contrast detection. We performed both in vitro and in vivo experiments to assess the performance of our proposed method in comparison with the current standard methods. A flow phantom study shows that HOSVD on AM pulsing exceeds the contrast-to-background ratio (CBR) of conventional AM and an SVD filter by 10 and 14 dB, respectively. In vivo porcine heart results also demonstrate that, compared to AM, HOSVD improves CBR in open-chest acquisition (up to 19 dB) and contrast ratio (CR) in closed-chest acquisition (3 dB).

Ultrasound-based shear wave elastography is a promising technique to non-invasively assess the dynamic stiffness variations of the heart. The technique is based on tracking the propagation of acoustically induced shear waves in the myocardium of which the propagation speed is linked to tissue stiffness. This measurement is repeated multiple times across the cardiac cycle to assess the natural variations in wave propagation speed. The interpretation of these measurements remains however complex, as factors such as loading and contractility affect wave propagation. We therefore applied transthoracic shear wave elastography in 13 pigs to investigate the dependencies of wave speed on pressure-volume derived indices of loading, myocardial stiffness, and contractility, while altering loading and inducing myocardial ischemia/reperfusion injury. Our results show that diastolic wave speed correlates to a pressure-volume derived index of operational myocardial stiffness (R = 0.75, p < 0.001), suggesting that both loading and intrinsic properties can affect diastolic wave speed. Additionally, the wave speed ratio, i.e. the ratio of systolic and diastolic speed, correlates to a pressure-volume derived index of contractility, i.e. preload-recruitable stroke work (R = 0.67, p < 0.001). Measuring wave speed ratio might thus provide a non-invasive index of contractility during ischemia/reperfusion injury., (© 2023. Springer Nature Limited.)

Background: The safety and efficacy of proton therapy is currently hampered by range uncertainties. The combination of ultrasound imaging with injectable radiation-sensitive superheated nanodroplets was recently proposed for in vivo range verification. The proton range can be estimated from the distribution of nanodroplet vaporization events, which is stochastically related to the stopping distribution of protons, as nanodroplets are vaporized by protons reaching their maximal LET at the end of their range., Purpose: Here, we aim to estimate the range estimation precision of this technique. As for any stochastic measurement, the precision will increase with the sample size, that is, the number of detected vaporizations. Thus, we first develop and validate a model to predict the number of vaporizations, which is then applied to estimate the range verification precision for a set of conditions (droplet size, droplet concentration, and proton beam parameters)., Methods: Starting from the thermal spike theory, we derived a model that predicts the expected number of droplet vaporizations in an irradiated sample as a function of the droplet size, concentration, and number of protons. The model was validated by irradiating phantoms consisting of size-sorted perfluorobutane droplets dispersed in an aqueous matrix. The number of protons was counted with an ionization chamber, and the droplet vaporizations were recorded and counted individually using high frame rate ultrasound imaging. After validation, the range estimate precision was determined for different conditions using a Monte Carlo algorithm., Results: A good agreement between theory and experiments was observed for the number of vaporizations, especially for large (5.8 ± 2.2 µm) and medium (3.5 ± 1.1 µm) sized droplets. The number of events was lower than expected in phantoms with small droplets (2.0 ± 0.7 µm), but still within the same order of magnitude. The inter-phantom variability was considerably larger (up to 30x) than predicted by the model. The validated model was then combined with Monte Carlo simulations, which predicted a theoretical range retrieval precision improving with the square-root of the number of vaporizations, and degrading at high beam energies due to range straggling. For single pencil beams with energies between 70 and 240 MeV, a range verification precision below 1% of the range required perfluorocarbon concentrations in the order of 0.3-2.4 µM., Conclusion: We proposed and experimentally validated a model to provide a quick estimate of the number of vaporizations for a given set of conditions (droplet size, droplet concentration, and proton beam parameters). From this model, promising range verification performances were predicted for realistic perfluorocarbon concentrations. These findings are an incentive to move towards preclinical studies, which are critical to assess the achievable droplet distribution in and around the tumor, and hence the in vivo range verification precision., (© 2023 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

High frame rate three-dimensional (3D) ultrasound imaging would offer excellent possibilities for the accurate assessment of carotid artery diseases. This calls for a matrix transducer with a large aperture and a vast number of elements. Such a matrix transducer should be interfaced with an application-specific integrated circuit (ASIC) for channel reduction. However, the fabrication of such a transducer integrated with one very large ASIC is very challenging and expensive. In this study, we develop a prototype matrix transducer mounted on top of multiple identical ASICs in a tiled configuration. The matrix was designed to have 7680 piezoelectric elements with a pitch of 300 μm × 150 μm integrated with an array of 8 × 1 tiled ASICs. The performance of the prototype is characterized by a series of measurements. The transducer exhibits a uniform behavior with the majority of the elements working within the -6 dB sensitivity range. In transmit, the individual elements show a center frequency of 7.5 MHz, a -6 dB bandwidth of 45%, and a transmit efficiency of 30 Pa/V at 200 mm. In receive, the dynamic range is 81 dB, and the minimum detectable pressure is 60 Pa per element. To demonstrate the imaging capabilities, we acquired 3D images using a commercial wire phantom.

Two-dimensional (2-D) arrays offer volumetric imaging capabilities without the need for probe translation or rotation. A sparse array with elements seeded in a tapering spiral pattern enables one-to-one connection to an ultrasound machine, thus allowing flexible transmission and reception strategies. To test the concept of sparse spiral array imaging, we have designed, realized, and characterized two prototype probes designed at 2.5-MHz low-frequency (LF) and 5-MHz high-frequency (HF) center frequencies. Both probes share the same electronic design, based on piezoelectric ceramics and rapid prototyping with printed circuit board substrates to wire the elements to external connectors. Different center frequencies were achieved by adjusting the piezoelectric layer thickness. The LF and HF prototype probes had 88% and 95% of working elements, producing peak pressures of 21 and 96 kPa/V when focused at 5 and 3 cm, respectively. The one-way -3-dB bandwidths were 26% and 32%. These results, together with experimental tests on tissue-mimicking phantoms, show that the probes are viable for volumetric imaging.

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is normally treated by RF ablation. Intracardiac echography (ICE) is widely employed during RF ablation procedures to guide the electrophysiologist in navigating the ablation catheter, although only 2-D probes are currently clinically used. A 3-D ICE catheter would not only improve visualization of the atrium and ablation catheter, but it might also provide the 3-D mapping of the electromechanical wave (EW) propagation pattern, which represents the mechanical response of cardiac tissue to electrical activity. The detection of this EW needs 3-D high-frame-rate imaging, which is generally only realizable in tradeoff with channel count and image quality. In this simulation-based study, we propose a high volume rate imaging scheme for a 3-D ICE probe design that employs 1-D micro-beamforming in the elevation direction. Such a probe can achieve a high frame rate while reducing the channel count sufficiently for realization in a 10-Fr catheter. To suppress the grating-lobe (GL) artifacts associated with micro-beamforming in the elevation direction, a limited number of fan-shaped beams with a wide azimuthal and narrow elevational opening angle are sequentially steered to insonify slices of the region of interest. An angular weighted averaging of reconstructed subvolumes further reduces the GL artifacts. We optimize the transmit beam divergence and central frequency based on the required image quality for EW imaging (EWI). Numerical simulation results show that a set of seven fan-shaped transmission beams can provide a frame rate of 1000 Hz and a sufficient spatial resolution to visualize the EW propagation on a large 3-D surface.

Suppressing tissue clutter is an essential step in blood flow estimation and visualization, even when using ultrasound contrast agents. Blind source separation (BSS)-based clutter filter for high-framerate ultrasound imaging has been reported to perform better in tissue clutter suppression than the conventional frequency-based wall filter and nonlinear contrast pulsing schemes. The most notable BSS technique, singular value decomposition (SVD) has shown compelling results in cases of slow tissue motion. However, its performance degrades when the tissue motion is faster than the blood flow speed, conditions that are likely to occur when imaging the small vessels, such as in the myocardium. Independent component analysis (ICA) is another BSS technique that has been implemented as a clutter filter in the spatiotemporal domain. Instead, we propose to implement ICA in the spatial domain where motion should have less impact. In this work, we propose a clutter filter with the combination of SVD and ICA to improve the contrast-to-background ratio (CBR) in cases where tissue velocity is significantly faster than the flow speed. In an in vitro study, the range of fast tissue motion velocity was 5-25 mm/s and the range of flow speed was 1-12 mm/s. Our results show that the combination of ICA and SVD yields 7-10 dB higher CBR than SVD alone, especially in the tissue high-velocity range. The improvement is crucial for cardiac imaging where relatively fast myocardial motions are expected.

Superheated nanodroplet (ND) vaporization by proton radiation was recently demonstrated, opening the door to ultrasound-based in vivo proton range verification. However, at body temperature and physiological pressures, perfluorobutane nanodroplets (PFB-NDs), which offer a good compromise between stability and radiation sensitivity, are not directly sensitive to primary protons. Instead, they are vaporized by infrequent secondary particles, which limits the precision for range verification. The radiation-induced vaporization threshold (i.e., sensitization threshold) can be reduced by lowering the pressure in the droplet such that ND vaporization by primary protons can occur. Here, we propose to use an acoustic field to modulate the pressure, intermittently lowering the proton sensitization threshold of PFB-NDs during the rarefactional phase of the ultrasound wave. Simultaneous proton irradiation and sonication with a 1.1 MHz focused transducer, using increasing peak negative pressures (PNPs), were applied on a dilution of PFB-NDs flowing in a tube, while vaporization was acoustically monitored with a linear array. Sensitization to primary protons was achieved at temperatures between [Formula: see text] and 40 °C using acoustic PNPs of relatively low amplitude (from 800 to 200 kPa, respectively), while sonication alone did not lead to ND vaporization at those PNPs. Sensitization was also measured at the clinically relevant body temperature (i.e., 37 °C) using a PNP of 400 kPa. These findings confirm that acoustic modulation lowers the sensitization threshold of superheated NDs, enabling a direct proton response at body temperature.

Background: Echo-Particle Image Velocimetry (echoPIV) tracks speckle patterns from ultrasound contrast agent(UCA), being less angle-sensitive than colour Doppler. High frame rate (HFR) echoPIV enables tracking of high velocity flow in the left ventricle (LV). We aimed to demonstrate the potential clinical use of HFR echoPIV and investigate the feasibility and accuracy in patients., Methods: Nineteen patients admitted for heart failure were included. HFR contrast images were acquired from an apical long axis view (ALAX), using a fully-programmable ultrasound system. A clinical UCA was continuously infused with a dedicated pump. Additionally, echocardiographic images were obtained using a clinical system, including LV contrast-enhanced images and pulsed-wave (PW) Doppler of the LV inflow and outflow in ALAX. 11 patients underwent CMR and 4 cardiac CT as clinically indicated. These CMR and CT images were used as reference. In 10 patients with good echoPIV tracking and reference imaging, the intracavitary flow was compared between echoPIV, conventional and UCA echocardiography., Results: EchoPIV tracking quality was good in 12/19 (63%), moderate in 2/19 (10%) and poor in 5/19 (26%) subjects. EchoPIV could determine inflow velocity in 17/19 (89%), and outflow in 14/19 (74%) patients. The correlation of echoPIV and PW Doppler was good for the inflow (R 2  = 0.77 to PW peak; R 2  = 0.80 PW mean velocity) and moderate for the outflow (R 2  = 0.54 to PW peak; R 2  = 0.44 to PW mean velocity), with a tendency for echoPIV to underestimate PW velocities. In selected patients, echoPIV was able in a single acquisition to demonstrate flow patterns which required multiple interrogations with classical echocardiography. Those flow patterns could also be linked to anatomical abnormalities as seen in CMR or CT., Conclusion: HFR echoPIV tracks multidirectional and complex flow patterns which are unapparent with conventional echocardiography, while having comparable feasibility. EchoPIV tends to underestimate flow velocities as compared to PW Doppler. It has the potential to provide in one acquisition all the functional information obtained by conventional imaging, overcoming the angle dependency of Doppler and low frame rate of classical contrast imaging., (© 2022. The Author(s).)

An in vivo range verification technology for proton beam cancer therapy, preferably in real-time and with submillimeter resolution, is desired to reduce the present uncertainty in dose localization. Acoustical imaging technologies exploiting possible local interactions between protons and microbubbles or nanodroplets might be an interesting option. Unfortunately, a theoretical model capable of characterising the acoustical field generated by an individual proton on nanometer and micrometer scales is still missing. In this work, such a model is presented. The proton acoustic field is generated by the adiabatic expansion of a region that is locally heated by a passing proton. To model the proton heat deposition, secondary electron production due to protons has been quantified using a semi-empirical model based on Rutherford's scattering theory, which reproduces experimentally obtained electronic stopping power values for protons in water within 10% over the full energy range. The electrons transfer energy into heat via electron-phonon coupling to atoms along the proton track. The resulting temperature increase is calculated using an inelastic thermal spike model. Heat deposition can be regarded as instantaneous, thus, stress confinement is ensured and acoustical initial conditions are set. The resulting thermoacoustic field in the nanometer and micrometer range from the single proton track is computed by solving the thermoacoustic wave equation using k-space Green's functions, yielding the characteristic amplitudes and frequencies present in the acoustic signal generated by a single proton in an aqueous medium. Wavefield expansion and asymptotic approximations are used to extend the spatial and temporal ranges of the proton acoustic field.

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 speed 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., (© 2022 Institute of Physics and Engineering in Medicine.)

The potential of proton therapy to improve the conformity of the delivered dose to the tumor volume is currently limited by range uncertainties. Injectable superheated nanodroplets have recently been proposed for ultrasound-based in vivo range verification, as these vaporize into echogenic microbubbles on proton irradiation. In previous studies, offline ultrasound images of phantoms with dispersed nanodroplets were acquired after irradiation, relating the induced vaporization profiles to the proton range. However, the aforementioned method did not enable the counting of individual vaporization events, and offline imaging cannot provide real-time feedback. In this study, we overcame these limitations using high-frame-rate ultrasound imaging with a linear array during proton irradiation of phantoms with dispersed perfluorobutane nanodroplets at 37°C and 50°C. Differential image analysis of subsequent frames allowed us to count individual vaporization events and to localize them with a resolution beyond the ultrasound diffraction limit, enabling spatial and temporal quantification of the interaction between ionizing radiation and nanodroplets. Vaporization maps were found to accurately correlate with the stopping distribution of protons (at 50°C) or secondary particles (at both temperatures). Furthermore, a linear relationship between the vaporization count and the number of incoming protons was observed. These results indicate the potential of real-time high-frame-rate contrast-enhanced ultrasound imaging for proton range verification and dosimetry., Competing Interests: Conflict of interest disclosure The authors declare no conflict of interest., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Intraosseous blood circulation is thought to have a critical role in bone growth and remodeling, fracture healing, and bone disorders. However, it is rarely considered in clinical practice because of the absence of a suitable noninvasive in vivo measurement technique. In this work, we assessed blood perfusion in tibial cortical bone simultaneously with blood flow in the superficial femoral artery with ultrasound imaging in five healthy volunteers. After suppression of stationary signal with singular-value-decomposition, pulsatile blood flow in cortical bone tissue is revealed, following the heart rate measured in the femoral artery. Using a method combining transverse oscillations and phase-based motion estimation, 2D vector flow was obtained in the cortex of the tibia. After spatial averaging over the cortex, the peak blood velocity along the long axis of the tibia was measured at four times larger than the peak blood velocity across the bone cortex. This suggests that blood flow in central (Haversian) canals is larger than in perforating (Volkmann's) canals, as expected from the intracortical vascular organization in humans. The peak blood velocity indicates a flow from the endosteum to the periosteum and from the heart to the foot for all subjects. Because aging and the development of bone disorders are thought to modify the direction and velocity of intracortical blood flow, their quantification is crucial. This work reports for the first time an in vivo quantification of the direction and velocity of blood flow in human cortical bone. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research., (© 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.)

Chemotherapy efficacy is often reduced by insufficient drug uptake in tumor cells. The combination of ultrasound and microbubbles (USMB) has been shown to improve drug delivery and to enhance the efficacy of several drugs in vitro and in vivo , through effects collectively known as sonopermeation. However, clinical translation of USMB therapy is hampered by the large variety of (non-clinical) US set-ups and US parameters that are used in these studies, which are not easily translated to clinical practice. In order to facilitate clinical translation, the aim of this study was to prove that USMB therapy using a clinical ultrasound system (Philips iU22) in combination with clinically approved microbubbles (SonoVue) leads to efficient in vitro sonopermeation. To this end, we measured the efficacy of USMB therapy for different US probes (S5-1, C5-1 and C9-4) and US parameters in FaDu cells. The US probe with the lowest central frequency (i.e. 1.6 MHz for S5-1) showed the highest USMB-induced intracellular uptake of the fluorescent dye SYTOX™ Green (SG). These SG uptake levels were comparable to or even higher than those obtained with a custom-built US system with optimized US parameters. Moreover, USMB therapy with both the clinical and the custom-built US system increased the cytotoxicity of the hydrophilic drug bleomycin. Our results demonstrate that a clinical US system can be used to perform USMB therapy as efficiently as a single-element transducer set-up with optimized US parameters. Therefore, future trials could be based on these clinical US systems, including validated US parameters, in order to accelerate successful translation of USMB therapy., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 de Maar, Rousou, van Elburg, Vos, Lajoinie, Bos, Moonen and Deckers.)

Volumetric ultrasound imaging of blood flow with microbubbles enables a more complete visualization of the microvasculature. Sparse arrays are ideal candidates to perform volumetric imaging at reduced manufacturing complexity and cable count. However, due to the small number of transducer elements, sparse arrays often come with high clutter levels, especially when wide beams are transmitted to increase the frame rate. In this study, we demonstrate with a prototype sparse array probe and a diverging wave transmission strategy, that a uniform transmission field can be achieved. With the implementation of a spatial coherence beamformer, the background clutter signal can be effectively suppressed, leading to a signal to background ratio improvement of 25 dB. With this approach, we demonstrate the volumetric visualization of single microbubbles in a tissue-mimicking phantom as well as vasculature mapping in a live chicken embryo chorioallantoic membrane.

Shear wave elastography (SWE) has the potential to determine cardiac tissue stiffness from non-invasive shear wave speed measurements, important, e.g., for predicting heart failure. Previous studies showed that waves traveling in the interventricular septum (IVS) may display Lamb-like dispersive behaviour, introducing a thickness-frequency dependency in the wave speed. However, the IVS tapers across its length, which complicates wave speed estimation by introducing an additional variable to account for. The goal of this work is to assess the impact of tapering thickness on SWE. The investigation is performed by combining in vitro experiments with acoustic radiation force (ARF) and 2D finite element simulations, to isolate the effect of the tapering curve on ARF-induced and natural waves in the heart. The experiments show a 11% deceleration during propagation from the thick to the thin end of an IVS-mimicking tapered phantom plate. The numerical analysis shows that neglecting the thickness variation in the wavenumber-frequency domain can introduce errors of more than 30% in the estimation of the shear modulus, and that the exact tapering curve, rather than the overall thickness reduction, determines the dispersive behaviour of the wave. These results suggest that septal geometry should be accounted for when deriving cardiac stiffness with SWE.

Collateral damage to healthy surrounding tissue during conventional radiotherapy increases when deviations from the treatment plan occur. Ultrasound contrast agents (UCAs) are a possible candidate for radiation dose monitoring. This study investigated the size distribution and acoustic response of two commercial formulations, SonoVue/Lumason and Definity/Luminity, as a function of dose on clinical megavoltage photon beam exposure (24 Gy). SonoVue samples exhibited a decrease in concentration of bubbles smaller than 7 µm, together with an increase in acoustic attenuation and a decrease in acoustic scattering. Definity samples did not exhibit a significant response to radiation, suggesting that the effect of megavoltage photons depends on the UCA formulation. For SonoVue, the influence of the megavoltage photon beam was especially apparent at the second harmonic frequency, and can be captured using pulse inversion and amplitude modulation (3.5-dB decrease for the maximum dose), which could eventually be used for dosimetry in a well-controlled environment., Competing Interests: Conflict of interest disclosure The authors declare no conflict of interest., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

High-frame-rate (HFR) echo-particle image velocimetry (echoPIV) is a promising tool for measuring intracardiac blood flow dynamics. In this study, we investigate the optimal ultrasound contrast agent (UCA: SonoVue) infusion rate and acoustic output to use for HFR echoPIV (PRF = 4900 Hz) in the left ventricle (LV) of patients. Three infusion rates (0.3, 0.6, and 1.2 ml/min) and five acoustic output amplitudes (by varying transmit voltage: 5, 10, 15, 20, and 30 V-corresponding to mechanical indices of 0.01, 0.02, 0.03, 0.04, and 0.06 at 60-mm depth) were tested in 20 patients admitted for symptoms of heart failure. We assess the accuracy of HFR echoPIV against pulsed-wave Doppler acquisitions obtained for mitral inflow and aortic outflow. In terms of image quality, the 1.2-ml/min infusion rate provided the highest contrast-to-background ratio (CBR) (3-dB improvement over 0.3 ml/min). The highest acoustic output tested resulted in the lowest CBR. Increased acoustic output also resulted in increased microbubble disruption. For the echoPIV results, the 1.2-ml/min infusion rate provided the best vector quality and accuracy; mid-range acoustic outputs (corresponding to 15-20-V transmit voltages) provided the best agreement with the pulsed-wave Doppler. Overall, the highest infusion rate (1.2 ml/min) and mid-range acoustic output amplitudes provided the best image quality and echoPIV results.

Background and Aim: Newborns needing major surgical intervention are at risk of brain injury and impaired neurodevelopment later in life. Disturbance of cerebral perfusion might be an underlying factor. This study investigates the feasibility of serial transfontanellar ultrasound measurements of the pial arteries during neonatal surgery, and whether perioperative changes in cerebral perfusion can be observed and related to changes in the perioperative management. Methods: In this prospective, observational feasibility study, neonates with congenital diaphragmatic hernia and esophageal atresia scheduled for surgical treatment within the first 28 days of life were eligible for inclusion. We performed transfontanellar directional power Doppler and pulsed wave Doppler ultrasound during major high-risk non-cardiac neonatal surgery. Pial arteries were of interest for the measurements. Extracted Doppler ultrasound parameters were: peak systolic velocity, end diastolic velocity, the resistivity index and pulsatility index. Results: In 10 out of 14 patients it was possible to perform perioperative measurements; the others failed for logistic and technical reasons. In 6 out of 10 patients, it was feasible to perform serial intraoperative transfontanellar ultrasound measurements with directional power Doppler and pulsed wave Doppler of the same pial artery during neonatal surgery. Median peak systolic velocity was ranging between 5.7 and 7.0 cm s -1 and end diastolic velocity between 1.9 and 3.2 cm s -1 . In patients with a vasoactive-inotropic score below 12 the trend of peak systolic velocity and end diastolic velocity corresponded with the mean arterial blood pressure trend. Conclusion: Perioperative transfontanellar ultrasound Doppler measurements of the pial arteries are feasible and provide new longitudinal data about perioperative cortical cerebral blood flow velocity. Trial Registration: https://www.trialregister.nl/trial/6972, identifier: NL6972., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Costerus, Kortenbout, Vos, Govaert, Tibboel, Wijnen, de Jong, Bosch and de Graaff.)

Cardiac function and vascular function are closely related to the flow of blood within. The flow velocities in these larger cavities easily reach 1 m/s, and generally complex spatiotemporal flow patterns are involved, especially in a non-physiologic state. Visualization of such flow patterns using ultrasound can be greatly enhanced by administration of contrast agents. Tracking the high-velocity complex flows is challenging with current clinical echographic tools, mostly because of limitations in signal-to-noise ratio; estimation of lateral velocities; and/or frame rate of the contrast-enhanced imaging mode. This review addresses the state of the art in 2-D high-frame-rate contrast-enhanced echography of ventricular and deep-vessel flow, from both technological and clinical perspectives. It concludes that current advanced ultrasound equipment is technologically ready for use in human contrast-enhanced studies, thus potentially leading to identification of the most clinically relevant flow parameters for quantifying cardiac and vascular function., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Natural and active shear wave elastography (SWE) are potential ultrasound-based techniques to non-invasively assess myocardial stiffness, which could improve current diagnosis of heart failure. This study aims to bridge the knowledge gap between both techniques and discuss their respective impacts on cardiac stiffness evaluation. We recorded the mechanical waves occurring after aortic and mitral valve closure (AVC, MVC) and those induced by acoustic radiation force throughout the cardiac cycle in four pigs after sternotomy. Natural SWE showed a higher feasibility than active SWE, which is an advantage for clinical application. Median propagation speeds of 2.5-4.0 m/s and 1.6-4.0 m/s were obtained after AVC and MVC, whereas ARF-based median speeds of 0.9-1.2 m/s and 2.1-3.8 m/s were reported for diastole and systole, respectively. The different wave characteristics in both methods, such as the frequency content, complicate the direct comparison of waves. Nevertheless, a good match was found in propagation speeds between natural and active SWE at the moment of valve closure, and the natural waves showed higher propagation speeds than in diastole. Furthermore, the results demonstrated that the natural waves occur in between diastole and systole identified with active SWE, and thus represent a myocardial stiffness in between relaxation and contraction.

Background: Shear waves are generated by the closure of the heart valves. Significant differences in shear wave velocity have been found recently between normal myocardium and disease models of diffusely increased muscle stiffness. In this study we correlate in vivo myocardial shear wave imaging (SWI) with presence of scarred tissue, as model for local increase of stiffness. Stiffness variation is hypothesized to appear as velocity variation., Methods: Ten healthy volunteers (group 1), 10 hypertrophic cardiomyopathy (HCM) patients without any cardiac intervention (group 2), and 10 HCM patients with prior septal reduction therapy (group 3) underwent high frame rate tissue Doppler echocardiography. The SW in the interventricular septum after aortic valve closure was mapped along two M-mode lines, in the inner and outer layer., Results: We compared SWI to 3D echocardiography and strain imaging. In groups 1 and 2, no change in velocity was detected. In group 3, 8/10 patients showed a variation in SW velocity. All three patients having transmural scar showed a simultaneous velocity variation in both layers. Out of six patients with endocardial scar, five showed variations in the inner layer., Conclusion: Local variations in stiffness, with myocardial remodeling post septal reduction therapy as model, can be detected by a local variation in the propagation velocity of naturally occurring shear waves.

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.

Microbubble translations driven by ultrasound-induced radiation forces can be beneficial for applications in ultrasound molecular imaging and drug delivery. Here, the effect of size range in microbubble populations on their translations is investigated experimentally and theoretically. The displacements within five distinct size-isolated microbubble populations are driven by a standard ultrasound-imaging probe at frequencies ranging from 3 to 7 MHz, and measured using the multi-gate spectral Doppler approach. Peak microbubble displacements, reaching up to 10 μm per pulse, are found to describe transient phenomena from the resonant proportion of each bubble population. The overall trend of the statistical behavior of the bubble displacements, quantified by the total number of identified displacements, reveals significant differences between the bubble populations as a function of the transmission frequency. A good agreement is found between the experiments and theory that includes a model parameter fit, which is further supported by separate measurements of individual microbubbles to characterize the viscoelasticity of their stabilizing lipid shell. These findings may help to tune the microbubble size distribution and ultrasound transmission parameters to optimize the radiation-force translations. They also demonstrate a simple technique to characterize the microbubble shell viscosity, the fitted model parameter, from freely floating microbubble populations using a standard ultrasound-imaging probe.

Shear wave elastography (SWE) might allow non-invasive assessment of cardiac stiffness by relating shear wave propagation speed to material properties. However, after aortic valve closure, when natural shear waves occur in the septal wall, the stiffness of the muscle decreases significantly, and the effects of such temporal variation of medium properties on shear wave propagation have not been investigated yet. The goal of this work is to fundamentally investigate these effects. To this aim, qualitative results were first obtained experimentally using a mechanical setup, and were then combined with quantitative results from finite difference simulations. The results show that the amplitude and period of the waves increase during propagation, proportional to the relaxation of the medium, and that reflected waves can originate from the temporal stiffness variation. These general results, applied to literature data on cardiac stiffness throughout the heart cycle, predict as a major effect a period increase of 20% in waves propagating during a healthy diastolic phase, whereas only a 10% increase would result from the impaired relaxation of an infarcted heart. Therefore, cardiac relaxation can affect the propagation of waves used for SWE measurements and might even provide direct information on the correct relaxation of a heart.

Left ventricular (LV) blood flow is an inherently complex time-varying 3-D phenomenon, where 2-D quantification often ignores the effect of out-of-plane motion. In this study, we describe high frame rate 4-D echocardiographic particle image velocimetry (echo-PIV) using a prototype matrix transesophageal transducer and a dynamic LV phantom for testing the accuracy of echo-PIV in the presence of complex flow patterns. Optical time-resolved tomographic PIV (tomo-PIV) was used as a reference standard for comparison. Echo-PIV and tomo-PIV agreed on the general profile of the LV flow patterns, but echo-PIV smoothed out the smaller flow structures. Echo-PIV also underestimated the flow rates at greater imaging depths, where the PIV kernel size and transducer point spread function were large relative to the velocity gradients. We demonstrate that 4-D echo-PIV could be performed in just four heart cycles, which would require only a short breath-hold, providing promising results. However, methods for resolving high velocity gradients in regions of poor spatial resolution are required before clinical translation., Competing Interests: Conflict of interest disclosure The authors declare no competing interests., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

This article quantitatively analyzes the impact of bit errors in digitized RF data on ultrasound image quality. The quality of B-mode images in both linear array and phased array imaging is evaluated by means of three objective image quality metrics: peak signal-to-noise ratio, structural similarity index, and contrast-to-noise ratio, when bit errors are introduced to the RF data with different bit-error rates (BERs). The effectiveness of coding schemes for forward error detection and correction to improve the image quality is also studied. The results show that ultrasound imaging is inherently resilient to high BER. The image quality suffers unnoticeable degradation for BER lower than 1E-6. Simple 1-bit parity coding with 9% added redundancy helps to retain similar image quality for BER up to 1E-4, and Hamming coding with 33.3% added redundancy allows the BER to increase to 1E-3. These results can serve as a guideline in the datalink design for ultrasound probes with in-probe receive digitization. With much more relaxed BER requirements than in typical datalinks, the design can be optimized by allowing fewer cables with higher data rate per cable or lower power consumption with the same cable count.

Due to the primary radiation force, microbubble displacement has been observed previously in the focal region of single-element and array ultrasound probes. This effect has been harnessed to increase the contact between the microbubbles and targeted endothelium for drug delivery and ultrasound molecular imaging. In this study, microbubble displacements associated with plane-wave (PW) transmission are thoroughly investigated and compared to those obtained in focused-wave (FW) transmission over a range of pulse repetition frequencies, burst lengths (BLs), peak negative pressures, and transmission frequencies. In PW mode, the displacements, depending upon the experimental conditions, are in some cases consistently higher (e.g., by 28%, when the longest BL was used at PRF = 4 kHz), and the axial displacements are spatially more uniform compared to FW mode. Statistical analysis on the measured displacements reveals a slightly different frequency dependence of statistical quantities compared to transient peak microbubble displacements, which may suggest the need to consider the size range within the tested microbubble population.

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., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

We apply a high frame rate (over 500 Hz) tissue Doppler method to measure the propagation velocity of naturally occurring shear waves (SW) generated by aortic and mitral valves closure. The aim of this work is to demonstrate clinical relevance. We included 45 healthy volunteers and 43 patients with hypertrophic cardiomyopathy (HCM). The mitral SW (4.68 ± 0.66 m/s) was consistently faster than the aortic (3.51 ± 0.38 m/s) in all volunteers (p < 0.0001). In HCM patients, SW velocity correlated with E/e' ratio (r = 0.346, p = 0.04 for aortic SW and r = 0.667, p = 0.04 for mitral SW). A subgroup of 20 volunteers were matched for age and gender to 20 HCM patients. In HCM, the mean velocity of 5.1 ± 0.7 m/s for the aortic SW (3.61 ± 0.46 m/s in matched volunteers, p < 0.0001) and 6.88 ± 1.12 m/s for the mitral SW(4.65 ± 0.77 m/s in matched volunteers, p < 0.0001). A threshold of 4 m/s for the aortic SW correctly classified pathologic myocardium with a sensitivity of 95% and specificity of 90%. Naturally occurring SW can be used to assess differences between normal and pathologic myocardium., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)