Clinical applications of shear wave dispersion imaging for breast lesions: a pictorial essay

Article information

Ultrasonography. 2023;42(4):589-599
Publication date (electronic) : 2023 August 7
doi : https://doi.org/10.14366/usg.23079
1Department of Radiology, Inha University Hospital and College of Medicine, Incheon, Korea
2Department of Radiology, Korea University Ansan Hospital, Korea University College of Medicine, Ansan, Korea
3Department of Pathology, Korea University Ansan Hospital, Korea University College of Medicine, Ansan, Korea
Correspondence to: Bo Kyoung Seo, MD, PhD, Department of Radiology, Korea University Ansan Hospital, Korea University College of Medicine, 123 Jeokgeum-ro, Danwon-gu, Ansan 15355, Korea Tel. +82-31-412-5229 Fax. +82-31-412-5224 E-mail: seoboky@korea.ac.kr, seoboky@gmail.com
*

These authors contributed equally to this work.

Received 2023 April 27; Revised 2023 August 6; Accepted 2023 August 7.

Abstract

Shear wave dispersion (SWD) imaging is a newly developed ultrasound technology designed to evaluate the dispersion slope of shear waves, which is related to tissue viscosity. This advanced imaging technique holds potential for distinguishing malignant lesions from benign lesions and normal breast tissue. The SWD slope, as determined by shear wave elastography (SWE), could offer crucial insights into the characterization of breast lesions. This article presents SWE and SWD images of both malignant and benign breast lesions, in addition to normal breast tissue.

Introduction

Shear wave elastography (SWE) is widely used to evaluate breast lesions [1,2]. The addition of SWE to B-mode ultrasound (US) can improve the diagnostic performance of US in differentiating benign and malignant breast lesions [3-5]. SWE is a US technique that quantifies tissue elasticity by using acoustic radiation forces to generate laterally propagating shear waves. It can measure the propagating shear wave velocity, which differs depending on tissue stiffness [1-5]. Most current SWE approaches use a linear elastic model that only quantifies elasticity, assuming that the tissue is purely elastic. However, human soft tissues show viscoelastic properties.

Shear wave dispersion (SWD) imaging is a newly developed US technique that evaluates the dispersion slope of shear waves, which is related to the viscosity of biological tissues [6-8]. Dispersion refers to the frequency dependence of both the speed and attenuation of shear waves in the viscous component. As the shear wave speed depends on the frequency of shear waves, the dispersion curve can be obtained from an analysis of shear wave speed according to frequency [9]. The dispersion slope is correlated with viscosity. Therefore, an analysis of SWD characteristics can measure tissue viscosity indirectly.

Several recent studies have evaluated the usefulness of SWD imaging in human and rat liver tissues [6-8,10,11]. SWD can be used for evaluating necroinflammation in the liver [6,7]. Lee et al. [7] investigated the role of SWD imaging in detecting allograft damage after liver transplantation. The SWD slope determined using SWE was related to both liver fibrosis and the degree of hepatic necroinflammatory changes, and the SWD slope showed better diagnostic performance than the SWE value in detecting allograft damage after liver transplantation. However, few studies have used SWD imaging to evaluate breast lesions [11,12]. Ormachea et al. [11] obtained SWD images to evaluate breast phantoms and two women (one with a fibroadenoma and one with dense breast tissue). Kumar et al. [12] found a significant difference in shear viscosity between malignant and benign breast masses, with a higher median value for the malignant masses.

This pictorial review aims to assess the feasibility of breast viscosity evaluation using SWD imaging in patients with benign and malignant breast lesions. All images were obtained from women who underwent breast imaging at Korea University Ansan Hospital.

SWD Imaging

In 1951, long before the advent of clinical elastography, Hans Oestreicher proposed that tissues behave as linear and isotropic viscoelastic media, with corresponding stiffness and viscosity [13]. Currently, SWE can measure the group velocity of a collection of shear waves with reasonable diagnostic performance. Tissue stiffness, measured with two-dimensional (2D) SWE, is a widely recommended parameter. However, stiffness and viscosity may provide distinct sets of clinical data because viscosity is a qualitatively different parameter.

