Ultrasonic Sound Speed Estimation for Liver Fat Quantification: A Review by the AIUM-RSNA QIBA Pulse-Echo Quantitative Ultrasound Initiative.

Non-alcoholic fatty liver disease (NAFLD) is a significant cause of diffuse liver disease, morbidity and mortality worldwide. Early and accurate diagnosis of NALFD is critical to identify patients at risk of disease progression. Liver biopsy is the current gold standard for diagnosis and prognosis. However, a non-invasive diagnostic tool is desired because of the high cost and risk of complications of tissue sampling. Medical ultrasound is a safe, inexpensive and widely available imaging tool for diagnosing NAFLD. Emerging sonographic tools to quantitatively estimate hepatic fat fraction, such as tissue sound speed estimation, are likely to improve diagnostic accuracy, precision and reproducibility compared with existing qualitative and semi-quantitative techniques. Various pulse-echo ultrasound speed of sound estimation methodologies have been investigated, and some have been recently commercialized. We review state-of-the-art in vivo speed of sound estimation techniques, including their advantages, limitations, technical sources of variability, biological confounders and existing commercial implementations. We report the expected range of hepatic speed of sound as a function of liver steatosis and fibrosis that may be encountered in clinical practice. Ongoing efforts seek to quantify sound speed measurement accuracy and precision to inform threshold development around meaningful differences in fat fraction and between sequential measurements.

Radiological Society of North America/Quantitative Imaging Biomarker Alliance Shear Wave Speed Bias Quantification in Elastic and Viscoelastic Phantoms.

OBJECTIVES: To quantify the bias of shear wave speed (SWS) measurements between different commercial ultrasonic shear elasticity systems and a magnetic resonance elastography (MRE) system in elastic and viscoelastic phantoms. METHODS: Two elastic phantoms, representing healthy through fibrotic liver, were measured with 5 different ultrasound platforms, and 3 viscoelastic phantoms, representing healthy through fibrotic liver tissue, were measured with 12 different ultrasound platforms. Measurements were performed with different systems at different sites, at 3 focal depths, and with different appraisers. The SWS bias across the systems was quantified as a function of the system, site, focal depth, and appraiser. A single MRE research system was also used to characterize these phantoms using discrete frequencies from 60 to 500 Hz. RESULTS: The SWS from different systems had mean difference 95% confidence intervals of ±0.145 m/s (±9.6%) across both elastic phantoms and ± 0.340 m/s (±15.3%) across the viscoelastic phantoms. The focal depth and appraiser were less significant sources of SWS variability than the system and site. Magnetic resonance elastography best matched the ultrasonic SWS in the viscoelastic phantoms using a 140 Hz source but had a - 0.27 ± 0.027-m/s (-12.2% ± 1.2%) bias when using the clinically implemented 60-Hz vibration source. CONCLUSIONS: Shear wave speed reconstruction across different manufacturer systems is more consistent in elastic than viscoelastic phantoms, with a mean difference bias of < ±10% in all cases. Magnetic resonance elastographic measurements in the elastic and viscoelastic phantoms best match the ultrasound systems with a 140-Hz excitation but have a significant negative bias operating at 60 Hz. This study establishes a foundation for meaningful comparison of SWS measurements made with different platforms.

Tumor Elastography and Its Association with Collagen and the Tumor Microenvironment.

Purpose: The tumor microenvironment presents with altered extracellular matrix (ECM) and stroma composition, which may affect treatment efficacy and contribute to tissue stiffness. Ultrasound (US) elastography can visualize and quantify tissue stiffness noninvasively. However, the contributions of ECM and stromal components to stiffness are poorly understood. We therefore set out to quantify ECM and stroma density and their relation to tumor stiffness.Experimental Design: A modified clinical ultrasound system was used to measure tumor stiffness and perfusion during tumor growth in preclinical tumor models. In vivo measurements were compared with collagen mass spectroscopy and automatic analysis of matrix and stromal markers derived from immunofluorescence images.Results: US elastography estimates of tumor stiffness were positively correlated with tumor volume in collagen and myofibroblast-rich tumors, while no correlations were found for tumors with low collagen and myofibroblast content. US elastography measurements were strongly correlated with ex vivo mechanical testing and mass spectroscopy-based measurements of total collagen and immature collagen crosslinks. Registration of ultrasound and confocal microscopy data showed strong correlations between blood vessel density and T-cell density in syngeneic tumors, while no correlations were found for genetic tumor models. In contrast to collagen density, which was positively correlated with stiffness, no significant correlations were observed for hyaluronic acid density. Finally, localized delivery of collagenase led to a significant reduction in tumor stiffness without changes in perfusion 24 hours after treatment.Conclusions: US elastography can be used as a potential biomarker to assess changes in the tumor microenvironment, particularly changes affecting the ECM. Clin Cancer Res; 24(18); 4455-67. ©2018 AACR.