SWD maps can be created using an imaging technique incorporated into a US system (Aplio i-series, Canon Medical Systems, Tokyo, Japan) [6-8]. This technique can estimate the dispersion slope of shear wave speed versus frequency to evaluate changes in tissue viscosity. SWD processing involves four main steps (Fig. 1). In the first step, the displacement induced by shear waves is obtained using a technique based on color Doppler scanning. Second, displacement at each location is transformed from the time to the frequency domain through the Fourier transform algorithm to estimate the phase changes in shear waves at different frequencies. Third, the shear wave speed is calculated using the phase-difference method between displacements at each frequency from the following equation:

Fig. 1.

Four main steps of shear wave dispersion map processing.

First, the displacement induced by the shear waves is obtained. Second, displacement at each location is transformed from the time to frequency domains by the fast Fourier transform to estimate phase changes in the shear waves at several frequencies. Third, the shear wave speed is calculated using a phase-difference method. Fourth, the gradient of the shear wave speed is calculated based on the distribution of the shear wave speed versus frequency. The calculated gradient values are then superimposed on the measurement locations to create a dispersion map.

2πfLø(nf)=oCs(nf)·(n=0,o1,o2,o3,o...)

where π is the circular constant, f is the frequency, ∆L is the displacement between the frequencies at each point, ∆ø is the difference between the phases at each point, Cs(nf) is the shear speed at each frequency, and ° defines each element in the equation [7]. Finally, the slope of the shear wave speed is calculated based on the distribution of shear wave speed versus frequency. The shear wave viscosity value can be estimated from the slope over the shear wave frequency bandwidth (Fig. 2). The slope of the dispersion slope is then superimposed on each measurement location.

Fig. 2.

Relationship between shear wave speed and frequency.

In perfectly elastic tissue, the shear wave speed is constant regardless of shear wave frequency. However, in viscoelastic tissue, shear wave speed varies depending on shear wave frequency. If the shear elasticity is fixed at 2.0 kPa and the shear viscosity varies from 0.1 to 0.5 Pa·s, the slope increases according to the shear viscosity level. The slope correlates with the viscosity coefficient.

SWD imaging can be activated automatically in SWE mode using the Aplio i-series (Canon Medical Systems). In SWE quad-view mode, it is possible to simultaneously view the shear wave elasticity (elasticity map), SWE arrival time contours (propagation map), grayscale US, and SWD slope (dispersion map) (Fig. 3). The dispersion slope value ([m/s]/kHz), elasticity value (kPa), and their standard deviations are displayed for quantitative assessment.

Fig. 3.

Quad-view mode of shear wave elastography.

In quad-view mode, the elasticity map, propagation map, grayscale image, and dispersion slope map are displayed simultaneously.

SWD Evaluation of Breast Lesions

A radiologist experienced in breast imaging (24 years of experience in breast US) performed 2D SWE and SWD imaging examinations using an Aplio i800 system with a linear probe (US frequency range, 5-18 MHz; PLI-1205BX, Canon Medical Systems). The SWD imaging technology for breast examination was a prototype. For quantitative measurements, the radiologist placed a 2-mm circular region of interest on the fibroglandular breast tissue in ultrasonographic negative patients and the areas with the highest SWE and SWD responses in patients with breast lesions [1]. Subcutaneous fat was selected as the reference area for both the elasticity and dispersion ratios. All SWE and SWD measurements were repeated at least three times per patient [14].

In normal breast tissue, tissue elasticity and viscosity tended to be similar to each other in the fibroglandular tissue, and were slightly higher than the corresponding values in the subcutaneous fat tissue (Fig. 4). In cases of benign breast lesions, the elasticity value was in the range of 10-20 kPa, while the dispersion slope value was approximately 10 (m/s)/kHz, which is similar to that of fibroglandular tissue (Figs. 5, 6). However, the elasticity values of malignant breast lesions notably exceeded 50 kPa, and the dispersion slope values were above 20 (m/s)/kHz in most malignant lesions (Figs. 7-14), which may help distinguish malignant lesions from both benign lesions and normal fibroglandular tissue. Several prior studies have evaluated the SWE characteristics among different histologic types or molecular subtypes of invasive breast cancers [15-18]. However, no studies have assessed these features in SWD imaging.