Identifying Clinically Significant Prostate Cancers using 3-D In Vivo Acoustic Radiation Force Impulse Imaging with Whole-Mount Histology Validation.

Overly aggressive prostate cancer (PCa) treatment adversely affects patients and places an unnecessary burden on our health care system. The inability to identify and grade clinically significant PCa lesions is a factor contributing to excessively aggressive PCa treatment, such as radical prostatectomy, instead of more focal, prostate-sparing procedures such as cryotherapy and high-dose radiation therapy. We have performed 3-D in vivo B-mode and acoustic radiation force impulse (ARFI) imaging using a mechanically rotated, side-fire endorectal imaging array to identify regions suspicious for PCa in 29 patients being treated with radical prostatectomies for biopsy-confirmed PCa. Whole-mount histopathology analyses were performed to identify regions of clinically significant/insignificant PCa lesions, atrophy and benign prostatic hyperplasia. Regions of suspicion for PCa were reader-identified in ARFI images based on boundary delineation, contrast, texture and location. These regions of suspicion were compared with histopathology identified lesions using a nearest-neighbor regional localization approach. Of all clinically significant lesions identified on histopathology, 71.4% were also identified using ARFI imaging, including 79.3% of posterior and 33.3% of anterior lesions. Among the ARFI-identified lesions, 79.3% corresponded to clinically significant PCa lesions, with these lesions having higher indices of suspicion than clinically insignificant PCa. ARFI imaging had greater sensitivity for posterior versus anterior lesions because of greater displacement signal-to-noise ratio and finer spatial sampling. Atrophy and benign prostatic hyperplasia can cause appreciable prostate anatomy distortion and heterogeneity that confounds ARFI PCa lesion identification; however, in general, ARFI regions of suspicion did not coincide with these benign pathologies.

Analyzing the Impact of Increasing Mechanical Index and Energy Deposition on Shear Wave Speed Reconstruction in Human Liver.

Shear wave elasticity imaging (SWEI) has found success in liver fibrosis staging. This work evaluates hepatic SWEI measurement success as a function of push pulse energy using two mechanical index (MI) values (1.6 and 2.2) over a range of pulse durations. Shear wave speed (SWS) was measured in the livers of 26 study subjects with known or potential chronic liver diseases. Each measurement consisted of eight SWEI sequences, each with different push energy configurations. The rate of successful SWS estimation was linearly proportional to the push energy. SWEI measurements with higher push energy were successful in patients for whom standard push energy levels failed. The findings also suggest that liver capsule depth could be used prospectively to identify patients who would benefit from elevated output. We conclude that there is clinical benefit to using elevated acoustic output for hepatic SWS measurement in patients with deeper livers.

Analysis of rapid multi-focal-zone ARFI imaging.

Acoustic radiation force impulse (ARFI) imaging has shown promise for visualizing structure and pathology within multiple organs; however, because the contrast depends on the push beam excitation width, image quality suffers outside of the region of excitation. Multi-focal-zone ARFI imaging has previously been used to extend the region of excitation (ROE), but the increased acquisition duration and acoustic exposure have limited its utility. Supersonic shear wave imaging has previously demonstrated that through technological improvements in ultrasound scanners and power supplies, it is possible to rapidly push at multiple locations before tracking displacements, facilitating extended depth of field shear wave sources. Similarly, ARFI imaging can utilize these same radiation force excitations to achieve tight pushing beams with a large depth of field. Finite element method simulations and experimental data are presented, demonstrating that single- and rapid multi-focal-zone ARFI have comparable image quality (less than 20% loss in contrast), but the multi-focal-zone approach has an extended axial region of excitation. Additionally, as compared with single-push sequences, the rapid multi-focalzone acquisitions improve the contrast-to-noise ratio by up to 40% in an example 4-mm-diameter lesion.

Single- and multiple-track-location shear wave and acoustic radiation force impulse imaging: matched comparison of contrast, contrast-to-noise ratio and resolution.