Fig. 4.

Shear wave dispersion imaging in normal breast tissue.

A circular 2-mm region of interest was placed on the fibroglandular tissue and subcutaneous fat. Homogeneous background echotexture (A) and heterogeneous background echotexture (B) are shown on breast ultrasonography.

Fig. 5.

A 50-year-old woman with a benign breast mass.

Grayscale ultrasound (US) shows an oval, microlobulated, hypoechoic mass at 1 o’clock in the left breast. The stiffness value on the elasticity map is 6.8 kPa, and the dispersion slope value on its map is 7.81 (m/s)/kHz. US-guided core-needle biopsy diagnosed this mass pathologically as sclerosing adenosis with fibrocystic change.

Fig. 6.

A 50-year-old woman with a benign non-mass lesion.

On grayscale ultrasound (US), a non-mass lesion is seen at 2 o’clock in the left breast. The stiffness value on the elasticity map is 6.0 kPa, and the dispersion slope value on its map is 6.05 (m/s)/kHz. US-guided coreneedle biopsy pathologically confirmed that the lesion was adenosis with microcalcifications and fibrocystic change.

Fig. 7.

A 57-year-old woman with ductal carcinoma in situ.

A grayscale ultrasound image shows an oval hypoechoic mass with partially indistinct margins at 3 o’clock in the right breast. The stiffness value on the elasticity map is 57.7 kPa, and the dispersion slope value on its map is 22.23 (m/s)/kHz. The mass was diagnosed as ductal carcinoma in situ by surgical excision.

Fig. 8.

A 48-year-old woman with microinvasive carcinoma.

On grayscale ultrasound (US), an indistinct hypoechoic non-mass lesion with internal calcifications is seen at 3 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 25.5 kPa, and the highest dispersion slope value on its map (T2) is 22.62 (m/s)/kHz. US-guided core-needle biopsy diagnosed the lesion as ductal carcinoma in situ, and surgical excision confirmed microinvasive carcinoma.

Fig. 9.

A 50-year-old woman with invasive lobular carcinoma, luminal A subtype.

A grayscale ultrasound image shows an irregular hypoechoic mass with angular margins at 12 o’clock in the left breast. The stiffness value on the elasticity map is 14.3 kPa, and the dispersion slope value on its map is 23.98 (m/s)/kHz. The breast cancer molecular subtype was luminal A (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index <5%).

Fig. 10.

A 49-year-old woman with invasive ductal carcinoma, luminal A subtype.

A grayscale ultrasound image demonstrates an irregular indistinct hypoechoic mass at 1 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 31.7 kPa, and the highest dispersion slope value on its map (T2) is 21.96 (m/s)/kHz. The breast cancer molecular subtype was luminal A (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index <10%).

Fig. 11.

A 67-year-old woman with invasive ductal carcinoma, luminal B subtype.

A grayscale ultrasound image shows a spiculated irregular hypoechoic mass at 10 o’clock in the right breast. The stiffness value on the elasticity map is 57.8 kPa, and the dispersion slope value on its map is 24.47 (m/s)/kHz. The breast cancer molecular subtype was luminal B (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index >14%).

Fig. 12.

A 62-year-old woman with mucinous carcinoma, luminal B subtype.

A grayscale ultrasound image shows an indistinct oval isoechoic mass at 7 o’clock direction of the right breast. The highest stiffness value on the elasticity map (T1) was 136.7 kPa, and the highest value of dispersion slope on its map was 48.10 (m/s)/kHz. The area showing low value on the elasticity map is correlated with the area of high value on dispersion map (yellow arrows), indicating high viscosity of the mass. Surgical pathology revealed a 1.8-cm mucinous carcinoma, luminal B subtype (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index >14%).

Fig. 13.

A 40-year-old woman with invasive ductal carcinoma, human epidermal growth factor receptor 2 (HER2)–enriched subtype.