Acoustic radiation force impulse imaging and shear wave elasticity imaging (SWEI) use the dynamic response of tissue to impulsive mechanical stimulus to characterize local elasticity. A variant of conventional, multiple-track-location SWEI, denoted single-track-location SWEI, offers the promise of creating speckle-free shear wave images. This work compares the three imaging modalities using a high push and track beam density combined acquisition sequence to image inclusions of different sizes and contrasts. Single-track-location SWEI is found to have a significantly higher contrast-to-noise ratio than multiple-track-location SWEI, allowing for operation at higher resolution. Acoustic radiation force impulse imaging and single-track-location SWEI perform similarly in the larger inclusions, with single-track-location SWEI providing better visualization of small targets ≤ 2.5 mm in diameter. The processing of each modality introduces different trade-offs between smoothness and resolution of edges and structures; these are discussed in detail.

Derivation and analysis of viscoelastic properties in human liver: impact of frequency on fibrosis and steatosis staging.

Commercially-available shear wave imaging systems measure group shear wave speed (SWS) and often report stiffness parameters applying purely elastic material models. Soft tissues, however, are viscoelastic, and higher-order material models are necessary to characterize the dispersion associated with broadband shear waves. In this paper, we describe a robust, model-based algorithm and use a linear dispersion model to perform shear wave dispersion analysis in traditionally difficult-to-image subjects. In a cohort of 135 non-alcoholic fatty liver disease patients, we compare the performance of group SWS with dispersion analysis-derived phase velocity c(200 Hz) and dispersion slope dc/df parameters to stage hepatic fibrosis and steatosis. Area under the ROC curve (AUROC) analysis demonstrates correlation between all parameters [group SWS, c(200 Hz), and, to a lesser extent dc/df ] and fibrosis stage, whereas no correlation was observed between steatosis stage and any of the material parameters. Interestingly, optimal AUROC threshold SWS values separating advanced liver fibrosis (≥F3) from mild-to-moderate fibrosis (≤F2) were shown to be frequency-dependent, and to increase from 1.8 to 3.3 m/s over the 0 to 400 Hz shear wave frequency range.

B-mode and acoustic radiation force impulse (ARFI) imaging of prostate zonal anatomy: comparison with 3T T2-weighted MR imaging.

Prostate cancer (PCa) is the most common non-cutaneous malignancy among men in the United States and the second leading cause of cancer-related death. Multi-parametric magnetic resonance imaging (mpMRI) has gained recent popularity to characterize PCa. Acoustic Radiation Force Impulse (ARFI) imaging has the potential to aid PCa diagnosis and management by using tissue stiffness to evaluate prostate zonal anatomy and lesions. MR and B-mode/ARFI in vivo imaging datasets were compared with one another and with gross pathology measurements made immediately after radical prostatectomy. Images were manually segmented in 3D Slicer to delineate the central gland (CG) and prostate capsule, and 3D models were rendered to evaluate zonal anatomy dimensions and volumes. Both imaging modalities showed good correlation between estimated organ volume and gross pathologic weights. Ultrasound and MR total prostate volumes were well correlated (R(2) = 0.77), but B-mode images yielded prostate volumes that were larger (16.82% ± 22.45%) than MR images, due to overestimation of the lateral dimension (18.4% ± 13.9%), with less significant differences in the other dimensions (7.4% ± 17.6%, anterior-to-posterior, and -10.8% ± 13.9%, apex-to-base). ARFI and MR CG volumes were also well correlated (R(2) = 0.85). CG volume differences were attributed to ARFI underestimation of the apex-to-base axis (-28.8% ± 9.4%) and ARFI overestimation of the lateral dimension (21.5% ± 14.3%). B-mode/ARFI imaging yielded prostate volumes and dimensions that were well correlated with MR T2-weighted image (T2WI) estimates, with biases in the lateral dimension due to poor contrast caused by extraprostatic fat. B-mode combined with ARFI imaging is a promising low-cost, portable, real-time modality that can complement mpMRI for PCa diagnosis, treatment planning, and management.

Functional neuroimaging using ultrasonic blood-brain barrier disruption and manganese-enhanced MRI.