A grayscale ultrasound image shows an irregular, microlobulated, hypoechoic mass at 3 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 151.1 kPa, and the highest dispersion slope value on its map (T2) is 50.09 (m/s)/kHz. The molecular subtype was the HER-enriched subtype (estrogen receptor-negative, progesterone receptor-negative, and human epidermal growth factor receptor 2-positive).

Fig. 14.

A 51-year-old woman with invasive ductal carcinoma, triple-negative subtype.

A grayscale ultrasound image shows an indistinct oval hypoechoic mass with a nonparallel orientation at 12 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 160.3 kPa, and the highest dispersion slope value on its map (T2) is 60.70 (m/s)/kHz. Surgical pathology revealed a 2.5-cm invasive ductal carcinoma of no specific type, which was confirmed as triple-negative (estrogen receptor–negative, progesterone receptor–negative, and human epidermal growth factor receptor 2–negative).

Conclusion

This pictorial essay shows the SWE and SWD images of benign and malignant lesions on breast US. The SWD slope determined by 2D SWE may provide more detailed and useful information regarding the characterization of breast lesions. We suggest that SWD imaging can help differentiate malignant from benign breast lesions. This possibility needs to be verified in further studies with larger numbers of participants.

Notes

Author Contributions

Conceptualization: Seo BK. Data acquisition: Seo BK. Data analysis or interpretation: Bae MS, Kim HY, Oh H, Seo BK. Drafting of the manuscript: Bae MS, Kim HY, Oh H, Seo BK. Critical revision of the manuscript: Bae MS, Seo BK. Approval of the final version of the manuscript: all authors.

No potential conflict of interest relevant to this article was reported.

Acknowledgements

This research was supported by a National Research Foundation of Korea (NRF) grant (No. NRF-2021R1A2C1010565) funded by the Korean government (Ministry of Science and ICT) and Canon Medical Systems.

References

1. Youk JH, Gweon HM, Son EJ. Shear-wave elastography in breast ultrasonography: the state of the art. Ultrasonography 2017;36:300–309.
2. Pillai A, Voruganti T, Barr R, Langdon J. Diagnostic accuracy of shear-wave elastography for breast lesion characterization in women: a systematic review and meta-analysis. J Am Coll Radiol 2022;19:625–634.
3. Park SY, Kang BJ. Combination of shear-wave elastography with ultrasonography for detection of breast cancer and reduction of unnecessary biopsies: a systematic review and meta-analysis. Ultrasonography 2021;40:318–332.
4. Berg WA, Cosgrove DO, Dore CJ, Schafer FK, Svensson WE, Hooley RJ, et al. Shear-wave elastography improves the specificity of breast US: the BE1 multinational study of 939 masses. Radiology 2012;262:435–449.
5. Lee SH, Chang JM, Kim WH, Bae MS, Seo M, Koo HR, et al. Added value of shear-wave elastography for evaluation of breast masses detected with screening US imaging. Radiology 2014;273:61–69.
6. Sugimoto K, Moriyasu F, Oshiro H, Takeuchi H, Yoshimasu Y, Kasai Y, et al. Clinical utilization of shear wave dispersion imaging in diffuse liver disease. Ultrasonography 2020;39:3–10.
7. Lee DH, Lee JY, Bae JS, Yi NJ, Lee KW, Suh KS, et al. Shear-wave dispersion slope from US shear-wave elastography: detection of allograft damage after liver transplantation. Radiology 2019;293:327–333.
8. Zhang XQ, Zheng RQ, Jin JY, Wang JF, Zhang T, Zeng J. US Shear-wave elastography dispersion for characterization of chronic liver disease. Radiology 2022;305:597–605.
9. Chen S, Urban MW, Pislaru C, Kinnick R, Greenleaf JF. Liver elasticity and viscosity quantification using shearwave dispersion ultrasound vibrometry (SDUV). Annu Int Conf IEEE Eng Med Biol Soc 2009;2009:2252–2255.
10. Sugimoto K, Moriyasu F, Oshiro H, Takeuchi H, Yoshimasu Y, Kasai Y, et al. Viscoelasticity measurement in rat livers using shear-wave US elastography. Ultrasound Med Biol 2018;44:2018–2024.
11. Ormachea J, Parker KJ, Barr RG. An initial study of complete 2D shear wave dispersion images using a reverberant shear wave field. Phys Med Biol 2019;64:145009.
12. Kumar V, Denis M, Gregory A, Bayat M, Mehrmohammadi M, Fazzio R, et al. Viscoelastic parameters as discriminators of breast masses: Initial human study results. PLoS One 2018;13e0205717.
13. Carstensen EL, Parker KJ. Oestreicher and elastography. J Acoust Soc Am 2015;138:2317–2325.
14. Barr RG, Nakashima K, Amy D, Cosgrove D, Farrokh A, Schafer F, et al. WFUMB guidelines and recommendations for clinical use of ultrasound elastography: Part 2: breast. Ultrasound Med Biol 2015;41:1148–1160.
15. Evans A, Sim YT, Thomson K, Jordan L, Purdie C, Vinnicombe SJ. Shear wave elastography of breast cancer: sensitivity according to histological type in a large cohort. Breast 2016;26:115–118.
16. Evans A, Whelehan P, Thomson K, McLean D, Brauer K, Purdie C, et al. Invasive breast cancer: relationship between shear-wave elastographic findings and histologic prognostic factors. Radiology 2012;263:673–677.
17. Chang JM, Park IA, Lee SH, Kim WH, Bae MS, Koo HR, et al. Stiffness of tumours measured by shear-wave elastography correlated with subtypes of breast cancer. Eur Radiol 2013;23:2450–2458.
18. Kim H, Lee J, Kang BJ, Kim SH. What shear wave elastography parameter best differentiates breast cancer and predicts its histologic aggressiveness? Ultrasonography 2021;40:265–273.

Article information Continued

Notes

Key point

Shear wave dispersion (SWD) imaging is a new ultrasound technique for assessing the dispersion slope of shear waves. This article reports the authors’ preliminary clinical experience with SWD imaging for breast lesions.

Fig. 1.

Four main steps of shear wave dispersion map processing.

First, the displacement induced by the shear waves is obtained. Second, displacement at each location is transformed from the time to frequency domains by the fast Fourier transform to estimate phase changes in the shear waves at several frequencies. Third, the shear wave speed is calculated using a phase-difference method. Fourth, the gradient of the shear wave speed is calculated based on the distribution of the shear wave speed versus frequency. The calculated gradient values are then superimposed on the measurement locations to create a dispersion map.

Fig. 2.

Relationship between shear wave speed and frequency.

In perfectly elastic tissue, the shear wave speed is constant regardless of shear wave frequency. However, in viscoelastic tissue, shear wave speed varies depending on shear wave frequency. If the shear elasticity is fixed at 2.0 kPa and the shear viscosity varies from 0.1 to 0.5 Pa·s, the slope increases according to the shear viscosity level. The slope correlates with the viscosity coefficient.

Fig. 3.

Quad-view mode of shear wave elastography.

In quad-view mode, the elasticity map, propagation map, grayscale image, and dispersion slope map are displayed simultaneously.

Fig. 4.

Shear wave dispersion imaging in normal breast tissue.

A circular 2-mm region of interest was placed on the fibroglandular tissue and subcutaneous fat. Homogeneous background echotexture (A) and heterogeneous background echotexture (B) are shown on breast ultrasonography.

Fig. 5.

A 50-year-old woman with a benign breast mass.

Grayscale ultrasound (US) shows an oval, microlobulated, hypoechoic mass at 1 o’clock in the left breast. The stiffness value on the elasticity map is 6.8 kPa, and the dispersion slope value on its map is 7.81 (m/s)/kHz. US-guided core-needle biopsy diagnosed this mass pathologically as sclerosing adenosis with fibrocystic change.

Fig. 6.

A 50-year-old woman with a benign non-mass lesion.