Although mice are the dominant model system for studying the genetic and molecular underpinnings of neuroscience, functional neuroimaging in mice remains technically challenging. One approach, Activation-Induced Manganese-enhanced MRI (AIM MRI), has been used successfully to map neuronal activity in rodents. In AIM MRI, Mn(2+) acts a calcium analog and accumulates in depolarized neurons. Because Mn(2+) shortens the T1 tissue property, regions of elevated neuronal activity will enhance in MRI. Furthermore, Mn(2+) clears slowly from the activated regions; therefore, stimulation can be performed outside the magnet prior to imaging, enabling greater experimental flexibility. However, because Mn(2+) does not readily cross the blood-brain barrier (BBB), the need to open the BBB has limited the use of AIM MRI, especially in mice. One tool for opening the BBB is ultrasound. Though potentially damaging, if ultrasound is administered in combination with gas-filled microbubbles (i.e., ultrasound contrast agents), the acoustic pressure required for BBB opening is considerably lower. This combination of ultrasound and microbubbles can be used to reliably open the BBB without causing tissue damage. Here, a method is presented for performing AIM MRI by using microbubbles and ultrasound to open the BBB. After an intravenous injection of perflutren microbubbles, an unfocused pulsed ultrasound beam is applied to the shaved mouse head for 3 minutes. For simplicity, we refer to this technique of BBB Opening with Microbubbles and UltraSound as BOMUS. Using BOMUS to open the BBB throughout both cerebral hemispheres, manganese is administered to the whole mouse brain. After experimental stimulation of the lightly sedated mice, AIM MRI is used to map the neuronal response. To demonstrate this approach, herein BOMUS and AIM MRI are used to map unilateral mechanical stimulation of the vibrissae in lightly sedated mice. Because BOMUS can open the BBB throughout both hemispheres, the unstimulated side of the brain is used to control for nonspecific background stimulation. The resultant 3D activation map agrees well with published representations of the vibrissae regions of the barrel field cortex. The ultrasonic opening of the BBB is fast, noninvasive, and reversible; and thus this approach is suitable for high-throughput and/or longitudinal studies in awake mice.

Acoustic radiation force impulse imaging of human prostates: initial in vivo demonstration.

Reliably detecting prostate cancer (PCa) has been a challenge for current imaging modalities. Acoustic radiation force impulse (ARFI) imaging is an elasticity imaging method that uses remotely generated, focused acoustic beams to probe tissue stiffness. A previous study on excised human prostates demonstrated ARFI images portray various prostatic structures and has the potential to guide prostate needle biopsy with improved sampling accuracy. The goal of this study is to demonstrate the feasibility of ARFI imaging to portray internal structures and PCa in the human prostate in vivo. Custom ARFI imaging sequences were designed and implemented using a modified Siemens Antares™ scanner with a three-dimensional (3-D) wobbler, end-firing, trans-cavity transducer, EV9F4. Nineteen patients were consented and imaged immediately preceding surgical prostatectomy. Pathologies and anatomic structures were identified in histologic slides by a pathologist blinded to ARFI data and were then registered with structures found in ARFI images. The results demonstrated that when PCa is visible, it generally appears as bilaterally asymmetric stiff structures; benign prostatic hyperplasia (BPH) appears heterogeneous with a nodular texture; the verumontanum and ejaculatory ducts appears softer compared with surrounding tissue, which form a unique 'V' shape; and the boundary of the transitional zone (TZ) forms a stiff rim separating the TZ from the peripheral zone (PZ). These characteristic appearances of prostatic structures are consistent with those found in our previous study of prostate ARFI imaging on excised human prostates. Compared with the matched B-mode images, ARFI images, in general, portray prostate structures with higher contrast. With the end-firing transducer used for this study, ARFI depth penetration was limited to 22 mm. Image contrast and resolution were decreased as compared with the previous ex vivo study due to the small transducer aperture. Even with these limitations, this study suggests ARFI imaging holds promise for guidance of targeted prostate needle biopsy and focal therapy, as well as aiding assessment of changes during watchful waiting/active surveillance.

Acoustic Radiation Force Impulse (ARFI) imaging-based needle visualization.

Ultrasound-guided needle placement is widely used in the clinical setting, particularly for central venous catheter placement, tissue biopsy and regional anesthesia. Difficulties with ultrasound guidance in these areas often result from steep needle insertion angles and spatial offsets between the imaging plane and the needle. Acoustic Radiation Force Impulse (ARFI) imaging leads to improved needle visualization because it uses a standard diagnostic scanner to perform radiation force based elasticity imaging, creating a displacement map that displays tissue stiffness variations. The needle visualization in ARFI images is independent of needle-insertion angle and also extends needle visibility out of plane. Although ARFI images portray needles well, they often do not contain the usual B-mode landmarks. Therefore, a three-step segmentation algorithm has been developed to identify a needle in an ARFI image and overlay the needle prediction on a coregistered B-mode image. The steps are: (1) contrast enhancement by median filtration and Laplacian operator filtration, (2) noise suppression through displacement estimate correlation coefficient thresholding and (3) smoothing by removal of outliers and best-fit line prediction. The algorithm was applied to data sets from horizontal 18, 21 and 25 gauge needles between 0-4 mm offset in elevation from the transducer imaging plane and to 18G needles on the transducer axis (in plane) between 10 degrees and 35 degrees from the horizontal. Needle tips were visualized within 2 mm of their actual position for both horizontal needle orientations up to 1.5 mm offset in elevation from the transducer imaging plane and on-axis angled needles between 10 degrees-35 degrees above the horizontal orientation. We conclude that segmented ARFI images overlaid on matched B-mode images hold promise for improved needle visibility in many clinical applications.