On grayscale ultrasound (US), a non-mass lesion is seen at 2 o’clock in the left breast. The stiffness value on the elasticity map is 6.0 kPa, and the dispersion slope value on its map is 6.05 (m/s)/kHz. US-guided coreneedle biopsy pathologically confirmed that the lesion was adenosis with microcalcifications and fibrocystic change.

Fig. 7.

A 57-year-old woman with ductal carcinoma in situ.

A grayscale ultrasound image shows an oval hypoechoic mass with partially indistinct margins at 3 o’clock in the right breast. The stiffness value on the elasticity map is 57.7 kPa, and the dispersion slope value on its map is 22.23 (m/s)/kHz. The mass was diagnosed as ductal carcinoma in situ by surgical excision.

Fig. 8.

A 48-year-old woman with microinvasive carcinoma.

On grayscale ultrasound (US), an indistinct hypoechoic non-mass lesion with internal calcifications is seen at 3 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 25.5 kPa, and the highest dispersion slope value on its map (T2) is 22.62 (m/s)/kHz. US-guided core-needle biopsy diagnosed the lesion as ductal carcinoma in situ, and surgical excision confirmed microinvasive carcinoma.

Fig. 9.

A 50-year-old woman with invasive lobular carcinoma, luminal A subtype.

A grayscale ultrasound image shows an irregular hypoechoic mass with angular margins at 12 o’clock in the left breast. The stiffness value on the elasticity map is 14.3 kPa, and the dispersion slope value on its map is 23.98 (m/s)/kHz. The breast cancer molecular subtype was luminal A (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index <5%).

Fig. 10.

A 49-year-old woman with invasive ductal carcinoma, luminal A subtype.

A grayscale ultrasound image demonstrates an irregular indistinct hypoechoic mass at 1 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 31.7 kPa, and the highest dispersion slope value on its map (T2) is 21.96 (m/s)/kHz. The breast cancer molecular subtype was luminal A (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index <10%).

Fig. 11.

A 67-year-old woman with invasive ductal carcinoma, luminal B subtype.

A grayscale ultrasound image shows a spiculated irregular hypoechoic mass at 10 o’clock in the right breast. The stiffness value on the elasticity map is 57.8 kPa, and the dispersion slope value on its map is 24.47 (m/s)/kHz. The breast cancer molecular subtype was luminal B (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index >14%).

Fig. 12.

A 62-year-old woman with mucinous carcinoma, luminal B subtype.

A grayscale ultrasound image shows an indistinct oval isoechoic mass at 7 o’clock direction of the right breast. The highest stiffness value on the elasticity map (T1) was 136.7 kPa, and the highest value of dispersion slope on its map was 48.10 (m/s)/kHz. The area showing low value on the elasticity map is correlated with the area of high value on dispersion map (yellow arrows), indicating high viscosity of the mass. Surgical pathology revealed a 1.8-cm mucinous carcinoma, luminal B subtype (estrogen receptor–positive, progesterone receptor–positive, human epidermal growth factor receptor 2–negative, and Ki-67 index >14%).

Fig. 13.

A 40-year-old woman with invasive ductal carcinoma, human epidermal growth factor receptor 2 (HER2)–enriched subtype.

A grayscale ultrasound image shows an irregular, microlobulated, hypoechoic mass at 3 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 151.1 kPa, and the highest dispersion slope value on its map (T2) is 50.09 (m/s)/kHz. The molecular subtype was the HER-enriched subtype (estrogen receptor-negative, progesterone receptor-negative, and human epidermal growth factor receptor 2-positive).

Fig. 14.

A 51-year-old woman with invasive ductal carcinoma, triple-negative subtype.

A grayscale ultrasound image shows an indistinct oval hypoechoic mass with a nonparallel orientation at 12 o’clock in the left breast. The highest stiffness value on the elasticity map (T1) is 160.3 kPa, and the highest dispersion slope value on its map (T2) is 60.70 (m/s)/kHz. Surgical pathology revealed a 2.5-cm invasive ductal carcinoma of no specific type, which was confirmed as triple-negative (estrogen receptor–negative, progesterone receptor–negative, and human epidermal growth factor receptor 2–negative).