GPU-based real-time small displacement estimation with ultrasound.

General purpose computing on graphics processing units (GPUs) has been previously shown to speed up computationally intensive data processing and image reconstruction algorithms for computed tomography (CT), magnetic resonance (MR), and ultrasound images. Although some algorithms in ultrasound have been converted to GPU processing, many investigative ultrasound research systems still use serial processing on a single CPU. One such ultrasound modality is acoustic radiation force impulse (ARFI) imaging, which investigates the mechanical properties of soft tissue. Traditionally, the raw data are processed offline to estimate the displacement of the tissue after the application of radiation force. It is highly advantageous to process the data in real-time to assess their quality and make modifications during a study. In this paper, we present algorithms for efficient GPU parallel processing of two widely used tools in ultrasound: cubic spline interpolation and Loupas' two-dimensional autocorrelator for displacement estimation. It is shown that a commercially available graphics card can be used for these computations, achieving speed increases up to 40x compared with single CPU processing. Thus, we conclude that the GPU-based data processing approach facilitates real-time (i.e., <1 second) display of ARFI data and is a promising approach for ultrasonic research systems.

Contrast-enhanced in vivo magnetic resonance microscopy of the mouse brain enabled by noninvasive opening of the blood-brain barrier with ultrasound.

The use of contrast agents for neuroimaging is limited by the blood-brain barrier (BBB), which restricts entry into the brain. To administer imaging agents to the brain of rats, intracarotid infusions of hypertonic mannitol have been used to open the BBB. However, this technically challenging approach is invasive, opens only a limited region of the BBB, and is difficult to extend to mice. In this work, the BBB was opened in mice, using unfocused ultrasound combined with an injection of microbubbles. This technique has several notable features: it (a) can be performed transcranially in mice; (b) takes only 3 min and uses only commercially available components; (c) opens the BBB throughout the brain; (d) causes no observed histologic damage or changes in behavior (with peak-negative acoustic pressures of 0.36 MPa); and (e) allows recovery of the BBB within 4 h. Using this technique, Gadopentetate Dimeglumine (Gd-DTPA) was administered to the mouse brain parenchyma, thereby shortening T(1) and enabling the acquisition of high-resolution (52 × 52 × 100 micrometers(3)) images in 51 min in vivo. By enabling the administration of both existing anatomic contrast agents and the newer molecular/sensing contrast agents, this technique may be useful for the study of mouse models of neurologic function and pathology with MRI.

A heterogeneous nonlinear attenuating full-wave model of ultrasound.

A full-wave equation that describes nonlinear propagation in a heterogeneous attenuating medium is solved numerically with finite differences in the time domain (FDTD). Three-dimensional solutions of the equation are verified with water tank measurements of a commercial diagnostic ultrasound transducer and are shown to be in excellent agreement in terms of the fundamental and harmonic acoustic fields and the power spectrum at the focus. The linear and nonlinear components of the algorithm are also verified independently. In the linear nonattenuating regime solutions match results from Field II, a well established software package used in transducer modeling, to within 0.3 dB. Nonlinear plane wave propagation is shown to closely match results from the Galerkin method up to 4 times the fundamental frequency. In addition to thermoviscous attenuation we present a numerical solution of the relaxation attenuation laws that allows modeling of arbitrary frequency dependent attenuation, such as that observed in tissue. A perfectly matched layer (PML) is implemented at the boundaries with a numerical implementation that allows the PML to be used with high-order discretizations. A -78 dB reduction in the reflected amplitude is demonstrated. The numerical algorithm is used to simulate a diagnostic ultrasound pulse propagating through a histologically measured representation of human abdominal wall with spatial variation in the speed of sound, attenuation, nonlinearity, and density. An ultrasound image is created in silico using the same physical and algorithmic process used in an ultrasound scanner: a series of pulses are transmitted through heterogeneous scattering tissue and the received echoes are used in a delay-and-sum beam-forming algorithm to generate a images. The resulting harmonic image exhibits characteristic improvement in lesion boundary definition and contrast when compared with the fundamental image. We demonstrate a mechanism of harmonic image quality improvement by showing that the harmonic point spread function is less sensitive to reverberation clutter.